Вычислить профиль отметки вдоль прямой линии

Вычислите профиль отметки вдоль прямой линии с помощью команды mapprofile функция. Расчет значений данных вдоль траектории является обычным применением при работе с геоданными с сеткой. Например, можно вычислить высоту рельефа местности вдоль участка, дороги или траектории полета. mapprofile функция делает это, основываясь на численных данных, определяющих набор ППМ, или определяя их в интерактивном режиме посредством графического ввода с дисплея карты. Значения, вычисленные для результирующего профиля, могут быть отображены на новом графике или возвращены в качестве выходных аргументов для дальнейшего анализа или отображения.

Загрузить данные о высоте и географические ячейки опорного объекта для Корейского полуострова.

load korea5c

Получение границ широты и долготы для данных отметки и настройка карты мира. Отобразите карту и примените карту цветов, соответствующую данным отметки.

latlim = korea5cR.LatitudeLimits;
lonlim = korea5cR.LongitudeLimits;
worldmap(latlim, lonlim)
meshm(korea5c,korea5cR,size(korea5c),korea5c)
demcmap(korea5c)

Определите конечные точки для прямолинейного перехода через область. Затем вычислите профиль отметки с помощью mapprofile. По умолчанию mapprofile использует билинейную интерполяцию вдоль большой дорожки круга.

plat = [40.5 30.7];
plon = [121.5 133.5];
[z,rng,lat,lon] = mapprofile(korea5c,korea5cR,plat,plon);

Постройте график трансекта в 3-D, чтобы он следовал за местностью.

plot3m(lat,lon,z,'w','LineWidth',2)

Постройте график диапазона пересечений и отметки на наборе декартовых осей.

figure
plot(rng,z,'r')