gpt4 book ai didi

python - OSMnx 获取干净交叉点节点的经纬度坐标

转载 作者:行者123 更新时间:2023-12-01 09:19:31 27 4
gpt4 key购买 nike

我正在使用 OSMnx 从 OpenStreetMaps 道路网络获取干净交叉路口。交叉点当前位于 (x,y) 坐标,但我想使用经纬度坐标绘制它们。

从 Jupiter 笔记本示例中,OSMnx Example #14 Clean Intersection Cluster Nodes ,我能够获取街道网络并调用 ox.clean_intersections 来生成干净的交叉路口。

import osmnx as ox, matplotlib.pyplot as plt, numpy as np
ox.config(use_cache=True, log_console=True)
%matplotlib inline

# get a street network and plot it with all edge intersections
address = '2700 Shattuck Ave, Berkeley, CA'
G = ox.graph_from_address(address, network_type='drive', distance=750)
G_proj = ox.project_graph(G)

# clean up the intersections and extract their xy coords
intersections = ox.clean_intersections(G_proj, tolerance=15, dead_ends=False)
points = np.array([point.xy for point in intersections])

我得到了交叉点的 Pandas 地理系列,如下所示:

0       POINT (564152.437121744 4189596.945341664)
1 POINT (564846.6779513165 4189615.534235776)
2 POINT (564571.2116373706 4189601.780093061)

由于干净的交叉点由 centroids of clusters of merged nodes 组成,它们不对应于具有 osm_id 的任何特定节点(具有经纬度坐标)。

如何将这些 (x,y) 点转换为经纬度坐标?

最佳答案

您将图表投影到米上,以使用合理的公差参数清理交叉点。现在您只需将清理后的交点质心投影回经纬度即可:

import geopandas as gpd
gdf = gpd.GeoDataFrame(geometry=intersections)
gdf.crs = G_proj.graph['crs']
ox.project_gdf(gdf, to_latlong=True)

关于python - OSMnx 获取干净交叉点节点的经纬度坐标,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/50913798/

27 4 0
Copyright 2021 - 2024 cfsdn All Rights Reserved 蜀ICP备2022000587号
广告合作:1813099741@qq.com 6ren.com