adopt polewards questinaire before processing

Author: BuddyVolly

ost/Project.py

@@ -652,10 +652,36 @@ def bursts_to_ards(
         # when outside SRTM
         center_lat = loads(self.aoi).centroid.y
         if float(center_lat) > 59 or float(center_lat) < -59:
-            logger.info('Scene is outside SRTM coverage. Snap will therefore '
-                        'use the Earth\'s geoid model.')
-            self.ard_parameters['single_ARD']['dem'][
-                'dem_name'] = 'Aster 1sec GDEM'
+
+            if 'SRTM' in self.ard_parameters['single_ARD']['dem']['dem_name']:
+                logger.info(
+                    'Scene is outside SRTM coverage. Snap will therefore use '
+                    'the Copernicus 30m Global DEM. '
+                )
+
+                self.ard_parameters['single_ARD']['dem'][
+                    'dem_name'] = 'Copernicus 30m Global DEM'
+
+            if self.ard_parameters['dem']['out_projection'] == 4326:
+
+                logger.info(
+                    'The scene\'s location is towards the poles. '
+                    'Consider to use a stereographic projection.'
+                )
+
+                epsg = input(
+                    'Type an alternative EPSG code for the projection of the '
+                    'output data or just press enter for keeping Lat/Lon '
+                    'coordinate system (e.g. 3413 for NSIDC Sea Ice Polar '
+                    'Stereographic North projection, or 3976 for '
+                    'NSIDC Sea Ice Polar Stereographic South projection'
+                )
+
+                if not epsg:
+                    epsg = 4326
+
+                self.ard_parameters['single_ARD']['dem'][
+                    'out_projection'] = int(epsg)
 
         # --------------------------------------------
         # 3 subset determination

ost/s1/s1scene.py

@@ -788,6 +788,7 @@ def create_ard(self, infile, out_dir, subset=None, overwrite=False):
         # when outside SRTM
         center_lat = self._get_center_lat(infile)
         if float(center_lat) > 59 or float(center_lat) < -59:
+
             if 'SRTM' in self.ard_parameters['single_ARD']['dem']['dem_name']:
 
                 logger.info(
@@ -798,24 +799,26 @@ def create_ard(self, infile, out_dir, subset=None, overwrite=False):
                 self.ard_parameters['single_ARD']['dem'][
                     'dem_name'] = 'Copernicus 30m Global DEM'
 
-            logger.info(
-                'The scene\'s location is towards the poles. '
-                'Consider to use a stereographic projection.'
-            )
+            if self.ard_parameters['dem']['out_projection'] == 4326:
 
-            epsg = input(
-                'Type an alternative EPSG code for the projection of the '
-                'output data or just press enter for keeping Lat/Lon '
-                'coordinate system (e.g. 3413 for NSIDC Sea Ice Polar '
-                'Stereographic North projection, or 3976 for NSIDC Sea Ice '
-                'Polar Stereographic South projection'
-            )
+                logger.info(
+                    'The scene\'s location is towards the poles. '
+                    'Consider to use a stereographic projection.'
+                )
+
+                epsg = input(
+                    'Type an alternative EPSG code for the projection of the '
+                    'output data or just press enter for keeping Lat/Lon '
+                    'coordinate system (e.g. 3413 for NSIDC Sea Ice Polar '
+                    'Stereographic North projection, or 3976 for '
+                    'NSIDC Sea Ice Polar Stereographic South projection'
+                )
 
-            if not epsg:
-                epsg = 4326
+                if not epsg:
+                    epsg = 4326
 
-            self.ard_parameters['single_ARD']['dem'][
-                'out_projection'] = int(epsg)
+                self.ard_parameters['single_ARD']['dem'][
+                    'out_projection'] = int(epsg)
 
         # --------------------------------------------
         # 3 Check ard parameters in case they have been updated,