mirror of
https://github.com/qgis/QGIS.git
synced 2025-10-30 00:07:09 -04:00
216 lines
9.5 KiB
Python
216 lines
9.5 KiB
Python
# -*- coding: utf-8 -*-
|
|
|
|
"""
|
|
***************************************************************************
|
|
RasterLayerCalculator.py
|
|
---------------------
|
|
Date : November 2016
|
|
Copyright : (C) 2016 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__ = 'November 2016'
|
|
__copyright__ = '(C) 2016, Victor Olaya'
|
|
|
|
# This will get replaced with a git SHA1 when you do a git archive
|
|
|
|
__revision__ = '$Format:%H$'
|
|
|
|
import os
|
|
import math
|
|
from processing.algs.qgis.QgisAlgorithm import QgisAlgorithm
|
|
from processing.algs.gdal.GdalUtils import GdalUtils
|
|
from qgis.core import (QgsProcessing,
|
|
QgsProcessingException,
|
|
QgsProcessingUtils,
|
|
QgsProcessingParameterCrs,
|
|
QgsProcessingParameterMultipleLayers,
|
|
QgsProcessingParameterNumber,
|
|
QgsProcessingParameterExtent,
|
|
QgsProcessingParameterRasterDestination,
|
|
QgsProcessingParameterRasterLayer,
|
|
QgsProcessingOutputRasterLayer,
|
|
QgsProcessingParameterString,
|
|
QgsCoordinateTransform)
|
|
from qgis.analysis import QgsRasterCalculator, QgsRasterCalculatorEntry
|
|
|
|
|
|
class RasterCalculator(QgisAlgorithm):
|
|
|
|
LAYERS = 'LAYERS'
|
|
EXTENT = 'EXTENT'
|
|
CELLSIZE = 'CELLSIZE'
|
|
EXPRESSION = 'EXPRESSION'
|
|
CRS = 'CRS'
|
|
OUTPUT = 'OUTPUT'
|
|
|
|
def group(self):
|
|
return self.tr('Raster analysis')
|
|
|
|
def groupId(self):
|
|
return 'rasteranalysis'
|
|
|
|
def __init__(self):
|
|
super().__init__()
|
|
|
|
def initAlgorithm(self, config=None):
|
|
class ParameterRasterCalculatorExpression(QgsProcessingParameterString):
|
|
|
|
def __init__(self, name='', description='', multiLine=False):
|
|
super().__init__(name, description, multiLine=multiLine)
|
|
self.setMetadata({
|
|
'widget_wrapper': 'processing.algs.qgis.ui.RasterCalculatorWidgets.ExpressionWidgetWrapper'
|
|
})
|
|
|
|
def type(self):
|
|
return 'raster_calc_expression'
|
|
|
|
def clone(self):
|
|
return ParameterRasterCalculatorExpression(self.name(), self.description(), self.multiLine())
|
|
|
|
def evaluateForModeler(self, value, model):
|
|
for i in list(model.inputs.values()):
|
|
param = i.param
|
|
if isinstance(param, QgsProcessingParameterRasterLayer):
|
|
new = "{}@".format(os.path.basename(param.value))
|
|
old = "{}@".format(param.name())
|
|
value = value.replace(old, new)
|
|
|
|
for alg in list(model.algs.values()):
|
|
for out in alg.algorithm.outputs:
|
|
if isinstance(out, QgsProcessingOutputRasterLayer):
|
|
if out.value:
|
|
new = "{}@".format(os.path.basename(out.value))
|
|
old = "{}:{}@".format(alg.modeler_name, out.name)
|
|
value = value.replace(old, new)
|
|
return value
|
|
|
|
self.addParameter(ParameterRasterCalculatorExpression(self.EXPRESSION, self.tr('Expression'),
|
|
multiLine=True))
|
|
self.addParameter(QgsProcessingParameterMultipleLayers(self.LAYERS,
|
|
self.tr('Reference layer(s) (used for automated extent, cellsize, and CRS)'),
|
|
layerType=QgsProcessing.TypeRaster,
|
|
optional=True))
|
|
self.addParameter(QgsProcessingParameterNumber(self.CELLSIZE,
|
|
self.tr('Cell size (use 0 or empty to set it automatically)'),
|
|
type=QgsProcessingParameterNumber.Double,
|
|
minValue=0.0, defaultValue=0.0, optional=True))
|
|
self.addParameter(QgsProcessingParameterExtent(self.EXTENT,
|
|
self.tr('Output extent'),
|
|
optional=True))
|
|
self.addParameter(QgsProcessingParameterCrs(self.CRS, 'Output CRS', optional=True))
|
|
self.addParameter(QgsProcessingParameterRasterDestination(self.OUTPUT, self.tr('Output')))
|
|
|
|
def name(self):
|
|
return 'rastercalculator'
|
|
|
|
def displayName(self):
|
|
return self.tr('Raster calculator')
|
|
|
|
def processAlgorithm(self, parameters, context, feedback):
|
|
expression = self.parameterAsString(parameters, self.EXPRESSION, context)
|
|
layers = self.parameterAsLayerList(parameters, self.LAYERS, context)
|
|
|
|
layersDict = {}
|
|
if layers:
|
|
layersDict = {os.path.basename(lyr.source().split(".")[0]): lyr for lyr in layers}
|
|
|
|
crs = self.parameterAsCrs(parameters, self.CRS, context)
|
|
if not layers and not crs.isValid():
|
|
raise QgsProcessingException(self.tr("No reference layer selected nor CRS provided"))
|
|
|
|
if not crs.isValid() and layers:
|
|
crs = list(layersDict.values())[0].crs()
|
|
|
|
bbox = self.parameterAsExtent(parameters, self.EXTENT, context)
|
|
if not layers and bbox.isNull():
|
|
raise QgsProcessingException(self.tr("No reference layer selected nor extent box provided"))
|
|
|
|
if not bbox.isNull():
|
|
bboxCrs = self.parameterAsExtentCrs(parameters, self.EXTENT, context)
|
|
if bboxCrs != crs:
|
|
transform = QgsCoordinateTransform(bboxCrs, crs, context.project())
|
|
bbox = transform.transformBoundingBox(bbox)
|
|
|
|
if bbox.isNull() and layers:
|
|
bbox = QgsProcessingUtils.combineLayerExtents(layers, crs)
|
|
|
|
cellsize = self.parameterAsDouble(parameters, self.CELLSIZE, context)
|
|
if not layers and cellsize == 0:
|
|
raise QgsProcessingException(self.tr("No reference layer selected nor cellsize value provided"))
|
|
|
|
def _cellsize(layer):
|
|
ext = layer.extent()
|
|
if layer.crs() != crs:
|
|
transform = QgsCoordinateTransform(layer.crs(), crs, context.project())
|
|
ext = transform.transformBoundingBox(ext)
|
|
return (ext.xMaximum() - ext.xMinimum()) / layer.width()
|
|
if cellsize == 0:
|
|
cellsize = min([_cellsize(lyr) for lyr in layersDict.values()])
|
|
|
|
for lyr in QgsProcessingUtils.compatibleRasterLayers(context.project()):
|
|
name = lyr.name()
|
|
if (name + "@") in expression:
|
|
layersDict[name] = lyr
|
|
|
|
entries = []
|
|
for name, lyr in layersDict.items():
|
|
for n in range(lyr.bandCount()):
|
|
ref = '{:s}@{:d}'.format(name, n + 1)
|
|
if ref in expression:
|
|
entry = QgsRasterCalculatorEntry()
|
|
entry.ref = ref
|
|
entry.raster = lyr
|
|
entry.bandNumber = n + 1
|
|
entries.append(entry)
|
|
|
|
output = self.parameterAsOutputLayer(parameters, self.OUTPUT, context)
|
|
|
|
width = math.floor((bbox.xMaximum() - bbox.xMinimum()) / cellsize)
|
|
height = math.floor((bbox.yMaximum() - bbox.yMinimum()) / cellsize)
|
|
driverName = GdalUtils.getFormatShortNameFromFilename(output)
|
|
calc = QgsRasterCalculator(expression,
|
|
output,
|
|
driverName,
|
|
bbox,
|
|
crs,
|
|
width,
|
|
height,
|
|
entries)
|
|
|
|
res = calc.processCalculation(feedback)
|
|
if res == QgsRasterCalculator.ParserError:
|
|
raise QgsProcessingException(self.tr("Error parsing formula"))
|
|
|
|
return {self.OUTPUT: output}
|
|
|
|
def processBeforeAddingToModeler(self, algorithm, model):
|
|
values = []
|
|
expression = algorithm.params[self.EXPRESSION]
|
|
for i in list(model.inputs.values()):
|
|
param = i.param
|
|
if isinstance(param, QgsProcessingParameterRasterLayer) and "{}@".format(param.name) in expression:
|
|
values.append(ValueFromInput(param.name()))
|
|
|
|
if algorithm.name:
|
|
dependent = model.getDependentAlgorithms(algorithm.name)
|
|
else:
|
|
dependent = []
|
|
for alg in list(model.algs.values()):
|
|
if alg.modeler_name not in dependent:
|
|
for out in alg.algorithm.outputs:
|
|
if (isinstance(out, QgsProcessingOutputRasterLayer) and
|
|
"{}:{}@".format(alg.modeler_name, out.name) in expression):
|
|
values.append(ValueFromOutput(alg.modeler_name, out.name))
|
|
|
|
algorithm.params[self.LAYERS] = values
|