Считывание векторных элементов и атрибутов из файла формы
S = shaperead(filename)filenameи возвращает N-by-1 массив структуры географических данных в спроецированных координатах карты (mapstruct). Структура географических данных объединяет информацию о геометрических и характерных атрибутах. shaperead поддерживает обычные типы форм 2-D: 'Point', 'Multipoint', 'PolyLine', и 'Polygon'.
S = shaperead(filename,Name,Value)S, как определено аргументами пары имя-значение. Использовать RecordNumbers, BoundingBox, и Selector для выбора элементов для чтения. При использовании более одного из этих параметров в одном вызове происходит пересечение записей, соответствующих отдельным спецификациям. Например, если указать значения для обоих RecordNumbers и BoundingBox, импортируются только те элементы с номерами записей, которые отображаются в списке и которые также имеют ограничивающие рамки, пересекающие указанную ограничивающую рамку.
Если параметры не указаны, shaperead возвращает запись для каждого элемента, отличного от NULL, и создает поле для каждого атрибута.
Прочитать весь файл формы под названием concord_hydro_line.shp, включая атрибуты в concord_hydro_line.dbf. shaperead функция возвращает mapstruct с координатными векторами x и y.
S = shaperead('concord_hydro_line.shp')S=237×1 struct array with fields:
Geometry
BoundingBox
X
Y
LENGTH
Укажите ограничительную рамку для ограничения возвращаемых данных shaperead . Кроме того, укажите имена атрибутов, которые требуется прочитать.
bbox = [2.08 9.11; 2.09 9.12] * 1e5; S = shaperead('concord_roads','BoundingBox',bbox,... 'Attributes',{'STREETNAME','CLASS'})
S=87×1 struct array with fields:
Geometry
BoundingBox
X
Y
STREETNAME
CLASS
Считывание дорожных данных только для участков дорог класса 4 и выше длиной не менее 200 метров. Обратите внимание на использование анонимной функции в селекторе.
S = shaperead('concord_roads.shp','Selector',... {@(v1,v2) (v1 >= 4) && (v2 >= 200),'CLASS','LENGTH'})
S=115×1 struct array with fields:
Geometry
BoundingBox
X
Y
STREETNAME
RT_NUMBER
CLASS
ADMIN_TYPE
LENGTH
Определите количество дорог каждого класса.
n = histcounts([S.CLASS],'BinLimits',[1 7],'BinMethod','integer')
n = 1×7
0 0 0 7 93 15 0
Отображение гистограммы количества дорог, приходящихся на каждую категорию длины.
figure histogram([S.LENGTH])
Укажите, что shapefile использует координаты широты и долготы, используя 'UseGeoCoords' пара имя-значение.
Например, вернуть информацию о файле формы в виде структуры. Проверьте, использует ли файл формы координаты широты и долготы, запросив CoordinateReferenceSystem поле. Shapefile использует координаты широты и долготы, если поле содержит geocrs объект.
info = shapeinfo('landareas.shp');
crs = info.CoordinateReferenceSystemcrs =
geocrs with properties:
Name: "WGS 84"
Datum: "World Geodetic System 1984"
Spheroid: [1x1 referenceEllipsoid]
PrimeMeridian: 0
AngleUnit: "degree"
Прочтите файл формы с помощью shaperead функция. Укажите, что в файле формы используются координаты широты и долготы с помощью 'UseGeoCoords' пара имя-значение.
S = shaperead('landareas.shp','UseGeoCoords',true)
S=537×1 struct array with fields:
Geometry
BoundingBox
Lon
Lat
Name
Обратите внимание, что shaperead функция возвращает структуру географических данных с полями широты и долготы (a geostruct).
Для получения информации о системе координат (CRS), связанной с файлом формы, используйте shapeinfo функция.