Sentinel-3 SLSTR L2 LST Data Access
Demonstrate how to access Sentinel-3 SLSTR L2 LST data in EOPF Zarr format
Table of Contents¶
Run this notebook interactively with all dependencies pre-installed
This notebook contain the most common data access procedures for the Sentinel-3 SLSTR L2 LST EOPF Zarr data collection.
For more examples on how to use the data in different application scenarios, please refer to the EOPF Sample Notebooks available at: https://
Discover Products using STAC¶
Discover Sentinel-3_SLSTR_L2_LST products using the EOPF STAC catalog at https://stac.core.eopf.eodc.eu.
The STAC Collection is https://stac.core.eopf.eodc.eu/collections/sentinel-3-slstr-l2-lst .
import pystac_client
import datetime
from datetime import timedelta
from xcube.core.store import new_data_store
import xarray as xr
xr.set_options(display_expand_attrs=False)<xarray.core.options.set_options at 0x7f754da8c5f0>STAC_URL = "https://stac.core.eopf.eodc.eu"
catalog = pystac_client.Client.open(STAC_URL)
# Define area of interest
bbox = [-72.7017, 61.1864, -28.7009, 75.3745]
# Temporal extent: today - 5 days to today
today = datetime.datetime.now(datetime.UTC).date()
start_date = today - timedelta(days=5)
temporal_extent = [str(start_date), str(today)]
print("Discovering Sentinel-3 products from STAC catalog...")
s3_items = list(
catalog.search(
collections=["sentinel-3-slstr-l2-lst"],
bbox=bbox,
datetime=temporal_extent,
).items()
)
print(f"Found {len(s3_items)} Sentinel-3 SLSTR L2 LST products")Discovering Sentinel-3 products from STAC catalog...
Found 140 Sentinel-3 SLSTR L2 LST products
Open a single product with xarray¶
This step will open a Sentinel-3 SLSTR L2 LST product as an Xarray.DataTree object.
stac_item = s3_items[-1]
zarr_url = stac_item.assets["product"].href
ds = xr.open_datatree(zarr_url, engine="zarr")
dsOpen a single product with xarray-eopf¶
This step will open a Sentinel-3 product as an Xarray.Dataset object.
For the package documentation please visit https://
ds = xr.open_dataset(zarr_url, engine="eopf-zarr")
dsOpen multiple products with xcube-eopf¶
This step will combine all the products into a datacube.
For the package documentation please visit https://
# Define area of interest
[-72.7017, 61.1864, -28.7009, 75.3745]
# Temporal extent: today - 5 day to today
today = datetime.datetime.now(datetime.UTC).date()
start_date = today - timedelta(days=5)
temporal_extent = [str(start_date), str(today)]
store = new_data_store("eopf-zarr")
ds = store.open_data(
data_id="sentinel-3-slstr-l2-lst",
bbox=bbox,
time_range=temporal_extent,
spatial_res=3000 / 111320, # meters → degrees (approx.)
crs="EPSG:4326",
interp_methods="nearest",
)
ds