将静脉坐标转换为 GPS

我正在使用从以下位置导入的真实街道网络打开街道地图用于静脉模拟，例如卢森堡情景来自拉拉·科德卡。现在，为了准备可视化（使用 Google Earth），我想将模拟中的车辆位置从 SUMO 或 OmNET 坐标导出到 GPS 坐标。

作为材料，我有用于生成场景的 OSM 文件，包括所有节点的 GPS 位置。我希望找到从模拟坐标到 GPS 坐标的简单映射，例如，通过了解边界框和模拟游乐场的角点的 GPS 坐标。

有没有一种简单的方法来进行这种转换，以及如何找到生成操场时 OSM 转换使用的实际角点？

转换过程如下：

1. 从 OmNET 访问位置信息

// Adapt your path to the mobility module here  
Veins::TraCIMobility* mobility =
  check_and_cast<Veins::TraCIMobility*>(
    getParentModule()->getSubmodule("veinsmobility"));

Veins::TraCICommandInterface* traci = mobility->getCommandInterface();

Coord currPos = mobility->getCurrentPosition();
std::pair<double, double> currLonLat = traci->getLonLat(currPos);

getLonLat()为我返回绝对二维坐标，因此需要一个转换步骤。

2. 寻找转变

SUMO 中的 .net.xml 文件包含所需的转换。这<location>标签包含属性netOffset and projParameters需要的。

对于卢森堡情景，这些是

netOffset="-285448.66,-5492398.13"
projParameter="+proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"

3. 反转转型

图书馆PROJ.4可以用来做反演。还提供 Python 接口（pyproj).

import pyproj

projection = pyproj.Proj(
  "+proj=utm +zone=32 +ellps=WGS84 +datum=WGS84 +units=m +no_defs")

# x, y obtained from OmNET
lon, lat = projection(x, y, inverse=True)

如果只有相对位置信息可用，则必须首先通过添加 netOffset 值来调整 x、y 值。

Edit

构建 SUMO 时只需要第一步--with-proj-gdal支持，结果getLonLat()将立即成为所需的格式。