如何在 osmnx 中将建筑物添加为图中的节点？

Question

我对 osmnx 不是很熟悉。我正在尝试创建包含节点的图 在每栋建筑的中心以及街道节点。以下代码将建筑物（不是节点）和街道节点一起显示：

import networkx as nx
import osmnx as ox
from matplotlib import *
import matplotlib.pyplot as plt
from descartes import PolygonPatch
from shapely.geometry import MultiPolygon
from shapely.geometry import Polygon
from IPython.display import Image
import geopandas as gpd
from IPython.display import IFrame


place = "Piedmont, California, US"
G = ox.graph_from_place(place)

tags = {'building': True}
buildings = ox.geometries_from_place(place, tags)



# or plot street network and the entities' footprints together
fig, ax = ox.plot_footprints(buildings, alpha=0.4, show=False)
fig, ax = ox.plot_graph(G, ax=ax, node_size=10, edge_color="w", edge_linewidth=0.7)

Answer 1

我可以通过提取每座建筑物的质心坐标来做到这一点 并在那里添加新节点：

import networkx as nx
import osmnx as ox
from matplotlib import *
import matplotlib.pyplot as plt
from descartes import PolygonPatch
from shapely.geometry import MultiPolygon
from shapely.geometry import Polygon
from shapely.geometry import Point
from IPython.display import Image
import geopandas as gpd
from IPython.display import IFrame

# Specify the name that is used to seach for the data
place_name = "kirchberg, luxembourg, luxembourg"

# Fetch OSM street network from the location
graph = ox.graph_from_place(place_name)

nodes= ox.graph_to_gdfs(graph, nodes=True, edges=False)
edges= ox.graph_to_gdfs(graph, edges=True, nodes=False)
#ox.plot_graph(graph)


for i in range(len(list(graph.nodes))):
    node_id = list(graph.nodes)[i]
#    print(f"lon: {graph.nodes[node_id]['x']}, lat: {graph.nodes[node_id]['y']}")



#print(graph.nodes)


buildings = ox.geometries_from_place(place_name, tags={
'building':True
}) # Retrieve buildings from the area:


#buildings = buildings.to_crs(4326)


centroids = buildings.centroid

#print(dir(centroids))

lons = list(centroids.x)

lats = list(centroids.y)


building_coords = [[lons[i],lats[i]] for i in range(len(lons))]

print(building_coords)

dictionary = dict()

for i in range(len(building_coords)):
     node_id = i+1
     dictionary[node_id] = {
    'x': lons[i],
    'y': lats[i]
     }

tmp_list = []
for item_key, item_value in dictionary.items() :
  tmp_list.append({
    'geometry' : Point(item_value['x'], item_value['y']),
    'osmid': item_key,
    'y' : item_value['y'],
      'x' : item_value['x'],
   })
my_nodes = gpd.GeoDataFrame(tmp_list)

# just plot the building nodes
nodes = my_nodes

# or plot building nodes and street nodes
#nodes = nodes.append(my_nodes, ignore_index = True)

graph2 = ox.graph_from_gdfs(nodes, edges)
#fig, ax = ox.plot_footprints(buildings, alpha=0.4, show=False)
#ox.plot_graph(graph2)
ox.plot_graph(graph2, node_size=10, edge_color="grey", node_color="green", edge_linewidth=0.7, bgcolor="white")

感谢 Debjit Bhowmick 的代码在此答案中。

Answer 2

这些节点然后用边连接到图吗？