Проектирование и отображение растровых данных

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

После преобразования растровых данных можно отобразить его на карте с помощью функций визуализации, таких как mapshow и geoshow. Использовать mapshow для проективных координат x-y и geoshow для нераскрытых координат широта-долгота.

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

Открыть Live Script

Чтобы спроектировать данные, которые связаны с географическим объектом растровой привязки, сначала создайте сетку координат широта-долгота для каждой точки растра. Затем проецируйте географические координаты на координаты x-y карты.

Например, импортируйте растровые данные о высоте как массив и объект-ссылку на географические ячейки. Получите координаты широта-долгота для каждой точки растра при помощи geographicGrid функция.

[Z,R] = readgeoraster('n39_w106_3arc_v2.dt1');
[lat,lon] = geographicGrid(R);

Теперь, когда у вас есть сетка, выберите проекцию карты, которая будет использоваться при проекции координат. В данном примере создайте projcrs объект для зоны 13 UTM в северном полушарии. Затем проецируйте координаты широта-долгота на координаты x-y.

p = projcrs(32613);
[x,y] = projfwd(p,lat,lon);

Отобразите проективный растр как поверхность путем вызова mapshow и определение координат x-y и массива повышений. Добавьте подписей по осям и примените палитру, соответствующую данным о повышении.

figure
mapshow(x,y,Z,'DisplayType','surface')
xlabel('x (meters)')
ylabel('y (meters)')
demcmap(Z)

Figure contains an axes. The axes contains an object of type surface.

Если географическая CRS координат широта-долгота не совпадает с географической CRS проектируемой CRS, то проектируемые координаты могут быть неточными. Вы можете найти географическую CRS projcrs объект или географический объект растровой привязки путем запроса их GeographicCRS свойства.

p.GeographicCRS.Name

ans = 
"WGS 84"

R.GeographicCRS.Name

ans = 
"WGS 84"

Файл DTED, используемый в этом примере, предоставлен Геологической службой США.

Депроектация растровых данных

Открыть Live Script

Чтобы отключить данные, которые связаны с объектом растровой привязки карты, сначала создайте сетку координат x-y для каждой точки растра. Затем постройте координаты карты x-y до географических координат.

Например, импортируйте изображение Бостона как массив и ссылку на объект камер карты. Получите информацию о проекции карты как о projcrs объект путем запроса ProjectedCRS свойство объекта ссылки.

[Z,R] = readgeoraster('boston.tif');
p = R.ProjectedCRS;

Получите координаты x-y для каждой точки растра при помощи worldGrid функция.

[x,y] = worldGrid(R);

Отключите координаты x-y в координаты широта-долгота при помощи projinv функция и установка projcrs объект и координатная сетка.

[lat,lon] = projinv(p,x,y);

Отображение неуправляемого изображения вызовом geoshow и определение координат широты-долготы и массива изображений. По умолчанию geoshow отображает координаты с помощью функции Plate Carr $\overset{´}{e}$ е проекция. Затем добавьте подписей по осям.

figure
geoshow(lat,lon,Z)
xlabel('Longitude (degrees)')
ylabel('Latitude (degrees)')

Документация