ArcGIS实战：如何将SHP文件高效转换为KML/KMZ格式并实现面要素标注

1次阅读

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

最近在参与一个应急指挥系统项目时，遇到一个典型问题：需要将行政区划 SHP 文件转换为 KML 格式供移动端使用，但转换后发现所有面要素的标注名称全部丢失。领导在平板电脑上查看时，只能看到一片空白区域，完全分不清哪个区域对应哪个乡镇。同样的问题也出现在环保监测系统中——排污许可证地块的编号在转换后消失，导致现场执法人员无法快速定位检查目标。

遇到这个问题后，我对比了三种主流解决方案：

GDAL 命令行工具：虽然轻量级，但需要额外安装库，且样式配置需要编写复杂的 XML
QGIS 图形界面：可视化操作友好，但批量处理需要手动点击，难以集成到自动化流程
ArcPy 脚本：作为 ArcGIS 原生 Python 模块，完美支持 ArcGIS 的符号系统，还能直接调用 ArcMap 中配置好的标注样式

最终选择 ArcPy 方案，因为它能：

与企业现有 ArcGIS 平台无缝集成
通过 Python 脚本实现批量自动化处理
完整保留原始数据的符号化和标注设置

import arcpy
from arcpy import env

# 设置工作空间
env.workspace = r"C:\GIS_Data\Admin"

# 关键转换函数
def shp_to_kmz(input_shp, output_kmz, label_field=None):
    try:
        # 坐标系处理（WGS84 的 WKT 定义）wgs84_wkt = """GEOGCS[\"GCS_WGS_1984\",DATUM[\"D_WGS_1984\",SPHEROID[\"WGS_1984\",6378137.0,298.257223563]],
        PRIMEM[\"Greenwich\",0.0],UNIT[\"Degree\",0.0174532925199433]]"""

        # 临时要素类（解决中文路径问题）temp_fc = arcpy.CreateUniqueName("temp_layer")
        arcpy.MakeFeatureLayer_management(input_shp, temp_fc)

        # 标注设置
        if label_field:
            arcpy.AddMessage(f"正在设置标注字段: {label_field}")
            arcpy.DisplayAnnotation_convert(temp_fc, "ANNOTATION_LAYER", 
                                          label_field, 
                                          "Arial 10")

        # 执行转换
        arcpy.FeatureClassToKML_conversion(
            temp_fc, 
            output_kmz,
            "0",  # 图层缩放范围
            False, # 是否包含要素属性
            "",    # 样本符号路径
            wgs84_wkt)

    except arcpy.ExecuteError:
        arcpy.AddError(arcpy.GetMessages(2))
    except Exception as e:
        arcpy.AddError(str(e))
    finally:
        # 清理临时数据
        if arcpy.Exists(temp_fc):
            arcpy.Delete_management(temp_fc)

# 在转换前检查并拆分多部件要素
def check_multipart(feature_class):
    with arcpy.da.SearchCursor(feature_class, ["OID@", "SHAPE@"]) as cursor:
        for row in cursor:
            if row[1].isMultipart:
                arcpy.AddWarning(f"要素 {row[0]} 是多部件面，建议使用 Multipart To Singlepart 工具拆分")

# 处理超长字段（KML 限制 255 字符）def truncate_fields(fc, max_length=250):
    fields = [f.name for f in arcpy.ListFields(fc) 
              if f.type in ["String", "Text"]]

    with arcpy.da.UpdateCursor(fc, fields) as cursor:
        for row in cursor:
            for i in range(len(row)):
                if row[i] and len(str(row[i])) > max_length:
                    row[i] = str(row[i])[:max_length] + "..."
            cursor.updateRow(row)

# 批量转换脚本
import os

def batch_convert_shp_to_kmz(input_folder, output_folder):
    # 创建输出目录
    if not os.path.exists(output_folder):
        os.makedirs(output_folder)

    # 遍历所有 SHP 文件
    for shp_file in arcpy.ListFiles("*.shp"):
        try:
            # 构建输出路径
            base_name = os.path.splitext(shp_file)[0]
            output_kmz = os.path.join(output_folder, f"{base_name}.kmz")

            # 预处理
            check_multipart(shp_file)
            truncate_fields(shp_file)

            # 执行转换（假设使用 NAME 字段作为标注）shp_to_kmz(shp_file, output_kmz, "NAME")

            arcpy.AddMessage(f"成功转换: {shp_file} -> {output_kmz}")
        except Exception as e:
            arcpy.AddError(f"转换失败 {shp_file}: {str(e)}")

# 使用示例
if __name__ == "__main__":
    batch_convert_shp_to_kmz("C:/input_shapes", "C:/output_kmz")