使用ArcPy批量统计Shapefile要素数量的高效方法

1. 项目背景与需求解析

在地理信息系统（GIS）数据处理工作中，我们经常需要处理大量分散存储的Shapefile文件（*.shp）。这些文件可能分布在不同的文件夹中，每个文件包含多个地理要素图层。面对这种情况，手动统计每个图层的要素数量不仅效率低下，而且容易出错。

最近我在处理一个县域基础地理信息项目时，就遇到了这样的挑战：需要统计颍上县200多个村级行政区划的建筑物轮廓数据，这些数据分散在15个不同的文件夹中，每个文件夹包含10-20个Shapefile文件。手动操作几乎不可能在短时间内完成，这正是arcpy脚本自动化处理的典型应用场景。

2. 技术方案设计

2.1 核心工具选型

ArcPy作为ArcGIS的Python站点包，提供了完整的GIS数据处理能力。本次任务主要用到以下核心模块：

• arcpy.da.Walk：递归遍历文件夹结构
• arcpy.Describe：获取图层属性信息
• arcpy.GetCount_management：统计要素数量

选择这些工具的主要考虑是：

1. arcpy.da.Walk比传统的os.walk能更好地识别GIS数据格式
2. Describe对象提供了完整的图层元数据访问接口
3. GetCount是经过优化的要素统计工具，比Python内置方法效率更高

2.2 处理流程设计

完整的处理流程分为四个阶段：

1. 文件遍历阶段：递归扫描目标文件夹，定位所有.shp文件
2. 数据加载阶段：将每个.shp文件作为临时图层加载
3. 统计分析阶段：获取图层要素数量和属性信息
4. 结果输出阶段：将统计结果写入CSV或Excel文件

python复制# 基础代码框架示例
import arcpy
import os

def batch_count_shp_features(root_folder, output_csv):
    # 初始化结果列表
    results = []
    
    # 第一阶段：文件遍历
    for dirpath, dirnames, filenames in arcpy.da.Walk(root_folder):
        for filename in filenames:
            if filename.lower().endswith('.shp'):
                shp_path = os.path.join(dirpath, filename)
                
                # 第二阶段：数据加载
                lyr_name = "temp_layer"
                arcpy.MakeFeatureLayer_management(shp_path, lyr_name)
                
                # 第三阶段：统计分析
                desc = arcpy.Describe(lyr_name)
                count = int(arcpy.GetCount_management(lyr_name).getOutput(0))
                
                # 记录结果
                results.append([shp_path, desc.name, count])
    
    # 第四阶段：结果输出
    with open(output_csv, 'w') as f:
        f.write("文件路径,图层名称,要素数量\n")
        for row in results:
            f.write(",".join(map(str, row)) + "\n")

3. 完整实现与优化

3.1 增强版实现代码

基础版本虽然能用，但在实际项目中还需要考虑更多细节。以下是经过生产环境验证的增强版实现：

python复制import arcpy
import os
import csv
from datetime import datetime

def batch_count_shp_features(root_folders, output_csv, coordinate_system=None):
    """
    批量统计多个文件夹中所有Shapefile的图层和要素数量
    
    参数:
        root_folders: