# -*- coding:utf-8 -*-
... import arcpy
... from arcpy import env #定义env
... from arcpy.sa import *
...
... arcpy.env.workspace = "H:/遥感数据/NDVI/" #定义工作空间及数据路径
... rasters = arcpy.ListRasters("*", "tif") #遍历工作空间中的tif格式数据
... inMask = "H:/遥感数据/NDVI/1.shp"#裁剪文件及路径,裁剪文件可以是tif文件
... for raster in rasters:
... outpath = "H:/遥感数据/OUT/" +str(raster) #输出存储路径以及输出文件名,这里使用的是读取到的数据名称
... outCJ = ExtractByMask(raster, inMask)
... outCJ.save(outpath) #输出存储裁剪的栅格数据
... print(str(raster)) #输出读取并裁剪的栅格数据名称
... print("over!!!!!!!!")
只需要修改输入空间路径,输出空间路径,以及裁剪的shp或者tif文件路径
打开方式打开【Arcmap】-【地理处理】-【Python】
使用方法将代码复制到窗口内,双击回车开始运行。