From 73d137b3222172ce64f3e26be1ec1669f2ff6fd9 Mon Sep 17 00:00:00 2001 From: "Morten W. Hansen" Date: Wed, 29 Jan 2025 09:39:58 +0100 Subject: [PATCH] add geolocation in ASAR mapper --- nansat/mappers/mapper_asar.py | 6 ++++-- nansat/nansat.py | 2 +- 2 files changed, 5 insertions(+), 3 deletions(-) diff --git a/nansat/mappers/mapper_asar.py b/nansat/mappers/mapper_asar.py index 60892e45..a33d2802 100644 --- 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 index afe9869c..20a30801 100644 --- 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