1. 项目背景与需求解析
在地理信息数据处理工作中，我们经常需要处理大量分散存储的Shapefile文件。最近我在处理一个县域基础地理信息项目时，就遇到了这样的需求：颍上县23个乡镇的规划数据分散存储在近百个不同文件夹中，每个文件夹包含若干.shp文件，需要快速统计每个图层的要素数量以及总数量。
传统的手动统计方式存在三个明显痛点：
- 文件数量庞大时（如超过50个.shp），ArcGIS桌面端逐个打开查看属性表的方式效率极低
- 不同文件夹中的同名.shp文件容易造成统计混淆
- 人工记录容易出错且无法形成结构化统计报表
2. 技术方案设计
2.1 核心工具选型
采用ArcPy作为解决方案的核心工具链，主要基于以下考虑：
- 原生支持ESRI格式数据读写，无需额外转换
- 提供完整的空间数据访问接口（如Describe函数）
- 可与ArcGIS Pro环境无缝集成
- 支持递归目录遍历和批量处理
2.2 关键函数说明
实现本需求主要依赖三个核心函数：
arcpy.da.Walk：递归遍历目录树，获取所有.shp文件路径arcpy.MakeFeatureLayer_management：创建临时图层对象arcpy.GetCount_management：获取要素数量统计
特别注意：MakeFeatureLayer创建的临时图层在脚本结束后会自动释放，不会占用额外存储空间
3. 完整实现代码
python复制import arcpy
import os
def batch_shp_stats(root_folder, output_csv):
"""
批量统计多个文件夹中shp文件的要素数量
:param root_folder: 根目录路径
:param output_csv: 输出统计结果CSV路径
"""
# 创建结果CSV文件并写入表头
with open(output_csv, 'w') as f:
f.write("文件路径,图层名称,要素数量\n")
# 设置工作空间环境
arcpy.env.workspace = root_folder
arcpy.env.overwriteOutput = True
# 递归遍历所有子目录
for dirpath, dirnames, filenames in arcpy.da.Walk(root_folder):
for filename in filenames:
if not filename.lower().endswith('.shp'):
continue
# 构建完整文件路径
full_path = os.path.join(dirpath, filename)
try:
# 创建临时图层
layer_name = "temp_layer"
arcpy.MakeFeatureLayer_management(full_path, layer_name)
# 获取要素数量
count = int(arcpy.GetCount_management(layer_name)[0])
# 写入统计结果
with open(output_csv, 'a') as f:
f.write(f"{full_path},{filename[:-4]},{count}\n")
except Exception as e:
print(f"处理文件 {full_path} 时出错: {str(e)}")
continue
if __name__ == "__main__":
# 使用示例
input_folder = r"D:\GIS数据\颍上县规划数据"
output_file = r"D:\统计结果.csv"
batch_shp_stats(input_folder, output_file)
4. 关键技术细节解析
4.1 路径处理注意事项
- 使用
os.path.join构建路径可确保跨平台兼容性 - 文件扩展名判断采用
.lower()处理大小写问题 - 输出CSV路径建议使用绝对路径
4.2 性能优化技巧
- 设置
arcpy.env.overwriteOutput = True避免重复确认 - 将结果写入操作放在with语句中确保文件正确关闭
- 使用try-except捕获单个文件处理异常，不影响整体流程
4.3 坐标系处理方案
当遇到CGCS2000等特定坐标系需求时，可在统计前添加坐标系检查：
python复制# 在GetCount之前添加坐标系检查
sr = arcpy.Describe(layer_name).spatialReference
if sr.name != "CGCS2000_3_Degree_GK_Zone_40":
print(f"警告：{filename} 使用非常用坐标系 {sr.name}")
5. 常见问题与解决方案
5.1 字段名称截断问题
当处理从地理数据库导出的shp时，可能会遇到字段名被截断为10个字符的情况。可通过以下方式获取原始字段信息：
python复制fields = [f.name for f in arcpy.ListFields(layer_name)]
print("完整字段列表:", fields)
5.2 要素数量统计异常
如果遇到统计数量为0的异常情况，建议按以下步骤排查：
- 检查文件是否被其他程序锁定
- 验证shp文件是否包含有效要素
- 检查是否有SQL查询条件意外过滤了所有要素
5.3 CAD转换数据问题
对于从CAD转换来的shp数据，需特别注意：
- 图层命名可能包含特殊字符
- 要素类型可能混合（点、线、面共存）
- 属性表结构可能非常规
6. 扩展应用场景
本脚本可轻松扩展支持以下需求：
- 按要素类型分类统计（点、线、面）
- 空间范围交叉分析
- 属性字段值分布统计
- 自动生成统计图表
例如添加要素类型统计：
python复制def get_geometry_type(layer):
desc = arcpy.Describe(layer)
return desc.shapeType
# 在统计循环中添加
geom_type = get_geometry_type(layer_name)
实际项目中，我将此脚本与QGIS结合使用，实现了：
- 自动生成各乡镇数据完整性报告
- 快速定位数据缺失区域
- 项目进度可视化看板