利用arcpy库实现按掩膜提取时边界出现异常值问题

一、介绍如何查询arcpy某个空间处理函数：

利用arcpy库实现按掩膜提取前，需先查找按掩膜提取的函数名称以及实现方法：
可以用ArcMap自带的帮助文档查看函数细节如下：

 

也可以利用ArcMap搜索框中项目描述查找如下：

 二、利用arcpy库实现按掩膜处理时遇到的问题

 首先查看ArcMap帮助文档中函数如下：

 按输入输出要求编写代码后，查看执行结果如下所示：

 执行得到的结果存在锯齿形边界，边界处存在部分负值，此方法无效。

三、解决这一问题

对比发现利用ArcMap搜索框中项目描述查找的函数与帮助文档中有所不同：

利用该 arcpy.gp.ExtractByMask_sa(in_raster, in_mask_data, out_raster)函数执行顺利，因此在利用arcpy库实现按掩膜提取时应当选用ExtractByMask_sa该函数。

最后实现这一过程编写完整的python代码如下：

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

#首先设置工作目录
# Set workspace and extent environments
arcpy.env.workspace = r"工作环境目录"
mask = r"待掩膜提取shp路径及名称"
# arcpy.env.extent = mask

# Set Snap Raster environment
dataset = r"要配准的栅格数据集路径及名称"
arcpy.env.snapRaster = dataset
spatial_ref = arcpy.Describe(dataset).spatialReference
arcpy.env.outputCoordinateSystem = spatial_ref

filePath = "待掩膜提取的文件目录"
outFilePath = "掩膜提取的结果存放目录"
fileList = os.listdir(filePath)

for file in fileList:
    if file.split('.')[-1] == "tif":
        comFilePath = os.path.join(filePath, file)
        outFileName = file.replace("chn", "Ytz")
        outFilePathName = os.path.join(outFilePath, outFileName)
        outExtractByMask = arcpy.gp.ExtractByMask_sa(comFilePath, mask, outFilePathName)
        print("Finished! -" + outFileName)

print("------ Finished all ------")