坡度坡向的提取算法.docVIP

a b c d e f g h i 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/dy], -[dz/dx]) 应为aspect = 57.29578 * atan2 (-[dz/dx],[dz/dy]) atan2（x_num,y_num）根据制定的X轴及Y轴坐标值，返回反正切值,(-pi,pi]。 if aspect 0 cell = 90.0 - aspect else if aspect 90.0 cell = 360.0 - aspect + 90.0 else cell = 90.0 - aspect 1.加载ArcTutorSpatial文件夹中的elevation2.利用Spatial AnalystSurface AnalysisSlope计算elevation数据的坡度，为避免求反正切函数，可选择Percent； 3.利用Spatial AnalystSurface AnalysisAspect计算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.利用InterpreterTopographic AnalysisSlope计算elevation1数据的坡度，为避免求反正切函数，可选择Percent； 3.利用InterpreterTopographic AnalysisAspect计算elevation数据的坡向； 4.利用绘图工具在数据视图中画一包含3*3个像元的窗口，将其高程数据依次输入Excel中； 5.按照坡度坡向的求取公式求坡度坡向； 6.验证Erdas的坡度坡向求取算法。