调用arcpy包批量进行矢量掩膜提取

使用一个polygon矢量提取某个文件夹中所有的tif格式栅格数据

（要确保先安装好arcpy包）

import arcpy
arcpy.checkoutextension("spatial")
arcpy.gp.overwriteoutput=1
arcpy.env.workspace = "e:\\data"  #所有栅格影像所在文件夹
rasters = arcpy.listrasters("*", "tif")
mask= "e:\\data\\bound.shp"  #用于提取的矢量掩膜
for raster in rasters:
    print(raster)
    out= "e:\\data\\images\\"+raster
    arcpy.gp.extractbymask_sa(raster, mask, out)
print("ok")

2019-09-02 23:53:36