ArcGIS批量合成文件夹图层数据的自动化解决方案

1次阅读

没有评论

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

在 GIS 数据处理工作中，经常需要处理分散在不同文件夹中的 Shapefile 或 FeatureClass 数据。传统的手动处理方法存在以下问题：

重复操作：每个文件夹都需要单独打开、选择图层、合并，操作繁琐
效率低下：当数据量超过 1000 个文件时，人工操作可能需要数小时甚至数天
易出错：手动操作容易遗漏文件或选择错误图层
一致性差：难以保证所有数据的坐标系统一和属性表结构一致

本方案基于 arcpy 模块实现，核心思路是：

使用 arcpy.ListDatasets 递归遍历文件夹结构
通过 arcpy.Merge_management 批量合并图层数据

与 GDAL/OGR 等替代方案相比，arcpy 具有以下优势：

原生支持 ArcGIS 数据格式
更好的拓扑关系保持能力
更完整的坐标系转换支持

以下是完整的 Python 脚本实现：

import arcpy
import os

def batch_merge_folders(input_folder, output_fc):
    """
    批量合并文件夹中的图层数据
    :param input_folder: 输入文件夹路径
    :param output_fc: 输出要素类路径
    """
    try:
        # 创建空要素类用于存储合并结果
        arcpy.CreateFeatureclass_management(os.path.dirname(output_fc), 
            os.path.basename(output_fc), 
            "POLYGON"  # 根据实际数据类型调整
        )

        # 递归遍历文件夹
        datasets = []
        for root, dirs, files in os.walk(input_folder):
            for file in files:
                if file.endswith('.shp'):  # 支持 Shapefile
                    datasets.append(os.path.join(root, file))
                elif file.endswith('.gdb'):  # 支持 File Geodatabase
                    arcpy.env.workspace = os.path.join(root, file)
                    for fc in arcpy.ListFeatureClasses():
                        datasets.append(os.path.join(root, file, fc))

        # 分块处理大文件
        chunk_size = 100  # 每批处理 100 个文件
        for i in range(0, len(datasets), chunk_size):
            chunk = datasets[i:i + chunk_size]
            # 坐标系一致性校验
            for fc in chunk:
                desc = arcpy.Describe(fc)
                if desc.spatialReference.name != "WGS_1984":  # 此处需处理 WGS84 与 CGCS2000 的投影转换
                    arcpy.Project_management(fc, fc + "_projected", arcpy.SpatialReference(4326))
                    chunk[chunk.index(fc)] = fc + "_projected"

            # 批量合并
            arcpy.Merge_management(chunk, output_fc)

            # 清理临时文件
            for fc in chunk:
                if "_projected" in fc:
                    arcpy.Delete_management(fc)

    except arcpy.ExecuteError as e:
        print(f"执行错误: {e}")
    except Exception as e:
        print(f"其他错误: {e}")

使用分块处理 (chunking) 技术，避免一次性加载所有数据
启用 arcpy.env.overwriteOutput = True，避免重复提示覆盖确认
在处理完成后及时删除临时文件，释放内存

# 合并后自动检查拓扑错误
try:
    arcpy.CheckGeometry_management(output_fc, "拓扑错误检查")
    arcpy.RepairGeometry_management(output_fc)
except:
    print("拓扑检查失败，请手动检查数据")