这个示例演示了如何计算常规数据网格的梯度、斜率和方向。梯度分量是在北和东方向每米距离的网格变量的变化。坡度定义为从一个网格单元格到它的八个相邻单元格之一的最陡上升或下降路径上每单位距离的高程变化,用反正切表示。如果网格包含以米为单位的高程,那么朝向和坡度就是地面从北向向上顺时针方向的角度。的gradientm
函数使用有限差分方法计算常规数据网格或地理引用数据网格的梯度。该功能返回北和东方向(即,北至南,东至西)的梯度分量,以及坡度和坡向。默认角度的单位是度。
构建一个100 × 100的网格使用山峰
函数,并为其构造一个引用矩阵。
datagrid = 500 *峰(100);Gridrv = [1000 00];
生成包含方向、坡度、向北的坡度和向东的坡度的网格。
[年级方面,斜坡,gradN] = gradientm (datagrid, gridrv);
在圆柱形等面积投影中映射表面数据。从原始的高度开始。
axesmeqacylMeshm (datagrid,gridrv) colormap (jet(64)) colorbar(“绿色”)标题(“峰:海拔”)轴广场
清除框架并显示坡度网格。
clma meshm(斜率,gridrv) colorbar (“绿色”);标题(“峰:坡”)
映射方面网格。
clma meshm方面,gridrv colorbar (“绿色”);标题(“峰:方面”)
在地图上标出向北的坡度。
clma meshm (gradN gridrv) colorbar (“绿色”);标题(峰:北梯度的)
在地图上标出向东的渐变。
clma meshm(年级,gridrv) colorbar (“绿色”);标题(“峰:东梯度”)