坡度坡向的提取算法

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.29578rise_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]。ifaspect0cell=90.0-aspectelseifaspect90.0cell=360.0-aspect+90.0elsecell=90.0-aspect1.加载ArcTutorSpatial文件夹中的elevation；2.利用SpatialAnalystSurfaceAnalysisSlope计算elevation数据的坡度，为避免求反正切函数，可选择Percent；3.利用SpatialAnalystSurfaceAnalysisAspect计算elevation数据的坡向；4.利用绘图工具在数据视图中画一包含3*3个像元的窗口，将其高程数据依次输入Excel中；5.按照坡度坡向的求取公式求坡度坡向；6.验证ArcGIS的坡度坡向求取算法。坡度（Percent）=rise_run*100；aspect=57.29578*atan2(-[dz/dx],[dz/dy])abcdefghiErdas提取坡度坡向的算法(向上为y轴正方向，向右为x轴正方向)三阶不带权差分坡度的计算slope(indegrees)=坡度百分比所以坡度百分比为1~200。坡向的提取Ifx=0andy=0thentheaspectisflat(codedto361degrees).Aspect*180/pi变为弧度1.将GRID数据转换为elevation1.img；2.利用InterpreterTopographicAnalysisSlope计算elevation1数据的坡度，为避免求反正切函数，可选择Percent；3.利用InterpreterTopographicAnalysisAspect计算elevation数据的坡向；4.利用绘图工具在数据视图中画一包含3*3个像元的窗口，将其高程数据依次输入Excel中；5.按照坡度坡向的求取公式求坡度坡向；6.验证Erdas的坡度坡向求取算法。