坡度坡向的提取算法

(向下为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 = A TAN ( √ ( [dz/dx]2 + [dz/dy]2 ) ) * 57.29578

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

if aspect < 0

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。

坡向的提取

If 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的坡度坡向求取算法。