On-Demand Earth System Data Cubes (ESDCs) in Python
cubo:distance_from_center
coordinate.get_all_items
by item_collection
.dask>=2023.7.0
, pandas>=2.0.3
, planetary_computer>=1.0.0
, pystac_client>=0.7.2
, stackstac>=0.4.4
, and xarray>=2023.6.0
.create()
cubo
is pretty straightforward, everything you need is in the create()
function:
da = cubo.create(
lat=4.31,
lon=-76.2,
collection="sentinel-2-l2a",
bands=["B02","B03","B04"],
start_date="2021-06-01",
end_date="2021-06-10",
edge_size=64,
resolution=10,
)
By default, cubo
uses Planetary Computer. But you can use another STAC provider endpoint if you want:
da = cubo.create(
lat=4.31,
lon=-76.2,
collection="sentinel-s2-l2a-cogs",
bands=["B05","B06","B07"],
start_date="2020-01-01",
end_date="2020-06-01",
edge_size=128,
resolution=20,
stac="https://earth-search.aws.element84.com/v0"
)
You can pass kwargs
to pystac_client.Client.search()
if required:
da = cubo.create(
lat=4.31,
lon=-76.2,
collection="sentinel-2-l2a",
bands=["B02","B03","B04"],
start_date="2021-01-01",
end_date="2021-06-10",
edge_size=64,
resolution=10,
query={"eo:cloud_cover": {"lt": 10}} # kwarg to pass
)