第一个:【计算坡度】
读取原始栅格影像:
rast_file=r'F:\arcpy_data\imagery\钦州湾.tif'
... rast_obj=arcpy.Raster(rast_file)
... rast_obj
【计算坡度】:
slp=arcpy.sa.Slope(rast_obj)
基于ArcGIS的python shell实现的页面如下所示:
输出结果:
第二个:【计算坡向】
代码:
aspect=arcpy.sa.Aspect(rast_obj)
结果输出:
第三个:【计算洪水淹没】
代码:
flood=rast_obj<0.5
结果输出:
4、保存栅格数据:
代码:
slp.save(r'F:\arcpy_data\slope.tif')
>>> flood.save(r'F:\arcpy_data\flood.tif')
>>> aspect.save(r'F:\arcpy_data\aspect.tif')
操作页面:
输出结果: