add debug + change reltive import and pin version xarray

Author: agrouaze

highleveltests/open_GRD_IW.py

@@ -0,0 +1,23 @@
+import pdb
+from safe_s1 import Sentinel1Reader, getconfig
+import time
+import logging
+logging.basicConfig(level=logging.DEBUG)
+logging.debug('start GRD test')
+conf = getconfig.get_config()
+subswath = conf['product_paths'][0]
+print(subswath)
+t0 = time.time()
+if 'GRD' in subswath:
+    sub_reader = Sentinel1Reader(subswath)
+else:
+    sub_reader = Sentinel1Reader('SENTINEL1_DS:'+subswath+':IW3')
+elapse_t = time.time()-t0
+
+dt = sub_reader.datatree
+print('out of the reader')
+print(dt)
+print('time to read the SAFE through nfs: %1.2f sec'%elapse_t)
+DN = sub_reader.load_digital_number(chunks={'pol':'VV','line':6000,'sample':8000})
+print('DN',DN)
+# pdb.set_trace()

highleveltests/open_SLC_IW.py

@@ -11,5 +11,5 @@
 dt = sub_reader.datatree
 print('out of the reader')
 print(dt)
-print('time to read the SAFE through S3: %1.2f sec'%elapse_t)
+print('time to read the SAFE through nfs: %1.2f sec'%elapse_t)
 pdb.set_trace()

pyproject.toml

@@ -5,7 +5,7 @@ license = {text = "MIT"}
 dependencies = [
     "geopandas",
     "numpy",
-    "xarray",
+    "xarray>=2024.10.0",
     "lxml",
     "rioxarray",
     "jmespath",

safe_s1/reader.py

@@ -8,9 +8,9 @@
 import yaml
 from affine import Affine
 from rioxarray import rioxarray
-
-from . import sentinel1_xml_mappings
-from .xml_parser import XmlParser
+import logging
+from safe_s1 import sentinel1_xml_mappings
+from safe_s1.xml_parser import XmlParser
 import xarray as xr
 import pandas as pd
 import warnings
@@ -19,6 +19,7 @@
 class Sentinel1Reader:
 
     def __init__(self, name, backend_kwargs=None):
+        logging.debug('input name: %s',name)
         if not isinstance(name, (str, os.PathLike)):
            raise ValueError(f"cannot deal with object of type {type(name)}: {name}")
         # gdal dataset name
@@ -28,13 +29,19 @@ def __init__(self, name, backend_kwargs=None):
         """Gdal dataset name"""
         name_parts = self.name.split(':')
         if len(name_parts) > 3:
+            logging.debug('windows case')
             # windows might have semicolon in path ('c:\...')
             name_parts[1] = ':'.join(name_parts[1:-1])
             del name_parts[2:-1]
         name_parts[1] = os.path.basename(name_parts[1])
         self.short_name = ':'.join(name_parts)
+        logging.debug('short_name : %s',self.short_name)
         """Like name, but without path"""
-        self.path = ':'.join(self.name.split(':')[1:-1])
+        if len(name_parts) == 2:
+            self.path = self.name.split(':')[1]
+        else:
+            self.path = ':'.join(self.name.split(':')[1:-1])
+        logging.debug('path: %s',self.path)
         # remove trailing slash in the safe path
         if self.path[-1]=='/':
             self.path = self.path.rstrip('/')