ArcGIS教程：按掩膜提取 (空间分析)

　　摘要

　　提取所定义掩膜区域内的栅格像元。

　　插图

　　



　　用法

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

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

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

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

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

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

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

　　代码实例

　　按掩膜提取 (ExtractByMask) 示例 1(Python 窗口)

　　本示例将从栅格内提取由输入面 shapefile 要素类所定义的掩膜内的像元。

　　import arcpy

　　from arcpy import env

　　from arcpy.sa import *

　　env.workspace = "C:/sapyexamples/data"

　　outExtractByMask = ExtractByMask("elevation", "mask.shp")

　　outExtractByMask.save("C:/sapyexamples/output/maskextract")

　　按掩膜提取 (ExtractByMask) 示例 2(独立脚本)

　　本示例将从栅格内提取由输入面 shapefile 要素类所定义的掩膜内的像元。

　　# Name: ExtractByMask_Ex_02.py

　　# Description: Extracts the cells of a raster that correspond with the areas

　　# defined by a mask.

　　# Requirements: Spatial Analyst Extension

　　# Import system modules

　　import arcpy

　　from arcpy import env

　　from arcpy.sa import *

　　# Set environment settings

　　env.workspace = "C:/sapyexamples/data"

　　# Set local variables

　　inRaster = "elevation"

　　inMaskData = "mask.shp"

　　# Check out the ArcGIS Spatial Analyst extension license

　　arcpy.CheckOutExtension("Spatial")

　　# Execute ExtractByMask

　　outExtractByMask = ExtractByMask(inRaster, inMaskData)

　　# Save the output

　　outExtractByMask.save("C:/sapyexamples/output/extractmask")