ArcGIS实战：批量转换SHP至KML并自动标注图斑编号的技术方案

1次阅读

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

在国土调查、林业普查等 GIS 项目中，我们经常需要将 SHP 格式的矢量数据批量转换为 KML 以便在 Google Earth 等平台展示。传统手动操作存在三大痛点：

效率低下：单个文件转换需重复点击 10 余次菜单，100 个文件需耗时约 2 小时
标注易错：图斑编号需手动复制粘贴，易出现遗漏或错位
信息丢失：转换过程中属性表字段经常无法完整保留

格式转换引擎 ：ArcPy 的ConvertFeaturesToKML 函数，支持保持拓扑结构
属性提取器 ：SearchCursor 遍历属性表，精准获取图斑编号字段
动态标注器 ：通过修改 KML 的<name> 标签实现编号自动注入

坐标系转换：强制输出为 WGS84（EPSG:4326）避免位置偏移
字段映射：建立 SHP 属性字段与 KML Placemark 的对应关系
批量处理：自动遍历文件夹内所有.shp 文件

import arcpy
import os
import logging

# 日志配置
logging.basicConfig(filename='shp2kml.log', level=logging.INFO)

def batch_convert(folder_path, output_folder, id_field='图斑编号'):
    """
    批量转换 SHP 到 KML 并标注编号
    :param folder_path: 输入文件夹路径
    :param output_folder: 输出文件夹路径
    :param id_field: 图斑编号字段名
    """
    try:
        # 创建输出目录
        if not os.path.exists(output_folder):
            os.makedirs(output_folder)

        # 遍历文件夹
        for shp_file in [f for f in os.listdir(folder_path) if f.endswith('.shp')]:
            input_path = os.path.join(folder_path, shp_file)
            output_name = os.path.splitext(shp_file)[0] + '.kmz'
            output_path = os.path.join(output_folder, output_name)

            # 动态设置图层名称（关键步骤）lyr_name = 'temp_lyr'
            arcpy.MakeFeatureLayer_management(input_path, lyr_name)

            # 使用属性标注要素
            with arcpy.da.SearchCursor(lyr_name, [id_field, 'SHAPE@']) as cursor:
                for row in cursor:
                    # 此处可扩展为更复杂的标注逻辑
                    arcpy.AddMessage(f"处理图斑：{row[0]}")

            # 执行 KML 转换
            arcpy.conversion.ConvertFeaturesToKML(
                lyr_name, 
                output_path,
                layer_output_scale=1,
                is_composite=False,
                label_field=id_field  # 关键参数：绑定标注字段
            )
            logging.info(f"成功转换：{shp_file}")

    except Exception as e:
        logging.error(f"转换失败：{str(e)}", exc_info=True)
        raise

if __name__ == '__main__':
    batch_convert(r"C:\input_shp", r"C:\output_kml")

from multiprocessing import Pool

def worker(args):
    """包装为可并行执行的函数"""
    shp_path, output_path, id_field = args
    # 此处插入单个文件转换逻辑

if __name__ == '__main__':
    file_list = [(f1, out1, 'ID'), (f2, out2, 'ID')]  # 构建参数列表
    with Pool(processes=4) as pool:  # 4 进程并行
        pool.map(worker, file_list)

显式释放游标：在 SearchCursor 使用后执行del cursor
禁用中间缓存：设置 arcpy.env.scratchWorkspace 到指定目录
分块处理：超大型文件建议按图幅拆分处理

# 解决方案：使用原始字符串 (raw string) 或双反斜杠
input_path = r"C:\ 国土数据 \ 林地.shp"  # 推荐方式
# 或
input_path = "C:\\ 国土数据 \\ 林地.shp"

简化几何：使用 arcpy.Generalize_edit() 减少节点数
强制 2D：arcpy.FeatureTo2D_management()去除 Z 值
分割处理：对超 5000 个顶点的面执行arcpy.Dissolve_management()

日期字段：预先使用 arcpy.ConvertTimeField_management() 标准化
浮点精度：设置 arcpy.env.outputCoordinateSystem 明确坐标系
文本截断：超过 255 字符的字段需提前截断

# Flask 示例：将转换功能封装为 API
from flask import Flask, send_file

app = Flask(__name__)

@app.route('/convert/<shp_name>')
def convert(shp_name):
    output = do_conversion(shp_name)
    return send_file(output, as_attachment=True)

维度	ArcPy 方案	FME 方案
处理速度	中等（单线程）	快（多线程优化）
开发成本	低（Python 生态）	高（需许可）
扩展性	灵活（可集成 AI 模型）	受限（依赖转换器）