GIS开发：如何将LineString道路网络转为NetworkX完全图？

你不需要强行把整个网络转成完全图（会生成O(n²)条边，节点多直接炸内存），以下是更合理的技术方案：

方案1：先修复路网连通性（优先推荐） #

NetworkXNoPath大多是因为原始路网本身断开（数据缺失、道路未衔接），先补全连通性比转完全图靠谱：

• 用momepy的close_gaps自动连接近距离的断开路段：

  import momepy
  # 假设gdf是你的LineString类型GeoDataFrame
  connected_gdf = momepy.close_gaps(gdf, threshold=10)  # threshold为允许的最大间隙（单位：米）
  # 重新转为带距离属性的NetworkX图
  G = momepy.gdf_to_nx(connected_gdf, approach='primal', length='length')
  

• 重新定位A、B的最近节点后，再调用nx.shortest_path(G, source=node_A, target=node_B, weight='length')即可。

方案2：仅为目标节点添加跨连通分量的边（按需生成） #

如果路网确实存在无法连通的区域，只给A、B所在的连通分量加一条直线距离边即可，不用全图改造：

1. 先确定A、B所在的连通分量：

   import networkx as nx
   components = list(nx.connected_components(G))
   comp_A = next(c for c in components if node_A in c)
   comp_B = next(c for c in components if node_B in c)
   

2. 若两个分量不同，计算节点地理距离并加边：

   from shapely.geometry import Point
   
   # 从节点属性中提取坐标（假设节点存了x、y属性）
   point_A = Point(G.nodes[node_A]['x'], G.nodes[node_A]['y'])
   point_B = Point(G.nodes[node_B]['x'], G.nodes[node_B]['y'])
   straight_dist = point_A.distance(point_B)
   
   # 添加跨分量边，权重设为直线距离（可乘以系数模拟绕行成本）
   G.add_edge(node_A, node_B, length=straight_dist)
   

3. 此时调用最短路径会自动使用这条跨分量边，同时保留原有路网的地理距离权重。

方案3：真转完全图（仅小数据集可用） #

如果节点数量极少（几十上百个），可以强行生成完全图：

nodes = list(G.nodes(data=True))
for i in range(len(nodes)):
    node_i, attrs_i = nodes[i]
    point_i = Point(attrs_i['x'], attrs_i['y'])
    for j in range(i+1, len(nodes)):
        node_j, attrs_j = nodes[j]
        point_j = Point(attrs_j['x'], attrs_j['y'])
        dist = point_i.distance(point_j)
        G.add_edge(node_i, node_j, length=dist)

注意：节点超过1000个时，边数会突破50万，内存和计算成本会急剧上升，不推荐使用。

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