ENH: Add TerraSAR-X management of default resolution and pixel size for RE data

Author: remi-braun

CHANGES.md

@@ -4,14 +4,15 @@
 
 - **ENH: Desambiguate condensed name of Custom stack in case of creation of several objects with the same datetime and same constellation and product type**
 - **ENH: Fix corrupted Maxar products with incoherent width between .IMD and .TIL files** [#242](https://github.com/sertit/eoreader/issues/242)
+- **ENH: Add TerraSAR-X management of default resolution and pixel size for RE data**
 - FIX: Fix regression when stacking with a custom nodata value with VHR data to be reprojected
 - FIX: Fix an unprecedented case with a PNEO having different name than usual (`DIM_PNEO3_STD_2025...` instead of `DIM_PNEO3_2025...`)
 - FIX: Fix WV Legion wrong band order in Multi Spectral 1 mode [#246](https://github.com/sertit/eoreader/issues/246)
 - FIX: Collocate Planet spectral bands with masks for some rare cases where it fails
 - FIX: Loosen the constraints on PlanetScope stack name as it may change, from `Analytic` to `composite`, etc. [#244](https://github.com/sertit/eoreader/issues/244)
 - FIX: Add a fallback in case of impossibleness of reading ICEYE `QUICKLOOK.kml` file
 - FIX: Manage the case of Maxar data with negative absolute calibration factor: don't compute the reflectance and leave it as is.
-- FIX: Fix DEM management with RPC orthorectification: handle correctly the vertical CRS (see dem notebook and `EOREADER_DEM_VCRS` environment variable). [#53](https://github.com/sertit/eoreader/issues/53)
+- FIX: Fix DEM management with RPC orthorectification: handle correctly the vertical CRS (see DEM notebook and `EOREADER_DEM_VCRS` environment variable). [#53](https://github.com/sertit/eoreader/issues/53)
 - CI: Filter some warnings in pytest
 
 ## 0.22.4 (2025-07-07)

eoreader/products/sar/tsx_product.py

@@ -196,15 +196,15 @@ def _set_pixel_size(self) -> None:
             - We force Spatially Enhanced Resolution (SE) as we keep SSC resolutions (as per the ESA Data Access Portfolio)
             - We use the pixel_size from SE products, with incidence angle corresponding to the SSC resolution (or the highest pixel_spacing)
         """
-        # TODO: Manage RE case? Not handled by Copernicus EMS, so be careful...
-
         # Read metadata
         try:
             root, _ = self.read_mtd()
             acq_info = root.find(".//acquisitionInfo")
             polarization = TsxPolarization.from_value(
                 acq_info.findtext(".//polarisationMode")
             )
+            res_variant = root.findtext(".//resolutionVariant")
+            is_re = res_variant == "RE"
         except (InvalidProductError, TypeError) as exc:
             raise InvalidProductError(
                 "acquisitionInfo or polarisationMode not found in metadata!"
@@ -214,29 +214,57 @@ def _set_pixel_size(self) -> None:
         def_pixel_size = None
         if self.sensor_mode == TsxSensorMode.HS:
             if polarization == TsxPolarization.S:
-                def_pixel_size = 0.5
-                def_res = 1.1
+                if is_re:
+                    def_pixel_size = 1.5
+                    def_res = 3.1
+                else:
+                    def_pixel_size = 0.5
+                    def_res = 1.4
             elif polarization == TsxPolarization.D:
-                def_pixel_size = 1.0
-                def_res = 2.2
+                if is_re:
+                    def_pixel_size = 2.0
+                    def_res = 4.4
+                else:
+                    def_pixel_size = 1.0
+                    def_res = 2.2
         elif self.sensor_mode == TsxSensorMode.SL:
             if polarization == TsxPolarization.S:
-                def_pixel_size = 0.75
-                def_res = 1.7
+                if is_re:
+                    def_pixel_size = 1.75
+                    def_res = 3.8
+                else:
+                    def_pixel_size = 0.75
+                    def_res = 1.7
             elif polarization == TsxPolarization.D:
-                def_pixel_size = 1.0
-                def_res = 3.4
+                if is_re:
+                    def_pixel_size = 2.5
+                    def_res = 5.5
+                else:
+                    def_pixel_size = 1.0
+                    def_res = 3.4
         elif self.sensor_mode == TsxSensorMode.ST:
             if polarization == TsxPolarization.S:
-                def_pixel_size = 0.2
-                def_res = 0.24
+                if is_re:
+                    def_pixel_size = 0.4
+                    def_res = 0.9
+                else:
+                    def_pixel_size = 0.2
+                    def_res = 0.24
         elif self.sensor_mode == TsxSensorMode.SM:
             if polarization == TsxPolarization.S:
-                def_pixel_size = 1.25
-                def_res = 3.3
+                if is_re:
+                    def_pixel_size = 3.25
+                    def_res = 7.0
+                else:
+                    def_pixel_size = 1.25
+                    def_res = 3.3
             elif polarization == TsxPolarization.D:
-                def_pixel_size = 3.0
-                def_res = 6.6
+                if is_re:
+                    def_pixel_size = 4.5
+                    def_res = 9.9
+                else:
+                    def_pixel_size = 3.0
+                    def_res = 6.6
         elif self.sensor_mode == TsxSensorMode.SC:
             # Read metadata
             try:
@@ -262,6 +290,12 @@ def _set_pixel_size(self) -> None:
                 f"Unknown sensor mode: {self.sensor_mode} or unknown polarization: {polarization}"
             )
 
+        # Overrides the pixel size for ortho EEC products
+        if self.product_type == TsxProductType.EEC:
+            with rasterio.open(
+                str(self.get_raw_band_paths()[self.get_default_band()])
+            ) as ds:
+                def_pixel_size = ds.res[0]
         self.pixel_size = def_pixel_size
         self.resolution = def_res