gradientm

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

Описание

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

пример

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

[aspect,slope,gradN,gradE] = gradientm(F,R,spheroid) использует заданный опорный сфероид вместо сфероида, содержавшегося в Spheroid свойство geocrs объект в GeographicCRS свойство R или сфероид 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'.

Примеры

свернуть все

Сгенерируйте демонстрационные данные о вертикальном изменении для области вокруг горной вершины с помощью географического объекта ссылки регистраций и 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 contains an axes object. The axes object with title Elevation Data contains an object of type surface.

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

figure
axesm('eqdcylin')
geoshow(gradN,R,'DisplayType','surface')
title('North Components of Gradient')
colorbar

Figure contains an axes object. The axes object with title North Components of Gradient contains an object of type surface.

figure
axesm('eqdcylin')
geoshow(gradE,R,'DisplayType','surface')
title('East Components of Gradient')
colorbar

Figure contains an axes object. The axes object with title East Components of Gradient contains an object of type surface.

Постройте углы наклона. Обратите внимание на то, что значение угла наклона является нулем на саммите.

figure
axesm('eqdcylin')
geoshow(slope,R,'DisplayType','surface')
title('Slope Angles')
colorbar

Figure contains an axes object. The axes object with title Slope Angles contains an object of type surface.

Постройте углы обзора. Угол обзора описывает направление горные поверхности наклона как азимут, измеренный по часовой стрелке от севера.

figure
axesm('eqdcylin')
geoshow(aspect,R,'DisplayType','surface')
title('Aspect Angles')
colorbar

Figure contains an axes object. The axes object with title Aspect Angles contains an object of type surface.

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

свернуть все

Сетка данных в виде числовой матрицы по крайней мере с двумя строками и двумя столбцами. Сетка данных может содержать 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
Для просмотра документации необходимо авторизоваться на сайте