Географические интерпретации определенных геолокацию сеток

Программное обеспечение Mapping Toolbox™ поддерживает три различных интерпретации определенных геолокацию сеток данных:

  • Во-первых, матрица карты, имеющая одинаковое число строк и столбцы как матрицы координаты широты и долготы, представляет значения данных о карте в соответствующих географических точках (центры ячеек данных).

  • Затем, матричное ссорение карты и один меньше столбца, чем географические координатные матрицы представляет значения данных о карте в области, сформированной четырьмя смежными широтами и долготами.

  • Наконец, если матрицы широты и долготы имеют меньшие размерности, чем матрица карты, можно интерпретировать их как описание более грубой координатной сетки или сетки ячеек широты и долготы, в которые деформированы блоки данных о карте.

В этом разделе рассматриваются первые две интерпретации определенных геолокацию сеток данных. Для получения дополнительной информации об использовании координатных сеток смотрите Сетку Карты.

Тип 1: значения, связанные с верхней левой декартовой координатой

Как пример первой интерпретации, рассмотрите матрицу карты 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 объекта подложки.

Тип 2: значения, в центре в четырех смежных координатах

Для второй интерпретации рассмотрите 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 упорядочены. В определенной геолокацию сетке данных порядок значений в двух координатных матрицах определяет расположение отображенных значений.

Преобразуйте регулярный к определенным геолокацию сеткам

При необходимости обычная сетка данных может быть преобразована в определенную геолокацию сетку данных. Это просто требует, чтобы пара матриц координат была вычислена при желаемом пространственном разрешении из обычной сетки. Сделайте это с функцией meshgrat, можно следующим образом:

load topo
[lat,lon] = meshgrat(topo,topolegend);
  Name              Size              Bytes  Class     Attributes

  lat             180x360            518400  double              
  lon             180x360            518400  double              
  topo            180x360            518400  double              
  topolegend        1x3                  24  double              
  topomap1         64x3                1536  double              
  topomap2        128x3                3072  double              

Преобразование определенного геолокацию к обычным сеткам

С другой стороны обычная сетка данных может также быть создана из определенной геолокацию сетки данных. Координаты и значения могут быть встроены в новую обычную сетку данных. Функцией, которая выполняет это преобразование, является geoloc2grid; это берет определенную геолокацию сетку данных и размер ячейки как входные параметры.