OSMnx nearest_nodes获取节点距GPS点过远问题排查

我尝试用OSMnx查找两地间的最短路径，调用shortest_path前需要获取GPS点附近的节点，但nearest_nodes函数返回的节点距离我的GPS坐标约5000米，远超预期。运行代码后得到的节点坐标为（33.8462657, -118.2883998），与原GPS点（33.8254852, -118.240463）差距明显。我已尝试修改图的network_type和simplify参数，但问题仍存在。

复现代码 #

import networkx as nx
import osmnx as ox

graph_area = 'Los Angeles, California, United States'
mode = 'all'

#G = ox.graph_from_place(graph_area, network_type=mode,simplify=False)
#G = ox.add_edge_speeds(G)
#G = ox.add_edge_travel_times(G)
#print("a")

#ox.save_graphml(G, "LA2.graphml") #
    
G = ox.load_graphml("LA2.graphml")

start_lat = 33.8254852
start_lon = -118.240463

start_node = ox.nearest_nodes(G, start_lon, start_lat, return_dist=True)

print(start_node)

print('lat:' + str(G.nodes[start_node[0]]['y']))
print('lon:' + str(G.nodes[start_node[0]]['x']))

原因分析 #

• 核心问题：你加载的LA2.graphml路网文件未覆盖目标GPS点。从坐标对比看，目标点位于洛杉矶东部，但返回的节点在更靠西的区域，说明当初保存的路网数据本身就不包含目标点所在的区域。
• 诱因：当初调用ox.graph_from_place时，'Los Angeles, California, United States'这个地名解析出的地理范围可能和你预期有偏差，OSMnx根据该地名获取的路网边界没有覆盖到目标GPS点；或者保存graphml时，路网数据生成不完整。

解决办法 #

• 重新生成覆盖目标点的路网
  取消代码中注释的生成路网部分，改用GPS点缓冲区来定义路网范围，确保覆盖目标点位：

  # 以目标GPS点为中心，生成10公里范围内的路网
  start_lat = 33.8254852
  start_lon = -118.240463
  G = ox.graph_from_point((start_lat, start_lon), dist=10000, network_type='all')
  G = ox.add_edge_speeds(G)
  G = ox.add_edge_travel_times(G)
  ox.save_graphml(G, "LA2.graphml")
  

  这样生成的路网会直接包含目标点周边的道路节点，再调用nearest_nodes就能得到距离最近的节点。

• 验证现有路网的覆盖范围
  先检查当前加载的路网是否包含目标点：

  # 获取路网的边界坐标（南纬、北纬、西经、东经）
  bbox = ox.utils_graph.get_bbox(G)
  print(f"路网边界：南纬{bbox[0]}, 北纬{bbox[1]}, 西经{bbox[2]}, 东经{bbox[3]}")
  # 判断目标点是否在路网范围内
  is_in_range = (bbox[0] <= start_lat <= bbox[1]) and (bbox[2] <= start_lon <= bbox[3])
  print(f"目标点是否在路网内：{is_in_range}")
  

  如果输出False，说明必须重新生成覆盖该点的路网。

• 确认参数顺序正确性
  注意ox.nearest_nodes的参数顺序是先经度（lon）、后纬度（lat），你的代码中ox.nearest_nodes(G, start_lon, start_lat)的写法是正确的，这部分无需调整，但如果参数写反会导致节点偏差，可作为排查点。

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