From 648e50fa7aefd66c1bce8b53e29fe41b6d2355a9 Mon Sep 17 00:00:00 2001
From: Scott Staniewicz <scott.j.staniewicz@jpl.nasa.gov>
Date: Tue, 4 Feb 2025 11:38:19 -0500
Subject: [PATCH] Add estimated geolocation accuracies for CEOS metadata (#261)

---
 src/disp_s1/product.py | 62 ++++++++++++++++++++++++++++++++++++++++++
 1 file changed, 62 insertions(+)

diff --git a/src/disp_s1/product.py b/src/disp_s1/product.py
index ff54b7a..c852e21 100644
--- a/src/disp_s1/product.py
+++ b/src/disp_s1/product.py
@@ -1263,6 +1263,17 @@ def _to_string(model: YamlModel):
             fillvalue=None,
             description="DOI of reference describing phase cosine similarity metric",
         )
+        _create_dataset(
+            group=metadata_group,
+            name="ceos_estimated_phase_quality_metric_algorithm",
+            dimensions=(),
+            data="Gaussian filtering",
+            fillvalue=None,
+            description=(
+                "Algorithm used on wrapped phase to create estimated phase"
+                " quality metric"
+            ),
+        )
         _create_dataset(
             group=metadata_group,
             name="ceos_phase_unwrapping_method",
@@ -1326,6 +1337,57 @@ def _to_string(model: YamlModel):
                 " DOI."
             ),
         )
+        # CEOS 4.3
+        _create_dataset(
+            group=metadata_group,
+            name="ceos_absolute_geolocation_ground_range_bias",
+            dimensions=(),
+            data=-0.06,
+            fillvalue=None,
+            attrs={"units": "meters"},
+            description=(
+                "Estimated mean absolute geolocation error for input geocoded SLCs, as"
+                " assessed via corner reflector analysis, in the ground range direction"
+            ),
+        )
+        _create_dataset(
+            group=metadata_group,
+            name="ceos_absolute_geolocation_ground_range_stddev",
+            dimensions=(),
+            data=0.38,
+            fillvalue=None,
+            attrs={"units": "meters"},
+            description=(
+                "Estimated standard deviation of absolute geolocation error for input"
+                " geocoded SLCs, as assessed via corner reflector analysis, in the"
+                " ground range direction"
+            ),
+        )
+        _create_dataset(
+            group=metadata_group,
+            name="ceos_absolute_geolocation_azimuth_bias",
+            dimensions=(),
+            data=-0.04,
+            fillvalue=None,
+            attrs={"units": "meters"},
+            description=(
+                "Estimated mean absolute geolocation error for input geocoded SLCs, as"
+                " assessed via corner reflector analysis, in the azimuth direction"
+            ),
+        )
+        _create_dataset(
+            group=metadata_group,
+            name="ceos_absolute_geolocation_azimuth_stddev",
+            dimensions=(),
+            data=0.46,
+            fillvalue=None,
+            attrs={"units": "meters"},
+            description=(
+                "Estimated standard deviation of absolute geolocation error for input"
+                " geocoded SLCs, as assessed via corner reflector analysis, in the"
+                " azimuth direction"
+            ),
+        )
 
 
 def _get_orbit_direction(cslc_filename: Filename) -> Literal["ascending", "descending"]: