|
| 1 | +# -*- coding: utf-8 -*- |
| 2 | + |
| 3 | +""" |
| 4 | +*************************************************************************** |
| 5 | + InterpolationDataWidget.py |
| 6 | + --------------------- |
| 7 | + Date : December 2016 |
| 8 | + Copyright : (C) 2016 by Alexander Bruy |
| 9 | + Email : alexander dot bruy at gmail dot com |
| 10 | +*************************************************************************** |
| 11 | +* * |
| 12 | +* This program is free software; you can redistribute it and/or modify * |
| 13 | +* it under the terms of the GNU General Public License as published by * |
| 14 | +* the Free Software Foundation; either version 2 of the License, or * |
| 15 | +* (at your option) any later version. * |
| 16 | +* * |
| 17 | +*************************************************************************** |
| 18 | +""" |
| 19 | + |
| 20 | +__author__ = 'Alexander Bruy' |
| 21 | +__date__ = 'December 2016' |
| 22 | +__copyright__ = '(C) 2016, Alexander Bruy' |
| 23 | + |
| 24 | +# This will get replaced with a git SHA1 when you do a git archive |
| 25 | + |
| 26 | +__revision__ = '$Format:%H$' |
| 27 | + |
| 28 | +import os |
| 29 | + |
| 30 | +from qgis.PyQt import uic |
| 31 | +from qgis.PyQt.QtCore import pyqtSignal |
| 32 | +from qgis.PyQt.QtWidgets import (QTreeWidgetItem, |
| 33 | + QComboBox |
| 34 | + ) |
| 35 | +from qgis.core import (QgsApplication, |
| 36 | + QgsMapLayer, |
| 37 | + QgsMapLayerProxyModel, |
| 38 | + QgsWkbTypes, |
| 39 | + QgsRectangle, |
| 40 | + QgsReferencedRectangle, |
| 41 | + QgsCoordinateReferenceSystem, |
| 42 | + QgsProcessingUtils, |
| 43 | + QgsProcessingParameterNumber, |
| 44 | + QgsProcessingParameterDefinition |
| 45 | + ) |
| 46 | +from qgis.core import QgsFieldProxyModel |
| 47 | +from qgis.analysis import QgsInterpolator |
| 48 | + |
| 49 | +from processing.gui.wrappers import WidgetWrapper, DIALOG_STANDARD |
| 50 | +from processing.tools import dataobjects |
| 51 | + |
| 52 | +pluginPath = os.path.dirname(__file__) |
| 53 | + |
| 54 | + |
| 55 | +class ParameterInterpolationData(QgsProcessingParameterDefinition): |
| 56 | + |
| 57 | + def __init__(self, name='', description=''): |
| 58 | + super().__init__(name, description) |
| 59 | + self.setMetadata({ |
| 60 | + 'widget_wrapper': 'processing.algs.qgis.ui.InterpolationWidgets.InterpolationDataWidgetWrapper' |
| 61 | + }) |
| 62 | + |
| 63 | + def type(self): |
| 64 | + return 'idw_interpolation_data' |
| 65 | + |
| 66 | + def clone(self): |
| 67 | + return ParameterInterpolationData(self.name(), self.description()) |
| 68 | + |
| 69 | + @staticmethod |
| 70 | + def parseValue(value): |
| 71 | + if value is None: |
| 72 | + return None |
| 73 | + |
| 74 | + if value == '': |
| 75 | + return None |
| 76 | + |
| 77 | + if isinstance(value, str): |
| 78 | + return value if value != '' else None |
| 79 | + else: |
| 80 | + return ParameterInterpolationData.dataToString(value) |
| 81 | + |
| 82 | + @staticmethod |
| 83 | + def dataToString(data): |
| 84 | + s = '' |
| 85 | + for c in data: |
| 86 | + s += '{}::~::{}::~::{:d}::~::{:d};'.format(c[0], |
| 87 | + c[1], |
| 88 | + c[2], |
| 89 | + c[3]) |
| 90 | + return s[:-1] |
| 91 | + |
| 92 | + |
| 93 | +WIDGET, BASE = uic.loadUiType(os.path.join(pluginPath, 'interpolationdatawidgetbase.ui')) |
| 94 | + |
| 95 | + |
| 96 | +class InterpolationDataWidget(BASE, WIDGET): |
| 97 | + |
| 98 | + hasChanged = pyqtSignal() |
| 99 | + |
| 100 | + def __init__(self): |
| 101 | + super(InterpolationDataWidget, self).__init__(None) |
| 102 | + self.setupUi(self) |
| 103 | + |
| 104 | + self.btnAdd.setIcon(QgsApplication.getThemeIcon('/symbologyAdd.svg')) |
| 105 | + self.btnRemove.setIcon(QgsApplication.getThemeIcon('/symbologyRemove.svg')) |
| 106 | + |
| 107 | + self.btnAdd.clicked.connect(self.addLayer) |
| 108 | + self.btnRemove.clicked.connect(self.removeLayer) |
| 109 | + |
| 110 | + self.cmbLayers.layerChanged.connect(self.layerChanged) |
| 111 | + self.cmbLayers.setFilters(QgsMapLayerProxyModel.VectorLayer) |
| 112 | + self.cmbFields.setFilters(QgsFieldProxyModel.Numeric) |
| 113 | + self.cmbFields.setLayer(self.cmbLayers.currentLayer()) |
| 114 | + |
| 115 | + def addLayer(self): |
| 116 | + layer = self.cmbLayers.currentLayer() |
| 117 | + |
| 118 | + attribute = '' |
| 119 | + if self.chkUseZCoordinate.isChecked(): |
| 120 | + attribute = 'Z_COORD' |
| 121 | + else: |
| 122 | + attribute = self.cmbFields.currentField() |
| 123 | + |
| 124 | + self._addLayerData(layer.name(), attribute) |
| 125 | + self.hasChanged.emit() |
| 126 | + |
| 127 | + def removeLayer(self): |
| 128 | + item = self.layersTree.currentItem() |
| 129 | + if not item: |
| 130 | + return |
| 131 | + self.layersTree.invisibleRootItem().removeChild(item) |
| 132 | + self.hasChanged.emit() |
| 133 | + |
| 134 | + def layerChanged(self, layer): |
| 135 | + self.chkUseZCoordinate.setEnabled(False) |
| 136 | + self.chkUseZCoordinate.setChecked(False) |
| 137 | + |
| 138 | + if layer is None or not layer.isValid(): |
| 139 | + return |
| 140 | + |
| 141 | + provider = layer.dataProvider() |
| 142 | + if not provider: |
| 143 | + return |
| 144 | + |
| 145 | + if QgsWkbTypes.hasZ(provider.wkbType()): |
| 146 | + self.chkUseZCoordinate.setEnabled(True) |
| 147 | + |
| 148 | + self.cmbFields.setLayer(layer) |
| 149 | + |
| 150 | + def _addLayerData(self, layerName, attribute): |
| 151 | + item = QTreeWidgetItem() |
| 152 | + item.setText(0, layerName) |
| 153 | + item.setText(1, attribute) |
| 154 | + self.layersTree.addTopLevelItem(item) |
| 155 | + |
| 156 | + comboBox = QComboBox() |
| 157 | + comboBox.addItem(self.tr('Points')) |
| 158 | + comboBox.addItem(self.tr('Structure lines')) |
| 159 | + comboBox.addItem(self.tr('Break lines')) |
| 160 | + comboBox.setCurrentIndex(0) |
| 161 | + self.layersTree.setItemWidget(item, 2, comboBox) |
| 162 | + |
| 163 | + def setValue(self, value): |
| 164 | + self.layersTree.clear() |
| 165 | + rows = value.split(';') |
| 166 | + for i, r in enumerate(rows): |
| 167 | + v = r.split('::~::') |
| 168 | + self.addLayerData(v[0], v[1]) |
| 169 | + |
| 170 | + comboBox = self.layersTree.itemWidget(self.layersTree.topLevelItem(i), 2) |
| 171 | + comboBox.setCurrentIndex(comboBox.findText(v[3])) |
| 172 | + self.hasChanged.emit() |
| 173 | + |
| 174 | + def value(self): |
| 175 | + layers = '' |
| 176 | + context = dataobjects.createContext() |
| 177 | + for i in range(self.layersTree.topLevelItemCount()): |
| 178 | + item = self.layersTree.topLevelItem(i) |
| 179 | + if item: |
| 180 | + layerName = item.text(0) |
| 181 | + layer = QgsProcessingUtils.mapLayerFromString(layerName, context) |
| 182 | + if not layer: |
| 183 | + continue |
| 184 | + |
| 185 | + provider = layer.dataProvider() |
| 186 | + if not provider: |
| 187 | + continue |
| 188 | + |
| 189 | + interpolationAttribute = item.text(1) |
| 190 | + interpolationSource = QgsInterpolator.ValueAttribute |
| 191 | + if interpolationAttribute == 'Z_COORD': |
| 192 | + interpolationSource = QgsInterpolator.ValueZ |
| 193 | + fieldIndex = -1 |
| 194 | + else: |
| 195 | + fieldIndex = layer.fields().indexFromName(interpolationAttribute) |
| 196 | + |
| 197 | + comboBox = self.layersTree.itemWidget(self.layersTree.topLevelItem(i), 2) |
| 198 | + inputTypeName = comboBox.currentText() |
| 199 | + if inputTypeName == self.tr('Points'): |
| 200 | + inputType = QgsInterpolator.SourcePoints |
| 201 | + elif inputTypeName == self.tr('Structure lines'): |
| 202 | + inputType = QgsInterpolator.SourceStructureLines |
| 203 | + else: |
| 204 | + inputType = QgsInterpolator.SourceBreakLines |
| 205 | + |
| 206 | + layers += '{}::~::{:d}::~::{:d}::~::{:d};'.format(layer.source(), |
| 207 | + interpolationSource, |
| 208 | + fieldIndex, |
| 209 | + inputType) |
| 210 | + return layers[:-1] |
| 211 | + |
| 212 | + |
| 213 | +class InterpolationDataWidgetWrapper(WidgetWrapper): |
| 214 | + |
| 215 | + def createWidget(self): |
| 216 | + widget = InterpolationDataWidget() |
| 217 | + widget.hasChanged.connect(lambda: self.widgetValueHasChanged.emit(self)) |
| 218 | + return widget |
| 219 | + |
| 220 | + def setValue(self, value): |
| 221 | + self.widget.setValue(value) |
| 222 | + |
| 223 | + def value(self): |
| 224 | + return self.widget.value() |
| 225 | + |
| 226 | + |
| 227 | +class ParameterPixelSize(QgsProcessingParameterNumber): |
| 228 | + |
| 229 | + def __init__(self, name='', description='', layersData=None, extent=None, minValue=None, default=None, optional=False): |
| 230 | + QgsProcessingParameterNumber.__init__(self, name, description, QgsProcessingParameterNumber.Double, default, optional, minValue) |
| 231 | + self.setMetadata({ |
| 232 | + 'widget_wrapper': 'processing.algs.qgis.ui.InterpolationWidgets.PixelSizeWidgetWrapper' |
| 233 | + }) |
| 234 | + |
| 235 | + self.layersData = layersData |
| 236 | + self.extent = extent |
| 237 | + self.layers = [] |
| 238 | + |
| 239 | + def clone(self): |
| 240 | + copy = ParameterPixelSize(self.name(), self.description(), self.layersData, self.extent, self.minimum(), self.defaultValue(), self.flags() & QgsProcessingParameterDefinition.FlagOptional) |
| 241 | + return copy |
| 242 | + |
| 243 | + |
| 244 | +WIDGET, BASE = uic.loadUiType(os.path.join(pluginPath, 'RasterResolutionWidget.ui')) |
| 245 | + |
| 246 | + |
| 247 | +class PixelSizeWidget(BASE, WIDGET): |
| 248 | + |
| 249 | + def __init__(self): |
| 250 | + super(PixelSizeWidget, self).__init__(None) |
| 251 | + self.setupUi(self) |
| 252 | + self.context = dataobjects.createContext() |
| 253 | + |
| 254 | + self.extent = QgsRectangle() |
| 255 | + self.layers = [] |
| 256 | + |
| 257 | + self.mCellXSpinBox.setShowClearButton(False) |
| 258 | + self.mCellYSpinBox.setShowClearButton(False) |
| 259 | + self.mRowsSpinBox.setShowClearButton(False) |
| 260 | + self.mColumnsSpinBox.setShowClearButton(False) |
| 261 | + |
| 262 | + self.mCellYSpinBox.valueChanged.connect(self.mCellXSpinBox.setValue) |
| 263 | + self.mCellXSpinBox.valueChanged.connect(self.pixelSizeChanged) |
| 264 | + self.mRowsSpinBox.valueChanged.connect(self.rowsChanged) |
| 265 | + self.mColumnsSpinBox.valueChanged.connect(self.columnsChanged) |
| 266 | + |
| 267 | + def setLayers(self, layersData): |
| 268 | + self.extent = QgsRectangle() |
| 269 | + self.layers = [] |
| 270 | + for row in layersData.split(';'): |
| 271 | + v = row.split('::~::') |
| 272 | + # need to keep a reference until interpolation is complete |
| 273 | + layer = QgsProcessingUtils.variantToSource(v[0], self.context) |
| 274 | + if layer: |
| 275 | + self.layers.append(layer) |
| 276 | + bbox = layer.sourceExtent() |
| 277 | + if self.extent.isEmpty(): |
| 278 | + self.extent = bbox |
| 279 | + else: |
| 280 | + self.extent.combineExtentWith(bbox) |
| 281 | + |
| 282 | + self.pixelSizeChanged() |
| 283 | + |
| 284 | + def setExtent(self, extent): |
| 285 | + if extent is not None: |
| 286 | + tokens = extent.split(' ')[0].split(',') |
| 287 | + ext = QgsRectangle(float(tokens[0]), float(tokens[2]), float(tokens[1]), float(tokens[3])) |
| 288 | + if len(tokens) > 1: |
| 289 | + self.extent = QgsReferencedRectangle(ext, QgsCoordinateReferenceSystem(tokens[1][1:-1])) |
| 290 | + else: |
| 291 | + self.extent = ext |
| 292 | + self.pixelSizeChanged() |
| 293 | + |
| 294 | + def pixelSizeChanged(self): |
| 295 | + cell_size = self.mCellXSpinBox.value() |
| 296 | + if cell_size <= 0: |
| 297 | + return |
| 298 | + |
| 299 | + self.mCellYSpinBox.blockSignals(True) |
| 300 | + self.mCellYSpinBox.setValue(cell_size) |
| 301 | + self.mCellYSpinBox.blockSignals(False) |
| 302 | + rows = max(round(self.extent.height() / cell_size) + 1, 1) |
| 303 | + cols = max(round(self.extent.width() / cell_size) + 1, 1) |
| 304 | + self.mRowsSpinBox.blockSignals(True) |
| 305 | + self.mRowsSpinBox.setValue(rows) |
| 306 | + self.mRowsSpinBox.blockSignals(False) |
| 307 | + self.mColumnsSpinBox.blockSignals(True) |
| 308 | + self.mColumnsSpinBox.setValue(cols) |
| 309 | + self.mColumnsSpinBox.blockSignals(False) |
| 310 | + |
| 311 | + def rowsChanged(self): |
| 312 | + rows = self.mRowsSpinBox.value() |
| 313 | + if rows <= 0: |
| 314 | + return |
| 315 | + cell_size = self.extent.height() / rows |
| 316 | + cols = max(round(self.extent.width() / cell_size) + 1, 1) |
| 317 | + self.mColumnsSpinBox.blockSignals(True) |
| 318 | + self.mColumnsSpinBox.setValue(cols) |
| 319 | + self.mColumnsSpinBox.blockSignals(False) |
| 320 | + for w in [self.mCellXSpinBox, self.mCellYSpinBox]: |
| 321 | + w.blockSignals(True) |
| 322 | + w.setValue(cell_size) |
| 323 | + w.blockSignals(False) |
| 324 | + |
| 325 | + def columnsChanged(self): |
| 326 | + cols = self.mColumnsSpinBox.value() |
| 327 | + if cols < 2: |
| 328 | + return |
| 329 | + cell_size = self.extent.width() / (cols - 1) |
| 330 | + rows = max(round(self.extent.height() / cell_size), 1) |
| 331 | + self.mRowsSpinBox.blockSignals(True) |
| 332 | + self.mRowsSpinBox.setValue(rows) |
| 333 | + self.mRowsSpinBox.blockSignals(False) |
| 334 | + for w in [self.mCellXSpinBox, self.mCellYSpinBox]: |
| 335 | + w.blockSignals(True) |
| 336 | + w.setValue(cell_size) |
| 337 | + w.blockSignals(False) |
| 338 | + |
| 339 | + def setValue(self, value): |
| 340 | + try: |
| 341 | + numeric_value = float(value) |
| 342 | + except: |
| 343 | + return False |
| 344 | + |
| 345 | + self.mCellXSpinBox.setValue(numeric_value) |
| 346 | + self.mCellYSpinBox.setValue(numeric_value) |
| 347 | + return True |
| 348 | + |
| 349 | + def value(self): |
| 350 | + return self.mCellXSpinBox.value() |
| 351 | + |
| 352 | + |
| 353 | +class PixelSizeWidgetWrapper(WidgetWrapper): |
| 354 | + |
| 355 | + def __init__(self, param, dialog, row=0, col=0, **kwargs): |
| 356 | + super().__init__(param, dialog, row, col, **kwargs) |
| 357 | + self.context = dataobjects.createContext() |
| 358 | + |
| 359 | + def _panel(self): |
| 360 | + return PixelSizeWidget() |
| 361 | + |
| 362 | + def createWidget(self): |
| 363 | + if self.dialogType == DIALOG_STANDARD: |
| 364 | + return self._panel() |
| 365 | + else: |
| 366 | + w = QgsDoubleSpinBox() |
| 367 | + w.setShowClearButton(False) |
| 368 | + w.setMinimum(0) |
| 369 | + w.setMaximum(99999999999) |
| 370 | + w.setDecimals(6) |
| 371 | + w.setToolTip(self.tr('Resolution of each pixel in output raster, in layer units')) |
| 372 | + return w |
| 373 | + |
| 374 | + def postInitialize(self, wrappers): |
| 375 | + if self.dialogType != DIALOG_STANDARD: |
| 376 | + return |
| 377 | + |
| 378 | + for wrapper in wrappers: |
| 379 | + if wrapper.parameterDefinition().name() == self.param.layersData: |
| 380 | + self.setLayers(wrapper.parameterValue()) |
| 381 | + wrapper.widgetValueHasChanged.connect(self.layersChanged) |
| 382 | + elif wrapper.parameterDefinition().name() == self.param.extent: |
| 383 | + self.setExtent(wrapper.parameterValue()) |
| 384 | + wrapper.widgetValueHasChanged.connect(self.extentChanged) |
| 385 | + |
| 386 | + def layersChanged(self, wrapper): |
| 387 | + self.setLayers(wrapper.parameterValue()) |
| 388 | + |
| 389 | + def setLayers(self, layersData): |
| 390 | + self.widget.setLayers(layersData) |
| 391 | + |
| 392 | + def extentChanged(self, wrapper): |
| 393 | + self.setExtent(wrapper.parameterValue()) |
| 394 | + |
| 395 | + def setExtent(self, extent): |
| 396 | + self.widget.setExtent(extent) |
| 397 | + |
| 398 | + def setValue(self, value): |
| 399 | + return self.widget.setValue(value) |
| 400 | + |
| 401 | + def value(self): |
| 402 | + return self.widget.value() |
0 commit comments