Skip to content
Snippets Groups Projects
Select Git revision
  • master default protected
1 result
Blame Permalink
3_process_glims.py 2.48 KiB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66 






# -*- coding: utf-8 -*-

"""
####################################################################
#                                                                  #
#                           AWI Basemap                            #
#       PROCESS GLOBAL LAND ICE MEASUREMENTS FROM SPACE DATA       #
#                                                                  #
####################################################################

This script is processing the Global Land Ice Measurements from
Space (GLIMS) input for the AWI Basemap. The GLIMS geometries are
fixed (with zero buffer) and the global dataset is clipped for the
Arctic and the Antarctic.

"""


#=============================================================================
#    SCRIPT INFO
#=============================================================================

__author__ = 'Simon Dreutter'
__version__ = '0.1'
__date__ = '2020-12-17'
__email__ = 'simon.dreutter@awi.de'
__status__ = 'Developement'


#=============================================================================
#    IMPORT
#=============================================================================

import os

# import AWI Basemap configuration
from config import *


#=============================================================================
#    PROCESS
#=============================================================================

# layer names
glims_in_name = os.path.splitext(os.path.basename(GLIMS_IN))[0]
glims_world_name = os.path.splitext(os.path.basename(GLIMS_WORLD))[0]
glims_arctic_name = os.path.splitext(os.path.basename(GLIMS_ARCTIC))[0]
glims_antarctic_name = os.path.splitext(os.path.basename(GLIMS_ANTARCTIC))[0]

# heading
print()
print('============================================================')
print('      GLOBAL LAND ICE MEASUREMENTS FROM SPACE / GLIMS       ')

msg = 'Fixing GLIMS geometries...'
cmd = f'ogr2ogr {GLIMS_WORLD} {GLIMS_IN} -dialect sqlite -sql "SELECT ST_Buffer(geometry, 0.0) AS geometry,* FROM \"{glims_in_name}\"" -nln "{glims_world_name}" -f "GPKG"'
run(msg,cmd)

msg = 'Clipping GLIMS for the Arctic...'
cmd = f'ogr2ogr -spat -180.0 90.0 180.0 {ARCTIC_EXTENT_LAT} -clipsrc spat_extent {GLIMS_ARCTIC} {GLIMS_WORLD} {glims_world_name} -nln "{glims_arctic_name}" -f "GPKG"'
run(msg,cmd)

msg = 'Clipping GLIMS for the Antarctic...'
cmd = f'ogr2ogr -spat -180.0 {ANTARCTIC_EXTENT_LAT} 180.0 -90.0 -clipsrc spat_extent {GLIMS_ANTARCTIC} {GLIMS_WORLD} {glims_world_name} -nln "{glims_antarctic_name}" -f "GPKG"'
run(msg,cmd)