以下代码,需要在arcgis系统的python编程窗口实现,需要替换栅格影像文件夹路径,矢量shp文件路径,及裁剪后的输出路径。
另外,需要注意,栅格格式如果不是tif,也需要修改为对应的栅格类型。
# -*- coding: utf-8 -*-
import arcpy
arcpy.CheckOutExtension("spatial")
arcpy.gp.overwriteOutput=0
arcpy.env.workspace = r'D:\1paper\1ETandP\dayquchunull' #被裁剪栅格影像所在文件夹
rasters = arcpy.ListRasters("*", 'tif') #栅格数据格式设置
mask= r'D:\1paper\1ETandP\liuyu\touy.shp' #裁剪范围shp文件路径
outpath=r'D:\1paper\1ETandP\dayquchunullcj' #影像裁剪结果输出文件夹
for raster in rasters:
print(raster)
out=outpath+'\\'+ raster.split(".")[1]+ ".tif" #重命名,将原文件以.为分割,提取哪一段为新名,0为.前的第1段字符
arcpy.gp.ExtractByMask_SA(raster, mask, out)
print(raster+" has done")
print("All done")