ArcPy 掩膜计算

### 使用 ArcPy 实现栅格掩膜计算

为了实现栅格掩膜计算，可以利用 `arcpy.sa.ExtractByMask` 函数。此函数允许提取输入栅格中的像元值，这些像元位于指定的掩膜要素内[^1]。

下面是一个简单的例子，展示如何设置工作空间并执行基于矢量图层作为掩膜的栅格裁剪操作：

import arcpy
from arcpy import env
from arcpy.sa import *

# 设置许可
arcpy.CheckOutExtension("Spatial")

# 定义工作区路径以及输入输出文件名
env.workspace = "C:/data"
inRaster = "landuse.tif"  # 输入栅格数据集名称
inMaskData = "mask.shp"   # 掩膜要素类或栅格数据集名称
outExtractByMask = "extracted_landuse.tif"

# 执行 Extract by Mask 工具
outExtract = ExtractByMask(inRaster, inMaskData)

# 将结果保存为 TIFF 文件
outExtract.save(outExtractByMask)

这段代码首先导入必要的模块，并设定好处理过程中需要用到的工作目录和其他参数；接着调用了 `ExtractByMask()` 方法来进行实际的数据提取操作；最后把得到的结果存储到了指定位置上[^2]。

需要注意的是，在运行上述脚本之前，请确保已经启用了 Spatial Analyst 扩展模块，并且拥有足够的权限访问所涉及的所有地理数据资源。