68 openaerialmap

Searching and visualizing OpenAerialMap imagery interactively

OpenAerialMap (OAM) provides openly licensed satellite and unmanned aerial vehicle (UAV) imagery. Currently, it has over 12,400+ images around the globe. This notebook demonstrates how to search and visualize OAM imagery interactively. You can download images automatically with one line of code.

Uncomment the following line to install leafmap if needed.

In [1]:

import leafmap

import leafmap

Search OAM images by a bounding box and a date range. The results can be returned as a GeoDataFrame or a list of image metadata.

In [2]:

bbox = [-79.6344, -0.9063, -77.3383, 1.0436]
start_date = '2016-04-01'
end_date = '2016-04-30'
gdf = leafmap.oam_search(
    bbox=bbox, start_date=start_date, end_date=end_date, return_gdf=True
)
print(f'Found {len(gdf)} images')

bbox = [-79.6344, -0.9063, -77.3383, 1.0436] start_date = '2016-04-01' end_date = '2016-04-30' gdf = leafmap.oam_search( bbox=bbox, start_date=start_date, end_date=end_date, return_gdf=True ) print(f'Found {len(gdf)} images')

Found 55 images

In [3]:

gdf.head()

gdf.head()

Out[3]:

	geometry	_id	uuid	__v	title	gsd	file_size	acquisition_start	acquisition_end	platform	provider	contact	sensor	thumbnail	tms	wmts	uploaded_at	meta_uri	bbox
0	POLYGON ((-78.31663 -0.26595, -78.31663 -0.080...	59e62b7e3d6412ef722096cf	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	0	Planet Labs PlanetScope Ecuador Earthquake	3	3018807	2016-04-20T00:00:00.000Z	2016-04-20T23:30:00.000Z	satellite	Planet Labs	Erika Reinhardt,erika@planet.com	PlanetScope	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	2016-04-25T00:00:00.000Z	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	[-78.3166259908633, -0.2659504218286762, -78.1...
1	POLYGON ((-78.48063 -0.50211, -78.48063 -0.318...	59e62b7e3d6412ef722096f4	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	0	Planet Labs PlanetScope Ecuador Earthquake	3	3873218	2016-04-20T00:00:00.000Z	2016-04-20T23:30:00.000Z	satellite	Planet Labs	Erika Reinhardt,erika@planet.com	PlanetScope	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	2016-04-25T00:00:00.000Z	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	[-78.48063444708373, -0.502107169583732, -78.2...
2	POLYGON ((-78.53432 -0.57968, -78.53432 -0.395...	59e62b7e3d6412ef722096ea	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	0	Planet Labs PlanetScope Ecuador Earthquake	3	3390335	2016-04-20T00:00:00.000Z	2016-04-20T23:30:00.000Z	satellite	Planet Labs	Erika Reinhardt,erika@planet.com	PlanetScope	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	2016-04-25T00:00:00.000Z	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	[-78.53432214477868, -0.5796813226530548, -78....
3	POLYGON ((-78.25983 -0.18810, -78.25983 -0.003...	59e62b7d3d6412ef7220969c	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	0	Planet Labs PlanetScope Ecuador Earthquake	3	3112239	2016-04-20T00:00:00.000Z	2016-04-20T23:30:00.000Z	satellite	Planet Labs	Erika Reinhardt,erika@planet.com	PlanetScope	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	2016-04-25T00:00:00.000Z	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	[-78.25983358913794, -0.18810305128235869, -78...
4	POLYGON ((-78.42670 -0.42386, -78.42670 -0.239...	59e62b7e3d6412ef722096e7	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	0	Planet Labs PlanetScope Ecuador Earthquake	3	3735755	2016-04-20T00:00:00.000Z	2016-04-20T23:30:00.000Z	satellite	Planet Labs	Erika Reinhardt,erika@planet.com	PlanetScope	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	http://tiles.openaerialmap.org/571e5c7fcd0663b...	2016-04-25T00:00:00.000Z	http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0...	[-78.42670295517972, -0.42386105287234965, -78...

The tile URLs are stored in the tms column of the GeoDataFrame.

In [4]:

tiles = gdf['tms'].tolist()
tiles[:5]

tiles = gdf['tms'].tolist() tiles[:5]

Out[4]:

['http://tiles.openaerialmap.org/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb35/{z}/{x}/{y}.png',
 'http://tiles.openaerialmap.org/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb41/{z}/{x}/{y}.png',
 'http://tiles.openaerialmap.org/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb3f/{z}/{x}/{y}.png',
 'http://tiles.openaerialmap.org/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb2c/{z}/{x}/{y}.png',
 'http://tiles.openaerialmap.org/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb40/{z}/{x}/{y}.png']

The image sources (downloadable URLs) are stored in the uuid column of the GeoDataFrame.

In [5]:

images = gdf['uuid'].tolist()
images[:5]

images = gdf['uuid'].tolist() images[:5]

Out[5]:

['http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb35.tif',
 'http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb41.tif',
 'http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb3f.tif',
 'http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb2c.tif',
 'http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb40.tif']

Download all images using the download_files() function.

In [6]:

leafmap.download_files(images[:2])

leafmap.download_files(images[:2])

Downloading...
From: http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb35.tif
To: /home/runner/work/leafmap/leafmap/docs/notebooks/571e60322b67227a79b4fb35.tif
100%|██████████| 3.02M/3.02M [00:00<00:00, 12.7MB/s]
Downloading...
From: http://oin-hotosm.s3.amazonaws.com/571e5c7fcd0663bb003c31f9/0/571e60322b67227a79b4fb41.tif
To: /home/runner/work/leafmap/leafmap/docs/notebooks/571e60322b67227a79b4fb41.tif
100%|██████████| 3.87M/3.87M [00:00<00:00, 12.4MB/s]

Add the image footprints to the map.

In [7]:

m = leafmap.Map()
m.add_gdf(gdf, layer_name='Footprints')
m

m = leafmap.Map() m.add_gdf(gdf, layer_name='Footprints') m

Out[7]:

Search OAM imagery within the current map extent or user drawn ROI.

In [8]:

m = leafmap.Map(center=[4.7955, -75.6899], zoom=15)
m.add_basemap('SATELLITE', show=False)

bbox = [-75.7138, 4.7826, -75.6659, 4.808332]
m.oam_search(bbox=bbox)

m

m = leafmap.Map(center=[4.7955, -75.6899], zoom=15) m.add_basemap('SATELLITE', show=False) bbox = [-75.7138, 4.7826, -75.6659, 4.808332] m.oam_search(bbox=bbox) m

Out[8]:

Search OAM imagery interactively using the interactive GUI.

In [9]:

m = leafmap.Map(center=[4.7955, -75.6899], zoom=15)
m

m = leafmap.Map(center=[4.7955, -75.6899], zoom=15) m

Out[9]:

---

Last update: 2023-03-04
Created: 2021-06-12