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

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

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

Растровые данные о проекте

Скрипт Open Live Script

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

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

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

Теперь, когда у вас есть своя сетка, выберите проекцию карты, чтобы использовать при проектировании координат. В данном примере создайте projcrs объект для зоны UTM 13 в северном полушарии. Затем спроектируйте координаты долготы широты к координатам 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 object. The axes object 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, используемый в этом примере, является любезностью Геологической службы США.

Растровые данные о непроекте

Скрипт Open 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 координаты отображений с помощью Пластины проекция Carrée. Затем добавьте подписи по осям.

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

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