15 openstreetmap

Downloading OpenStreetMap data with a single line of code

Uncomment the following line to install leafmap if needed.

In [1]:

# !pip install geopandas osmnx leafmap

# !pip install geopandas osmnx leafmap

In [2]:

import leafmap

import leafmap

Add OSM data of place(s) by name or ID to the map. Note that the leafmap custom layer control does not support GeoJSON, we need to use the ipyleaflet built-in layer control.

In [3]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.add_osm_from_geocode("New York City", layer_name='NYC')
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.add_osm_from_geocode("New York City", layer_name='NYC') m

Out[3]:

Make this Notebook Trusted to load map: File -> Trust Notebook

In [4]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.add_osm_from_geocode("Chicago, Illinois", layer_name='Chicago, IL')
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.add_osm_from_geocode("Chicago, Illinois", layer_name='Chicago, IL') m

Out[4]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Add OSM entities within boundaries of geocodable place(s) to the map.

In [5]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
place = "Bunker Hill, Los Angeles, California"
tags = {"building": True}
m.add_osm_from_place(place, tags, layer_name="Los Angeles, CA")
m

m = leafmap.Map(toolbar_control=False, layers_control=True) place = "Bunker Hill, Los Angeles, California" tags = {"building": True} m.add_osm_from_place(place, tags, layer_name="Los Angeles, CA") m

Out[5]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Show OSM feature tags.

https://wiki.openstreetmap.org/wiki/Map_features

In [6]:

# leafmap.osm_tags_list()

# leafmap.osm_tags_list()

Add OSM entities within some distance N, S, E, W of address to the map.

In [7]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.add_osm_from_address(
    address="New York City", tags={"amenity": "bar"}, dist=1500, layer_name="NYC bars"
)
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.add_osm_from_address( address="New York City", tags={"amenity": "bar"}, dist=1500, layer_name="NYC bars" ) m

Out[7]:

Make this Notebook Trusted to load map: File -> Trust Notebook

In [8]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.add_osm_from_address(
    address="New York City",
    tags={"landuse": ["retail", "commercial"], "building": True},
    dist=1000,
    layer_name="NYC buildings",
)
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.add_osm_from_address( address="New York City", tags={"landuse": ["retail", "commercial"], "building": True}, dist=1000, layer_name="NYC buildings", ) m

Out[8]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Add OSM entities within a N, S, E, W bounding box to the map.

In [9]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
north, south, east, west = 40.7551, 40.7454, -73.9738, -73.9965
m.add_osm_from_bbox(
    north, south, east, west, tags={"amenity": "bar"}, layer_name="NYC bars"
)
m

m = leafmap.Map(toolbar_control=False, layers_control=True) north, south, east, west = 40.7551, 40.7454, -73.9738, -73.9965 m.add_osm_from_bbox( north, south, east, west, tags={"amenity": "bar"}, layer_name="NYC bars" ) m

Out[9]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Add OSM entities within some distance N, S, E, W of a point to the map.

In [10]:

m = leafmap.Map(
    center=[46.7808, -96.0156], zoom=12, toolbar_control=False, layers_control=True
)
m.add_osm_from_point(
    center_point=(46.7808, -96.0156),
    tags={"natural": "water"},
    dist=10000,
    layer_name="Lakes",
)
m

m = leafmap.Map( center=[46.7808, -96.0156], zoom=12, toolbar_control=False, layers_control=True ) m.add_osm_from_point( center_point=(46.7808, -96.0156), tags={"natural": "water"}, dist=10000, layer_name="Lakes", ) m

Out[10]:

Make this Notebook Trusted to load map: File -> Trust Notebook

In [11]:

m = leafmap.Map(
    center=[39.9170, 116.3908], zoom=15, toolbar_control=False, layers_control=True
)
m.add_osm_from_point(
    center_point=(39.9170, 116.3908),
    tags={"building": True, "natural": "water"},
    dist=1000,
    layer_name="Beijing",
)
m

m = leafmap.Map( center=[39.9170, 116.3908], zoom=15, toolbar_control=False, layers_control=True ) m.add_osm_from_point( center_point=(39.9170, 116.3908), tags={"building": True, "natural": "water"}, dist=1000, layer_name="Beijing", ) m

/home/runner/.local/lib/python3.9/site-packages/osmnx/geometries.py:805: ShapelyDeprecationWarning: __len__ for multi-part geometries is deprecated and will be removed in Shapely 2.0. Check the length of the `geoms` property instead to get the  number of parts of a multi-part geometry.
  for merged_outer_linestring in list(merged_outer_linestrings):
/home/runner/.local/lib/python3.9/site-packages/osmnx/geometries.py:805: ShapelyDeprecationWarning: Iteration over multi-part geometries is deprecated and will be removed in Shapely 2.0. Use the `geoms` property to access the constituent parts of a multi-part geometry.
  for merged_outer_linestring in list(merged_outer_linestrings):

Out[11]:

Make this Notebook Trusted to load map: File -> Trust Notebook

Add OSM entities within the current map view to the map.

In [12]:

m = leafmap.Map(toolbar_control=False, layers_control=True)
m.set_center(-73.9854, 40.7500, 16)
m

m = leafmap.Map(toolbar_control=False, layers_control=True) m.set_center(-73.9854, 40.7500, 16) m

Out[12]:

Make this Notebook Trusted to load map: File -> Trust Notebook

The add_osm_from_view() is only available for the ipyleaflet plotting backend.

In [13]:

# m.add_osm_from_view(tags={"amenity": "bar", "building": True}, layer_name="New York")

# m.add_osm_from_view(tags={"amenity": "bar", "building": True}, layer_name="New York")

Create a GeoPandas GeoDataFrame from place.

In [14]:

gdf = leafmap.osm_gdf_from_place("New York City", tags={"amenity": "bar"})
gdf

gdf = leafmap.osm_gdf_from_place("New York City", tags={"amenity": "bar"}) gdf

/home/runner/.local/lib/python3.9/site-packages/osmnx/geocoder.py:110: FutureWarning: The frame.append method is deprecated and will be removed from pandas in a future version. Use pandas.concat instead.
  gdf = gdf.append(_geocode_query_to_gdf(q, wr, by_osmid))

Out[14]:

		addr:city	addr:country	addr:housenumber	addr:postcode	addr:state	addr:street	amenity	ele	gnis:county_name	gnis:feature_id	...	disused:amenity	screen	drink:cocktails	building:colour	building:levels	roof:material	roof:shape	roof:colour	ways	type
element_type	osmid
node	368043598	New York	US	53	10014	NY	Christopher Street	bar	7	New York	2082946	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	419366609	NaN	NaN	NaN	NaN	NaN	NaN	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	419367969	NaN	NaN	NaN	NaN	NaN	NaN	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	484078905	Brooklyn	NaN	388	11211	NY	Union Avenue	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	617740843	NaN	NaN	NaN	NaN	NaN	NaN	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...	...
way	620691053	NaN	NaN	480	11231	NaN	Union Street	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	black	NaN	NaN	NaN	#999999	NaN	NaN
	637354701	NaN	NaN	308	11238	NaN	Flatbush Avenue	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	664632831	NaN	NaN	1840	10128	NaN	2nd Avenue	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
	665897391	New York	NaN	54	10001	NaN	West 31st Street	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN
relation	12554650	NaN	NaN	NaN	NaN	NaN	NaN	bar	NaN	NaN	NaN	...	NaN	NaN	NaN	NaN	NaN	NaN	NaN	NaN	[278963046, 278963089]	multipolygon

1027 rows × 169 columns

---

Last update: 2022-03-14