solving map projection issue

Author: BuddyVolly

ost/Project.py

@@ -1,6 +1,14 @@
 #!/usr/bin/env python
 # -*- coding: utf-8 -*-
 
+# ------ bug of rasterio --------
+import os
+if "GDAL_DATA" in list(os.environ.keys()):
+    del os.environ["GDAL_DATA"]
+if "PROJ_LIB" in list(os.environ.keys()):
+    del os.environ["PROJ_LIB"]
+# ------ bug of rasterio --------
+
 import json
 import urllib.request
 import urllib.parse

ost/generic/common_wrappers.py

@@ -3,7 +3,6 @@
 
 import logging
 from retrying import retry
-
 from rasterio.crs import CRS
 
 from ost.helpers import helpers as h
@@ -189,14 +188,15 @@ def terrain_correction(infile, outfile, logfile, config_dict):
     cpus = config_dict['snap_cpu_parallelism']
 
     # auto projections of snap
-    # if 42001 <= dem_dict['out_projection'] <= 97002:
-    #    projection = f"AUTO:{dem_dict['out_projection']}"
+    if 42001 <= dem_dict['out_projection'] <= 97002:
+        projection = f"AUTO:{dem_dict['out_projection']}"
     # epsg codes
-    # else:
-    #    projection = f"EPSG:{dem_dict['out_projection']}"
+    elif int(dem_dict['out_projection']) == 4326:
+        projection = CRS.from_epsg(4326).to_wkt()
+    else:
+        projection = f"EPSG:{dem_dict['out_projection']}"
 
-    projection = CRS.from_epsg(dem_dict['out_projection']).to_wkt()
-    #projection = 'GEOGCRS["WGS 84",DATUM["World Geodetic System 1984",ELLIPSOID["WGS 84",6378137,298.257223563,LENGTHUNIT["metre",1]]],PRIMEM["Greenwich",0,ANGLEUNIT["degree",0.0174532925199433]],CS[ellipsoidal,2],AXIS["geodetic latitude (Lat)",north,ORDER[1],ANGLEUNIT["degree",0.0174532925199433]],AXIS["geodetic longitude (Lon)",east,ORDER[2],ANGLEUNIT["degree",0.0174532925199433]],USAGE[SCOPE["unknown"],AREA["World"],BBOX[-90,-180,90,180]],ID["EPSG",4326]]'
+    #
     logger.debug('Geocoding product.')
 
     if ard['geocoding'] == 'terrain':
@@ -211,7 +211,7 @@ def terrain_correction(infile, outfile, logfile, config_dict):
             f"-PimgResamplingMethod=\'{dem_dict['image_resampling']}\' "
             f"-PpixelSpacingInMeter={ard['resolution']} "
             f"-PalignToStandardGrid=true "
-            #f"-PmapProjection=\'{projection}\' "
+            f"-PmapProjection=\'{projection}\' "
             f"-t \'{str(outfile)}\' \'{str(infile)}\' "
         )
     elif ard['geocoding'] == 'ellipsoid':
@@ -226,7 +226,7 @@ def terrain_correction(infile, outfile, logfile, config_dict):
             f"-PimgResamplingMethod=\'{dem_dict['image_resampling']}\' "
             f"-PpixelSpacingInMeter={ard['resolution']} "
             f"-PalignToStandardGrid=true "
-            #f"-PmapProjection=\'{projection}\' "
+            f"-PmapProjection=\'{projection}\' "
             f"-t \'{str(outfile)}\' \'{str(infile)}\' "
         )
     else:
@@ -276,6 +276,8 @@ def ls_mask(infile, outfile, logfile, config_dict):
     if 42001 <= dem_dict['out_projection'] <= 97002:
         projection = f"AUTO:{dem_dict['out_projection']}"
     # epsg codes
+    elif int(dem_dict['out_projection']) == 4326:
+        projection = CRS.from_epsg(4326).to_wkt()
     else:
         projection = f"EPSG:{dem_dict['out_projection']}"
 

ost/helpers/settings.py

@@ -240,7 +240,7 @@ def get_gpt():
                            'choices': [i for i in range(3, 27, 2)]},
     'scale_size': {'type': int, 'choices': range(0, 2)},
     'create_ls_mask': {'type': bool},
-    'dem_name': {'type': str, 'choices': ['SRTM 1Sec HGT', 'SRTM 3Sec', 'Aster 1sec GDEM', 'GETASSE30', 'External DEM']},
+    'dem_name': {'type': str, 'choices': ['Copernicus 30m Global DEM', 'Copernicus 90m Global DEM', 'SRTM 1Sec HGT', 'SRTM 3Sec', 'Aster 1sec GDEM', 'GETASSE30', 'External DEM']},
     'dem_file': {'type': str},
     'dem_nodata': {'type': int, 'choices': range(0, 66000)},
     'dem_resampling': {'type': str,

ost/s1/s1scene.py

@@ -12,7 +12,14 @@
 based on a ARD product type.
 """
 
+# ------ bug of rasterio --------
 import os
+if "GDAL_DATA" in list(os.environ.keys()):
+    del os.environ["GDAL_DATA"]
+if "PROJ_LIB" in list(os.environ.keys()):
+    del os.environ["PROJ_LIB"]
+# ------ bug of rasterio --------
+
 import sys
 import json
 import logging