Программное обеспечение Mapping Toolbox™ поддерживает три различных интерпретации геоориентированных сетей данных:
Во-первых, матрица отображения, имеющая такое же количество строк и столбцов, как матрицы координат широты и долготы, представляет значения данных отображения в соответствующих географических точках (центрах ячеек данных).
Затем матрица карты, имеющая на одну строку меньше и на один столбец меньше, чем матрицы географических координат, представляет значения данных карты в области, образованной четырьмя соседними широтами и долготами.
Наконец, если матрицы широты и долготы имеют меньшие размеры, чем матрица карты, их можно интерпретировать как описание более грубой решетки или сетки ячеек широты и долготы, в которые деформируются блоки данных карты.
В этом разделе рассматриваются первые две интерпретации геоориентированных сетей данных. Дополнительные сведения об использовании гратикулов см. в разделе Сетка карты.
В качестве примера первой интерпретации рассмотрим матрицу карты 4 на 4, размер ячейки которой равен 30 на 30 градусов, вместе с соответствующими матрицами широты и долготы 4 на 4:
Z = [ ... 1 2 3 4; ... 5 6 7 8; ... 9 10 11 12; ... 13 14 15 16]; lat = [ ... 30 30 30 30; ... 0 0 0 0; ... -30 -30 -30 -30; ... -60 -60 -60 -60]; lon = [ ... 0 30 60 90;... 0 30 60 90;... 0 30 60 90;... 0 30 60 90];
Отобразить геоориентированную сетку данных со значениями map показаны в соответствующих широтах и долготах:
figure('Color','white','Colormap',autumn(64))
axesm('pcarree','Grid','on','Frame','on',...
'PLineLocation',30,'PLabelLocation',30)
mlabel; plabel; axis off; tightmap
h = geoshow(lat,lon,Z,'DisplayType','surface');
set(h,'ZData',zeros(size(Z)))
ht = textm(lat(:),lon(:),num2str(Z(:)), ...
'Color','blue','FontSize',14);
colorbar('southoutside')
Обратите внимание, что отображаются только 9 из 16 ячеек. Значение, отображаемое для каждой ячейки, является значением в левом верхнем углу этой ячейки, координаты которой задаются соответствующим lat и lon элементы. По соглашению, последняя строка и столбец матрицы карты не отображаются, хотя они существуют в CData свойства объекта поверхности.
Для второй интерпретации рассмотрим матрицу карты 3 на 3 с той же самой lat и lon переменные:
delete(h)
delete(ht)
Z3by3 = [ ...
1 2 3; ...
4 5 6; ...
7 8 9];
h = geoshow(lat,lon,Z3by3,'DisplayType','texturemap');
tlat = [ ...
15 15 15; ...
-15 -15 -15; ...
-45 -45 -45];
tlon = [ ...
15 45 75; ...
15 45 75; ...
15 45 74];
textm(tlat(:),tlon(:),num2str(Z3by3(:)), ...
'Color','blue','FontSize',14)
Отображение графика поверхности матрицы карты со значениями map показано в центре ассоциированных ячеек:
Все данные карты отображаются для этой географической сетки данных. Значение каждой ячейки является значением в центре ячейки, а широты и долготы в матрицах координат являются границами ячеек.
Возможно, вы заметили, что первая строка матрицы отображается в верхней части карты, тогда как для обычной сетки данных было верно обратное: первая строка соответствовала нижней части карты. Это различие полностью связано с тем, как lat и lon матрицы упорядочены. В геоориентированной сетке данных порядок значений в двух матрицах координат определяет расположение отображаемых значений.
При необходимости регулярная сетка данных может быть преобразована в геоориентированную сетку данных. Это просто требует, чтобы пара матриц координат была вычислена с желаемым пространственным разрешением из регулярной сетки. Например, загрузить растровые данные отметки и объект ссылки на географические ячейки. Затем создайте сетку широты-долготы для данных, вызвав geographicGrid функция.
load topo60c
[lat,lon] = geographicGrid(topo60cR);
И наоборот, регулярная сетка данных может быть также построена из геоориентированной сетки данных. Координаты и значения могут быть встроены в новую регулярную сетку данных. Функция, выполняющая это преобразование: geoloc2grid; в качестве входных данных используется геоориентированная сетка данных и размер ячейки.