使用Python脚本批量裁切栅格

对栅格的裁切，我们通常使用裁切（数据管理-栅格-栅格处理）或按掩膜提取（空间分析-提取分析）来裁切，裁切的矢量要素通常是一个要素图层或Shape文件。如果要进行批量处理，可以使用ToolBox中的批量处理工具。但是，有时我们选择的裁切矢量要素不是一个要素图层，而是要素图层中的一个要素，如按标准图幅裁切栅格影像，以前的做法是将要素一个一个的导出来，再进行批量处理。快速导出方法：ArcGIS按字段属性分割文件。

现在，可以使用Python脚本来批量裁切，在网上找了位大神的原始代码学习一下：

#-*- encoding:UTF-8 -*-

# Tool Name: 遍历要素裁剪栅格

# Source Name: extractByMask.py

# Version: ArcGIS 10.1

# Author: GIScloud

# Created: 2012/11/1

import arcpy

import string

from arcpy.sa import *

try:

    raster = arcpy.GetParameterAsText(0) #要裁剪的栅格

    clip_feat = arcpy.GetParameterAsText(1) #裁剪要素类

    field = arcpy.GetParameterAsText(2) #命名字段

    outworkspace = arcpy.GetParameterAsText(3) #命名字段裁剪后输出目录

    for row in arcpy.SearchCursor(clip_feat):

        mask=row.getValue("Shape")

        outPath=outworkspace+"\\"+str(row.getValue(field))

        outExtractByMask = ExtractByMask(raster,mask)

        outExtractByMask.save(outPath)

except arcpy.ExecuteError:

    print arcpy.GetMessages()

　　

将脚本导入ToolBox中（导入方法：ArcGIS使用Python脚本工具），执行报错：UnicodeEncodeError: 'ascii' codec can't encode characters in position 32-34: ordinal not in range(128)问题，使用的是ulipad编译器。

看样子是编码的问题，查了一下解决方案：在开头加上代码设置默认编辑

import sys

reload(sys)

sys.setdefaultencoding( "utf-8" )

　　

• 能够实现部分裁切了，但最后一个出了问题，看消息输出的文件名中怎么多了一个句号？

查看属性表，没有问题。将输出字段"北川羌族自治县"改为"北川羌族自治"能够正确输出，如果改成"北川羌族自治北川羌族自治"，又报相同的错误。说明应该就是输出名称长度的问题。

还有一个问题，输出的不同栅格边界接连有点小问题，如果避免这种情况，只能更改系统变量，输出更小的单元格大小。

参考：http://bbs.esrichina-bj.cn/ESRI/viewthread.php?tid=121642

http://blog.sina.com.cn/s/blog_64a3795a01018vyp.html