Add suport for S3 SLSTR #11

Author: mario-winkler

CHANGELOG.md

@@ -4,6 +4,14 @@ All notable changes to this project will be documented in this file.
 
 The format is based on [Keep a Changelog](https://keepachangelog.com/en/1.1.0/) and this project adheres to [Semantic Versioning](https://semver.org/).
 
+## [Unreleased]
+
+### Added
+
+- Add support for the following Sentinel 3 OLCI products: `S03OLCEFR`, `S03OLCERR`, `S03OLCLFR`, `S03OLCLRR` [#9](https://github.com/EOPF-Sample-Service/eopf-stac/issues/9)
+- Add support for the following Sentinel 3 SLSTR products: `S03SLSRBT`, `S03SLSLST` [#11](https://github.com/EOPF-Sample-Service/eopf-stac/issues/11)
+- Add tests [#10](https://github.com/EOPF-Sample-Service/eopf-stac/issues/10)
+
 ## [0.8.0] - 2025-04-15
 
 ### Added

src/eopf_stac/common/constants.py

@@ -24,13 +24,35 @@
 SUPPORTED_PRODUCT_TYPES_S2 = ["S02MSIL1C", "S02MSIL2A"]
 
 SUPPORTED_S3_OLCI_L1_PRODUCT_TYPES = ["S03OLCEFR", "S03OLCERR"]
-SUPPORTED_S3_OLCI_L2_PRODUCT_TYPES = ["S03OLCLFR"]
-SUPPORTED_PRODUCT_TYPES_S3 = SUPPORTED_S3_OLCI_L1_PRODUCT_TYPES + SUPPORTED_S3_OLCI_L2_PRODUCT_TYPES
+SUPPORTED_S3_OLCI_L2_PRODUCT_TYPES = ["S03OLCLFR", "S03OLCLRR"]
+SUPPORTED_S3_SLSTR_L1_PRODUCT_TYPES = ["S03SLSRBT"]
+SUPPORTED_S3_SLSTR_L2_LST_PRODUCT_TYPE = ["S03SLSLST"]
+# SUPPORTED_S3_SLSTR_L2_FRP_PRODUCT_TYPE = ["S03SLSFRP"]  # conversion error
+SUPPORTED_PRODUCT_TYPES_S3 = (
+    SUPPORTED_S3_OLCI_L1_PRODUCT_TYPES
+    + SUPPORTED_S3_OLCI_L2_PRODUCT_TYPES
+    + SUPPORTED_S3_SLSTR_L1_PRODUCT_TYPES
+    + SUPPORTED_S3_SLSTR_L2_LST_PRODUCT_TYPE
+)
 
 # Other Sentinen-3 product types to support
-# SRAL
-# SLSTR
-# SYN
+SUPPORTED_S3_SRAL_L1_PRODUCT_TYPES = ["S03AHRL1B"]  # sentinel-3-sra-l1b
+SUPPORTED_S3_SRAL_L2_PRODUCT_TYPES = ["S03AHRL2H"]  # sentinel-3-sra-l2-lan-hy
+SUPPORTED_S3_SYN_L2_PRODUCT_TYPES = [
+    "S03SYNAOD",
+    "S03SYNSDR",
+    "S03SYNV10",
+    "S03SYNVG1",
+    "S03SYNVGK",
+    "S03SYNVGP",
+]  # sentinel-3-syn-l2-aod, sentinel-3-syn-l2, sentinel-3-syn-l2-v10, sentinel-3-syn-l2-vg1, ?, sentinel-3-syn-l2-vgp
+
+
+# other SRAL listed in [1]
+#   - S03AHRL1A (SR_1_SRA_A_), S03ALRL1A (SR_1_SRA_A_), S03ALRL1B (SR_1_SRA_BS), S03ALRL2H (SR_2_LAN_HY)
+# other SYN listed in [1]: S03SYNMIS
+
+# [1] https://cpm.pages.eopf.copernicus.eu/eopf-cpm/main/PSFD/3-product-types-naming-rules.html
 
 PRODUCT_TYPE_TO_COLLECTION: Final[dict] = {
     "S01SIWGRH": "sentinel-1-l1-grd",
@@ -49,6 +71,10 @@
     "S03OLCEFR": "sentinel-3-olci-l1-efr",
     "S03OLCERR": "sentinel-3-olci-l1-err",
     "S03OLCLFR": "sentinel-3-olci-l2-lfr",
+    "S03OLCLRR": "sentinel-3-olci-l2-lrr",
+    "S03SLSRBT": "sentinel-3-slstr-l1-rbt",
+    "S03SLSFRP": "sentinel-3-slstr-l2-frp",
+    "S03SLSLST": "sentinel-3-slstr-l2-lst",
 }
 
 MEDIA_TYPE_ZARR = "application/vnd+zarr"

src/eopf_stac/sentinel3/constants.py

@@ -6,7 +6,7 @@
 from pystac.extensions.sat import OrbitState
 from pystac.item_assets import ItemAssetDefinition
 from pystac.utils import str_to_datetime
-from stactools.sentinel3.constants import SENTINEL_OLCI_BANDS
+from stactools.sentinel3.constants import SENTINEL_OLCI_BANDS, SENTINEL_SLSTR_BANDS
 
 from eopf_stac.common.constants import (
     DATASET_ASSET_EXTRA_FIELDS,
@@ -205,6 +205,124 @@ def get_olci_bands(band_keys: list[str] | None = None) -> list[dict]:
     PRODUCT_METADATA_ASSET_KEY: PRODUCT_METADATA_PATH,
 }
 
+# SENTINEL_SLSTR_BANDS
+# SLSTR_BANDS_TO_RESOLUTIONS
+
+
+def get_slstr_bands(band_keys: list[str] | None = None) -> list[dict]:
+    bands = []
+    if band_keys is None:
+        for _, band in SENTINEL_SLSTR_BANDS.items():
+            bands.append(band.to_dict())
+    else:
+        for key in band_keys:
+            bands.append(SENTINEL_SLSTR_BANDS[key].to_dict())
+
+    return bands
+
+
+SLSTR_L1_ASSETS: dict[str, ItemAssetDefinition] = {
+    "radiance_an": ItemAssetDefinition.create(
+        title="TOA radiance - stripe A, nadir view",
+        media_type=pystac.MediaType.ZARR,
+        description=("Dataset of the TOA radiances for the 500m grid, stripe A, nadir view"),
+        roles=[ROLE_DATA, ROLE_DATASET],
+        extra_fields={"bands": get_slstr_bands(["S01", "S02", "S03", "S04", "S05", "S06"]), "gsd": 500},
+    ),
+    "radiance_ao": ItemAssetDefinition.create(
+        title="TOA radiance - stripe A, oblique view",
+        media_type=pystac.MediaType.ZARR,
+        description=("Dataset of the TOA radiances for the 500m grid, stripe A, oblique view"),
+        roles=[ROLE_DATA, ROLE_DATASET],
+        extra_fields={"bands": get_slstr_bands(["S01", "S02", "S03", "S04", "S05", "S06"]), "gsd": 500},
+    ),
+    "radiance_bn": ItemAssetDefinition.create(
+        title="TOA radiance - stripe B, nadir view",
+        media_type=pystac.MediaType.ZARR,
+        description=("Dataset of the TOA radiances for the 500m grid, stripe B, nadir view"),
+        roles=[ROLE_DATA, ROLE_DATASET],
+        extra_fields={"bands": get_slstr_bands(["S04", "S05", "S06"]), "gsd": 500},
+    ),
+    "radiance_bo": ItemAssetDefinition.create(
+        title="TOA radiance - stripe B, oblique view",
+        media_type=pystac.MediaType.ZARR,
+        description=("Dataset of the TOA radiances for the 500m grid, stripe B, oblique view"),
+        roles=[ROLE_DATA, ROLE_DATASET],
+        extra_fields={"bands": get_slstr_bands(["S04", "S05", "S06"]), "gsd": 500},
+    ),
+    "BT_in": ItemAssetDefinition.create(
+        title="TOA brightness temperature - TIR, nadir view",
+        media_type=pystac.MediaType.ZARR,
+        description=("Dataset of the TOA brightness temperature for channels S7-S9 and F2 in the 1km grid, nadir view"),
+        roles=[ROLE_DATA, ROLE_DATASET],
+        extra_fields={"bands": get_slstr_bands(["S07", "S08", "S09", "S11"]), "gsd": 1000},
+    ),
+    "BT_io": ItemAssetDefinition.create(
+        title="TOA brightness temperature - TIR, oblique view",
+        media_type=pystac.MediaType.ZARR,
+        description=("Dataset of the TOA brightness temperature for channels S7-S9 and F2, 1km grid, oblique view"),
+        roles=[ROLE_DATA, ROLE_DATASET],
+        extra_fields={"bands": get_slstr_bands(["S07", "S08", "S09", "S11"]), "gsd": 1000},
+    ),
+    "BT_fn": ItemAssetDefinition.create(
+        title="TOA brightness temperature - F1, nadir view",
+        media_type=pystac.MediaType.ZARR,
+        description=("Dataset of the TOA brightness temperature for the F1 channel, 1km grid, nadir view"),
+        roles=[ROLE_DATA, ROLE_DATASET],
+        extra_fields={"bands": get_slstr_bands(["S10"]), "gsd": 1000},
+    ),
+    "BT_fo": ItemAssetDefinition.create(
+        title="TOA brightness temperature - F1, oblique view",
+        media_type=pystac.MediaType.ZARR,
+        description=("Dataset of the TOA brightness temperature for the F1 channel, 1km grid, oblique view"),
+        roles=[ROLE_DATA, ROLE_DATASET],
+        extra_fields={"bands": get_slstr_bands(["S10"]), "gsd": 1000},
+    ),
+    PRODUCT_ASSET_KEY: get_item_asset_product(),
+    PRODUCT_METADATA_ASSET_KEY: get_item_asset_metadata(),
+}
+
+SLSTR_L1_ASSETS_KEY_TO_PATH: dict[str:str] = {
+    "radiance_an": "measurements/anadir",
+    "radiance_ao": "measurements/aoblique",
+    "radiance_bn": "measurements/bnadir",
+    "radiance_bo": "measurements/boblique",
+    "BT_in": "measurements/inadir",
+    "BT_io": "measurements/ioblique",
+    "BT_fn": "measurements/fnadir",
+    "BT_fo": "measurements/foblique",
+    PRODUCT_ASSET_KEY: "",
+    PRODUCT_METADATA_ASSET_KEY: PRODUCT_METADATA_PATH,
+}
+
+SLSTR_L2_LST_ASSETS: dict[str, ItemAssetDefinition] = {
+    "lst": ItemAssetDefinition.create(
+        title="Land Surface Temperature (LST)",
+        media_type=pystac.MediaType.ZARR,
+        description=("Gridded Land Surface Temperature generated on the wide 1 km measurement grid"),
+        roles=[ROLE_DATA, ROLE_DATASET],
+        extra_fields={"bands": get_slstr_bands(["S07", "S08", "S09"]), "gsd": 1000},
+    ),
+    PRODUCT_ASSET_KEY: get_item_asset_product(),
+    PRODUCT_METADATA_ASSET_KEY: get_item_asset_metadata(),
+}
+
+SLSTR_L2_LST_ASSETS_KEY_TO_PATH: dict[str:str] = {
+    "lst": "measurements",
+    PRODUCT_ASSET_KEY: "",
+    PRODUCT_METADATA_ASSET_KEY: PRODUCT_METADATA_PATH,
+}
+
+SLSTR_L2_FRP_ASSETS: dict[str, ItemAssetDefinition] = {
+    PRODUCT_ASSET_KEY: get_item_asset_product(),
+    PRODUCT_METADATA_ASSET_KEY: get_item_asset_metadata(),
+}
+
+SLSTR_L2_FRP_ASSETS_KEY_TO_PATH: dict[str:str] = {
+    PRODUCT_ASSET_KEY: "",
+    PRODUCT_METADATA_ASSET_KEY: PRODUCT_METADATA_PATH,
+}
+
 # -- Collection metadata
 
 S3_OLCI_L1_EFR = {
@@ -218,7 +336,7 @@ def get_olci_bands(band_keys: list[str] | None = None) -> list[dict]:
     "product_type": "S03OLCEFR",
     "processing_level": "L1",
     "instrument": "olci",
-    "gsd": 300,
+    "gsd": [300],
     "item_assets": {**OLCI_L1_ASSETS},
 }
 
@@ -233,7 +351,7 @@ def get_olci_bands(band_keys: list[str] | None = None) -> list[dict]:
     "product_type": "S03OLCERR",
     "processing_level": "L1",
     "instrument": "olci",
-    "gsd": 1200,
+    "gsd": [1200],
     "item_assets": {**OLCI_L1_ASSETS},
 }
 
@@ -246,21 +364,10 @@ def get_olci_bands(band_keys: list[str] | None = None) -> list[dict]:
     "product_type": "S03OLCLFR",
     "processing_level": "L2",
     "instrument": "olci",
-    "gsd": 300,
+    "gsd": [300],
     "item_assets": {**OLCI_L2_ASSETS},
 }
 
-
-# S03SLSRBT / SL_1_RBT
-S3_SLSTR_L1_RBT = {}
-
-# S03SLSLST / SL_2_LST
-S3_SLSTR_L2_LST = {}
-
-# TBD: SRAL, SYN
-
-# Conversion not supported by CPM; no mapping
-
 S3_OLCI_L2_LRR = {
     "id": "sentinel-3-olci-l2-lrr",
     "title": "Sentinel-3 OLCI Level-2 LRR",
@@ -270,10 +377,53 @@ def get_olci_bands(band_keys: list[str] | None = None) -> list[dict]:
     "product_type": "S03OLCLRR",
     "processing_level": "L2",
     "instrument": "olci",
-    "gsd": 1200,
+    "gsd": [1200],
     "item_assets": {**OLCI_L2_ASSETS},
 }
 
+S3_SLSTR_L1_RBT = {
+    "id": "sentinel-3-slstr-l1-rbt",
+    "title": "Sentinel-3 SLSTR Level-1 RBT",
+    "description": (
+        "The Sentinel-3 SLSTR Level-1B RBT product provides radiances and brightness temperatures for each pixel "
+        "in a regular image grid for each view and SLSTR channel. In addition, it also contains annotations data "
+        "associated with each image pixels."
+    ),
+    "product_type": "S03SLSRBT",
+    "processing_level": "L1",
+    "instrument": "slstr",
+    "gsd": [500, 1000],
+    "item_assets": {**SLSTR_L1_ASSETS},
+}
+
+S3_SLSTR_L2_LST = {
+    "id": "sentinel-3-slstr-l2-lst",
+    "title": "Sentinel-3 SLSTR Level-2 LST",
+    "description": "The Sentinel-3 SLSTR Level-2 LST product provides land surface temperature.",
+    "product_type": "S03SLSLST",
+    "processing_level": "L2",
+    "instrument": "slstr",
+    "gsd": [500, 1000],
+    "item_assets": {**SLSTR_L2_LST_ASSETS},
+}
+
+S3_SLSTR_L2_FRP = {
+    "id": "sentinel-3-slstr-l2-frp",
+    "title": "Sentinel-3 SLSTR Level-2 FRP",
+    "description": (
+        "The Sentinel-3 SLSTR Level-2 FRP product provides global (over land and water) fire radiative power."
+    ),
+    "product_type": "S03SLSFRP",
+    "processing_level": "L2",
+    "instrument": "slstr",
+    "gsd": [500, 1000],
+    "item_assets": {**SLSTR_L2_FRP_ASSETS},
+}
+
+# TBD: SRAL, SYN
+
+# Conversion not supported by CPM; no mapping
+
 S3_OLCI_L2_WFR = {
     "id": "sentinel-3-olci-l2-wfr",
     "title": "Sentinel-3 OLCI Level-2 WFR",

src/eopf_stac/sentinel3/stac.py

@@ -11,6 +11,8 @@
     SENTINEL_LICENSE,
     SUPPORTED_S3_OLCI_L1_PRODUCT_TYPES,
     SUPPORTED_S3_OLCI_L2_PRODUCT_TYPES,
+    SUPPORTED_S3_SLSTR_L1_PRODUCT_TYPES,
+    SUPPORTED_S3_SLSTR_L2_LST_PRODUCT_TYPE,
     THUMBNAIL_ASSET,
 )
 from eopf_stac.common.stac import (
@@ -31,6 +33,10 @@
     OLCI_L2_ASSETS,
     OLCI_L2_ASSETS_KEY_TO_PATH,
     SENTINEL3_METADATA,
+    SLSTR_L1_ASSETS,
+    SLSTR_L1_ASSETS_KEY_TO_PATH,
+    SLSTR_L2_LST_ASSETS,
+    SLSTR_L2_LST_ASSETS_KEY_TO_PATH,
 )
 
 logger = logging.getLogger(__name__)
@@ -42,7 +48,7 @@ def create_collection(collection_metadata: dict, thumbnail_href: str) -> pystac.
         "constellation": [mission_metadata.get("constellation")],
         "platform": mission_metadata.get("platforms"),
         "instruments": [collection_metadata.get("instrument")],
-        "gsd": [collection_metadata.get("gsd")],
+        "gsd": collection_metadata.get("gsd"),
         "processing:level": [collection_metadata.get("processing_level")],
         "product:type": [collection_metadata.get("product_type")],
     }
@@ -152,6 +158,12 @@ def create_item(metadata: dict, product_type: str, asset_href_prefix: str) -> py
     elif product_type in SUPPORTED_S3_OLCI_L2_PRODUCT_TYPES:
         asset_defintions = OLCI_L2_ASSETS
         asset_path_lookups = OLCI_L2_ASSETS_KEY_TO_PATH
+    elif product_type in SUPPORTED_S3_SLSTR_L1_PRODUCT_TYPES:
+        asset_defintions = SLSTR_L1_ASSETS
+        asset_path_lookups = SLSTR_L1_ASSETS_KEY_TO_PATH
+    elif product_type in SUPPORTED_S3_SLSTR_L2_LST_PRODUCT_TYPE:
+        asset_defintions = SLSTR_L2_LST_ASSETS
+        asset_path_lookups = SLSTR_L2_LST_ASSETS_KEY_TO_PATH
     else:
         raise ValueError(f"Unsupported Sentinel-3 product type '{product_type}'")