switching to Copernicus Global 30m DEM as default

Author: BuddyVolly

ost/generic/common_wrappers.py

@@ -300,6 +300,8 @@ def ls_mask(infile, outfile, logfile, config_dict):
         f'-Poutput=\'{str(outfile)}\''
     )
 
+    print(command)
+
     # run command and get return code
     return_code = h.run_command(command, logfile)
 

ost/graphs/ard_json/grd.ceos.json

@@ -27,7 +27,7 @@
             "create_ls_mask": true,
             "apply_ls_mask": false,
             "dem": {
-                "dem_name": "SRTM 1Sec HGT",
+                "dem_name": "Copernicus 30m Global DEM",
                 "dem_file": "",
                 "dem_nodata": 0,
                 "dem_resampling": "BILINEAR_INTERPOLATION",

ost/graphs/ard_json/grd.earth_engine.json

@@ -27,7 +27,7 @@
             "create_ls_mask": false,
             "apply_ls_mask": false,
             "dem": {
-                "dem_name": "SRTM 1Sec HGT",
+                "dem_name": "Copernicus 30m Global DEM",
                 "dem_file": "",
                 "dem_nodata": 0,
                 "dem_resampling": "BILINEAR_INTERPOLATION",

ost/graphs/ard_json/grd.ost_gtc.json

@@ -27,7 +27,7 @@
             "create_ls_mask": false,
             "apply_ls_mask": false,
             "dem": {
-                "dem_name": "SRTM 1Sec HGT",
+                "dem_name": "Copernicus 30m Global DEM",
                 "dem_file": "",
                 "dem_nodata": 0,
                 "dem_resampling": "BILINEAR_INTERPOLATION",

ost/graphs/ard_json/grd.ost_rtc.json

@@ -27,7 +27,7 @@
             "create_ls_mask": true,
             "apply_ls_mask": false,
             "dem": {
-                "dem_name": "SRTM 1Sec HGT",
+                "dem_name": "Copernicus 30m Global DEM",
                 "dem_file": "",
                 "dem_nodata": 0,
                 "dem_resampling": "BILINEAR_INTERPOLATION",

ost/graphs/ard_json/slc.ost_gtc.json

@@ -38,7 +38,7 @@
             },
             "create_ls_mask": true,
             "dem": {
-                "dem_name": "SRTM 1Sec HGT",
+                "dem_name": "Copernicus 30m Global DEM",
                 "dem_file": " ",
                 "dem_nodata": 0,
                 "dem_resampling": "BILINEAR_INTERPOLATION",

ost/graphs/ard_json/slc.ost_minimal.json

@@ -38,7 +38,7 @@
             },
             "create_ls_mask": false,
             "dem": {
-                "dem_name": "SRTM 1Sec HGT",
+                "dem_name": "Copernicus 30m Global DEM",
                 "dem_file": "",
                 "dem_nodata": 0,
                 "dem_resampling": "BILINEAR_INTERPOLATION",

ost/graphs/ard_json/slc.ost_rtc.json

@@ -39,7 +39,7 @@
             "create_ls_mask": true,
             "apply_ls_mask": false,
             "dem": {
-                "dem_name": "SRTM 1Sec HGT",
+                "dem_name": "Copernicus 30m Global DEM",
                 "dem_file": "",
                 "dem_nodata": 0,
                 "dem_resampling": "BILINEAR_INTERPOLATION",

ost/s1/s1scene.py

@@ -788,23 +788,32 @@ 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(
+                    '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'
+
             logger.info(
-                'Scene is outside SRTM coverage. Snap will therefore use '
-                'the GETASSE30 DEM. Also consider to use a stereographic '
-                'projection.'
+                'The scene\'s location is towards the poles. '
+                'Consider to use a stereographic projection.'
             )
+
             epsg = input(
-                'Please type the EPSG you want to project 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'
+                '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'][
-                'dem_name'] = 'GETASSE30'
             self.ard_parameters['single_ARD']['dem'][
                 'out_projection'] = int(epsg)