В дополнение к обычным сеткам данных тулбокс обеспечивает другой формат для геоданных: определенные геолокацию сетки данных. Эти многомерные наборы данных могут быть отображены, и их значения и координаты могут быть запрошены, но к сожалению большая часть функциональности, поддерживающей обычные сетки данных, не доступна для определенных геолокацию сеток данных.
Регулярное покрытие сеток данных простые, регулярные четырехугольники, то есть, географически прямоугольный и выровненный с параллелями и меридианами. Определенные геолокацию сетки данных, в дополнение к этим прямоугольным ориентациям, могут иметь другие формы также.
Чтобы задать определенную геолокацию сетку данных, необходимо задать три переменные: матрица индексов или значений, сопоставленных с сопоставленной областью, матрица, дающая координаты широты ячейки ячейкой и матрицу, дающую координаты долготы ячейки ячейкой.
Загрузите MAT-файл, содержащий определенную геолокацию сетку данных неправильной формы под названием mapmtx
.
load mapmtx
Просмотрите переменные, созданные из этого MAT-файла. Две определенных геолокацию сетки данных находятся в этом наборе данных, каждый требующий трех переменных. Значения содержатся в map1
соответствуйте координатам широты и долготы, соответственно, в lt1
и lg1
. Заметьте, что все три матрицы одного размера. Точно так же map2
, lt2
, и lg2
вместе сформируйте вторую определенную геолокацию сетку данных. Эти наборы данных были извлечены из topo60c
сетка данных. Ни одна из этих карт не обычна, потому что их столбцы не запускаются с севера на юг.
whos
Name Size Bytes Class Attributes description 1x54 108 char lg1 50x50 20000 double lg2 50x50 20000 double lt1 50x50 20000 double lt2 50x50 20000 double map1 50x50 20000 double map2 50x50 20000 double source 1x43 86 char
Отобразите сетки один за другим, чтобы видеть их географию.
close all axesm mercator gridm on framem on h1 = surfm(lt1,lg1,map1); h2 = surfm(lt2,lg2,map2);
Показ береговых линий поможет ориентировать вас на эти скошенные сетки. Заметьте, что никакая сетка не является регулярным прямоугольником. Каждый походит на ромб географически, другой как трапецоид. Трапецоид отображен в двух частях, потому что он пересекает ребро карты. Эти формы могут считаться географической организацией данных, как прямоугольники для обычных сеток данных. Но, так же, как для обычных сеток данных, эта организационная логика не означает, что отображения этих карт являются обязательно определенной формой.
load coastlines plotm(coastlat,coastlon,'r')
Теперь измените представление в поликоническую проекцию с источником в 0°N, 90°E. Когда поликоническая проекция ограничивается областью значений на 150 ° в долготе, те фрагменты карт за пределами этой области автоматически обрезаются.
setm(gca,'MapProjection','polycon','Origin',[0 90 0])