在 OSmnx 中将图形投影到 epsg:4326 的最快方法

Question

我使用 OSmnx 的 osmnx.graph_project 首先将我的图形投影到 UTM 坐标系上，以便使用 consolidate_intersections 函数，然后 ox.project_graph(G_proj_ConsolidateIntersection, to_crs='epsg:4326') 将其投影到纬度-经度坐标系上，以便使用 ox.get_nearest_nodes .最后一个投影需要很长时间（5 分钟），而第一个投影在几秒钟内完成。有没有更好更快的投影到 epsg:4326 的方法？我有looked at this answer 但是当我执行这行gdf = gpd.GeoDataFrame(geometry=intersections) 时，我收到以下错误：

ValueError: Unknown column

这是一个极端的例子（我正在处理的图表较小，但我希望它能说明我的观点）：

import osmnx as ox
G = ox.graph_from_place('Edmonton, Canada', network_type='drive', simplify=False)
# 39.1 s

G_proj = ox.project_graph(G)
# 15.9 s

G_proj_ConsolidateIntersection = ox.consolidate_intersections(G_proj)
# 14 min 45 s

G_proj_ConsolidateIntersection_LatLon = ox.project_graph(G_proj_ConsolidateIntersection, to_crs='epsg:4326')
# 44 min 35 s

gpd.GeoDataFrame(geometry=G_proj_ConsolidateIntersection) 也会引发 Unknown column 错误。

Answer 1

两件事。首先，确保您使用的是最新版本的 OSMnx，因为在之前的版本中有显着的性能增强。其次，如果您先简化图表，这一切都会快得多。

但更重要的是，合并交叉点实际上没有意义除非您已经简化了图表，因为现实世界中的“交叉点”的概念在未简化（扩展）中是没有意义的图形。只有当图被简化时，模型的节点才会映射到现实世界的交叉点和死胡同的概念。有关详细信息，请参阅文档和相关文献。

这是一个比较原始时序（左）和我的简化时序（右）的代码 sn-p：

import osmnx as ox
ox.config(use_cache=True, log_console=True)

G = ox.graph_from_place('Edmonton, Canada', network_type='drive', simplify=True)
# 39.1 s vs 18.2 s

G_proj = ox.project_graph(G)
# 15.9 s vs 8.9 s

G_proj_con = ox.consolidate_intersections(G_proj)
# 14 m 45 s vs 24.2 s

G_con = ox.project_graph(G_proj_con, to_crs='epsg:4326')
# 44 m 35 s vs 6.4 s

编辑：我刚刚重新运行了上面共享的代码 sn-p，但更改了 simplify=False，我使用最新版本的时间仍然比您看到的要快得多。例如，上面的最后两行代码分别在 6m 20s 和 12s 内完成（与您的 14m 45s 和 44m 35s 时间相比）。