Add a more complete and modern STAC Item maker

Author: alexamici

xarray_sentinel/esa_safe.py

@@ -15,6 +15,7 @@
     "safe": "http://www.esa.int/safe/sentinel-1.0",
     "s1": "http://www.esa.int/safe/sentinel-1.0/sentinel-1",
     "s1sarl1": "http://www.esa.int/safe/sentinel-1.0/sentinel-1/sar/level-1",
+    "gml": "http://www.opengis.net/gml",
 }
 
 SENTINEL1_FOLDER = importlib.resources.files("xarray_sentinel") / "resources/sentinel1"
@@ -176,6 +177,98 @@ def parse_manifest_sentinel1(
     return attributes, files
 
 
+def make_sentinel1_stac_item(
+    item_id: str,
+    manifest_path: PathOrFileType,
+    annotation: PathOrFileType,
+    namespaces=SENTINEL1_NAMESPACES,
+) -> dict[str, Any]:
+    manifest = ElementTree.parse(manifest_path).getroot()
+
+    product_information = parse_tag(annotation, ".//productInformation")
+    image_information = parse_tag(annotation, ".//imageInformation")
+
+    coordinates = [
+        [float(v) for v in token.split(",")]
+        for token in findtext(manifest, ".//gml:coordinates").split()
+    ]
+    coordinates += coordinates[:1]
+    product_timeliness = findtext(manifest, ".//s1sarl1:productTimelinessCategory")
+    product_timeliness_map = {
+        "Fast-24h": {
+            "product:timeliness_category": "STC",
+            "product:timeliness": "PT24H",
+        }
+    }
+
+    stac_item = {
+        "type": "Feature",
+        "stac_version": "1.1.0",
+        "stac_extensions": [
+            "https://stac-extensions.github.io/product/v0.1.0/schema.json",
+            "https://stac-extensions.github.io/processing/v1.2.0/schema.json",
+            "https://stac-extensions.github.io/sat/v1.0.0/schema.json",
+            "https://stac-extensions.github.io/view/v1.0.0/schema.json",
+            "https://stac-extensions.github.io/sar/v1.2.0/schema.json",
+            "https://stac-extensions.github.io/eopf/v1.0.0/schema.json",
+        ],
+        "id": item_id,
+        "properties": {
+            "datetime": None,
+            "start_datetime": findtext(manifest, ".//safe:startTime") + "Z",
+            "end_datetime": findtext(manifest, ".//safe:stopTime") + "Z",
+            "created": manifest.find(
+                ".//safe:processing", namespaces=namespaces
+            ).attrib["stop"]
+            + "Z",
+            "platform": "sentinel-1"
+            + findtext(manifest, ".//safe:platform/safe:number").lower(),
+            "instruments": ["sar"],
+            "constellation": "sentinel-1",
+            "product:type": "S01SSMSLC",
+            "product:timeliness_category": product_timeliness_map[product_timeliness][
+                "product:timeliness_category"
+            ],
+            "product:timeliness": product_timeliness_map[product_timeliness][
+                "product:timeliness"
+            ],
+            "processing:software": manifest.find(
+                ".//safe:software", namespaces=namespaces
+            ).attrib,
+            "sat:platform_international_designator": findtext(
+                manifest, ".//safe:nssdcIdentifier"
+            ),
+            "sat:absolute_orbit": int(findall(manifest, ".//safe:orbitNumber")[0]),
+            "sat:relative_orbit": int(
+                findall(manifest, ".//safe:relativeOrbitNumber")[0]
+            ),
+            "sat:orbit_state": findtext(manifest, ".//s1:pass").lower(),
+            "sat:anx_datetime": findtext(manifest, ".//s1:ascendingNodeTime") + "Z",
+            "view:incidence_angle": image_information["incidenceAngleMidSwath"],
+            "sar:polarizations": findall(
+                manifest, ".//s1sarl1:transmitterReceiverPolarisation"
+            ),
+            "sar:instrument_mode": findtext(
+                manifest, ".//s1sarl1:instrumentMode/s1sarl1:mode"
+            ),
+            "sar:frequency_band": "C",
+            "sar:center_frequency": product_information["radarFrequency"] / 1e9,
+            "sar:pixel_spacing_range": image_information["rangePixelSpacing"],
+            "sar:pixel_spacing_azimuth": image_information["azimuthPixelSpacing"],
+            "sar:observation_direction": "right",
+            "sar:beam_ids": findall(
+                manifest, ".//s1sarl1:instrumentMode/s1sarl1:swath"
+            ),
+            "eopf:datatake_id": int(findtext(manifest, ".//s1sarl1:missionDataTakeID")),
+        },
+        "geometry": {"type": "Polygon", "coordinates": [coordinates]},
+        "links": [],
+        "assets": {},
+    }
+    return stac_item
+
+
+# DEPRECATED
 def make_stac_item(attrs: Mapping[str, Any]) -> dict[str, Any]:
     assert attrs["family_name"] == "SENTINEL-1"