Python批量裁剪MODIS栅格影像数据——利用shp文件裁剪tif影像

任务描述： 利用同一个shp文件裁剪多幅栅格影像数据
数据准备： 文件夹(工作空间)内含有需要批量裁剪的tif影像
工具准备： Python arcpy 环境
*** 作： 替换代码中的路径即可

# -*- coding: utf-8 -*-
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "D:/Year_MEAN"  # 待裁剪tif影像存放路径
outputpath = "D:/Year_MEAN"  # 输出结果影像存放路径

# Set local variables
rasters = arcpy.ListRasters("*", "tif")
inMaskData = "D:/data/YellowRiver_nineprovince.shp"  # shp文件存放路径及名称

for raster in rasters:
    # 去掉后缀和文件名中的'.'
    i = 0
    file_name = raster.split('.')[i]
    i = i + 1
    print(file_name)

    # Check out the ArcGIS Spatial Analyst extension license
    arcpy.CheckOutExtension("Spatial")

    # Execute ExtractByMask
    outExtractByMask = ExtractByMask(raster, inMaskData)

    # Save the output
    outExtractByMask.save(outputpath + '/' + file_name + '_clip.tif')

print('finished!')

此处mark一个可以多级目录批量裁剪的方法 留作备用：ArcGIS批量裁剪栅格数据（ArcPy方法）