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
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
|
# -*- coding: utf-8 -*-
"""
***************************************************************************
proximity.py
---------------------
Date : August 2012
Copyright : (C) 2012 by Victor Olaya
Email : volayaf at gmail dot com
***************************************************************************
* *
* 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__ = 'Victor Olaya'
__date__ = 'August 2012'
__copyright__ = '(C) 2012, Victor Olaya'
# This will get replaced with a git SHA1 when you do a git archive
__revision__ = '$Format:%H$'
from processing.algs.gdal.GdalAlgorithm import GdalAlgorithm
from processing.parameters.ParameterRaster import ParameterRaster
from processing.parameters.ParameterString import ParameterString
from processing.parameters.ParameterSelection import ParameterSelection
from processing.parameters.ParameterNumber import ParameterNumber
from processing.outputs.OutputRaster import OutputRaster
from processing.tools.system import *
from processing.algs.gdal.GdalUtils import GdalUtils
class proximity(GdalAlgorithm):
INPUT = 'INPUT'
VALUES = 'VALUES'
UNITS = 'UNITS'
MAX_DIST = 'MAX_DIST'
NODATA = 'NODATA'
BUF_VAL = 'BUF_VAL'
OUTPUT = 'OUTPUT'
DISTUNITS = ['GEO', 'PIXEL']
def commandLineName(self):
return "gdalogr:proximity"
def defineCharacteristics(self):
self.name = 'Proximity (raster distance)'
self.group = '[GDAL] Analysis'
self.addParameter(ParameterRaster(self.INPUT, 'Input layer', False))
self.addParameter(ParameterString(self.VALUES, 'Values', ''))
self.addParameter(ParameterSelection(self.UNITS, 'Dist units',
self.DISTUNITS, 0))
self.addParameter(ParameterNumber(self.MAX_DIST,
'Max dist (negative value to ignore)', -1, 9999, -1))
self.addParameter(ParameterNumber(self.NODATA,
'No data (negative value to ignore)', -1, 9999, -1))
self.addParameter(ParameterNumber(self.BUF_VAL,
'Fixed buf val (negative value to ignore)', -1,
9999, -1))
self.addOutput(OutputRaster(self.OUTPUT, 'Output layer'))
def processAlgorithm(self, progress):
output = self.getOutputValue(self.OUTPUT)
arguments = []
arguments.append(self.getParameterValue(self.INPUT))
arguments.append(output)
arguments.append('-of')
arguments.append(GdalUtils.getFormatShortNameFromFilename(output))
arguments.append('-distunits')
arguments.append(self.DISTUNITS[self.getParameterValue(self.UNITS)])
values = self.getParameterValue(self.VALUES)
if len(values) > 0:
arguments.append('-values')
arguments.append(values)
values = str(self.getParameterValue(self.MAX_DIST))
if values < 0:
arguments.append('-maxdist')
arguments.append(values)
values = str(self.getParameterValue(self.NODATA))
if values < 0:
arguments.append('-nodata')
arguments.append(values)
values = str(self.getParameterValue(self.BUF_VAL))
if values < 0:
arguments.append('-fixed-buf-val')
arguments.append(values)
commands = []
if isWindows():
commands = ['cmd.exe', '/C ', 'gdal_proximity.bat',
GdalUtils.escapeAndJoin(arguments)]
else:
commands = ['gdal_proximity.py',
GdalUtils.escapeAndJoin(arguments)]
GdalUtils.runGdal(commands, progress)
|
