ArcPy按掩膜提取ExtractByMask的边界问题

问题描述

按掩膜提取是GIS数据处理中常用的功能之一。按照官方的示例，使用ArcPy的ExtractByMask方法可实现提取，代码如下：

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")

但实际提取结果，边界出现明显锯齿和异常值。

使用工具箱工具截切的默认效果如下：

解决方案

这应该是ArcPy的Bug，也有解决方案，可以使用arcpy.gp.ExtractByMask_sa方法来平替：

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

env.workspace = "C:/sapyexamples/data"
arcpy.gp.ExtractByMask_sa("elevation", "mask.shp","maskextract")