diff --git a/process/1_process_add_backup1.py b/process/1_process_add_backup1.py deleted file mode 100644 index 6617b7d42a86c6705430e0c57f7edb415d413afb..0000000000000000000000000000000000000000 --- a/process/1_process_add_backup1.py +++ /dev/null @@ -1,77 +0,0 @@ -# -*- coding: utf-8 -*- - -""" -#################################################################### -# # -# AWI Basemap # -# PROCESS ANTARCTIC DIGITAL DATANASE DATA # -# # -#################################################################### - -This script is processing the Antarctic Digital Database (ADD) input -for the AWI Basemap. Input layers are the ADD coastlines and the ADD -rock outcrops. Both layers are reprojected to a global CRS -(EPSG_WORLD), coastlines are separated by shelf ice and ice sheet, -and the rock outcrop geometries are fixed. - -""" - - -#============================================================================= -# SCRIPT INFO -#============================================================================= - -__author__ = 'Simon Dreutter' -__version__ = '0.1' -__date__ = '2020-12-17' -__email__ = 'simon.dreutter@awi.de' -__status__ = 'Developement' - - -#============================================================================= -# IMPORT -#============================================================================= - -# import AWI Basemap configuration -from config import * - - -#============================================================================= -# PROCESS -#============================================================================= - -EPSG_ADD = 'EPSG:3031' - -msg = 'Reprojecting and extracting land mask from ADD coastlines...' -cmd = f'ogr2ogr -s_srs {EPSG_ADD} -t_srs {EPSG_WORLD} -nln "add_coastline_land" -where "\"surface\" = \'land\'" -wrapdateline -f GPKG {ADD_COASTLINE_LAND} {ADD_COASTLINE_IN}' -run(msg,cmd) - -msg = 'Reprojecting and extracting ice mask from ADD coastlines...' -cmd = f'ogr2ogr -s_srs {EPSG_ADD} -t_srs {EPSG_WORLD} -nln "add_coastline_ice" -where "\"surface\" = \'ice shelf\' OR \"surface\" = \'rumple\' OR \"surface\" = \'ocean\' OR \"surface\" = \'ice tongue\'" -wrapdateline -f GPKG {ADD_SHELF_ICE} {ADD_COASTLINE_IN}' -run(msg,cmd) - -msg = 'Reprojecting ADD rock outcrops...' -cmd = f'ogr2ogr -s_srs {EPSG_ADD} -t_srs {EPSG_WORLD} -nln "add_rock_outcrop" -wrapdateline -f GPKG {ADD_ROCK_OUTCROP} {ADD_ROCK_OUTCROP_IN}' -run(msg,cmd) - -msg = 'Fixing ADD rock outcrops geometries...' -cmd = f'ogr2ogr -dialect sqlite -sql "SELECT ST_Buffer(geom, 0.0) AS geom,* FROM \"add_rock_outcrop\"" -nln "add_rock_outcrop_fixed" -f GPKG {ADD_ROCK_OUTCROP_FIXED} {ADD_ROCK_OUTCROP}' -run(msg,cmd) - -print() -print('ADD done!') -print() -print('Now do not forget to manually cut the rock outcrops from the ADD land layer!') -print() -input('Press Enter to exit!') -print() - - -#============================================================================= -# TODO -#============================================================================= - -""" -- find a scripted solution for 'Difference' -- once found, add cleanup -""" \ No newline at end of file