ArcGIS空间数据处理：如何高效合并两个单独面为单一数据集

1次阅读

共计 3557 个字符，预计需要花费 9 分钟才能阅读完成。

在 GIS 日常工作中，合并多个面要素的需求非常普遍。比如需要将相邻的行政区划合并为更大的管理单元，或者把零散的土地利用图斑整合为连片区域。传统的手动操作方式是通过 ArcMap 的 Editor 工具逐个选择面要素进行合并，这种方式存在几个明显问题：

操作繁琐耗时，特别是面对成百上千个面要素时
容易遗漏某些面要素，导致数据不完整
属性信息处理困难，难以保持数据一致性
无法形成可重复使用的流程

ArcGIS 平台提供了三种主要的自动化合并方案：

ArcPy 脚本
优点：灵活性最高，可定制所有处理细节，适合批量处理
缺点：需要 Python 编程基础，调试成本较高
ModelBuilder
优点：可视化搭建流程，适合简单合并需求
缺点：处理复杂逻辑时模型会变得臃肿
内置地理处理工具
优点：开箱即用，操作简单
缺点：功能相对固定，无法处理特殊需求

对于专业 GIS 开发人员，我们推荐使用 ArcPy 方案，它能在保证灵活性的同时实现高效处理。

使用 arcpy.Merge_management() 是最直接的合并方式：

import arcpy

# 设置工作空间
arcpy.env.workspace = r"C:\data\gis_db.gdb"

# 输入要素列表
input_features = ["parcels_A", "parcels_B"]

# 输出合并后的要素类
output_feature = "merged_parcels"

# 执行合并
arcpy.Merge_management(input_features, output_feature)

实际工作中经常会遇到面要素边界不完全重合的情况，需要添加拓扑检查：

try:
    # 先检查并修复几何错误
    arcpy.CheckGeometry_management(input_features, "geometry_errors")
    arcpy.RepairGeometry_management(input_features)

    # 执行合并
    arcpy.Merge_management(input_features, output_feature)

except arcpy.ExecuteError as e:
    print(f"合并过程中出现错误: {e}")
    # 记录错误日志
    with open("merge_errors.log", "a") as log_file:
        log_file.write(f"{arcpy.GetMessages()}\n")

当输入要素的属性结构不一致时，需要特别处理字段映射：

# 创建字段映射对象
field_mappings = arcpy.FieldMappings()

# 添加需要保留的字段
for field in ["LANDUSE", "AREA", "OWNER"]:
    fmap = arcpy.FieldMap()
    fmap.addInputField(input_features[0], field)
    field_mappings.addFieldMap(fmap)

# 执行带字段映射的合并
arcpy.Merge_management(input_features, output_feature, field_mappings)

下面是一个包含完整功能的 Python 脚本：

# -*- coding: utf-8 -*-
"""
面要素合并工具
功能：合并多个面要素类，处理拓扑错误和属性映射
"""
import arcpy
import os
import time

def merge_features(input_list, output_fc, keep_fields=None):
    """
    合并多个面要素类
    :param input_list: 输入要素路径列表
    :param output_fc: 输出要素路径
    :param keep_fields: 需要保留的字段列表
    """
    # 设置进度条
    arcpy.SetProgressor("step", "正在合并面要素...", 0, 4)

    try:
        # 步骤 1：检查输入数据
        arcpy.SetProgressorLabel("检查输入数据...")
        for fc in input_list:
            if not arcpy.Exists(fc):
                raise ValueError(f"输入要素 {fc} 不存在")
        arcpy.SetProgressorPosition(1)

        # 步骤 2：处理几何错误
        arcpy.SetProgressorLabel("检查并修复几何错误...")
        for fc in input_list:
            arcpy.CheckGeometry_management(fc, "in_memory/geom_errors")
            arcpy.RepairGeometry_management(fc)
        arcpy.SetProgressorPosition(2)

        # 步骤 3：设置字段映射
        arcpy.SetProgressorLabel("配置字段映射...")
        field_mappings = arcpy.FieldMappings()

        if keep_fields:
            for field in keep_fields:
                fmap = arcpy.FieldMap()
                for fc in input_list:
                    if field in [f.name for f in arcpy.ListFields(fc)]:
                        fmap.addInputField(fc, field)
                if fmap.inputFieldCount > 0:
                    field_mappings.addFieldMap(fmap)
        arcpy.SetProgressorPosition(3)

        # 步骤 4：执行合并
        arcpy.SetProgressorLabel("执行合并操作...")
        arcpy.Merge_management(input_list, output_fc, field_mappings)
        arcpy.SetProgressorPosition(4)

        print(f"成功合并要素到 {output_fc}")

    except Exception as e:
        print(f"处理失败: {str(e)}")
        raise

    finally:
        # 清理临时数据
        if arcpy.Exists("in_memory/geom_errors"):
            arcpy.Delete_management("in_memory/geom_errors")

# 使用示例
if __name__ == "__main__":
    # 输入参数
    input_features = [
        r"C:\data\gis_db.gdb\parcels_east",
        r"C:\data\gis_db.gdb\parcels_west"
    ]
    output_feature = r"C:\data\gis_db.gdb\all_parcels"
    fields_to_keep = ["LANDUSE", "AREA", "OWNER"]

    # 记录开始时间
    start_time = time.time()

    # 执行合并
    merge_features(input_features, output_feature, fields_to_keep)

    # 计算耗时
    print(f"处理完成，耗时 {time.time()-start_time:.2f} 秒")

当处理大型数据集时，可以考虑以下优化策略：