dvcrn GIS 坐标转换 Skill
快速转换点、线、面数据的坐标系：WGS84（经纬度）↔ Web Mercator（EPSG:3857）
快速开始
转换 GeoJSON 到墨卡托
python3 ~/.openclaw/workspace/skills/gis-coord-transform/scripts/transform_coords.py input.geojson -t mercator -o output.geojson --pretty
转换墨卡托回经纬度
python3 ~/.openclaw/workspace/skills/gis-coord-transform/scripts/transform_coords.py input.geojson -t wgs84 -o output.geojson --pretty
支持的几何类型
| 类型 | 说明 |
|---|---|
| Point | 点 |
| LineString | 线 |
| Polygon | 面 |
| MultiPoint | 多点 |
| MultiLineString | 多线 |
| MultiPolygon | 多面 |
| Feature | GeoJSON 要素 |
| FeatureCollection | 要素集合 |
命令行参数
| 参数 | 说明 |
|---|---|
input |
输入 GeoJSON 文件路径 |
-o, --output |
输出文件路径（默认输出到 stdout） |
-t, --to |
目标坐标系：mercator 或 wgs84 |
--pretty |
美化输出 JSON |
使用示例
示例 1：转换单个点
// input.geojson
{
"type": "Point",
"coordinates": [116.4074, 39.9042]
}
python3 scripts/transform_coords.py input.geojson -t mercator
输出（墨卡托坐标）：
{
"type": "Point",
"coordinates": [12958396.7, 4865942.3]
}
示例 2：转换 FeatureCollection
python3 scripts/transform_coords.py cities.geojson -t mercator -o cities_mercator.geojson --pretty
示例 3：批量转换
for file in *.geojson; do
python3 scripts/transform_coords.py "$file" -t mercator -o "${file%.geojson}_mercator.geojson"
done
依赖安装
pip3 install pyproj
坐标系说明
| 坐标系 | EPSG 代码 | 用途 |
|---|---|---|
| WGS84 | EPSG:4326 | GPS、经纬度、地理坐标 |
| Web Mercator | EPSG:3857 | Web 地图（Google Maps、OpenStreetMap、Mapbox 等） |
常见问题
为什么转换后坐标数字变大了？
Web Mercator 使用米为单位，而 WGS84 使用度。墨卡托坐标通常在百万级别。
转换精度如何？
使用 pyproj 库，精度非常高，适合专业 GIS 应用。
支持其他坐标系吗？
当前版本仅支持 WGS84 ↔ Web Mercator。如需其他坐标系，可扩展 pyproj.Transformer。
扩展建议
需要添加其他坐标系转换时，修改脚本中的 Transformer 定义：
# 例如：WGS84 → CGCS2000 (EPSG:4490)
transformer = Transformer.from_crs("EPSG:4326", "EPSG:4490", always_xy=True)
相关文件
scripts/transform_coords.py- 坐标转换主脚本