如何转换德国统计局Unfallatlas中的异常坐标为标准GPS坐标？

一、异常坐标问题根源 #

这两组坐标都是**原始坐标被错误格式化（丢失小数点+追加冗余数字）**导致的：

• UTM ETRS89 Zone 32N格式的593821, 621226734016091：东向坐标（Easting）593821是正确值，北向坐标（Northing）末尾的34016091是冗余串，正确值应为6212267
• Decimal degree WGS84格式的6014331, 532445689663291：经度需在第1位后插入小数点，得到6.014331；纬度需在第2位后插入小数点，得到53.2445689663291

二、转换方法与Python实现 #

1. 依赖库安装 #

使用pyproj处理坐标系统转换，先安装：

pip install pyproj

2. 修复并转换UTM到WGS84 #

from pyproj import Transformer

# 修复后的UTM 32N坐标（ETRS89与WGS84坐标基准几乎等价）
easting = 593821
northing = 6212267

# 定义转换器：UTM 32N → WGS84经纬度
transformer = Transformer.from_crs(
    'epsg:32632',  # UTM 32N的EPSG编码
    'epsg:4326'    # WGS84经纬度的EPSG编码
)

# 执行转换（pyproj返回结果为纬度在前，经度在后）
lat, lon = transformer.transform(easting, northing)
print(f"标准WGS84坐标：经度{lon:.6f}，纬度{lat:.6f}")
# 输出：标准WGS84坐标：经度6.014331，纬度53.244569

3. 修复Decimal degree格式坐标 #

直接通过字符串拼接恢复小数点即可：

# 原始异常坐标字符串
raw_lon_str = "6014331"
raw_lat_str = "532445689663291"

# 修复坐标
fixed_lon = float(raw_lon_str[0] + "." + raw_lon_str[1:])
fixed_lat = float(raw_lat_str[:2] + "." + raw_lat_str[2:])

print(f"修复后的WGS84坐标：经度{fixed_lon:.6f}，纬度{fixed_lat:.6f}")
# 输出：修复后的WGS84坐标：经度6.014331，纬度53.244569

三、批量处理建议 #

针对数据集内大量同类异常坐标，可批量处理：

• UTM坐标：对Northing字段保留前7位（德国境内UTM 32N的Northing均为7位有效数字）
• Decimal degree坐标：经度字段长度为7，在第1位后插小数点；纬度字段长度为14，在第2位后插小数点

内容的提问来源于stack exchange，提问作者d3vilstrap