arcgis中坡向计算工作原理说明

用于识别出从每个像元到其相邻像元方向上值的变化率最大的下坡方向。坡向可以被视为坡度方向。输出栅格中各像元的值可指示出各像元位置处表面的朝向的罗盘方向。将按照顺时针方向进行测量,角度范围介于 0(正北)到 360(仍是正北)之间,即完整的圆。不具有下坡方向的平坦区域将赋值为 -1。

坡向数据集中每个像元的值都可指示出该像元的坡度朝向。

从概念上讲,坡向工具将根据要处理的像元或中心像元周围一个 3 x 3 的像元邻域的 z 值拟合出一个平面。该平面的朝向就是要处理的像元的坡向。

下图显示的是输入高程数据集和输出坡向栅格。

1. 坡向算法

移动的 3 x 3 窗口会访问输入栅格中的每个像元,而每次位于窗口中心的像元的坡向值将通过一种将纳入八个相邻像元值的算法进行计算。这些像元使用字母 ai 进行标识,其中 e 表示当前正在计算坡向的像元。

像元 e 在 x 方向上的变化率将通过以下算法进行计算:

复制代码
[dz/dx] = ((c + 2f + i) - (a + 2d + g)) / 8

像元 e 在 y 方向上的变化率将通过以下算法进行计算:

复制代码
[dz/dy] = ((g + 2h + i) - (a + 2b + c)) / 8

代入像元 e 在 x 方向和 y 方向上的变化率,坡向将通过以下算法进行计算:

复制代码
aspect = 57.29578 * atan2 ([dz/dy], -[dz/dx])

然后,坡向值将根据以下规则转换为罗盘方向值(0 到 360 度):

2. 坡向计算示例

示例中,将计算移动窗口内中心像元的坡向值。

中心像元 e 在 x 方向上的变化率为:

dz/dx\] =((*c* + 2*f* + *i* )-(*a* + 2*d* + *g*))/8=((85 + 170 + 84))-(101 + 202 + 101))/8=-8.125 像元 *e* 在 y 方向上的变化率为: \[dz/dy\] = ((*g* + 2*h* + *i* )-(*a* + 2*b* + *c*) / 8 =((101 + 182 + 84)-(101 + 184 + 85))/ 8 = -0.375 坡向计算如下: a*spect* = 57.29578 \* atan2 (\[dz/dy\], -\[dz/dx\])= 57.29578 \* atan2 (-0.375, 8.125) = -2.64 由于计算得出的值小于零,则根据最终规则得出: c*ell* = 90.0 - *aspect* = 90 - (-2.64)= 90 + 2.64 = 92.64 中心像元 *e* 的值 92.64 表明它的坡向为朝东。 ![](https://file.jishuzhan.net/article/1782115455889051650/7023949beab525331a63405f88afb214.webp) **3. 参考文献** Burrough, P. A., and McDonell, R. A., 1998. Principles of Geographical Information Systems (Oxford University Press, New York), 190 pp.

相关推荐
嘘嘘出差1 小时前
ArcGIS切片方案记录bundle文件
arcgis
摆烂老大1 天前
SWAT| 水文 | SWAT模型(三):土壤数据库制备
arcgis·水文·swat模型
鸿业远图科技2 天前
分式注记种表达方式arcgis
python·arcgis
城市数据匠2 天前
31【干货】Arcgis属性表常用查询表达式实战大全
arcgis·gis·cad·国土空间规划·自然资源局
新中地GIS开发老师2 天前
【Cesium入门教程】第七课:Primitive图元
arcgis·信息可视化·gis开发·webgis·地理信息系统·地理信息科学
王孝点3 天前
geoserver发布arcgis瓦片地图服务(最新版本)
arcgis
没有梦想的咸鱼185-1037-16635 天前
【大语言模型ChatGPT4/4o 】“AI大模型+”多技术融合:赋能自然科学暨ChatGPT在地学、GIS、气象、农业、生态与环境领域中的应用
人工智能·python·机器学习·arcgis·语言模型·chatgpt·数据分析
新中地GIS开发老师5 天前
【Cesium入门教程】第五课:数据源
arcgis·遥感·gis开发·webgis·地理信息科学
没有梦想的咸鱼185-1037-16636 天前
【大模型ChatGPT+ArcGIS】数据处理、空间分析、可视化及多案例综合应用
人工智能·arcgis·chatgpt·数据分析
玩大数据的龙威7 天前
【ArcGIS技巧】用地块生成界址点去重、顺时针编号挂接DKBM属性
arcgis