add geolocation in ASAR mapper

mortenwh · Jan 29, 2025 · 73d137b · 73d137b
1 parent 3cd073a
commit 73d137b
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 3 deletions.
diff --git a/nansat/mappers/mapper_asar.py b/nansat/mappers/mapper_asar.py
@@ -34,7 +34,7 @@ class Mapper(VRT, Envisat):
             http://envisat.esa.int/handbooks/asar/CNTR6-6-9.htm#eph.asar.asardf.asarrec.ASAR_Geo_Grid_ADSR
     '''
 
-    def __init__(self, filename, gdalDataset, gdalMetadata, **kwargs):
+    def __init__(self, filename, gdalDataset, gdalMetadata, step=1, **kwargs):
 
         '''
         Parameters
@@ -235,7 +235,7 @@ def __init__(self, filename, gdalDataset, gdalMetadata, **kwargs):
         # add bands with metadata and corresponding values to the empty VRT
         self.create_bands(metaDict)
 
-        # Add oribit and look information to metadata domain
+        # Add orbit and look information to metadata domain
         # ASAR is always right-looking
         self.dataset.SetMetadataItem('ANTENNA_POINTING', 'RIGHT')
         self.dataset.SetMetadataItem('ORBIT_DIRECTION',
@@ -293,6 +293,8 @@ def __init__(self, filename, gdalDataset, gdalMetadata, **kwargs):
         self.dataset.SetMetadataItem('instrument', json.dumps(mm))
         self.dataset.SetMetadataItem('platform', json.dumps(ee))
 
+        self.add_geolocation_from_ads(gdalDataset)
+
         # Set projection
         if self.dataset.GetProjection() == "":
             self.dataset.SetProjection(NSR().wkt)
diff --git a/nansat/nansat.py b/nansat/nansat.py
@@ -1241,7 +1241,7 @@ def get_transect(self, points, bands,
         Parameters
         ----------
         points : 2xN list or array, N (number of points) >= 1
-            coordinates [[x1, x2, y2], [y1, y2, y3]]
+            coordinates [[x1, x2, x3], [y1, y2, y3]]
         bands : list of int or string
             elements of the list are band number or band Name
         lonlat : bool