Use OSMnx to plot street network over place shape¶

Author: Geoff Boeing

This example uses Portland, Maine - a city with several islands within its municipal boundaries. Thus, we set retain_all=True when getting the network so that we keep all the graph components, not just the largest connected component.

In [ ]:

import matplotlib.pyplot as plt
import osmnx as ox

%matplotlib inline
ox.__version__

In [ ]:

# get the place boundaries
place = "Portland, Maine"
gdf = ox.geocode_to_gdf(place)

# get the street network, with retain_all=True to retain all the disconnected islands' networks
G = ox.graph_from_place(place, network_type="drive", retain_all=True)

In [ ]:

# plot the network, but do not show it or close it yet
fig, ax = ox.plot_graph(
    G,
    show=False,
    close=False,
    bgcolor="#333333",
    edge_color="w",
    edge_linewidth=0.3,
    node_size=0,
)

# to this matplotlib axis, add the place shape(s)
gdf.plot(ax=ax, fc="k", ec="#666666", lw=1, alpha=1, zorder=-1)

# optionally set up the axes extents
margin = 0.02
west, south, east, north = gdf.unary_union.bounds
margin_ns = (north - south) * margin
margin_ew = (east - west) * margin
ax.set_ylim((south - margin_ns, north + margin_ns))
ax.set_xlim((west - margin_ew, east + margin_ew))
plt.show()

Notice this municipal boundary is an administrative boundary, not a physical boundary, so it represents jurisdictional bounds, not individual physical features like islands. To retrieve individual islands' geometries, use the features module to search for features matching certain OSM tags:

In [ ]:

islands = ox.features_from_place(place, tags={"place": ["island", "islet"]})
islands.shape

In [ ]: