Visualize SpatioTemporal Asset Catalog (STAC)

This notebook demonstrates how to visualize SpatioTemporal Asset Catalog (STAC) items using the TiTiler demo endpoint (titiler.xyz). Please be gentle with the demo endpoint.

Uncomment the following line to install leafmap if needed.

In [1]:

# %pip install "leafmap[maplibre]"

# %pip install "leafmap[maplibre]"

In [2]:

import leafmap.maplibregl as leafmap

import leafmap.maplibregl as leafmap

To run this notebook, you will need an API key from MapTiler. Once you have the API key, you can uncomment the following code block and replace YOUR_API_KEY with your actual API key. Then, run the code block code to set the API key as an environment variable.

In [3]:

# import os
# os.environ["MAPTILER_KEY"] = "YOUR_API_KEY"

# import os # os.environ["MAPTILER_KEY"] = "YOUR_API_KEY"

In [4]:

m = leafmap.Map(style="streets")
url = "https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/S5_11055_6057_20070622.json"
m.add_stac_layer(url, bands=["pan"], name="Panchromatic", vmin=0, vmax=150)
m.add_stac_layer(url, bands=["B4", "B3", "B2"], name="RGB", vmin=0, vmax=150)
m

m = leafmap.Map(style="streets") url = "https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/S5_11055_6057_20070622.json" m.add_stac_layer(url, bands=["pan"], name="Panchromatic", vmin=0, vmax=150) m.add_stac_layer(url, bands=["B4", "B3", "B2"], name="RGB", vmin=0, vmax=150) m

Failed to retrieve the MapTiler style. Defaulting to OpenFreeMap 'liberty' style.

WARNING:rasterio._env:CPLE_AppDefined in HTTP response code on https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/s5_11055_6057_20070622_p10_1_lcc00_cog.tif.aux.xml: 403

WARNING:rasterio._env:CPLE_AppDefined in HTTP response code on https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/s5_11055_6057_20070622_p10_1_lcc00_cog.aux: 403

WARNING:rasterio._env:CPLE_AppDefined in HTTP response code on https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/s5_11055_6057_20070622_p10_1_lcc00_cog.AUX: 403

WARNING:rasterio._env:CPLE_AppDefined in HTTP response code on https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/s5_11055_6057_20070622_p10_1_lcc00_cog.tif.aux: 403

WARNING:rasterio._env:CPLE_AppDefined in HTTP response code on https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/s5_11055_6057_20070622_p10_1_lcc00_cog.tif.AUX: 403

WARNING:rasterio._env:CPLE_AppDefined in HTTP response code on https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/s5_11055_6057_20070622_p10_1_lcc00_cog.xml: 403

WARNING:rasterio._env:CPLE_AppDefined in HTTP response code on https://canada-spot-ortho.s3.amazonaws.com/canada_spot_orthoimages/canada_spot5_orthoimages/S5_2007/S5_11055_6057_20070622/s5_11055_6057_20070622_p10_1_lcc00_cog.XML: 403

In [5]:

m.layer_interact()

m.layer_interact()

In [6]:

collection = "landsat-8-c2-l2"
item = "LC08_L2SP_047027_20201204_02_T1"

collection = "landsat-8-c2-l2" item = "LC08_L2SP_047027_20201204_02_T1"

In [7]:

leafmap.stac_assets(collection=collection, item=item, titiler_endpoint="pc")

leafmap.stac_assets(collection=collection, item=item, titiler_endpoint="pc")

Out[7]:

['SR_B1',
 'SR_B2',
 'SR_B3',
 'SR_B4',
 'SR_B5',
 'SR_B6',
 'SR_B7',
 'ST_QA',
 'ST_B10',
 'ST_DRAD',
 'ST_EMIS',
 'ST_EMSD',
 'ST_TRAD',
 'ST_URAD',
 'QA_PIXEL',
 'ST_ATRAN',
 'ST_CDIST',
 'QA_RADSAT',
 'SR_QA_AEROSOL']

In [8]:

m = leafmap.Map(style="streets")
m.add_stac_layer(
    collection=collection,
    item=item,
    assets=["SR_B5", "SR_B4", "SR_B3"],
    name="Color infrared",
)
m

m = leafmap.Map(style="streets") m.add_stac_layer( collection=collection, item=item, assets=["SR_B5", "SR_B4", "SR_B3"], name="Color infrared", ) m

Failed to retrieve the MapTiler style. Defaulting to OpenFreeMap 'liberty' style.