基于Python的批量裁剪遥感影像

需求

现有一些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!")
最后编辑于 ：
©著作权归作者所有,转载或内容合作请联系作者
平台声明：文章内容（如有图片或视频亦包括在内）由作者上传并发布，文章内容仅代表作者本人观点，简书系信息发布平台，仅提供信息存储服务。