主要内容

gradientm

数据网格的梯度、斜率和方向

描述

常规数据网格

例子

方面gradN年级) = gradientm (FR返回常规数据网格的方向角、倾斜角以及梯度的北和东组件F关于地理参考R.默认情况下,gradientm定位引用的经纬度坐标R使用包含在球体财产的geocrs对象GeographicCRS的属性R.如果GeographicCRS的属性R是空的,然后geopeaks使用GRS80。

方面gradN年级) = gradientm (FR球体使用指定的引用球体而不是包含在球体财产的geocrs对象GeographicCRS的属性R或椭球体GRS80。

定位数据网格

方面gradN年级) = gradientm (纬度F为已定位的数据网格返回相同的值F定义的经度-纬度网格纬度.默认情况下,纬度和经度以度数表示。默认的参考球体是GRS80。

方面gradN年级) = gradientm (纬度F球体使用指定的参考球体而不是GRS80。

方面gradN年级) = gradientm (纬度F球体angleUnit指定纬度和经度的单位为“度”(默认)或“弧度”

例子

全部折叠

使用地理贴子参考对象和geopeaks函数。为此,首先通过指定纬度和经度限制以及高程数据网格的大小,为该区域创建一个引用对象。接下来,为使用的区域生成高程数据geopeaks

Latlim = [10 45];Lonlim = [60 100];Size = [100 100];R = georefpostings (latlim、lonlim、大小);F = geopeaks (R);

计算数据的纵横角、斜角和梯度分量。

[年级方面,斜坡,gradN] = gradientm (F, R);

通过绘制数据来可视化结果。首先,用等距圆柱投影绘制高程数据。为此,创建一组地图轴并指定使用的投影axesm.绘制数据作为表面使用geoshow,调整轴的长宽比使用daspect.显示轴的三维视图使用视图

图axesm (“eqdcylin”);geoshow (F R“DisplayType”“表面”([1 1 5])高程数据的) colorbar视图(3)

图中包含一个轴对象。标题为Elevation Data的axes对象包含一个类型为surface的对象。

然后,用相同的投影绘制梯度分量。注意,在顶点处,北分量和东分量的值都是零。

图axesm (“eqdcylin”) geoshow (gradN, R,“DisplayType”“表面”)标题(“梯度的北分量”) colorbar

图中包含一个轴对象。标题为North Components of Gradient的轴对象包含一个类型为surface的对象。

图axesm (“eqdcylin”) geoshow(年级、R,“DisplayType”“表面”)标题(“梯度的东分量”) colorbar

图中包含一个轴对象。标题为East Components of Gradient的轴对象包含一个类型为surface的对象。

画出斜率角。注意,坡角在峰顶处的值为零。

图axesm (“eqdcylin”) geoshow(斜率,R,“DisplayType”“表面”)标题(“坡角”) colorbar

图中包含一个轴对象。标题为“斜角”的轴对象包含一个类型为surface的对象。

画出侧面角。坡向角用从北顺时针方向测量的方位角来描述山坡的朝向。

图axesm (“eqdcylin”) geoshow(方面,R,“DisplayType”“表面”)标题(的角度方面) colorbar

图中包含一个轴对象。标题为Aspect Angles的axes对象包含一个类型为surface的对象。

输入参数

全部折叠

数据网格,指定为至少有两行两列的数字矩阵。数据网格可能包含值。F是与地理栅格参考对象相关联的常规数据网格,或与经纬度网格相关联的地理引用数据网格。

如果F是常规的数据网格和R是引用对象吗大小(F)一定和R.RasterSize.如果F是一个地理定位的数据网格吗大小(F)一定和大小(lat)大小(朗)

数据类型:|

包含的地理空间引用信息的地理引用F,指定为以下值之一:

  • GeographicCellsReferenceGeographicPostingsReference对象,R.RasterSize是一样的大小(F)

  • 将数据网格的行和列索引与地理坐标关联起来的3 × 2数字矩阵[color = # 0000ff] [color = # 0000ff] [color = # 0000ff.矩阵必须定义一个非旋转和非倾斜的关系,其中数据网格的每一列沿着子午线,每一行沿着平行线。

  • 1乘3的带有元素的数字向量[s nlat wlon],在那里年代为每度的数据网格样本数,nlat数据网格的最北端纬度是否以度数和wlon是最西端的经度。

有关引用向量和矩阵的更多信息,请参见地理栅格数据

参考球体,指定为referenceEllipsoid对象,oblateSpheroid对象,referenceSphere对象,或形式向量(semimajorAxis偏心)

有关参考球体的更多信息,请参阅参考球状体

例子:球体= referenceEllipsoid (wgs84Ellipsoid);

纬度,指定为至少有两行两列的数字矩阵。默认情况下,以度数指定纬度。要使用以弧度为单位的值,请指定angleUnit参数为“弧度”

纬度必须和F

数据类型:|

经度,指定为至少有两行两列的数字矩阵。默认情况下,以度数指定经度。要使用以弧度为单位的值,请指定angleUnit参数为“弧度”

必须和纬度F

数据类型:|

角度单位,指定为“度”(默认)或“弧度”

输出参数

全部折叠

方面角度,返回的矩阵大小相同F.纵横角是指F递减最快,用从北顺时针方向测量的方位角表示。

在默认情况下,方面的角度是角度。若要返回以弧度为单位的值,请指定纬度弧度和angleUnit作为“弧度”

如果梯度的两个分量都为零,那么纵横角度返回为

数据类型:

斜率角度,以相同大小的矩阵返回F.为了使坡度角具有物理意义,数据网格必须指定高程,其距离单位必须与参考椭球的长度单位相匹配。否则,斜率角就是梯度大小的反正切。

默认情况下,斜率以度数表示。若要返回以弧度为单位的值,请指定纬度弧度和angleUnit作为“弧度”

数据类型:

梯度的北分量,返回为相同大小的矩阵F.梯度的北分量是变化量R在北方向上的每一单位距离,其中距离单位与参考球体的长度单位相匹配。

数据类型:

梯度的东分量,返回为相同大小的矩阵F.梯度的东分量是变化量R在东侧的距离单位,其中距离单位与参考椭球的长度单位相匹配。

数据类型:

之前介绍过的R2006a