Exponenta Banner

Соответствуйте данным с координатной сеткой к координатной сетке

Объекты подложки проектов тулбокса способом, подобным традиционным методам создания карты. Картограф сначала размечает сетку меридианов и находит что-либо подобное, вызвал координатную сетку. Каждая ячейка координатной сетки является географическим четырехугольником. Картограф вычисляет или интерполирует соответствующие x-y местоположения для каждой вершины в сетке координатной сетки и чертит спроектированную координатную сетку путем соединения точек. Наконец, картограф чертит данные о карте, от руки, пытаясь составлять форму ячеек координатной сетки, которые обычно изменяют форму через карту. Точно так же тулбокс вычисляет x-y местоположения четырех вершин каждой ячейки координатной сетки и деформирует или выбирает матричные данные, чтобы соответствовать получившемуся четырехугольнику.

В отображении сеток данных с помощью тулбокса, как в традиционной картографии, чем более прекрасный mesh (аналогичный использованию координатной сетки с большим количеством меридианов и параллелей), тем большая точность спроектированное отображение карты будет иметь, за счет большего усилия и время. Координатная сетка в печатной карте походит на интервал элементов сетки в обычной сетке данных, представление Mapping Toolbox™ которой является двухэлементными векторами формы [ number-of-parallels, number-of-meridians ]. Координатная сетка для определенных геолокацию сеток данных подобна; это - размер матриц координаты широты и долготы: number-of-parallels = mrows-1 и number-of-meridians = ncols-1. Однако, потому что определенные геолокацию сетки данных имеют произвольные угловые местоположения ячейки, интервал может отличаться, и таким образом их координатная сетка не является регулярной mesh.

Соответствуйте данным с координатной сеткой к прекрасным и крупным координатным сеткам

Скрипт Open Live Script

Этот пример показывает, как соответствовать данным с координатной сеткой к прекрасным и крупным координатным сеткам. Выбором координатной сетки является баланс скорости по точности с точки зрения расположения сетки на карте. Как правило, нет никакого смысла к определению mesh, более прекрасной, чем разрешение данных (в этом примере, 180 360 ячейки сетки). На практике это целесообразно использовать крупные координатные сетки для задач разработки и прекрасные координатные сетки для производства конечного изображения.

Обратите внимание на то, что, независимо от разрешения координатной сетки, данные о сетке неизменны. В этом случае сетка данных 180 360 матрица topo, и независимо от того, где это расположено, значения данных неизменны.

Загрузите сетку данных. Это загружает несколько переменных в рабочую область включая пределы широты (topolatlim), и долгота ограничивает (topolonlim).

load topo

Создайте ссылочный объект для сетки данных topo. 

topoR = georefcells(topolatlim,topolonlim,size(topo))
topoR = 
  GeographicCellsReference with properties:

             LatitudeLimits: [-90 90]
            LongitudeLimits: [0 360]
                 RasterSize: [180 360]
       RasterInterpretation: 'cells'
           ColumnsStartFrom: 'south'
              RowsStartFrom: 'west'
       CellExtentInLatitude: 1
      CellExtentInLongitude: 1
     RasterExtentInLatitude: 180
    RasterExtentInLongitude: 360
           XIntrinsicLimits: [0.5 360.5]
           YIntrinsicLimits: [0.5 180.5]
       CoordinateSystemType: 'geographic'
                  AngleUnit: 'degree'

Настройте проекцию Робинсона, задайте крупное (10 20) координатная сетка ячейки и отобразите данные, сопоставленные с координатной сеткой с помощью палитры DEM. Заметьте, что для этой крупной координатной сетки, ребра карты не появляются как плавные кривые.

figure
axesm robinson
spacing = [10 20];
h = meshm(topo,topoR,spacing);
demcmap(topo)

Теперь сбросьте координатную сетку, с помощью функции setm, чтобы сделать его менее крупным, [50 100]. (Можно также сбросить координатную сетку с помощью функции meshgrat.) Замечают, что зубчатый эффект ребер теперь незначителен.

setm(h,'MeshGrat',[50 100])

Сбросьте координатную сетку снова, на этот раз к очень прекрасной сетке с помощью функции setm. Заметьте, что результат, кажется, не немного лучше, чем исходное отображение со значением по умолчанию [50 100] координатная сетка, но это взяло намного дольше, чтобы произвести. Создание более точной mesh является компромиссом разрешения по сравнению со временем и использованием памяти. 

setm(h,'MeshGrat',[200 400])

Документация Mapping Toolbox
Поддержка