需求

现有一些tif格式的东亚范围MODIS EVI数据，希望将其裁剪为中国范围。

简介

掩膜裁剪（提取分析/Extract by Mask）

arcpy.gp.ExtractByMask_sa (in_raster, in_mask_data,out)

• in_raster:提取像元的输入栅格
• in_mask_data:用于定义提取区域的输入掩膜数据。可以是栅格，也可以是要素数据集。
• out:输出栅格

代码

import arcpy
from arcpy import env
env.workspace = "I:\\MODIS\\EVI\\2001_2017"
inMaskData = "E:\\GPP\\DEM\\DEM_1km\\china_dem1km";
rasterlist = arcpy.ListRasters("*", "tif")   
for raster in rasterlist:
   print str(raster)
   out = "I:\\MODIS\\EVI\\cj_2006_17\\"+raster
   arcpy.gp.ExtractByMask_sa(raster, inMaskData,out)
print("OK!")