方法总结类：栅格数据批量裁剪--Python脚本

原创于 2025-10-15 22:22:51 发布 · 138 阅读

5 ·

CC 4.0 BY-SA版权

文章标签：

#python

这个脚本是2015年处理MODIS影像时，在网上查了很多类似帖子后，总结的一个技术贴。该贴的首发在气象家园上。

下面这个就是栅格数据批量处理的脚本：

import arcpy
import os
from arcpy import env
from arcpy.sa import *
arcpy.CheckOutExtension("Spatial")
env.workspace = r'K:\\modissnow\\MOD10A2\\exracttif'
outMaskDate = 'K:\\modisrslt\\result\\arcgistemp.gdb'
for fc in arcpy.ListRasters():
    output = os.path.join(outMaskDate,fc.split('.')[0])
    for mask in arcpy.ListFeatureClasses():        
            outMask=ExtractByMask(fc,mask)
            outMask.save("K:\\modisrslt\\result\\"+fc)

上面这个是针对大批量的栅格数据裁剪的，这样一个是速度快，还可以避免手动裁剪过程中的错误。
注：把矢量边界和待裁剪的栅格数据放到一个文件下。

效果图：

裁剪前：

裁剪后：

另外，如果是需要处理的影像少的话，就可以用arcgis中自带的Batch工具进行裁剪，可以参考这个博文（http://blog.sina.com.cn/s/blog_67bd825f0101d6qp.html），讲的很详细。

在这个过程中，参考了好几篇博文(下面的链接)，感谢这些博主的无私奉献，再次表示感谢。
1.http://blog.youkuaiyun.com/lucky51222/article/details/39100109
2.http://www.gisall.com/forum.php?mod=viewthread&tid=132112&extra=&page=1

另外，帖子中关于Python的内容，若有不足之处，还请大家多多指导。

附件是栅格数据批量裁剪的脚本。