В дополнение к обычным сеткам данных, инструментарий предоставляет другой формат для геоданных: геолокационные сетки данных. Эти многомерные наборы данных могут быть отображены, и их значения и координаты могут быть запрошены, но, к сожалению, большая часть функциональных возможностей, поддерживающих регулярные сетки данных, не доступна для геоориентированных сетей данных.
Регулярные сетки данных охватывают простые, правильные четырехугольники, то есть географически прямоугольные и выровненные с параллелями и меридианами. Кроме этих прямоугольных ориентаций, геоцентрированные сетки данных могут иметь и другие формы.
Чтобы определить геоориентированную сетку данных, необходимо определить три переменные: матрицу индексов или значений, связанных с отображенной областью, матрицу, дающую координаты широты ячейки за ячейкой, и матрицу, дающую координаты долготы ячейки за ячейкой.
Загрузка 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])
