Arcgis — 坡向批处理

• 理论：

  • 获取DEM影像上各个像素的坡向值，即每个像元到相邻像元数值变化率最大的下坡方向

  • 坡向，即坡度方向，最后输出的是0-359.9的数值，从正北方向顺时针方向算起：
    

  • 求解坡向的算法即将需要处理的像元3*3领域内的高程值拟合为一个平面，平面的朝向即为坡向值：
    

• 代码实现：

语法：Aspect_3d(in_raster, out_raster)，批处理则多了个循环读影像。

# -*- coding: cp936 -*-
def main():
    # Import system modules
    import arcpy
    arcpy.env.workspace="F:\\lj_0220\\nac_dem\\all_5m_dem" #待求坡向图的DEM文件夹路径
    rasters=arcpy.ListRasters("*","tif")
    for raster in rasters:        
        outAspect = "F:\\lj_0220\\nac_dem\\all_5m_aspect\\"+raster.strip(".tif")+"_"+"aspect.tif"
        arcpy.CheckOutExtension("3D")
        arcpy.Aspect_3d(raster, outAspect)
        print(raster+" has done!")
    print("All done")
main()

• 代码中对应的地方改成自己的文件夹路径即可。

• 参考链接：

arcmap坡向计算原理