坐标转换技能 (coord-transform)

坐标系转换工具，支持 GeoJSON 和 WKT 格式，支持中国常用坐标系。

核心功能

• - 格式互转：WKT ↔ GeoJSON 互转
• 坐标转换：支持各种坐标系之间的转换
• 中国坐标系：完整支持 CGCS2000、WGS84、火星坐标系(GCJ-02)、百度坐标系(BD-09)

支持的坐标系

依赖

bash
apt-get install python3-pyproj

使用方式

基本语法

python3 scripts/coord_convert.py --input <文件或字符串> --from <源坐标系> --to <目标坐标系> [--format <输出格式>] [--output <输出文件>]

参数说明

使用示例

1. 坐标系统转换

CGCS2000 (4490) → WGS84 (4326)：
bash
python3 scripts/coord_convert.py --wkt POINT(113 23) --from 4490 --to 4326

CGCS2000 (4490) → 投影坐标系 (4547)：
bash
python3 scripts/coord_convert.py --wkt POLYGON((112.58 26.56,112.60 26.56,112.60 26.58,112.58 26.58,112.58 26.56)) --from 4490 --to 4547

2. WKT ↔ GeoJSON 格式互转

WKT 转 GeoJSON：
bash
python3 scripts/coord_convert.py --wkt POINT(113 23) --from 4326 --to 4490 --format geojson

GeoJSON 转 WKT：
bash
python3 scripts/coord_convert.py --input data.geojson --from 4490 --to 4326 --format wkt

3. 火星坐标系转换

WGS84 → 火星坐标系 (GCJ-02)：
bash
python3 scripts/coord_convert.py --wkt POINT(113 23) --from 4326 --to gcj02

火星坐标系 → WGS84：
bash
python3 scripts/coord_convert.py --wkt POINT(113 23) --from gcj02 --to 4326

4. 百度坐标系转换

GCJ-02 → 百度坐标系 (BD-09)：
bash
python3 scripts/coord_convert.py --wkt POINT(113 23) --from gcj02 --to bd09

BD-09 → WGS84：
bash
python3 scripts/coord_convert.py --wkt POINT(113 23) --from bd09 --to 4326

5. 完整流程示例

从百度坐标系的 GeoJSON 文件转换到 CGCS2000 投影坐标系：
bash
python3 scripts/coordconvert.py --input baidudata.geojson --from bd09 --to 4547 --output result.geojson

工作流程

1. 1. 解析输入（判断是 WKT 还是 GeoJSON）
2. 解析源坐标系和目标坐标系
3. 如果需要格式转换，先进行格式转换
4. 使用合适的转换器（pyproj 或中国坐标转换算法）转换坐标
5. 输出转换后的结果

输出格式

• - GeoJSON 输入 → 默认输出 GeoJSON
• WKT 输入 → 默认输出 WKT
• 可用 --format 参数强制指定输出格式

坐标系转换关系图

WGS84 (4326) ←→ GCJ-02 (火星坐标系) ←→ BD-09 (百度坐标系)
↑ ↑ ↑
└────────────────┼────────────────────┘
↓
CGCS2000 (4490)
↓
CGCS2000 分带 (4547-4556)

注意事项

1. 1. 火星坐标系和百度坐标系是中国特有的加密坐标系，无法直接用 pyproj 转换，使用内置算法
2. WGS84 和 CGCS2000 在中国境内差异很小（通常 <1米），实际使用中可互换
3. 投影坐标系（如 4547）是平面坐标，单位为米，适用于工程测量
4. 地理坐标系（如 4326）是经纬度，单位为度