使用arcpy的方式
import arcpy
from arcpy import env
#柵格文件路徑
env.workspace = "D:/研究數(shù)據/MODIS數(shù)據下載/論文數(shù)據集/LST2003_2016"
#獲取文件下所有tif格式文件
rasterList = arcpy.ListRasters("*","TIF")
#裁剪的shapefile文件
mask = "D:/hxl/HHHXJDY.shp"
#輸出路徑
output_path = "D:/LST_MAX_2003-2016/"
#循環(huán)讀取文件名
for raster in rasterList:
? ? print raster
? ? out = output_path + raster? #對生成文件進行命名
? ? arcpy.Clip_management(raster,"#",out,mask,"0","ClippingGeometry")
? ? print out+" clip has done!"