exponenta event banner

Создание глобальной сетки данных

Представьте себе чрезвычайно грубую карту мира, в которой каждая клетка представляет 60 °. Такая матрица карты будет 3 на 6.

  1. Создайте сетку 3 на 6:

    miniZ = [1 2 3 4 5 6; 7 8 9 10 11 12; 13 14 15 16 17 18];
    
  2. Теперь создайте ссылочный объект:

    miniR = georasterref('RasterSize', size(miniZ), ...
       'Latlim', [-90 90], 'Lonlim', [-180 180])
    

    Выходные данные выглядят следующим образом:

    miniR = 
    
      GeographicCellsReference with properties:
    
                 LatitudeLimits: [-90 90]
                LongitudeLimits: [-180 180]
                     RasterSize: [3 6]
           RasterInterpretation: 'cells'
               ColumnsStartFrom: 'south'
                  RowsStartFrom: 'west'
           CellExtentInLatitude: 60
          CellExtentInLongitude: 60
         RasterExtentInLatitude: 180
        RasterExtentInLongitude: 360
               XIntrinsicLimits: [0.5 6.5]
               YIntrinsicLimits: [0.5 3.5]
           CoordinateSystemType: 'geographic'
                      AngleUnit: 'degree'
  3. Настройте эквидистантную цилиндрическую проекцию карты:

    figure('Color','white')
    ax = axesm('MapProjection', 'eqdcylin');
    axis off
    setm(ax,'GLineStyle','-', 'Grid','on','Frame','on')
    
  4. Нарисуйте решетку с параллельными и меридианальными метками с интервалами 60 °:

    setm(ax, 'MlabelLocation', 60, 'PlabelLocation',[-30 30],...
        'MLabelParallel','north', 'MeridianLabel','on',...
        'ParallelLabel','on','MlineLocation',60,...
        'PlineLocation',[-30 30])
    
  5. Отображение данных с помощью geoshow и отображение с цветным пандусом и легендой:

    geoshow(miniZ, miniR, 'DisplayType', 'texturemap');
    colormap('autumn')
    colorbar
    

Обратите внимание, что первая строка матрицы отображается внизу карты, а последняя - вверху.