需求
现有一些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!")
网友评论