Градиент, наклон и аспект сетки данных
[ возвращает угол обзора, угол наклона и северные и восточные компоненты градиента для обычной сетки данных aspect,slope,gradN,gradE] = gradientm(F,R)F относительно географического ссылочного R. По умолчанию, gradientm определяет местоположение координат широты и долготы, на которые ссылается R использование сфероида содержится в Spheroid свойство geocrs объект в GeographicCRS свойство R. Если GeographicCRS свойство R пусто, затем geopeaks использование GRS80.
Сгенерируйте демонстрационные данные о вертикальном изменении для области вокруг горной вершины с помощью географического объекта ссылки регистраций и geopeaks функция. Для этого сначала создайте ссылочный объект для области путем определения пределов широты и долготы и размера сетки данных о вертикальном изменении. Затем сгенерируйте данные о вертикальном изменении для области с помощью geopeaks.
latlim = [10 45]; lonlim = [60 100]; size = [100 100]; R = georefpostings(latlim,lonlim,size); F = geopeaks(R);
Вычислите углы обзора, углы наклона и компоненты градиента данных.
[aspect,slope,gradN,gradE] = gradientm(F,R);
Визуализируйте результаты путем отображения данных на графике. Во-первых, отобразите данные о вертикальном изменении на графике с помощью равноотстоящей цилиндрической проекции. Для этого создайте набор осей карты и задайте проекцию с помощью axesm. Отобразите данные на графике как поверхность с помощью geoshow, и настройте соотношение сторон осей с помощью daspect. Отобразите 3-D представление осей с помощью view.
figure axesm('eqdcylin'); geoshow(F,R,'DisplayType','surface') daspect([1 1 5]) title('Elevation Data') colorbar view(3)
Затем постройте компоненты градиента с помощью той же проекции. Обратите внимание на то, что и северные и восточные значения компонента являются нулем на саммите.
figure axesm('eqdcylin') geoshow(gradN,R,'DisplayType','surface') title('North Components of Gradient') colorbar
figure axesm('eqdcylin') geoshow(gradE,R,'DisplayType','surface') title('East Components of Gradient') colorbar
Постройте углы наклона. Обратите внимание на то, что значение угла наклона является нулем на саммите.
figure axesm('eqdcylin') geoshow(slope,R,'DisplayType','surface') title('Slope Angles') colorbar
Постройте углы обзора. Угол обзора описывает направление горные поверхности наклона как азимут, измеренный по часовой стрелке от севера.
figure axesm('eqdcylin') geoshow(aspect,R,'DisplayType','surface') title('Aspect Angles') colorbar
