t.rast.in.sos

t.rast.in.sos/Makefile

@@ -0,0 +1,7 @@
+MODULE_TOPDIR = ../..
+
+PGM = t.rast.in.sos
+
+include $(MODULE_TOPDIR)/include/Make/Script.make
+
+default: script

t.rast.in.sos/t.rast.in.sos.html

@@ -0,0 +1,46 @@
+<h2>DESCRIPTION</h2>
+
+This module imports data from SOS server to GRASS GIS.
+
+<h2>DESCRIPTION</h2>
+
+<em>t.rast.in.sos</em> imports data from SOS server to GRASS GIS. It
+creates new spatio-temporal map for each observed property of each procedure
+(registered from vector maps created as intermediates).
+
+<h2>EXAMPLES</h2>
+
+<h3>Standard usage:</h3>
+
+<div class="code"><pre>
+    t.rast.in.sos url=http://localhost/istsos/demo? output=workshop \
+    offering=workshop procedure=BELLINZONA \
+    observed_properties=meteo:air:humidity:relative \
+    event_time='2015-05-03 18:28:01+01:00/2015-05-03 19:33:01+01:00'
+</pre></div>
+
+<h3>Info about services:</h3>
+
+<div class="code"><pre>
+    > t.rast.in.sos url=http://localhost/istsos/demo? offering=workshop -p -v
+     Observed properties of workshop offering:
+     urn:ogc:def:parameter:x-istsos:1.0:meteo:air:humidity:relative
+     urn:ogc:def:parameter:x-istsos:1.0:meteo:air:rainfall
+     urn:ogc:def:parameter:x-istsos:1.0:meteo:air:temperature
+     urn:ogc:def:parameter:x-istsos:1.0:meteo:air:wind:velocity
+     Procedures of workshop offering:
+     BELLINZONA
+     V_LUGANO
+</pre></div>
+
+<h2>SEE ALSO</h2>
+
+<em>
+    <a href="v.in.sos.html">v.in.sos</a>
+    <a href="r.in.sos.html">r.in.sos</a>
+    <a href="t.vect.in.sos.html">t.vect.in.sos</a>
+</em>
+
+<h2>AUTHORS</h2>
+
+Ondrej Pesek under the supervision of Luca Delucchi. 

t.rast.in.sos/t.rast.in.sos.py

@@ -0,0 +1,246 @@
+#!/usr/bin/env python
+#
+############################################################################
+#
+# MODULE:	    t.rast.in.sos
+# AUTHOR(S):	Ondrej Pesek <pesej.ondrek@gmail.com>
+# PURPOSE:	    Import data from SOS server as a vector layer to GRASS GIS
+# COPYRIGHT:	(C) 2017 Ondrej Pesek and the GRASS Development Team
+#
+#		This program is free software under the GNU General
+#		Public License (>=v2). Read the file COPYING that
+#		comes with GRASS for details.
+#
+#############################################################################
+
+#%module
+#% description: Import data from SOS server as space temporal maps to GRASS.
+#% keyword: raster
+#% keyword: import
+#% keyword: SOS
+#%end
+#%flag
+#% key: v
+#% description: Print observed properties for given url and offering
+#% guisection: SOS description
+#%end
+#%flag
+#% key: o
+#% description: Print offerings for given url
+#% guisection: SOS description
+#%end
+#%flag
+#% key: p
+#% description: Print procedures for given url and offering
+#% guisection: SOS description
+#%end
+#%flag
+#% key: t
+#% description: Print begin and end timestamps for given url and offering
+#% guisection: SOS description
+#%end
+#%flag
+#% key: g
+#% description: Print informations in shell script style
+#% guisection: SOS description
+#%end
+#%option
+#% key: url
+#% type: string
+#% description: Base URL starting with 'http' and ending in '?'
+#% required: yes
+#%end
+#%option
+#% key: output
+#% type: string
+#% description: Prefix for output maps
+#% required: yes
+#% guisection: Request
+#%end
+#%option
+#% key: offering
+#% type: string
+#% description: A collection of sensor used to conveniently group them up
+#% required: no
+#% multiple: yes
+#% guisection: Request
+#%end
+#%option
+#% key: response_format
+#% type: string
+#% options: text/xml;subtype="om/1.0.0", application/json
+#% description: Format of data output
+#% answer: text/xml;subtype="om/1.0.0"
+#% required: no
+#% guisection: Request
+#%end
+#%option
+#% key: observed_properties
+#% type: string
+#% description: The phenomena that are observed
+#% required: no
+#% guisection: Request
+#% multiple: yes
+#%end
+#%option
+#% key: procedure
+#% type: string
+#% description: Who provide the observations
+#% required: no
+#% guisection: Request
+#%end
+#%option
+#% key: event_time
+#% type: string
+#% label: Timestamp of first/timestamp of last requested observation
+#% description: Exmaple: 2015-06-01T00:00:00+0200/2015-06-03T00:00:00+0200
+#% required: no
+#% guisection: Request
+#%end
+#%option
+#% key: version
+#% type: string
+#% description: Version of SOS server
+#% guisection: Request
+#% options: 1.0.0, 2.0.0
+#% answer: 1.0.0
+#%end
+#%option
+#% key: username
+#% type: string
+#% description: Username with access to server
+#% guisection: User
+#%end
+#%option
+#% key: password
+#% type: string
+#% description: Password according to username
+#% guisection: User
+#%end
+#%rules
+#% requires_all: -v, offering, url
+#% requires_all: -p, offering, url
+#% requires_all: -t, offering, url
+#% requires: -o, url
+#%end
+
+
+import sys
+import json
+from sqlite3 import OperationalError
+try:
+    from owslib.sos import SensorObservationService
+    from grass.script import parser, run_command, overwrite
+    from grass.script import core as grass
+    from grass.script import vector
+    from grass.pygrass.vector import VectorTopo
+    from grass.pygrass.vector.geometry import Point
+    from grass.pygrass.vector.table import Link
+    import grass.temporal as tgis
+except ImportError as e:
+    sys.stderr.write('Error importing internal libs. '
+                     'Did you run the script from GRASS GIS?\n')
+    raise(e)
+
+sys.path.append('/home/ondrej/workspace/GRASS-GIS-SOS-tools/format_conversions')
+# TODO: Incorporate format conversions into OWSLib and don't use absolute path
+from xml2geojson import xml2geojson
+from json2geojson import json2geojson
+
+
+def cleanup():
+    pass
+
+
+def main():
+
+    fl = 'f'
+    for f, val in flags.iteritems():
+        if val is True:
+            fl += f
+
+    run_command('r.in.sos', flags=fl, **options)
+    if any(value is True and key in [
+        'o', 'v', 'p', 't'] for key, value in flags.iteritems()):
+        return 0
+
+    service = SensorObservationService(options['url'],
+                                       version=options['version'])
+
+    for offering in options['offering'].split(','):
+        procedure, observed_properties, event_time = handle_not_given_options(
+            service, offering)
+        for observedProperty in observed_properties.split(','):
+            mapName = '{}_{}_{}'.format(options['output'], offering,
+                                        observedProperty)
+            if ':' in mapName:
+                mapName = '_'.join(mapName.split(':'))
+            if '-' in mapName:
+                mapName = '_'.join(mapName.split('-'))
+            if '.' in mapName:
+                mapName = '_'.join(mapName.split('.'))
+
+            mapsListFile = get_maps(mapName)
+            create_temporal(mapsListFile, mapName)
+
+    return 0
+
+
+def get_maps(mapName):
+    tmpFile = grass.tempfile()
+    run_command('g.list', type='raster',
+                pattern='{}_*'.format(mapName),
+                output=tmpFile)
+
+    return tmpFile
+
+
+def create_temporal(mapsListFile, mapName):
+
+    run_command('t.create',
+                output=mapName,
+                type='strds',
+                title='Dataset for offering {} and observed '
+                      'property {}'.format(mapName.split('_')[1],
+                                           '_'.join(mapName.split('_')[2:])),
+                description='Raster space time dataset')
+
+    with open(mapsListFile, 'r') as maps:
+        for rasterMap in maps.readlines():
+            a = rasterMap.split('t')[-1]
+            mapTimestamp = '{}-{}-{} {}:{}'.format(a[0:4], a[4:6], a[6:8],
+                                                   a[9:11], a[11:13])
+            run_command('t.register',
+                        type='raster',
+                        input=mapName,
+                        maps='{}'.format(rasterMap.strip()),
+                        start=mapTimestamp)
+
+
+def handle_not_given_options(service, offering=None):
+    # DUPLICATED: Also in v.in.sos
+    if options['procedure'] == '':
+        procedure = None
+    else:
+        procedure = options['procedure']
+
+    if options['observed_properties'] == '':
+        observed_properties = ''
+        for observed_property in service[offering].observed_properties:
+            observed_properties += '{},'.format(observed_property)
+        observed_properties = observed_properties[:-1]
+    else:
+        observed_properties = options['observed_properties']
+
+    if options['event_time'] == '':
+        event_time = '{}/{}'.format(service[offering].begin_position,
+                                    service[offering].end_position)
+    else:
+        event_time = options['event_time']
+
+    return procedure, observed_properties, event_time
+
+
+if __name__ == "__main__":
+    options, flags = parser()
+    main()