尝试使用OSMnx

Question

我正在尝试使用OSMnx在同一图中绘制多个基础设施网络（例如街道，铁路和建筑物），但实际上没有成功。

这是我的一次尝试：

import osmnx as ox
dist = 2000
point = (41.877092, -87.628)
north, south, east, west = ox.bbox_from_point(point, distance=dist)
bbox_proj = ox.bbox_from_point(point, dist, project_utm=True)
streets = ox.core.osm_net_download(
    north=north,
    south=south,
    east=east,
    west=west,
    infrastructure='way["highway"]'
    )
railways = ox.core.osm_net_download(
    north=north,
    south=south,
    east=east,
    west=west,
    infrastructure='way["railway"]'
    )
buildings = ox.core.osm_net_download(
    north=north,
    south=south,
    east=east,
    west=west,
    infrastructure='way["building"]'
    )
streets[0]['elements'] = streets[0]['elements'] + railways[0]['elements'] + buildings[0]['elements']
net = streets
G = ox.core.create_graph(net)
G = ox.truncate_graph_bbox(G, north, south, east, west, truncate_by_edge=True)
G = ox.project_graph(G)
_, _ = ox.plot.plot_graph(G, bbox=bbox_proj, fig_height=10, node_size=0, edge_color='black', edge_linewidth=0.5, save=True)

此代码的结果仅绘制了前2个基础设施，街道和轨道，而不是建筑物：

ox.plot_figure_ground的结果（仅绘制街道基础设施）：

尽管没有绘制，但正在下载建筑物数据（我知道plot_buildings，但我不想要彩色建筑，只需要行）。

在此之前，我试图找到一种方法在infrastructure参数中一次添加多个过滤器。类似的东西：

nets = ox.core.osm_net_download(
    north=north,
    south=south,
    east=east,
    west=west,
    infrastructure='way["highway"],way["railway"],way["buildings"]'
    )

但不确定这是否可行。

任何方式可以在一个图中绘制两个以上的图形，如果可以的话，在OSMnx中简洁地绘制它？

Answer 1

尝试将import { p } from './logging'; 替换为console.log。正确的OSM密钥是building，最后没有“s”。事实上，大多数OSM键都是单数。

Answer 2

是的，这是可能的，并且这样做会很简洁：

import matplotlib.pyplot as plt
import networkx as nx
import osmnx as ox
ox.config(use_cache=True, log_console=True)

# get graphs of different infrastructure types, then combine
place = 'Berkeley, California, USA'
G1 = ox.graph_from_place(place, custom_filter='["highway"~"residential"]')
G2 = ox.graph_from_place(place, custom_filter='["railway"~"rail"]')
G = nx.compose(G1, G2)

# get building footprints
fp = ox.footprints_from_place(place)

# plot highway edges in yellow, railway edges in red
ec = ['y' if 'highway' in d else 'r' for _, _, _, d in G.edges(keys=True, data=True)]
fig, ax = ox.plot_graph(G, bgcolor='k', edge_color=ec,
                        node_size=0, edge_linewidth=0.5,
                        show=False, close=False)

# add building footprints in 50% opacity white
fp.plot(ax=ax, color='w', alpha=0.5)
plt.show()

另请参阅https://stackoverflow.com/a/62883614/7321942和https://stackoverflow.com/a/62720802/7321942