final touch of v0.1

process/4_process_gebco.py

@@ -42,9 +42,6 @@ from config import *
 #    PROCESS
 #=============================================================================
 
-# number of pixels for interpolation (fill NoData at poles)
-interpolation = 1 / RESOLUTION_M * 5000         # TEMPORARY
-
 # heading
 print()
 print('============================================================')
@@ -82,27 +79,3 @@ cleanup(GEBCO_ARCTIC_CLIP_VRT)
 cleanup(GEBCO_ARCTIC_FULL_VRT)
 cleanup(GEBCO_ANTARCTIC_CLIP_VRT)
 cleanup(GEBCO_ANTARCTIC_FULL_VRT)
-
-
-#=============================================================================
-#    NOTES / TODO
-#=============================================================================
-
-"""
-The current solution for creating the polar stereographic GEBCO grids (VRT)
-produces holes on the north and south pole, which are then interpolated. This
-likely comes from creating polar stereographic grids with such large extents
-(low latitudes). Everything larger than down to 70°N / up to 70°S creates a
-hole.
-The interpolation (gdal_fillnodata), however, extrapolates
-the grid around the minmax latitude circle.
-Another solution needs to be found. Normal clip and warp, however, produces
-issues:
-
-- Nearest Neighbor: raster artifacts
-- Cubic (and others): dateline blur
-
-What to do?
-
-- somehow get rid of interpolation
-"""

process/config.py

@@ -144,12 +144,6 @@ ANTARCTIC_EXTENT_LAT= -30.0
 # These are file paths to temporary files that will be created (and cleaned
 # up) during the process. The varibales care used throughout the different
 # scripts.
-# The varibales marked as "DEBUGGING" are for the GEBCO process which needs to
-# be tweaked to avoid the gap at the south and north poles. This is likely due
-# to projecting to polar stereographic with such low latitudes (large
-# extents). Once a nicer solution is found, the variables will be adjusted.
-# This should, however, not influence the general process, but is just an
-# internal change.
 
 # ADD:
 # ADD land with cut rock outcrops
@@ -185,18 +179,14 @@ GEBCO_WORLD = os.path.join(PATH_GEBCO,'GEBCO_world.tif')
 GEBCO_WORLD_HILLSHADE = os.path.join(PATH_GEBCO,'GEBCO_world_hillshade.tif')
 
 # GEBCO Arctic
-#GEBCO_ARCTIC_CLIP = os.path.join(PATH_GEBCO,'GEBCO_arctic_clip.tif') # DEBUGGING
-GEBCO_ARCTIC_CLIP_VRT = os.path.join(PATH_GEBCO,'GEBCO_arctic_clip.vrt') # DEBUGGING
-#GEBCO_ARCTIC_FULL = os.path.join(PATH_GEBCO,'GEBCO_arctic_full.tif') # DEBUGGING
-GEBCO_ARCTIC_FULL_VRT = os.path.join(PATH_GEBCO,'GEBCO_arctic_full.vrt') # DEBUGGING
+GEBCO_ARCTIC_CLIP_VRT = os.path.join(PATH_GEBCO,'GEBCO_arctic_clip.vrt')
+GEBCO_ARCTIC_FULL_VRT = os.path.join(PATH_GEBCO,'GEBCO_arctic_full.vrt')
 GEBCO_ARCTIC = os.path.join(PATH_GEBCO,'GEBCO_arctic.tif')
 GEBCO_ARCTIC_HILLSHADE = os.path.join(PATH_GEBCO,'GEBCO_arctic_hillshade.tif')
 
 # GEBCO Antarctic
-#GEBCO_ANTARCTIC_CLIP = os.path.join(PATH_GEBCO,'GEBCO_antarctic_clip.tif') # DEBUGGING
-GEBCO_ANTARCTIC_CLIP_VRT = os.path.join(PATH_GEBCO,'GEBCO_antarctic_clip.vrt') # DEBUGGING
-#GEBCO_ANTARCTIC_FULL = os.path.join(PATH_GEBCO,'GEBCO_antarctic_full.tif') # DEBUGGING
-GEBCO_ANTARCTIC_FULL_VRT = os.path.join(PATH_GEBCO,'GEBCO_antarctic_full.vrt') # DEBUGGING
+GEBCO_ANTARCTIC_CLIP_VRT = os.path.join(PATH_GEBCO,'GEBCO_antarctic_clip.vrt')
+GEBCO_ANTARCTIC_FULL_VRT = os.path.join(PATH_GEBCO,'GEBCO_antarctic_full.vrt')
 GEBCO_ANTARCTIC = os.path.join(PATH_GEBCO,'GEBCO_antarctic.tif')
 GEBCO_ANTARCTIC_HILLSHADE = os.path.join(PATH_GEBCO,'GEBCO_antarctic_hillshade.tif')