Access Sentinel-2 assets

from shapely import box
from tilebox.datasets import Client, field

datasets = Client()
collection = datasets.dataset("open_data.aws_earth.sentinel2").collection("L2A")

scenes = collection.query(
    temporal_extent=("2025-10-01", "2025-11-01"),
    spatial_extent=box(-106.0, 38.0, -105.9, 38.1),
    filter=field("cloud_cover") < 10,
)

datapoint = scenes.isel(time=0)
print(datapoint.stac_id.item())

from tilebox.datasets.assets import AssetCollection

assets = AssetCollection.from_datapoint(datapoint)

for key, asset in assets.items():
    print(key, asset.media_type)

red = assets["red"]

import asyncio

from tilebox.storage.aio import Client
from tilebox.storage.geotiff import window_from_bounds

async def read_area():
    storage = Client()
    geotiff = await storage.open_geotiff(red)
    window = window_from_bounds(
        geotiff,
        (-106.0, 38.0, -105.9, 38.1),
        crs="EPSG:4326",
    )
    return await geotiff.read(window=window)

pixels = asyncio.run(read_area())
print(pixels.shape)

async def download_red_band():
    storage = Client()
    return await storage.download(red, "data/sentinel-2-red.tif")

path = asyncio.run(download_red_band())
print(path)