Ссылочные растровые регистрации, чтобы сопоставить координаты
Объект растровой привязки регистраций карты содержит информацию, которая описывает отношение между плоской системой координат карты и внутренней системой координат, которая привязывается на столбцах и строках 2D сетки, на которую пространственно ссылаются, выборок точки (или регистрации).
Как правило, растр регулярно производится в плоском мировом x - и мировом y - координаты системы карты, такой, что внутренний x и мировые оси x выравниваются и внутренний y, и мировые оси y выравниваются. Когда это верно, отношение между этими двумя системами прямолинейно. В более общем плане, и намного более редко, их отношение является аффинным. Аффинное отношение допускает возможное вращение (и скос). Или в случае, прямолинейном или в аффинном, демонстрационный интервал от строки до строки не должен равняться демонстрационному интервалу от столбца до столбца.
Это изображение иллюстрирует 9 11 плоский растр регистраций. Регистрации расположены с интервалами на расстоянии в 1 метр в x и y, покрывают значения x в области значений [40,50] в метрах и покрывают значения y в области значений [20,28] в метрах. Контур растра регистраций является прямоугольником с ребрами, которые проходят через наиболее удаленные регистрации.
Чтобы сослаться на растр ячеек к плоским координатам карты, используйте MapCellsReference объект вместо этого.
Можно использовать эти функции, чтобы создать MapPostingsReference объекты:
readgeoraster — Импортируйте растровый файл данных как массив и ссылочный объект.
maprefpostings — Создайте объект растровой привязки карты путем определения x и пределов мира y растра и или размер растра или интервал регистраций.
maprasterref — Преобразуйте файл привязки в объект растровой привязки карты.
refmatToMapRasterReference — Преобразуйте матрицу привязки в объект растровой привязки карты.
contains | Определите, сопоставляют ли географический или растр, содержит точки |
firstCornerX | Возвратите мировой x - координата растрового индекса (1,1) карты |
firstCornerY | Возвратите мировой y - координата растрового индекса (1,1) карты |
intrinsicToWorld | Преобразуйте внутренний параметр к плоским мировым координатам |
sizesMatch | Определите, сопоставляют ли географический или растровый объект, и изображение или растр совместимы с размером |
worldFileMatrix | Возвратите параметры файла привязки для преобразования |
worldGrid | Мировые координаты растровых элементов |
worldToDiscrete | Преобразуйте плоский мир к дискретным координатам |
worldToIntrinsic | Преобразуйте плоский мир к внутренним координатам |
Импортируйте данные о вертикальном изменении как массив и MapPostingsReference объект. Затем отобразите данные по карте.
В данном примере задайте файл, чтобы импортировать путем извлечения файла GZ во временную директорию. Затем определите имя первого извлеченного файла к переменной.
filenames = gunzip('sanfranciscos.dem.gz',tempdir);
filename1 = filenames{1};Импортируйте данные при помощи readgeoraster функция. Подготовьте данные к графическому выводу путем определения выходного типа как 'double' и заменяя недостающие данные на NaN значения.
[Z,R] = readgeoraster(filename1,'OutputType','double'); info = georasterinfo(filename1); Z = standardizeMissing(Z,info.MissingDataIndicator);
Просмотрите свойства MapPostingsReference объект.
R
R =
MapPostingsReference with properties:
XWorldLimits: [544050 555150]
YWorldLimits: [4164120 4178070]
RasterSize: [466 371]
RasterInterpretation: 'postings'
ColumnsStartFrom: 'north'
RowsStartFrom: 'west'
SampleSpacingInWorldX: 30
SampleSpacingInWorldY: 30
RasterExtentInWorldX: 11100
RasterExtentInWorldY: 13950
XIntrinsicLimits: [1 371]
YIntrinsicLimits: [1 466]
TransformationType: 'rectilinear'
CoordinateSystemType: 'planar'
ProjectedCRS: [1x1 projcrs]
Отобразите данные о вертикальном изменении как поверхность на карте. Примените палитру, которая подходит для данных о вертикальном изменении.
mapshow(Z,R,'DisplayType','surface') demcmap(Z)
