diff --git a/process/4_process_gebco.py b/process/4_process_gebco.py index 15fd5c71b94c1c3d5fefff152ffa9c63ea0c358b..f3e25080cfb6106b03dac12fadf020a9bdc22c9a 100644 --- a/process/4_process_gebco.py +++ b/process/4_process_gebco.py @@ -56,12 +56,20 @@ print('====================================================================') print(f'{"GENERAL BATHYMETRIC CHART OF THE OCEANS / GEBCO":^68}') # World + +# DEBUG +if RESOLUTION_DEG == (1 / 60 / 60 * 15): + tr = '' +else: + tr = f' -tr {RESOLUTION_DEG} {RESOLUTION_DEG}' +#/DEBUG + # ice surface msg = 'Resampling GEBCO ice surface...' cmd = (f'gdalwarp' f' -s_srs EPSG:{EPSG_WORLD}' f' -t_srs EPSG:{EPSG_WORLD}' - f' -tr {RESOLUTION_DEG} {RESOLUTION_DEG}' + f'{tr}' f' -r cubic' f' -of GTiff' f' {WARP_OPTIONS}' @@ -100,7 +108,7 @@ msg = 'Resampling GEBCO sub ice...' cmd = (f'gdalwarp' f' -s_srs EPSG:{EPSG_WORLD}' f' -t_srs EPSG:{EPSG_WORLD}' - f' -tr {RESOLUTION_DEG} {RESOLUTION_DEG}' + f'{tr}' f' -r cubic' f' -of GTiff' f' {WARP_OPTIONS}' diff --git a/process/config.py b/process/config.py index 063176fc108fdcdbf47d804790cd976ba03b438f..f1413631b476bf7b184c912bc6a9abeaf7708cfd 100644 --- a/process/config.py +++ b/process/config.py @@ -194,8 +194,8 @@ GLIMS_IN = PATH_GLIMS / 'glims_polygons.shp' # warping in the lower latitudes. # Final AWI Basemap (World) resolution -RESOLUTION_DEG = 0.004166666666666666 # 15 arc seconds (original GEBCO resolution) -# RESOLUTION_DEG = 0.008333333333333333 # 30 arc seconds +RESOLUTION_DEG = 1 / 60 / 60 * 15 # 15 arc seconds (original GEBCO resolution) +# RESOLUTION_DEG = 1 / 60 / 60 * 30 # 30 arc seconds # Final AWI Basemap (Arctic and Antarctic) resolution RESOLUTION_M = 500.0 # approximate pendant to 15 arc seconds