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

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

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

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

Скрипт Open Live Script

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

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

[Z,R] = readgeoraster('n39_w106_3arc_v2.dt1');
[lat,lon] = meshgrat(Z,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)

Если географический 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 для каждой точки в растре. Для этого сначала создайте сетку внутренних координат, которая совпадает с размером растра при помощи meshgrid функция. Внутренние y-координаты увеличиваются от строки до строки, и внутренние x-координаты увеличиваются от столбца до столбца. Затем преобразуйте внутренние координаты к координатам x-y при помощи intrinsicToWorld функция.

s = size(Z);
[xIntrinsic,yIntrinsic] = meshgrid(1:s(2),1:s(1));
[x,y] = intrinsicToWorld(R,xIntrinsic,yIntrinsic);

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

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

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

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

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