ArcGIS Python使用及批量裁剪示例

版权声明：本文为博主原创文章，欢迎署名转载。 https://blog.csdn.net/dou3516/article/details/82762789

ArcGIS中的工具均提供Python函数，方便支持批量处理。以掩膜裁剪（提取分析/Extract by Mask）为例，

可以查找到帮助文档“按掩膜提取”（源1）（源2），路径“工具参考/工具/Spatial Analyst 工具箱/按掩膜提取”函数为：Extract by Mask

Extract by Mask的基本语法是：

ExtractByMask (in_raster, in_mask_data)

in_raster：提取像元的输入栅格。

in_mask_data：用于定义提取区域的输入掩膜数据。它可以是栅格，也可以是要素数据集。

代码示例

以下以批量读取裁剪某文件夹下的tif数据为例(特别注意Python2脚本循环语句格式的空格不能用Tab，否则会报错Indent；

另外在使用不同的矢量裁剪同一影像时，为保证裁剪后范围一致需设置范围为输入栅格范围。例如下面的：

rd = arcpy.sa.Raster(inRaster)
# Set the extent environment as the raster, very important for clip with different vector
arcpy.env.extent = rd.extent)：

# Import system modules
import arcpy
from arcpy import env
from arcpy.sa import *

# Set environment settings
env.workspace = "D:/img"
rasterList = arcpy.ListRasters("*","tif")

#输出路径
output_path = "D:/imgMask/"
# mask shp
inMaskData = "D:/range.shp"

for raster in rasterList:
    print raster
    # Set local variables
    inRaster = raster
    rd = arcpy.sa.Raster(inRaster)
    # Set the extent environment as the raster, very important for clip with different vector
    arcpy.env.extent = rd.extent
    # Execute ExtractByMask
    outExtractByMask = ExtractByMask(inRaster, inMaskData)
    # Save the output
    out = output_path + inRaster #对生成文件进行命名
    outExtractByMask1.save(out)

代码执行

方式一：将代码保存为单独的py文件，cmd或其它IDE执行。cmd执行命令（如何独立运行脚本）：

C:\Python27\ArcGIS10.3\python.exe my_script.py

方式二：ArcGIS安装的Python IDE里执行