Exponenta Banner

Преобразуйте векторные данные в растровый формат

Можно преобразовать данные о векторе долготы широты в сетку в любом разрешении, вы принимаете решение сделать растровую базовую карту или слой сетки. Определенные инструменты Mapping Toolbox™ GUI помогают вам сделать часть этого, но можно также выполнить преобразования векторной графики в растровую из командной строки. Основной функцией для gridding векторных данных является vec2mtx, который выделяет строки сетке любого размера, на который вы указываете, отмечая строки 1 с и незанятые ячейки сетки с 0s. Сетка содержит, удваивается, но если вы хотите логическую сетку (см. Сетки Данных как Логические переменные), бросьте результат быть логическим массивом. Видеть пример, представление Creating Data Grids from Vector Data.

Если векторные данные состоят из полигонов (закрашенные фигуры), основы с координатной сеткой являются всей пустотой. Можно дифференцировать их использующий функцию encodem, вызывая его с массивом строк, столбцов, и отобрать значения, чтобы произвести новую сетку, содержащую многоугольные области, заполненные со значениями seed, чтобы заменить двоичные значения, сгенерированные vec2mtx. Видеть пример, представление Rasterize Polygons Interactively.

Создание сеток данных от векторных данных

Скрипт Open Live Script

Этот пример показывает, как преобразовать векторные данные в растровые данные с помощью функции vec2mtx. Пример использует данные о закрашенной фигуре для Индианы из файла форм usastatehi. Для получения дополнительной информации смотрите, Преобразовывают Векторные Данные в Растровый Формат.

Используйте shaperead, чтобы получить данные о закрашенной фигуре для контура.

indiana = shaperead('usastatehi.shp',...
    'UseGeoCoords', true,...
    'Selector', {@(name)strcmpi('Indiana',name), 'Name'});
inLat = indiana.Lat;
inLon = indiana.Lon;

Преобразуйте векторы в обычную сетку данных с помощью vec2mtx. Установите плотность сетки быть 40 ячейками на степень. Растеризируйте контур и сгенерируйте вектор ссылки для него.

gridDensity = 40;
[inGrid, inRefVec] = vec2mtx(inLat, inLon, gridDensity);
whos
  Name               Size              Bytes  Class     Attributes

  gridDensity        1x1                   8  double              
  inGrid           165x137            180840  double              
  inLat              1x626              5008  double              
  inLon              1x626              5008  double              
  inRefVec           1x3                  24  double              
  indiana            1x1               11324  struct

Сделайте карту сетки данных в контрастировании цветов.

figure
axesm eqdcyl
meshm(inGrid, inRefVec)
colormap jet(4)

Настройте пределы карты.

[latlim, lonlim] = limitm(inGrid, inRefVec);
setm(gca, 'Flatlimit', latlim, 'FlonLimit', lonlim)
tightmap

Задайте точку seed и отберите значение. Чтобы заполнить (перекод) внутреннюю часть Индианы, вам нужна точка seed (который должен быть идентифицирован строкой и столбцом), и значение seed (чтобы быть выделенным всем ячейкам в полигоне). Выберите среднюю строку и столбец сетки и выберите индексное значение 3, чтобы идентифицировать территорию при вызове encodem, чтобы сгенерировать новую сетку. Последний аргумент (1) идентифицирует код для граничных ячеек, где заполнение должно остановиться.

inPt = round([size(inGrid)/2, 3]);
inGrid3 = encodem(inGrid, inPt,1);

Очистите и перерисуйте карту с помощью заполненной сетки.

meshm(inGrid3, inRefVec)

Постройте исходные векторы на сетке, чтобы видеть, как хорошо данные были растеризированы. Можно использовать инструмент Zoom на окне рисунка, чтобы исследовать результаты gridding более тесно.

plotm(inLat, inLon,'k')

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