数据存放在E盘,[data]文件夹,其中[mask.shp]为矢量边界。输出结果存放在E盘[data_extract] 文件夹中。
在ArcGIS中实现。
import arcpy arcpy.CheckOutExtension("spatial") arcpy.gp.overwriteOutput=1 arcpy.env.workspace = "E:\\data" rasters = arcpy.ListRasters("*", "tif") mask= "E:\\data\\mask.shp" for raster in rasters: print(raster) out= "E:\\data_extract\\"+"_"+raster[0:18] arcpy.gp.ExtractByMask_sa(raster, mask, out) print("OK")