arcpy/python循环批量裁剪tiff栅格文件

**

1. arcgis利用一个矢量批量裁剪多个栅格，可参考

**批量裁剪栅格文件

2.arcgis利用多个矢量要裁剪多个栅格文件，裁剪后自动生成文件夹保存。
例如：有全国每省的矢量和全国2000-2020年的植被覆盖图，现在要分省统计进行分析。不同省shp文件可以裁剪全国影像后，裁剪好的影像在路径中用省市名字保存。

# -*- coding: UTF-8 -*-
import arcpy
import os
import glob
import arcpy
from arcpy.sa import *

arcpy.CheckOutExtension("ImageAnalyst")  # 检查许可
arcpy.CheckOutExtension("spatial")

from arcpy.sa import *
arcpy.env.workspace = "D:\NMNEW\soil\soil_TG"
#定义工作空间及数据路径，路径下有不同区域矢量文件和被裁剪的全部TIFF文件；例如河北省、山东省各省shp,全国多年份 tiff文件（2000-2020）
rasters = arcpy.ListRasters("*", "tif") #遍历工作空间中的tif格式数据
inMasks = arcpy.ListFeatureClasses()#遍历工作空间中的shp格式数据
for inMask in inMasks:
  for raster in rasters:
    outpath = r"D:\NMNEW\soil\TG"
    if not os.path.exists(outpath + os.sep + str(inMask).replace('.','_')):#如果储存路径下没有以矢量文件命名的文件夹
      os.mkdir(outpath + os.sep + str(inMask).replace('.','_'))#生成以矢量文件命名的文件夹（注意.shp中的点要替换）
    outCJ = ExtractByMask(raster, inMask)#批量裁剪文件
    print (outpath + os.sep + str(inMask).replace('.','_')+os.sep+str(raster))
    outCJ.save(outpath + os.sep + str(inMask).replace('.','_') +os.sep+str(raster))  # 输出存储裁剪的栅格数据，存储到新建文件夹里
    print(str(raster))    #输出读取并裁剪的栅格数据名称
print("over!!!!!!!!")