File: r_support.py

package info (click to toggle)
qgis 2.18.28%2Bdfsg-2
  • links: PTS, VCS
  • area: main
  • in suites: buster
  • size: 1,007,948 kB
  • sloc: cpp: 671,774; python: 158,539; xml: 35,690; ansic: 8,346; sh: 1,766; perl: 1,669; sql: 999; yacc: 836; lex: 461; makefile: 292
file content (66 lines) | stat: -rw-r--r-- 2,799 bytes parent folder | download
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 -*-

"""
***************************************************************************
    r_support.py
    ------------
    Date                 : February 2016
    Copyright            : (C) 2016 by Médéric Ribreux
    Email                : medspx at medspx dot fr
***************************************************************************
*                                                                         *
*   This program is free software; you can redistribute it and/or modify  *
*   it under the terms of the GNU General Public License as published by  *
*   the Free Software Foundation; either version 2 of the License, or     *
*   (at your option) any later version.                                   *
*                                                                         *
***************************************************************************
"""

__author__ = 'Médéric Ribreux'
__date__ = 'February 2016'
__copyright__ = '(C) 2016, Médéric Ribreux'

# This will get replaced with a git SHA1 when you do a git archive

__revision__ = '$Format:%H$'

import os


def processCommand(alg):
    # We temporary remove the output
    out = alg.getOutputFromName('output')
    mapParam = alg.getParameterValue('map')
    alg.exportedLayers[out.value] = alg.exportedLayers[mapParam]
    alg.removeOutputFromName('output')
    timestamp = alg.getParameterValue('timestamp')
    if timestamp:
        command = "r.timestamp map={} date='{}'".format(alg.exportedLayers[mapParam],
                                                        timestamp)
        alg.commands.append(command)
        alg.parameters.remove(alg.getParameterFromName('timestamp'))
    alg.processCommand()

    # We re-add the new output
    alg.addOutput(out)


def processOutputs(alg):
    tags = {'title': 'TIFFTAG_DOCUMENTNAME', 'description': 'TIFFTAG_IMAGEDESCRIPTION',
            'creator': 'TIFFTAG_ARTIST', 'timestamp': 'TIFFTAG_DATETIME',
            'source1': 'TIFFTAG_COPYRIGHT', 'units': 'TIFFTAG_RESOLUTIONUNIT',
            'source2': 'GRASS_SOURCE2', 'comments': 'GRASS_HISTORY', 'vdatum': 'GRASS_VDATUM'}
    awk = "awk -F '=' '"
    for support, tag in tags.items():
        awk = '{0} /{1}=".+"/{{ print \"{2}=\"substr($0,{3},length($0) - {3})"," }}'.format(
            awk, support, tag, len(support) + 3)

    # Output results ('from' table and output table)
    out = alg.getOutputValue('output')
    command = u"SDF=$(r.info -e map={} | {}')".format(alg.exportedLayers[out], awk)
    alg.commands.append(command)
    command = u"r.out.gdal --overwrite -c createopt=\"TFW=YES,COMPRESS=LZW\" input={} output=\"{}\" metaopt=\"${{SDF%%,}}\"".format(
        alg.exportedLayers[out], out)
    alg.commands.append(command)
    alg.outputCommands.append(command)