坡度坡向的提取算法

abcdefghi))*57.29578ArcGIS提取坡度坡向的算法))*57.29578（向下为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(V([dz/dx】2+[dz/dy】2)slope_degrees=ATAN(V([dz/dx]2+[dz/dy]2)rise_run=V([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]。ifaspect<0cell=90.0-aspectelseifaspect>90.0cell=360.0-aspect+90.0elsecell=90.0-aspect加载ArcTutor>Spatial文件夹中的elevation；利用SpatialAnalyst>SurfaceAnalysis>Slope计算elevation数据的坡度，为避免求反正切函数，可选择Percent；利用SpatialAnalyst>SurfaceAnalysis>Aspect计算elevation数据的坡向；利用绘图工具在数据视图中画一包含3*3个像元的窗口，将其高程数据依次输入Excel中；按照坡度坡向的求取公式求坡度坡向；验证ArcGIS的坡度坡向求取算法。坡度(Percent)=rise_run*100；aspect=57.29578*atan2(-[dz/dx],[dz/dy])

Erdas提取坡度坡向的算法A.vj=c-lja7iErdas提取坡度坡向的算法A.vj=c-lja7i=a—gAa't=f-dAy2=h-gg二〔-gA^3=c-!坡度的计算蜃—禺+g+对$_〔期十4十轨）（向上为y轴正方向，向右为x轴正方向）三阶不带权差分abcd日JghI,-uA.ISO

taniAz.)x slope(indegrees)= 八坡度百分比IfA;<1.percentslope=[A：x100)IfA：>lhpercentslope=200-所以坡度百分比为1~200。坡向的提取(AA'i十十4JA.i= = —3(Ay}+A^+A^3)Ay= 才二 —Ifvx=0andy=0thentheaspectisflat(codedto361degrees).Aspect*180/pi变为弧度将GRID数据转换为elevation1.img；利用Interpreter>TopographicAnalysis>Slope计算elevation1数据的坡度，为避免求反正切函数，可选择Percent；利用Interpreter>TopographicAnalysis>Aspect计算elevation数据的坡向；4.利用绘图工具在数