Fix gcps type for newer rioxarray versions

mraspaud · Jun 25, 2024 · ad353f0 · ad353f0
1 parent df4a1ba
commit ad353f0
Showing 1 changed file with 11 additions and 2 deletions.
diff --git a/satpy/readers/sar_c_safe.py b/satpy/readers/sar_c_safe.py
@@ -35,6 +35,7 @@
 """
 
 import functools
+import json
 import logging
 import warnings
 from collections import defaultdict
@@ -664,7 +665,7 @@ def get_gcps(self):
            gcp_coords (tuple): longitude and latitude 1d arrays
 
         """
-        gcps = self._data.coords["spatial_ref"].attrs["gcps"]
+        gcps = get_gcps_from_array(self._data)
         crs = self._data.rio.crs
 
         gcp_list = [(feature["properties"]["row"], feature["properties"]["col"], *feature["geometry"]["coordinates"])
@@ -726,7 +727,7 @@ def load(self, dataset_keys, **kwargs):
                     if key["name"] not in ["longitude", "latitude"]:
                         lonlats = self.load([DataID(self._id_keys, name="longitude", polarization=key["polarization"]),
                                              DataID(self._id_keys, name="latitude", polarization=key["polarization"])])
-                        gcps = val.coords["spatial_ref"].attrs["gcps"]
+                        gcps = get_gcps_from_array(val)
                         from pyresample.future.geometry import SwathDefinition
                         val.attrs["area"] = SwathDefinition(lonlats["longitude"], lonlats["latitude"],
                                                             attrs=dict(gcps=gcps))
@@ -797,3 +798,11 @@ def _create_measurement_handlers(self, calibrators, denoisers):
                                                          filetype_info=None)
 
         return measurement_handlers
+
+
+def get_gcps_from_array(val):
+    """Get the gcps from the spatial_ref coordinate as a geojson dict."""
+    gcps = val.coords["spatial_ref"].attrs["gcps"]
+    if isinstance(gcps, str):
+        gcps = json.loads(gcps)
+    return gcps