File: r_stats_quantile_rast.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,431 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_stats_quantile_rast.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$'

from processing.core.parameters import getParameterFromString
import os


def processCommand(alg):
    # We create the output sequence according to percentiles number
    base = alg.getParameterValue('base')
    quantiles = alg.getParameterValue('quantiles') - 1
    outputs = []
    for i in range(0, int(quantiles)):
        outputs.append('output_{}'.format(i))

    output = getParameterFromString('ParameterString|output|Output Rasters|None|False|True')
    output.value = ','.join(outputs)
    alg.addParameter(output)

    output_dir = alg.getOutputFromName('output_dir')
    alg.removeOutputFromName('output_dir')

    # Launch the algorithm
    alg.processCommand()

    # We re-add the previous output
    alg.addOutput(output_dir)


def processOutputs(alg):
    # We need to export each of the output
    output_dir = alg.getOutputValue('output_dir')
    outputParam = alg.getParameterFromName('output')
    outputs = outputParam.value.split(',')
    alg.parameters.remove(outputParam)
    for output in outputs:
        command = u"r.out.gdal -c createopt=\"TFW=YES,COMPRESS=LZW\" input={} output=\"{}\" --overwrite".format(
            output,
            os.path.join(output_dir, output + '.tif')
        )
        alg.commands.append(command)
        alg.outputCommands.append(command)