Access Copernicus data

from shapely import Polygon
from tilebox.datasets import Client

area = Polygon(
    [
        (-109.05, 37.0),
        (-102.05, 37.0),
        (-102.05, 41.0),
        (-109.05, 41.0),
        (-109.05, 37.0),
    ]
)

client = Client()
collection = client.dataset("open_data.copernicus.sentinel2_msi").collection("S2A_S2MSI2A")

scenes = collection.query(
    temporal_extent=("2025-10-01", "2025-11-01"),
    spatial_extent=area,
    show_progress=True,
)

selected = scenes.where(scenes.cloud_cover < 10, drop=True).isel(time=0)
print(selected.granule_name.item())

from pathlib import Path

from tilebox.storage import CopernicusStorageClient

storage = CopernicusStorageClient(
    access_key="YOUR_COPERNICUS_ACCESS_KEY",
    secret_access_key="YOUR_COPERNICUS_SECRET_ACCESS_KEY",
    cache_directory=Path("./data"),
)

product_path = storage.download(selected)

print(f"Downloaded {product_path.name} to {product_path}")
print("Contents:")
for path in product_path.iterdir():
    print(f"- {path.relative_to(product_path)}")

Downloaded S2A_MSIL2A_20251002T180751_N0511_R084_T13TEE_20251002T225842.SAFE to data/Sentinel-2/MSI/L2A/2025/10/02/S2A_MSIL2A_20251002T180751_N0511_R084_T13TEE_20251002T225842.SAFE
Contents:
- manifest.safe
- GRANULE
- INSPIRE.xml
- MTD_MSIL2A.xml
- DATASTRIP
- HTML
- rep_info
- S2A_MSIL2A_20251002T180751_N0511_R084_T13TEE_20251002T225842-ql.jpg

objects = storage.list_objects(selected)

wanted_bands = ["B02_10m", "B03_10m", "B04_10m", "B08_10m"]
band_objects = [
    obj for obj in objects
    if any(band in obj for band in wanted_bands)
]

for obj in band_objects:
    print(obj)

downloaded_files = storage.download_objects(selected, band_objects)
print(downloaded_files)

storage.quicklook(selected)