OPTIM: Use raw band path over default band path to speed up the computation of image information (extent, footprint, CRS, transform) of SAR ortho products

Author: remi-braun

CHANGES.md

@@ -13,6 +13,7 @@
 - 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)
+- OPTIM: Use raw band path over default band path to speed up the computation of image information (extent, footprint, CRS, transform) of SAR ortho products
 - CI: Filter some warnings in pytest
 
 ## 0.22.4 (2025-07-07)

eoreader/products/sar/sar_product.py

@@ -28,7 +28,8 @@
 import rasterio
 import rioxarray
 import xarray as xr
-from rasterio import crs
+from affine import Affine
+from rasterio import CRS, crs
 from rasterio.enums import Resampling
 from rasterio.windows import Window
 from sertit import AnyPath, geometry, misc, path, rasters, snap, strings, types, vectors
@@ -287,15 +288,46 @@ def footprint(self) -> gpd.GeoDataFrame:
         Returns:
             gpd.GeoDataFrame: Footprint as a GeoDataFrame
         """
+        if self.is_ortho:
+            default_band_path = self.get_raw_band_paths()[self.get_default_band()]
+        else:
+            default_band_path = self.get_default_band_path()
+
         # Processed by SNAP: the nodata is set -> use get_footprint instead of vectorize
         downsampled_band = utils.read(
-            self.get_default_band_path(),
+            default_band_path,
             pixel_size=max(
                 self.resolution * 10, float(os.environ.get(SAR_DEF_PIXEL_SIZE, 0))
             ),
         )
         return rasters.get_footprint(downsampled_band)
 
+    @cache
+    def default_transform(self, **kwargs) -> (Affine, int, int, CRS):
+        """
+        Returns default transform data of the default band (UTM),
+        as the :code:`rasterio.warp.calculate_default_transform` does:
+        - transform
+        - width
+        - height
+        - crs
+
+        Args:
+            kwargs: Additional arguments
+        Returns:
+            Affine, int, int, CRS: transform, width, height, CRS
+
+        """
+        if self.is_ortho:
+            default_band_path = self.get_raw_band_paths(**kwargs)[
+                self.get_default_band()
+            ]
+        else:
+            default_band_path = self.get_default_band_path(**kwargs)
+
+        with rasterio.open(str(default_band_path)) as dst:
+            return dst.transform, dst.width, dst.height, dst.crs
+
     def get_default_band(self) -> BandNames:
         """
         Get default band:
@@ -398,11 +430,16 @@ def extent(self) -> gpd.GeoDataFrame:
         Returns:
             gpd.GeoDataFrame: Extent in UTM
         """
-        # Get WGS84 extent
-        extent_wgs84 = self.wgs84_extent()
+        if self.is_ortho:
+            return rasters.get_extent(
+                self.get_raw_band_paths()[self.get_default_band()]
+            ).to_crs(self.crs())
+        else:
+            # Get WGS84 extent
+            extent_wgs84 = self.wgs84_extent()
 
-        # Convert to UTM
-        return extent_wgs84.to_crs(self.crs())
+            # Convert to UTM
+            return extent_wgs84.to_crs(self.crs())
 
     @cache
     def crs(self) -> crs.CRS:
@@ -420,11 +457,17 @@ def crs(self) -> crs.CRS:
         Returns:
             crs.CRS: CRS object
         """
-        # Get WGS84 extent
-        extent_wgs84 = self.wgs84_extent()
+        if self.is_ortho:
+            with rasterio.open(
+                str(self.get_raw_band_paths()[self.get_default_band()])
+            ) as ds:
+                return ds.crs
+        else:
+            # Get WGS84 extent
+            extent_wgs84 = self.wgs84_extent()
 
-        # Estimate UTM from extent
-        return extent_wgs84.estimate_utm_crs()
+            # Estimate UTM from extent
+            return extent_wgs84.estimate_utm_crs()
 
     @abstractmethod
     def _set_sensor_mode(self) -> None:

eoreader/products/sar/tsx_product.py

@@ -26,8 +26,7 @@
 import geopandas as gpd
 import rasterio
 from lxml import etree
-from rasterio import crs
-from sertit import path, rasters, vectors
+from sertit import path, vectors
 from sertit.misc import ListEnum
 from sertit.types import AnyPathStrType
 
@@ -346,51 +345,6 @@ def _post_init(self, **kwargs) -> None:
         # Post init done by the super class
         super()._post_init(**kwargs)
 
-    @cache
-    def extent(self) -> gpd.GeoDataFrame:
-        """
-        Get UTM extent of the tile
-
-        .. code-block:: python
-
-            >>> from eoreader.reader import Reader
-            >>> path = r"S1A_IW_GRDH_1SDV_20191215T060906_20191215T060931_030355_0378F7_3696.zip"
-            >>> prod = Reader().open(path)
-            >>> prod.utm_extent()
-                                   Name  ...                                           geometry
-            0  Sentinel-1 Image Overlay  ...  POLYGON ((817914.501 4684349.823, 555708.624 4...
-            [1 rows x 12 columns]
-
-        Returns:
-            gpd.GeoDataFrame: Extent in UTM
-        """
-        if self.product_type == TsxProductType.EEC:
-            return rasters.get_extent(self.get_default_band_path()).to_crs(self.crs())
-        else:
-            return super().extent()
-
-    @cache
-    def crs(self) -> crs.CRS:
-        """
-        Get UTM projection
-
-        .. code-block:: python
-
-            >>> from eoreader.reader import Reader
-            >>> path = r"S1A_IW_GRDH_1SDV_20191215T060906_20191215T060931_030355_0378F7_3696.zip"
-            >>> prod = Reader().open(path)
-            >>> prod.utm_crs()
-            CRS.from_epsg(32630)
-
-        Returns:
-            crs.CRS: CRS object
-        """
-        if self.product_type == TsxProductType.EEC:
-            with rasterio.open(str(self.get_default_band_path())) as ds:
-                return ds.crs
-        else:
-            return super().crs()
-
     @cache
     def wgs84_extent(self) -> gpd.GeoDataFrame:
         """