Вычислите градиент, наклон и аспект сетки данных
[ASPECT, SLOPE, gradN, gradE] = gradientm(Z,
R)
[...] = gradientm(lat, lon, Z)
[...] = gradientm(..., ellipsoid)
[...] = gradientm(lat, lon, Z, ellipsoid, units)
[ASPECT, SLOPE, gradN, gradE] = gradientm(Z,
R) вычисляет наклон, аспект и северные и восточные компоненты градиента для обычной сетки данных Z относительно ссылочного R. Если сетка содержит вертикальные изменения в метрах, получившийся аспект и наклон находятся в модулях степеней по часовой стрелке от севера и от горизонтали. Северные и восточные компоненты градиента являются изменением в переменной карты на метр расстояния в северных и восточных направлениях. Расчет использует конечные разности в переменной карты на наземном эллипсоиде по умолчанию.
R может быть географический объект растровой привязки, вектор ссылки или матрица привязки.
Если R географический объект растровой привязки, его RasterSize свойство должно быть сопоставимо с size(Z).
Если R вектор ссылки, это должно быть 1 3 с элементами:
[cells/degree northern_latitude_limit western_longitude_limit]
Если R матрица привязки, это должно быть 3 2 и преобразовать растровые индексы строки и столбца к или от географических координат согласно:
[lon lat] = [row col 1] * R
Если R матрица привязки, она должна задать (невращательный, нескошенный) отношение, в котором каждый столбец сетки данных падает вдоль меридиана и каждой строки падения вдоль параллели. Интерполяция ближайшего соседа используется по умолчанию. NaN возвращен для точек вне пределов сетки или для который lat или lon содержите NaN. Все углы находятся в модулях степеней.
[...] = gradientm(lat, lon, Z) делает расчет для определенной геолокацию сетки данных. lat и lon, широты и долготы точек геолокации, в градусах.
[...] = gradientm(..., ellipsoid) использует ссылочный эллипсоид, заданный входом ellipsoid, который может быть referenceSphere, referenceEllipsoid, или oblateSpheroid объект или вектор формы [semimajor_axis eccentricity]. Если карта содержит вертикальные изменения в тех же единицах длины как большая полуось эллипсоида, наклон и аспект находятся в модулях степеней. Эта форма вызова является самой полезной для расчетов на телах кроме земли.
[...] = gradientm(lat, lon, Z, ellipsoid, units) задает угловые модули входных параметров широты и долготы. Значением по умолчанию является 'degrees'. Поскольку вертикальное изменение сопоставляет в тех же модулях как ellipsoid(1), получившийся наклон и аспект находятся в заданных модулях. Компоненты градиента являются изменением в переменной карты на единицу длины, с помощью той же единицы длины в качестве большой полуоси эллипсоида.
Вычислите и отобразите наклон для 30 арксекунд (10 км) данные о вертикальном изменении Кореи. Наклоны в Японском море являются до 12 градусов в этом разрешении сетки.
load korea [aspect,slope,gradN,gradE] = gradientm(map,refvec); worldmap(slope,refvec) geoshow(slope,refvec,'DisplayType','texturemap') cmap = cool(10); demcmap('inc',slope,1,[],cmap) colorbar latlim = getm(gca,'maplatlimit'); lonlim = getm(gca,'maplonlimit'); land = shaperead('landareas',... 'UseGeoCoords',true,'BoundingBox',[lonlim' latlim']); geoshow(land,'FaceColor','none')
set(gca,'Visible','off')
Крупные цифровые модели вертикального изменения могут значительно недооценить локальный наклон. Для предыдущей карты точки вертикального изменения разделяются приблизительно на 10 километров. Ландшафт между двумя смежными точками моделируется как линейное изменение, в то время как фактический ландшафт может варьироваться намного более резко по такому расстоянию.
