如何将WGS84经纬度坐标转换为ONE模拟器兼容格式？

核心逻辑 #

你手里的是WGS84经纬度（EPSG:4326），ONE模拟器需要的是局部米制平面坐标——本质是把地理经纬度通过投影转换为平面坐标，再做偏移得到相对位置。曼哈顿区域属于UTM 18N投影带（EPSG:32618），这是转换的关键依据。

Python实现步骤 #

1. 安装依赖库 #

先安装处理坐标转换的pyproj库：

pip install pyproj

2. 转换代码实现 #

以下代码可读取WKT文件，将经纬度转换为ONE兼容的平面坐标，再输出新的WKT：

from pyproj import Transformer
import re

def convert_wkt_coords(input_wkt, origin_utm=None):
    # 定义投影转换器：WGS84经纬度 -> UTM 18N
    transformer = Transformer.from_crs("EPSG:4326", "EPSG:32618", always_xy=True)
    
    # 提取WKT中的坐标对
    coord_pairs = re.findall(r'(\d+\.\d+)\s+(-?\d+\.\d+)', input_wkt)
    converted_coords = []
    
    for lat, lon in coord_pairs:
        # 转换为UTM坐标（x:东向米，y:北向米）
        utm_x, utm_y = transformer.transform(float(lon), float(lat))
        
        # 若指定原点，计算相对坐标（匹配ONE示例的数值范围）
        if origin_utm:
            rel_x = utm_x - origin_utm[0]
            rel_y = utm_y - origin_utm[1]
            converted_coords.append(f"{rel_x:.3f} {rel_y:.3f}")
        else:
            converted_coords.append(f"{utm_x:.3f} {utm_y:.3f}")
    
    # 重构WKT字符串
    wkt_type = re.match(r'(\w+)\s+\(', input_wkt).group(1)
    return f"{wkt_type} ({', '.join(converted_coords)})"

# 示例用法
if __name__ == "__main__":
    # 原始WKT内容
    input_wkt = "LINESTRING (40.619605 -73.9584046, 40.6205286 -73.9585793, 40.6206127 -73.9585952)"
    
    # 可选：设置曼哈顿区域的UTM原点（示例用曼哈顿西南角附近的UTM值，可自行替换）
    manhattan_origin_utm = (583700.0, 4498000.0)
    
    converted_wkt = convert_wkt_coords(input_wkt, manhattan_origin_utm)
    print("转换后的WKT：")
    print(converted_wkt)
    
    # 保存到文件
    with open("converted_manhattan.wkt", "w") as f:
        f.write(converted_wkt)

3. 关键说明 #

• 投影带选择：曼哈顿属于UTM 18N带，对应EPSG代码32618，这是转换的核心，不能选错。
• 原点偏移：ONE示例的数值是相对坐标，需找曼哈顿区域内的参考点作为原点，用其UTM坐标做偏移，就能得到和示例量级一致的数值。你可以通过坐标转换工具查询特定点的UTM坐标，替换代码中的manhattan_origin_utm值。
• 坐标顺序：WGS84是纬度在前、经度在后，但pyproj要求输入经度在前、纬度在后，代码中已做对应处理，无需额外调整。

验证方法 #

转换后可对比ONE模拟器的示例坐标范围，微调原点偏移值，直到生成的坐标数值量级匹配（比如几千的范围）。

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