坡度坡向的提取算法

ArcGIS 提取坡度坡向的算法

(向下为y 轴正方向，向右为x 轴正方向) 三阶反距离平方权差分

[dz/dx] = ((c + 2f + i) - (a + 2d + g) / (8 * x_cell_size)

[dz/dy] = ((g + 2h + i) - (a + 2b + c)) / (8 * y_cell_size)

slope_radians = ATAN ( √ ( [dz/dx]2 + [dz/dy]2 ) )

slope_degrees = ATAN ( √ ( [dz/dx]2 + [dz/dy]2 ) ) * 57.29578

rise_run = √ ( [dz/dx]2 + [dz/dy]2 ]

应为aspect = 57.29578 * atan2 (-[dz/dx],[dz/dy]) atan2（x_num,y_num）根据制定的X 轴及Y 轴坐标值，返回反正切值,(-pi,pi]。

if aspect

cell = 90.0 - aspect else if aspect > 90.0

cell = 360.0 - aspect + 90.0 else

cell = 90.0 - aspect

1. 加载ArcTutor>Spatial文件夹中的elevation ；

2. 利用Spatial Analyst>Surface Analysis>Slope计算elevation 数据的坡度，为避免求反正切函数，可选择Percent ；

3. 利用Spatial Analyst>Surface Analysis>Aspect计算elevation 数据的坡向； 4. 利用绘图工具在数据视图中画一包含3*3个像元的窗口，将其高程数据依次输入Excel 中；

5. 按照坡度坡向的求取公式求坡度坡向； 6. 验证ArcGIS 的坡度坡向求取算法。坡度（Percent ）=rise_run*100；

aspect = 57.29578 * atan2 (-[dz/dx],[dz/dy])

Erdas 提取坡度坡向的算法

(向上为y 轴正方向，向右为x 轴正方向) 三阶不带权差分

坡度的计算

slope (in degrees) = 坡度百分比

所以坡度百分比为1~200。

坡向的提取

x=0 and y=0 then the aspect is flat (coded to 361 degrees).

Aspect*180/pi 变为弧度

1. 将GRID 数据转换为elevation1.img ；

2. 利用Interpreter>Topographic Analysis>Slope计算elevation1数据的坡度，为避免求反正切函数，可选择Percent ；

3. 利用Interpreter>Topographic Analysis>Aspect计算elevation 数据的坡向； 4. 利用绘图工具在数据视图中画一包含3*3个像元的窗口，将其高程数据依次输入Excel 中；

5. 按照坡度坡向的求取公式求坡度坡向； 6. 验证Erdas 的坡度坡向求取算法。