gradientm

Градиент, наклон и аспект сетки данных

свернуть все на странице

Синтаксис

[aspect,slope,gradN,gradE] = gradientm(F,R)
[aspect,slope,gradN,gradE] = gradientm(F,R,spheroid)
[aspect,slope,gradN,gradE] = gradientm(lat,lon,F)
[aspect,slope,gradN,gradE] = gradientm(lat,lon,F,spheroid)
[aspect,slope,gradN,gradE] = gradientm(lat,lon,F,spheroid,angleUnit)

Описание

Обычные сетки данных

пример

[aspect,slope,gradN,gradE] = gradientm(F,R) возвращает угол обзора, угол наклона и северные и восточные компоненты градиента для обычной сетки данных F относительно географического ссылочного R. По умолчанию, gradientm на использование GRS80 как ссылочный сфероид для значений широты и долготы ссылается R.

[aspect,slope,gradN,gradE] = gradientm(F,R,spheroid) использует заданный ссылочный сфероид вместо GRS80.

Определенные геолокацию сетки данных

[aspect,slope,gradN,gradE] = gradientm(lat,lon,F) возвращает те же значения для определенной геолокацию сетки данных F относительно mesh долготы широты, заданной lat и lon. По умолчанию широта и долгота в градусах. Ссылочный сфероид по умолчанию является GRS80.

[aspect,slope,gradN,gradE] = gradientm(lat,lon,F,spheroid) использует заданный ссылочный сфероид вместо GRS80.

[aspect,slope,gradN,gradE] = gradientm(lat,lon,F,spheroid,angleUnit) задает модули для широты и долготы как 'degrees' (значение по умолчанию) или 'radians'.

Примеры

свернуть все

Скрипт Open Live Script

Сгенерируйте демонстрационные данные о вертикальном изменении для области вокруг горной вершины с помощью географического объекта ссылки регистраций и 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

Входные параметры

свернуть все

Сетка данных в виде числовой матрицы по крайней мере с двумя строками и двумя столбцами. Сетка данных может содержать NaN значения. F или обычная сетка данных, сопоставленная с географическим объектом растровой привязки или сетка геосправочных данных относительно mesh долготы широты.

Если F обычная сетка данных и R ссылочный объект, затем size(F) должен совпасть с R.RasterSize. Если F определенная геолокацию сетка данных, затем size(F) должен совпасть с size(lat) и size(lon).

Типы данных: single | double

Географическая ссылка, которая содержит геопространственную информацию о ссылке для FВ виде одного из этих значений:

  • GeographicCellsReference или GeographicPostingsReference объект, где R.RasterSize совпадает с size(F).

  • 3 2 числовая матрица, которая сопоставляет индексы строки и столбца сетки данных с географическими координатами, такими что [lon lat] = [row col 1] * R. Матрица должна задать невращательное и нескошенное отношение, в котором каждый столбец сетки данных падает вдоль меридиана и каждой строки падения вдоль параллели.

  • 1 3 числовой вектор с элементами [s nlat wlon], где s количество выборок сетки данных на степень, nlat самая северная широта сетки данных в градусах и wlon самая западная долгота в градусах.

Для получения дополнительной информации о ссылке на векторы и матрицы, смотрите Растровые Данные, на Которые геоссылаются.

Ссылочный сфероид в виде referenceEllipsoid объект, oblateSpheroid объект, referenceSphere объект или вектор формы [semimajorAxis eccentricity].

Для получения дополнительной информации о ссылочных сфероидах, смотрите Ссылочные Сфероиды.

Пример: spheroid = referenceEllipsoid(wgs84Ellipsoid);

Широты в виде числовой матрицы по крайней мере с двумя строками и двумя столбцами. По умолчанию задайте широты в градусах. Чтобы использовать значения в радианах, задайте angleUnit аргумент как 'radians'.

lat должен быть одного размера с lon и F.

Типы данных: single | double

Долготы в виде числовой матрицы по крайней мере с двумя строками и двумя столбцами. По умолчанию задайте долготы в градусах. Чтобы использовать значения в радианах, задайте angleUnit аргумент как 'radians'.

lon должен быть одного размера с lat и F.

Типы данных: single | double

Угловые модули в виде 'degrees' (значение по умолчанию) или 'radians'.

Выходные аргументы

свернуть все

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

По умолчанию углы обзора в градусах. К возвращаемым значениям в радианах задайте lat и lon в радианах и angleUnit как 'radians'.

Если оба компонента градиента являются нулем, то угол обзора возвращен как NaN.

Типы данных: double

Углы наклона, возвращенные как матрица одного размера с F. Для углов наклона, чтобы иметь физический смысл, сетка данных должна задать вертикальное изменение, и его единица расстояния должна совпадать с единицей длины ссылочного сфероида. В противном случае угол наклона является арктангенсом величины градиента.

По умолчанию углы наклона в градусах. К возвращаемым значениям в радианах задайте lat и lon в радианах и angleUnit как 'radians'.

Типы данных: double

Северные компоненты градиента, возвращенного как матрица одного размера с F. Северный компонент градиента является изменением в R на единицу расстояния в северном направлении, где единица расстояния совпадает с единицей длины ссылочного сфероида.

Типы данных: double

Восточные компоненты градиента, возвращенного как матрица одного размера с F. Восточный компонент градиента является изменением в R на единицу расстояния в восточном направлении, где единица расстояния совпадает с единицей длины ссылочного сфероида.

Типы данных: double

Смотрите также

Функции

Объекты

Темы

Представлено до R2006a

Документация Mapping Toolbox
Поддержка
Для просмотра документации необходимо авторизоваться на сайте