Simple Image Visualization

Visualize a true color Landsat 8 image.

from descarteslabs.catalog import Product, properties as p
from descarteslabs.geo import DLTile
from descarteslabs.utils import display

Create a tile around Pisa, Italy.

tile = DLTile.from_latlon(43.7230, 10.3966, resolution=20.0, tilesize=1024, pad=0)

Use the Catalog API to search for imagery available over the area of interest.

search = (
    Product.get("usgs:landsat:oli-tirs:c2:l2:v0")
    .images()
    .intersects(tile)
    .filter("2022-04-01" <= p.acquired < "2022-05-01")
    .filter(p.cloud_fraction < 0.001)
    .sort("acquired")
    .limit(1)
)
images = search.collect()

Pick just one image to raster and display.

image = images[0]

# Load the data as an ndarray
arr = image.ndarray("red green blue", geocontext=images.geocontext)

# Display the image
display(arr, size=5, title=image.id)
usgs:landsat:oli-tirs:c2:l2:v0:LC08_L2SP_193030_20220429_20220504_02_T1

Total running time of the script: ( 0 minutes 2.553 seconds)

Gallery generated by Sphinx-Gallery