Uncomment the following line to install leafmap if needed.

In [ ]:

# !pip install leafmap

In [ ]:

import leafmap

Add a STAC item via an HTTP URL

In [ ]:

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"

In [ ]:

leafmap.stac_assets(url)

In [ ]:

leafmap.stac_bounds(url)

In [ ]:

leafmap.stac_center(url)

In [ ]:

# leafmap.stac_info(url)

In [ ]:

# leafmap.stac_stats(url)

In [ ]:

m = leafmap.Map()
m.add_stac_layer(url, bands=["B3", "B2", "B1"])
m

Add a Microsoft Planetry Computer STAC item. The titiler endpoint can set in one of the ways below:

os.environ["TITILER_ENDPOINT"] = "planetary-computer"
titiler_endpoint="pc"
titiler_endpoint="planetary-computer"

In [ ]:

# import os
# os.environ["TITILER_ENDPOINT"] = "planetary-computer"

In [ ]:

collection = "landsat-8-c2-l2"

In [ ]:

item = "LC08_L2SP_047027_20201204_02_T1"

In [ ]:

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

In [ ]:

leafmap.stac_bounds(collection=collection, item=item)

In [ ]:

leafmap.stac_info(collection=collection, item=item, assets="SR_B7")

In [ ]:

leafmap.stac_stats(collection=collection, item=item, assets="SR_B7")

Color infrared composite.

In [ ]:

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

False color composite.

In [ ]:

m = leafmap.Map()
m.add_stac_layer(
    collection=collection, item=item, assets="SR_B7,SR_B5,SR_B4", name="False color"
)
m

Calculate NDVI.

In [ ]:

m = leafmap.Map()
m.add_stac_layer(
    collection=collection,
    item=item,
    expression="(SR_B5-SR_B4)/(SR_B5+SR_B4)",
    rescale="-1,1",
    name="NDVI",
)
m

Calculate NDVI and add a colormap. See available colormaps at https://planetarycomputer.microsoft.com/docs/reference/data/

In [ ]:

m = leafmap.Map()
m.add_stac_layer(
    collection=collection, item=item, assets="SR_B5,SR_B4,SR_B3", name="Color infrared"
)
m.add_stac_layer(
    collection=collection,
    item=item,
    expression="(SR_B5-SR_B4)/(SR_B5+SR_B4)",
    rescale="-1,1",
    colormap_name="greens",
    name="NDVI Green",
)
m