在土地资源管理、城市规划、生态保护等 GIS 工作中，你是否经常遇到这样的场景：需要用多个图层依次对同一目标要素进行擦除操作？比如在更新土地利用现状时，用不同年度的建设用地范围图层逐步擦除旧图斑；或是在处理生态红线与永久基本农田冲突时，依次用多个管控图层剔除不符合要求的地块。

一、什么是梯次擦除？

目标要素：某区域的土地利用现状图斑

擦除图层 1：生态保护红线范围

擦除图层 2：永久基本农田范围

擦除图层 3：规划道路红线

二、核心应用领域

三、为什么说它便捷高效？

1. 少量代码，强大功能

defgradual_erase(target_feature,erase_features,output_name,tolerance):

1. 复制初始数据，避免修改原始数据

initial_copy="in_memory\\initial_{}".format(timestamp)

arcpy.management.CopyFeatures(target_feature,initial_copy)

current_result=initial_copy

2. 循环执行擦除操作

fori,erase_featureinenumerate(erase_features,1):

temp_output="in_memory\\erase_{}_{}".format(timestamp,i)

arcpy.analysis.Erase(current_result,erase_feature,temp_output,tolerance)

current_result=temp_output

3. 保存最终结果并清理临时文件

arcpy.management.CopyFeatures(current_result,output_name)

2. 全流程自动化

自动验证输入要素是否存在、输出目录是否有效

自动生成唯一的临时文件，避免数据冲突

自动清理内存中的临时数据，释放系统资源

自动将结果添加到当前地图，方便后续分析

3. 操作简单，上手快

在 ArcMap 中添加脚本工具

设置参数：输入源要素、擦除要素列表、输出要素路径、容差（可选）

点击 “确定”，等待工具自动完成所有擦除操作

四、代码背后的设计巧思

内存工作空间

：所有中间数据都存储在in_memory工作空间中，大幅提升处理速度，同时避免磁盘空间占用。

进度条反馈

：通过 ArcGIS 的进度条功能，实时显示擦除进度，让你清晰掌握工具执行状态。

完善的错误处理

：对输入参数、数据存在性、执行过程进行全面校验，遇到错误时给出明确提示，避免数据损坏。

兼容性优化

：同时支持 Python 2 和 Python 3 环境，兼容不同版本的 ArcGIS Desktop。

五、结语