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

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

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

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')