需求
将原始多幅影像批量生成坡度和坡向信息,原始多幅拼接DEM:
去看原文
过程
#坡度代码
import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.CheckOutExtension(“spatial”)
arcpy.env.workspace = r"C:\Users\Desktop\test\input_dems" #dem所在文件夹
rasters = arcpy.ListRasters("*", “tif”) #dem格式
for ras in rasters:
print(ras)
out_slope=Slope(ras)
out_slope.save(“C:\Users\Desktop\test\output_slope\”+ras ) #坡度输出的路径
print(“All done”)
得到批量生成的坡度图:
#坡向代码
import arcpy
from arcpy import env
from arcpy.sa import *
arcpy.CheckOutExtension(“spatial”)
arcpy.env.workspace = r"C:\Users\Desktop\test\input_dems" #dem所在文件夹
rasters = arcpy.ListRasters("*", “tif”) #dem格式
for ras in rasters:
print(ras)
out_slope=Aspect(ras)
out_slope.save(“C:\Users\Desktop\test\output_aspect\”+ras ) #坡向输出的路径
print(“All done”)
原文有源码,更多内容,请关注地学数据处理分析。