KML/SHP/CAD格式转换 V1.0

12 种 GIS 格式任意双向转换 | DWG 丢进去直接出 KML/SHP | 坐标系统全自动识别

坤图_GIS

一、快速开始（3 步上手）

Step 1: 安装依赖

pip install geopandas pyproj ezdxf shapely pyarrow

Step 2: 开始转换

# DWG → KML（全自动，无需配置坐标系）
python scripts/gis_converter.py 我的图纸.dwg 输出.kml

# DWG → SHP（自动识别 CGCS2000 分带）
python scripts/gis_converter.py 地形图.dwg 成果.shp

# 其他格式互转
python scripts/gis_converter.py 边界.shp 数据.geojson

Step 3: 查看结果

转换完成后自动输出日志，包含：

• 自动识别的坐标系及置信度
• 要素数量和几何类型统计
• DWG→DXF 转换耗时（如有）

如果自动识别不确定，工具会主动提示：

检测到坐标位于湖北神农架区域，推测为 CGCS2000 3°带 第37带 (EPSG:4546)，是否正确？[Y/n]

直接回车确认即可，或输入 n 手动指定。

二、支持格式全矩阵

图例：GPKG = GeoPackage（OGC标准，SHP推荐替代）, WPMZ = DJI-WPMZ, Huace = 华测GNSS KML

三、完整命令参考

3.1 DWG/CAD 转 GIS（零配置推荐）

# 单文件转换
python scripts/gis_converter.py 地形图.dwg 成果.kml
python scripts/gis_converter.py 地形图.dwg 成果.shp
python scripts/gis_converter.py 地形图.dwg 成果.dxf       # DWG→DXF(版本转换)
python scripts/gis_converter.py 地形图.dwg 成果.geojson
python scripts/gis_converter.py 地形图.dwg 成果.gpkg       # GeoPackage

# 手动指定坐标系（覆盖自动检测）
python scripts/gis_converter.py 地形图.dxf 成果.shp --source-crs EPSG:3857

DWG 自动转换原理：

1. 检查系统是否安装 ODA File Converter（免费命令行工具）
2. 未安装则自动下载（~28MB）并静默安装
3. 在临时目录执行 DWG→DXF 转换
4. 后续自动识别坐标系 → 标准 GIS 转换

首次转换 DWG 需联网下载 ODA 引擎，之后无需重复下载。

3.2 奥维地图格式

# OVKML 自动 GCJ-02→WGS84 纠偏
python scripts/gis_converter.py 奥维收藏夹.ovkml 成果.shp

# OVJSN（奥维JSON）读写
python scripts/gis_converter.py 数据.ovjsn 成果.geojson
python scripts/gis_converter.py 成果.geojson 奥维格式.ovjsn

# 禁用自动纠偏
python scripts/gis_converter.py 数据.ovkml 成果.shp --no-gcj02-correct

3.3 大疆无人机航线

python scripts/gis_converter.py 测区边界.geojson 航线任务.wpmz \
  --drone "M350 RTK" \
  --payload "Zenmuse P1" \
  --altitude 80 \
  --speed 10 \
  --action takePhoto

DJI 支持设备：

航线参数说明：

3.4 华测 LandStar 格式

python scripts/gis_converter.py 地块.shp 华测导出.kml \
  --name-field "地块编号" \
  --desc-field "备注"

3.5 GeoPackage / GeoParquet（现代化格式）

# SHP → GeoPackage（UTF-8原生，无字段名限制）
python scripts/gis_converter.py 数据.shp 数据.gpkg

# GeoParquet（列式压缩 10:1~20:1，DuckDB兼容）
python scripts/gis_converter.py 大数据.shp 大数据.parquet

3.6 Python API 调用

from gis_converter import convert

# DWG 零配置转换（自动DWG→DXF + 自动识别坐标系）
convert("地形图.dwg", "成果.kml")
convert("地形图.dwg", "成果.shp")

# DXF 自动识别坐标系
convert("管线图.dxf", "成果.shp")

# 手动指定坐标系
convert("图纸.dxf", "成果.shp", source_crs="EPSG:4490")

# 生成大疆航线
convert("测区.geojson", "任务.wpmz",
        drone_model="M350 RTK",
        payload_model="Zenmuse H20",
        altitude=80, speed=10,
        action_type="takePhoto")

# 北京54→WGS84（自动诊断偏移量）
convert("老数据.shp", "成果.kml", source_crs="EPSG:4214")

四、坐标系处理规则

自动坐标系识别算法

适用于中国 CGCS2000 Gauss-Kruger 投影坐标：

1. 地理坐标检测：X 在 ±180°、Y 在 ±90° → EPSG:4326
2. Web Mercator 检测：投影坐标超出 GK 范围 → EPSG:3857
3. 中国 GK 分带匹配：
   • 从北坐标(Y)估算纬度 ≈ Y / 111,000
   • 根据纬度带确定区域参考经度（按中国行政区划）
   • 计算经度偏差 → 匹配最近 CGCS2000 3°/6° 分带

支持的分带：CGCS2000 3° 带 (EPSG:4534~4553，CM 75°135°E) 和 6° 带 (EPSG:45504562)

五、常见问题与解决方案

Q1: DWG 首次转换很慢

原因：首次需下载 ODA File Converter（~28MB） 解决：等待下载完成，之后同类转换秒级完成

Q2: SHP 中文乱码

原因：Shapefile 编码为 GBK/GB2312 但读取用 UTF-8 解决：python gis_converter.py input.shp output.kml --encoding gbk

Q3: SHP 字段名变短/截断

原因：DBF 格式限制，字段名不能超过 10 个英文字符 解决：使用 GeoPackage 替代（无字段名限制）： python gis_converter.py data.shp data.gpkg

Q4: 坐标系识别错误

原因：数据不在中国境内（境外投影坐标系），或为 CUSTOM 用户自定义投影 解决：手动指定： python gis_converter.py input.dxf output.shp --source-crs EPSG:XXXX

Q5: 北京54/西安80 转 WGS84 偏移 100~200 米

原因：pyproj 默认变换不使用精确的七参数/格网改正 解决：工具已自动诊断并输出偏移警告。如需精确转换，提供七参数后配置 pyproj 管线

Q6: GeoPackage 文件被锁定

原因：SQLite 单写入者限制，QGIS 等软件同时打开 解决：关闭所有打开该 gpkg 的软件后重试

Q7: DXF 属性丢失

原因：DXF 格式不支持属性表，仅保留几何图形和图层名 解决：转换时保留 layer 字段，可用 --name-field layer 按图层命名

Q8: OVOBJ 二进制无法解析

原因：奥维私有二进制格式，无公开解析规范 解决：通过奥维桌面端导出为 OVKML/OVJSN 后再转换

六、适用场景

七、依赖环境

八、已知限制

1. OVOBJ 二进制无法直接解析，需通过奥维桌面端导出为 OVKML/OVJSN
2. CASS XDATA 不支持（测绘行业专用 SOUTH-XDATA 编码）
3. BJ54/西安80 精确转换需七参数（已内置诊断提醒）
4. CRS 自动检测仅覆盖中国 CGCS2000 3°/6°带 + WGS84 + Web Mercator
5. DWG 首次需联网下载 ODA File Converter
6. SHP 字段名超过 10 字符自动截断（建议用 GeoPackage 替代）
7. 高纬度（>50°N）坐标系识别精度下降

九、日志与排错

每次转换自动生成日志文件：logs/gis_convert_YYYYMMDD_HHMMSS.log

日志包含：

• 输入/输出文件信息
• DWG→DXF 转换详情（如有）
• 自动坐标系识别结果及置信度
• 坐标转换参数
• 要素数量和几何类型统计
• 错误堆栈（如有异常）

十、格式规范参考

详见 references/formats.md

\x3C!-- wm:坤图_GIS -->