python批量化裁剪遥感影像

import arcpy

import glob

import os

arcpy.CheckOutExtension('Spatial')

inws=r"E:\MODIS\MODIS_Daily_Cloudfree_2002-2015\MODIS_Daily_Cloudfree_2002-2015\2015"

outws=r"E:\MODIS\MODIS_Daily_Cloudfree_2002-2015\MODIS_Indus_2015"

mask=r"E:\UIB\UIB_Basin_Khan\UIB_Subbasins_Pro.shp"

rasters = glob.glob(os.path.join(inws, "*.tif"))

for ras in rasters:

outname = os.path.join(outws, os.path.basename(ras).split(".")[0] +

"_clp.tif")

out_extract = arcpy.sa.ExtractByMask(ras, mask)

out_extract.save(outname) ​​​​