Arcpy 实现批量按掩膜提取

Arcpy 实现批量按掩膜提取，具体基于像元的原理如上图所示。

 

用法

• 按掩膜提取工具与设置掩膜环境大体相同，但二者的区别在于输入掩膜只适用于即时实例，而在环境中设置的掩膜可在发生更改或被禁用之前应用于所有工具。

• 当输入为多波段栅格时，将输出一个新的多波段栅格。对输入多波段栅格中的每一单个波段都会进行相应分析。

  默认输出格式是Esri Grid 堆栈。请注意，Esri Grid 堆栈的名称不能以数字开头、不能使用空格，并且长度不能超过 9 个字符。

• 如果输入内容为根据多波段栅格（超过三个波段）创建的图层，提取操作将只考虑图层加载（符号化）的波段。因此，输出的多波段栅格只会具有三个波段，分别对应于输入图层显示中使用的多波段栅格。

• 如果输入掩膜是栅格，则会将非 NoData 输入像元位置的值复制到输出栅格。可用来创建掩膜栅格的工具包括条件函数、条件测试，以及提取分析工具集中的其他工具。

• 当为输入栅格掩膜指定多波段栅格时，将只在运算中使用第一个波段。

• 如果输入栅格为整型，那么输出栅格也为整型。如果输入栅格为浮点型，则输出栅格也为浮点型。 

• 摘自ArcGIS帮助文档

python代码实现：

# author:zhoulong_GISER
# -*- coding: utf-8 -*-
# blog:https://blog.csdn.net/qq_33356563
# E-mail:zhoulong_giser@163.com

import arcpy,os,glob
from arcpy import env

arcpy.CheckOutExtension("Spatial")
filepath=r"E:\biyesheji\t_1488599756442_0531\data"#输入文件文件夹
env.workspace = filepath
os.chdir(filepath)
Rasters = glob.glob("*.tif")
outfile=r'E:\biyesheji\t_1488599756442_0531\out'#输出文件文件夹
inMaskData = r"E:\biyesheji\t_1488599756442_0531\filter\EYM1.tif"#过滤文件文
for Raster in Rasters:
    inRaster = Raster
    arcpy.gp.ExtractByMask_sa(Raster, inMaskData, outfile+"\\"+Raster.split(".")[0]+'.img')

效果图示：

工具箱：

地址：https://download.csdn.net/download/qq_33356563/11020148