Newer
Older
# -*- coding: utf-8 -*-
"""
/***************************************************************************

Benjamin Jakimow
committed
EO Time Series Viewer
-------------------
begin : 2015-08-20
git sha : $Format:%H$
copyright : (C) 2017 by HU-Berlin
email : benjamin.jakimow@geo.hu-berlin.de
***************************************************************************/
/***************************************************************************
* *
* 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. *
* *
***************************************************************************/
"""
import os, sys, re, fnmatch, collections, copy, traceback, bisect

benjamin.jakimow@geo.hu-berlin.de
committed
from qgis.core import *
from qgis.core import QgsContrastEnhancement, QgsRasterShader, QgsColorRampShader, QgsProject, QgsCoordinateReferenceSystem, \
QgsRasterLayer, QgsVectorLayer, QgsMapLayer, QgsMapLayerProxyModel, QgsColorRamp, QgsSingleBandPseudoColorRenderer
from qgis.gui import *
from qgis.gui import QgsDockWidget, QgsMapCanvas, QgsMapTool, QgsCollapsibleGroupBox
from PyQt5.QtXml import *
from PyQt5.QtCore import *
from PyQt5.QtGui import *

benjamin.jakimow@geo.hu-berlin.de
committed
import numpy as np
from timeseriesviewer.utils import *
from timeseriesviewer.timeseries import SensorInstrument, TimeSeriesDatum, TimeSeries
from timeseriesviewer.utils import loadUI

benjamin.jakimow@geo.hu-berlin.de
committed
from timeseriesviewer.ui.mapviewscrollarea import MapViewScrollArea
from timeseriesviewer.mapcanvas import MapCanvas
#assert os.path.isfile(dummyPath)
#lyr = QgsRasterLayer(dummyPath)
#assert lyr.isValid()
DUMMY_RASTERINTERFACE = QgsSingleBandGrayRenderer(None, 0)
class MapViewUI(QFrame, loadUI('mapviewdefinition.ui')):
def __init__(self, parent=None):
super(MapViewUI, self).__init__(parent)
self.setupUi(self)
self.mSensors = collections.OrderedDict()
m = QMenu(self.btnToggleCrosshair)
m.addAction(self.actionSetCrosshairStyle)
#a = m.addAction('Set Crosshair Style')
self.btnToggleCrosshair.setMenu(m)
from timeseriesviewer.main import TimeSeriesViewer
tsv = TimeSeriesViewer.instance()
self.mVectorSourceModel = self.cbQgsVectorLayer.model().sourceModel()
QgsProject.instance().layersAdded.connect(self.mVectorSourceModel.addLayers)
QgsProject.instance().layersRemoved.connect(self.mVectorSourceModel.removeLayers)
self.mStore.layersAdded.connect(self.mVectorSourceModel.addLayers)
self.mStore.layersRemoved.connect(self.mVectorSourceModel.removeLayers)
#connect the QActions with the QgsCollapsibleGroupBoxes

Benjamin Jakimow
committed
self.gbVectorRendering.toggled.connect(self.actionToggleVectorVisibility.setChecked)
self.gbRasterRendering.toggled.connect(self.actionToggleRasterVisibility.setChecked)
self.btnToggleCrosshair.setDefaultAction(self.actionToggleCrosshairVisibility)
self.btnToggleMapViewVisibility.setDefaultAction(self.actionToggleMapViewHidden)

benjamin.jakimow@geo.hu-berlin.de
committed
self.btnSetVectorStyle.setDefaultAction(self.actionSetVectorStyle)

Benjamin Jakimow
committed
def addSensor(self, sensor:SensorInstrument):
"""
Registers a new SensorInstrument to this map view. Initializes widgets to handle sensor specific visualization properties.
:param sensor: SensorInstrument
"""

Benjamin Jakimow
committed
w = MapViewRenderSettings(sensor, parent=self)

Benjamin Jakimow
committed
w.collapsedStateChanged.connect(self.onSensorBoxCollapsed)
l = self.gbRasterRendering.layout()
i = l.count()-1
while i > 0 and not isinstance(l.itemAt(i), QWidget):
i -= 1
l.insertWidget(i, w, stretch=0, alignment=Qt.AlignTop)
self.mSensors[sensor] = w
return w
def removeSensor(self, sensor):
assert isinstance(sensor, SensorInstrument)
sensorSettings = self.mSensors.pop(sensor)
#l = self.renderSettingsLayout
l = self.gbRasterRendering.layout()
l.removeWidget(sensorSettings)
sensorSettings.close()

Benjamin Jakimow
committed
def onSensorBoxCollapsed(self, b:bool):
l = self.gbRasterRendering.layout()
for i in range(l.count()):
item = l.itemAt(i)

Benjamin Jakimow
committed
s = ""
class RendererWidgetModifications(object):
def __init__(self, *args):

Benjamin Jakimow
committed
self.initWidgetNames()

Benjamin Jakimow
committed
s = ""
gridLayoutOld = self.layout().children()[0]
self.gridLayout = QGridLayout()
while gridLayoutOld.count() > 0:
w = gridLayoutOld.takeAt(0)
w = w.widget()
gridLayoutOld.removeWidget(w)
w.setVisible(False)
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
l = self.layout()
l.removeItem(gridLayoutOld)
if isinstance(l, QBoxLayout):
l.insertItem(0, self.gridLayout)
self.layout().addStretch()
elif isinstance(l, QGridLayout):
l.addItem(self.gridLayout, 0, 0)
minMaxWidget = self.minMaxWidget()
if isinstance(minMaxWidget, QWidget):
minMaxWidget.layout().itemAt(0).widget().collapsedStateChanged.connect(self.onCollapsed)
def initWidgetNames(self, parent=None):
"""
Create a python variables to access QObjects which are child of parent
:param parent: QObject, self by default
"""
if parent is None:
parent = self
for c in parent.children():
setattr(parent, c.objectName(), c)
def onCollapsed(self, b):
hint = self.sizeHint()
self.parent().adjustSize()
# self.parent().setFixedSize(hint)
self.parent().parent().adjustSize()
def connectSliderWithBandComboBox(self, slider, combobox):
"""
Connects a band-selection slider with a band-selection combobox
:param widget: QgsRasterRendererWidget
:param slider: QSlider to show the band number
:param combobox: QComboBox to show the band name
:return:
"""
assert isinstance(self, QgsRasterRendererWidget)
assert isinstance(slider, QSlider)
assert isinstance(combobox, QComboBox)
# init the slider
lyr = self.rasterLayer()
if lyr.isValid():
nb = lyr.dataProvider().bandCount()
else:
ds = gdal.Open(lyr.source())
if isinstance(ds, gdal.Dataset):
nb = ds.RasterCount
else:
nb = 1
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
slider.setTickPosition(QSlider.TicksAbove)
slider.valueChanged.connect(combobox.setCurrentIndex)
slider.setMinimum(1)
slider.setMaximum(nb)
intervals = [1, 2, 5, 10, 25, 50]
for interval in intervals:
if nb / interval < 10:
break
slider.setTickInterval(interval)
slider.setPageStep(interval)
def onBandValueChanged(self, idx, slider):
assert isinstance(self, QgsRasterRendererWidget)
assert isinstance(idx, int)
assert isinstance(slider, QSlider)
# i = slider.value()
slider.blockSignals(True)
slider.setValue(idx)
slider.blockSignals(False)
# self.minMaxWidget().setBands(myBands)
# self.widgetChanged.emit()
if self.comboBoxWithNotSetItem(combobox):
combobox.currentIndexChanged[int].connect(lambda idx: onBandValueChanged(self, idx, slider))
else:
combobox.currentIndexChanged[int].connect(lambda idx: onBandValueChanged(self, idx + 1, slider))
s = ""
def comboBoxWithNotSetItem(self, cb)->bool:
assert isinstance(cb, QComboBox)
data = cb.itemData(0, role=Qt.DisplayRole)
return re.search(r'^(not set|none|nonetype)$', str(data).strip(), re.I) is not None
def setLayoutItemVisibility(self, grid, isVisible):
assert isinstance(self, QgsRasterRendererWidget)
for i in range(grid.count()):
item = grid.itemAt(i)
if isinstance(item, QLayout):
s = ""
elif isinstance(item, QWidgetItem):
item.widget().setVisible(isVisible)
item.widget().setParent(self)
else:
s = ""
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
def setBandSelection(self, key):
key = key.upper()
if key == 'DEFAULT':
bandIndices = defaultBands(self.rasterLayer())
else:
colors = re.split('[ ,;:]', key)
bandIndices = [bandClosestToWavelength(self.rasterLayer(), c) for c in colors]
n = min(len(bandIndices), len(self.mBandComboBoxes))
for i in range(n):
cb = self.mBandComboBoxes[i]
bandIndex = bandIndices[i]
if self.comboBoxWithNotSetItem(cb):
cb.setCurrentIndex(bandIndex+1)
else:
cb.setCurrentIndex(bandIndex)
def fixBandNames(self, comboBox):
"""
Changes the QGIS default bandnames ("Band 001") to more meaningfull information including gdal.Dataset.Descriptions.
:param widget:
:param comboBox:
"""
nb = self.rasterLayer().bandCount()
assert isinstance(self, QgsRasterRendererWidget)
assert isinstance(comboBox, QComboBox)
#comboBox.clear()
m = comboBox.model()
assert isinstance(m, QStandardItemModel)
bandNames = displayBandNames(self.rasterLayer())
b = 1 if nb < comboBox.count() else 0
for i in range(nb):
item = m.item(i+b,0)
assert isinstance(item, QStandardItem)
item.setData(bandNames[i], Qt.DisplayRole)
item.setData('Band {} "{}"'.format(i+1, bandNames[i]), Qt.ToolTipRole)
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
def displayBandNames(provider_or_dataset, bands=None):
results = None
if isinstance(provider_or_dataset, QgsRasterLayer):
return displayBandNames(provider_or_dataset.dataProvider())
elif isinstance(provider_or_dataset, QgsRasterDataProvider):
if provider_or_dataset.name() == 'gdal':
ds = gdal.Open(provider_or_dataset.dataSourceUri())
results = displayBandNames(ds, bands=bands)
else:
# same as in QgsRasterRendererWidget::displayBandName
results = []
if bands is None:
bands = range(1, provider_or_dataset.bandCount() + 1)
for band in bands:
result = provider_or_dataset.generateBandName(band)
colorInterp ='{}'.format(provider_or_dataset.colorInterpretationName(band))
if colorInterp != 'Undefined':
result += '({})'.format(colorInterp)
results.append(result)
elif isinstance(provider_or_dataset, gdal.Dataset):
results = []
if bands is None:
bands = range(1, provider_or_dataset.RasterCount+1)
for band in bands:
b = provider_or_dataset.GetRasterBand(band)
descr = b.GetDescription()
if len(descr) == 0:
descr = 'Band {}'.format(band)
results.append(descr)
return results
class SingleBandPseudoColorRendererWidget(QgsSingleBandPseudoColorRendererWidget, RendererWidgetModifications):
@staticmethod
def create(layer, extent):
return SingleBandPseudoColorRendererWidget(layer, extent)
def __init__(self, layer, extent):
super(SingleBandPseudoColorRendererWidget, self).__init__(layer, extent)

Benjamin Jakimow
committed
self.gridLayout = self.layout().children()[0]
assert isinstance(self.gridLayout, QGridLayout)
for i in range(self.gridLayout.count()):
w = self.gridLayout.itemAt(i)
w = w.widget()
if isinstance(w, QWidget):
setattr(self, w.objectName(), w)

Benjamin Jakimow
committed
toReplace = [self.mBandComboBox, self.mMinLabel, self.mMaxLabel, self.mMinLineEdit, self.mMaxLineEdit]
for w in toReplace:
self.gridLayout.removeWidget(w)
w.setVisible(False)

Benjamin Jakimow
committed
self.mBandSlider = QSlider(Qt.Horizontal, self)
self.mBandComboBoxes.append(self.mBandComboBox)
self.fixBandNames(self.mBandComboBox)
self.connectSliderWithBandComboBox(self.mBandSlider, self.mBandComboBox)

Benjamin Jakimow
committed
self.mBtnBar = QFrame(self)

Benjamin Jakimow
committed
grid = QGridLayout(self)
grid.addWidget(self.mBtnBar, 0, 0, 1, 4, Qt.AlignLeft)
grid.addWidget(self.mBandSlider, 1, 0, 1, 2)
grid.addWidget(self.mBandComboBox, 1, 2, 1, 2)
grid.addWidget(self.mMinLabel, 2, 0)
grid.addWidget(self.mMinLineEdit, 2, 1)
grid.addWidget(self.mMaxLabel, 2, 2)
grid.addWidget(self.mMaxLineEdit, 2, 3)

Benjamin Jakimow
committed
grid.setColumnStretch(0, 0)
grid.setColumnStretch(1, 2)
grid.setColumnStretch(2, 0)
grid.setColumnStretch(3, 2)
grid.setSpacing(2)

Benjamin Jakimow
committed
self.gridLayout.addItem(grid, 0, 1, 2, 4)
self.gridLayout.setSpacing(2)
self.setLayoutItemVisibility(grid, True)
def initActionButtons(self):
wl, wlu = parseWavelength(self.rasterLayer())
self.wavelengths = wl
self.wavelengthUnit = wlu

Benjamin Jakimow
committed
self.mBtnBar.setLayout(QHBoxLayout(self.mBtnBar))
self.mBtnBar.layout().addStretch()
self.mBtnBar.layout().setContentsMargins(0, 0, 0, 0)
self.mBtnBar.layout().setSpacing(2)
self.actionSetDefault = QAction('Default', self)
self.actionSetRed = QAction('R', self)
self.actionSetGreen = QAction('G', self)
self.actionSetBlue = QAction('B', self)
self.actionSetNIR = QAction('nIR', self)
self.actionSetSWIR = QAction('swIR', self)
self.actionSetDefault.triggered.connect(lambda: self.setBandSelection('default'))
self.actionSetRed.triggered.connect(lambda: self.setBandSelection('R'))
self.actionSetGreen.triggered.connect(lambda: self.setBandSelection('G'))
self.actionSetBlue.triggered.connect(lambda: self.setBandSelection('B'))
self.actionSetNIR.triggered.connect(lambda: self.setBandSelection('nIR'))
self.actionSetSWIR.triggered.connect(lambda: self.setBandSelection('swIR'))
def addBtnAction(action):

Benjamin Jakimow
committed
btn = QToolButton(self.mBtnBar)
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
btn.setDefaultAction(action)
self.mBtnBar.layout().addWidget(btn)
self.insertAction(None, action)
return btn
self.btnDefault = addBtnAction(self.actionSetDefault)
self.btnRed = addBtnAction(self.actionSetRed)
self.btnGreen = addBtnAction(self.actionSetGreen)
self.btnBlue = addBtnAction(self.actionSetRed)
self.btnNIR = addBtnAction(self.actionSetNIR)
self.btnSWIR = addBtnAction(self.actionSetSWIR)
b = self.wavelengths is not None
for a in [self.actionSetRed, self.actionSetGreen, self.actionSetBlue, self.actionSetNIR, self.actionSetSWIR]:
a.setEnabled(b)
def displayBandNames(provider_or_dataset, bands=None):
results = None
if isinstance(provider_or_dataset, QgsRasterLayer):
return displayBandNames(provider_or_dataset.dataProvider())
elif isinstance(provider_or_dataset, QgsRasterDataProvider):
if provider_or_dataset.name() == 'gdal':
ds = gdal.Open(provider_or_dataset.dataSourceUri())
results = displayBandNames(ds, bands=bands)
else:
# same as in QgsRasterRendererWidget::displayBandName
results = []
if bands is None:
bands = range(1, provider_or_dataset.bandCount() + 1)
for band in bands:
result = provider_or_dataset.generateBandName(band)
colorInterp ='{}'.format(provider_or_dataset.colorInterpretationName(band))
if colorInterp != 'Undefined':
result += '({})'.format(colorInterp)
results.append(result)
elif isinstance(provider_or_dataset, gdal.Dataset):
results = []
if bands is None:
bands = range(1, provider_or_dataset.RasterCount+1)
for band in bands:
b = provider_or_dataset.GetRasterBand(band)
descr = b.GetDescription()
if len(descr) == 0:
descr = 'Band {}'.format(band)
results.append(descr)
return results

Benjamin Jakimow
committed
class SingleBandGrayRendererWidget(QgsSingleBandGrayRendererWidget, RendererWidgetModifications):
@staticmethod
def create(layer, extent):
return SingleBandGrayRendererWidget(layer, extent)
def __init__(self, layer, extent):
super(SingleBandGrayRendererWidget, self).__init__(layer, extent)
self.modifyGridLayout()

Benjamin Jakimow
committed
self.mGrayBandSlider = QSlider(Qt.Horizontal, self)
self.mBandComboBoxes.append(self.mGrayBandComboBox)
self.fixBandNames(self.mGrayBandComboBox)
self.connectSliderWithBandComboBox(self.mGrayBandSlider, self.mGrayBandComboBox)

Benjamin Jakimow
committed
self.mBtnBar = QFrame(self)
self.initActionButtons()
self.gridLayout.addWidget(self.mGrayBandLabel, 0, 0)
self.gridLayout.addWidget(self.mBtnBar, 0, 1, 1, 4, Qt.AlignLeft)
self.gridLayout.addWidget(self.mGrayBandSlider, 1, 1, 1, 2)

Benjamin Jakimow
committed
self.gridLayout.addWidget(self.mGrayBandComboBox, 1, 3, 1, 2)
self.gridLayout.addWidget(self.label, 2, 0)
self.gridLayout.addWidget(self.mGradientComboBox, 2, 1, 1, 4)
self.gridLayout.addWidget(self.mMinLabel, 3, 1)
self.gridLayout.addWidget(self.mMinLineEdit, 3, 2)
self.gridLayout.addWidget(self.mMaxLabel, 3, 3)
self.gridLayout.addWidget(self.mMaxLineEdit, 3, 4)
self.gridLayout.addWidget(self.mContrastEnhancementLabel, 4, 0)

Benjamin Jakimow
committed
self.gridLayout.addWidget(self.mContrastEnhancementComboBox, 4, 1, 1, 4)
self.gridLayout.setSpacing(2)
self.setLayoutItemVisibility(self.gridLayout, True)
self.mDefaultRenderer = layer.renderer()
self.setFromRenderer(self.mDefaultRenderer)
def initActionButtons(self):
wl, wlu = parseWavelength(self.rasterLayer())
self.wavelengths = wl
self.wavelengthUnit = wlu

Benjamin Jakimow
committed
self.mBtnBar.setLayout(QHBoxLayout(self))
self.mBtnBar.layout().addStretch()
self.mBtnBar.layout().setContentsMargins(0, 0, 0, 0)
self.mBtnBar.layout().setSpacing(2)
self.actionSetDefault = QAction('Default', self)
self.actionSetRed = QAction('R', self)
self.actionSetGreen = QAction('G', self)
self.actionSetBlue = QAction('B', self)
self.actionSetNIR = QAction('nIR', self)
self.actionSetSWIR = QAction('swIR', self)
self.actionSetDefault.triggered.connect(lambda: self.setBandSelection('default'))
self.actionSetRed.triggered.connect(lambda: self.setBandSelection('R'))
self.actionSetGreen.triggered.connect(lambda: self.setBandSelection('G'))
self.actionSetBlue.triggered.connect(lambda: self.setBandSelection('B'))
self.actionSetNIR.triggered.connect(lambda: self.setBandSelection('nIR'))
self.actionSetSWIR.triggered.connect(lambda: self.setBandSelection('swIR'))
def addBtnAction(action):

Benjamin Jakimow
committed
btn = QToolButton(self.mBtnBar)
btn.setDefaultAction(action)
self.mBtnBar.layout().addWidget(btn)
self.insertAction(None, action)
return btn
self.btnDefault = addBtnAction(self.actionSetDefault)
self.btnBlue = addBtnAction(self.actionSetBlue)
self.btnGreen = addBtnAction(self.actionSetGreen)
self.btnRed = addBtnAction(self.actionSetRed)
self.btnNIR = addBtnAction(self.actionSetNIR)
self.btnSWIR = addBtnAction(self.actionSetSWIR)
b = self.wavelengths is not None
for a in [self.actionSetRed, self.actionSetGreen, self.actionSetBlue, self.actionSetNIR, self.actionSetSWIR]:
a.setEnabled(b)
class PalettedRendererWidget(QgsPalettedRendererWidget, RendererWidgetModifications):
@staticmethod
def create(layer, extent):
return PalettedRendererWidget(layer, extent)
def __init__(self, layer, extent):
super(PalettedRendererWidget, self).__init__(layer, extent)
#self.modifyGridLayout()
self.fixBandNames(self.mBandComboBox)
self.mTreeView.setMinimumSize(QSize(10,10))
s = ""
class MultiBandColorRendererWidget(QgsMultiBandColorRendererWidget, RendererWidgetModifications):
@staticmethod
def create(layer, extent):
return MultiBandColorRendererWidget(layer, extent)
def __init__(self, layer, extent):
super(MultiBandColorRendererWidget, self).__init__(layer, extent)
self.modifyGridLayout()

Benjamin Jakimow
committed
self.mRedBandSlider = QSlider(Qt.Horizontal, self)
self.mGreenBandSlider = QSlider(Qt.Horizontal, self)
self.mBlueBandSlider = QSlider(Qt.Horizontal, self)
self.mBandComboBoxes.extend([self.mRedBandComboBox, self.mGreenBandComboBox, self.mBlueBandComboBox])
self.mSliders = [self.mRedBandSlider, self.mGreenBandSlider, self.mBlueBandSlider]
for cbox, slider in zip(self.mBandComboBoxes, self.mSliders):
self.connectSliderWithBandComboBox(slider, cbox)
self.fixBandNames(self.mRedBandComboBox)
self.fixBandNames(self.mGreenBandComboBox)
self.fixBandNames(self.mBlueBandComboBox)

Benjamin Jakimow
committed
self.mBtnBar = QFrame(self)
self.mBtnBar.setLayout(QHBoxLayout(self.mBtnBar))
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
self.initActionButtons()
self.mBtnBar.layout().addStretch()
self.mBtnBar.layout().setContentsMargins(0, 0, 0, 0)
self.mBtnBar.layout().setSpacing(2)
#self.gridLayout.deleteLater()
# self.gridLayout = newGrid
self.gridLayout.addWidget(self.mBtnBar, 0, 1, 1, 3)
self.gridLayout.addWidget(self.mRedBandLabel, 1, 0)
self.gridLayout.addWidget(self.mRedBandSlider, 1, 1)
self.gridLayout.addWidget(self.mRedBandComboBox, 1, 2)
self.gridLayout.addWidget(self.mRedMinLineEdit, 1, 3)
self.gridLayout.addWidget(self.mRedMaxLineEdit, 1, 4)
self.gridLayout.addWidget(self.mGreenBandLabel, 2, 0)
self.gridLayout.addWidget(self.mGreenBandSlider, 2, 1)
self.gridLayout.addWidget(self.mGreenBandComboBox, 2, 2)
self.gridLayout.addWidget(self.mGreenMinLineEdit, 2, 3)
self.gridLayout.addWidget(self.mGreenMaxLineEdit, 2, 4)
self.gridLayout.addWidget(self.mBlueBandLabel, 3, 0)
self.gridLayout.addWidget(self.mBlueBandSlider, 3, 1)
self.gridLayout.addWidget(self.mBlueBandComboBox, 3, 2)
self.gridLayout.addWidget(self.mBlueMinLineEdit, 3, 3)
self.gridLayout.addWidget(self.mBlueMaxLineEdit, 3, 4)
self.gridLayout.addWidget(self.mContrastEnhancementAlgorithmLabel, 4, 0, 1, 2)
self.gridLayout.addWidget(self.mContrastEnhancementAlgorithmComboBox, 4, 2, 1, 3)
self.setLayoutItemVisibility(self.gridLayout, True)
self.mRedBandLabel.setText('R')
self.mGreenBandLabel.setText('G')
self.mBlueBandLabel.setText('B')
self.mDefaultRenderer = layer.renderer()
def initActionButtons(self):
wl, wlu = parseWavelength(self.rasterLayer())
self.wavelengths = wl
self.wavelengthUnit = wlu
self.actionSetDefault = QAction('Default', self)
self.actionSetTrueColor = QAction('RGB', self)
self.actionSetCIR = QAction('nIR', self)
self.actionSet453 = QAction('swIR', self)
self.actionSetDefault.triggered.connect(lambda: self.setBandSelection('default'))
self.actionSetTrueColor.triggered.connect(lambda: self.setBandSelection('R,G,B'))
self.actionSetCIR.triggered.connect(lambda: self.setBandSelection('nIR,R,G'))
self.actionSet453.triggered.connect(lambda: self.setBandSelection('nIR,swIR,R'))
def addBtnAction(action):

Benjamin Jakimow
committed
btn = QToolButton(self.mBtnBar)
btn.setDefaultAction(action)
self.mBtnBar.layout().addWidget(btn)
self.insertAction(None, action)
return btn
self.btnDefault = addBtnAction(self.actionSetDefault)
self.btnTrueColor = addBtnAction(self.actionSetTrueColor)
self.btnCIR = addBtnAction(self.actionSetCIR)
self.btn453 = addBtnAction(self.actionSet453)
b = self.wavelengths is not None
for a in [self.actionSetCIR, self.actionSet453, self.actionSetTrueColor]:
a.setEnabled(b)

benjamin.jakimow@geo.hu-berlin.de
committed
class MapView(QObject):
sigRemoveMapView = pyqtSignal(object)
sigMapViewVisibility = pyqtSignal(bool)
#sigVectorVisibility = pyqtSignal(bool)
#sigRasterVisibility = pyqtSignal(bool)

benjamin.jakimow@geo.hu-berlin.de
committed
sigTitleChanged = pyqtSignal(str)

benjamin.jakimow@geo.hu-berlin.de
committed
sigSensorRendererChanged = pyqtSignal(SensorInstrument, QgsRasterRenderer)

Benjamin Jakimow
committed

benjamin.jakimow@geo.hu-berlin.de
committed
sigVectorLayerChanged = pyqtSignal()
sigShowProfiles = pyqtSignal(SpatialPoint, MapCanvas, str)

benjamin.jakimow@geo.hu-berlin.de
committed

benjamin.jakimow@geo.hu-berlin.de
committed
def __init__(self, mapViewCollectionDock, name='Map View', recommended_bands=None, parent=None):

benjamin.jakimow@geo.hu-berlin.de
committed
super(MapView, self).__init__()

benjamin.jakimow@geo.hu-berlin.de
committed
assert isinstance(mapViewCollectionDock, MapViewCollectionDock)

benjamin.jakimow@geo.hu-berlin.de
committed
self.ui = MapViewUI(mapViewCollectionDock.stackedWidget)

Benjamin Jakimow
committed
self.ui.cbQgsVectorLayer.setFilters(QgsMapLayerProxyModel.VectorLayer)
self.ui.cbQgsVectorLayer.layerChanged.connect(self.setVectorLayer)
self.ui.tbName.textChanged.connect(self.sigTitleChanged.emit)
from timeseriesviewer.crosshair import getCrosshairStyle
self.ui.actionSetCrosshairStyle.triggered.connect(

Benjamin Jakimow
committed
lambda : self.onCrosshairChanged(getCrosshairStyle(

Benjamin Jakimow
committed
crosshairStyle=self.crosshairStyle()))

benjamin.jakimow@geo.hu-berlin.de
committed
self.mapViewCollection = mapViewCollectionDock
self.mSensorViews = collections.OrderedDict()

benjamin.jakimow@geo.hu-berlin.de
committed
self.mVectorLayer = None

benjamin.jakimow@geo.hu-berlin.de
committed
self.setVectorLayer(None)
self.ui.actionToggleVectorVisibility.toggled.connect(self.setVectorVisibility)
self.ui.actionToggleRasterVisibility.toggled.connect(self.setRasterVisibility)

Benjamin Jakimow
committed
self.ui.actionToggleCrosshairVisibility.toggled.connect(self.onCrosshairChanged)
self.ui.actionToggleMapViewHidden.toggled.connect(lambda b: self.setIsVisible(not b))
self.ui.actionToggleVectorVisibility.setChecked(False)
self.ui.actionToggleRasterVisibility.setChecked(True)

benjamin.jakimow@geo.hu-berlin.de
committed
self.ui.actionSetVectorStyle.triggered.connect(self.setVectorLayerStyle)
self.addSensor(sensor)

Benjamin Jakimow
committed
def setIsVisible(self, b: bool):
"""
Sets the map view visibility
:param b: bool
"""
for mapCanvas in self.mapCanvases():
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.setVisible(b)
if self.ui.actionToggleMapViewHidden.isChecked() == b:
self.ui.actionToggleMapViewHidden.setChecked(not b)
if changed:
self.sigMapViewVisibility.emit(b)

Benjamin Jakimow
committed
def isVisible(self)->bool:
"""
Returns the map view visibility
:return: bool
"""
return not self.ui.actionToggleMapViewHidden.isChecked()

Benjamin Jakimow
committed
def mapCanvases(self)->list:
"""
Returns the MapCanvases related to this map view
:return: [list-of-MapCanvases]
"""
m = []
for sensor, sensorView in self.mSensorViews.items():
m.extend(sensorView.mapCanvases())
return m

benjamin.jakimow@geo.hu-berlin.de
committed
def setVectorLayerStyle(self, *args):
if isinstance(self.mVectorLayer, QgsVectorLayer):
d = QgsRendererPropertiesDialog(self.mVectorLayer, QgsStyle.defaultStyle())

Benjamin Jakimow
committed
mc = self.mapCanvases()
if len(mc) > 0:
d.setMapCanvas(mc[0])

benjamin.jakimow@geo.hu-berlin.de
committed
d.exec_()

benjamin.jakimow@geo.hu-berlin.de
committed

Benjamin Jakimow
committed
def vectorLayerRenderer(self)->QgsFeatureRenderer:
if isinstance(self.mVectorLayer, QgsVectorLayer):

benjamin.jakimow@geo.hu-berlin.de
committed
return self.mVectorLayer.renderer()
return None

benjamin.jakimow@geo.hu-berlin.de
committed
def setVectorLayerRenderer(self, renderer):

benjamin.jakimow@geo.hu-berlin.de
committed
if isinstance(renderer, QgsFeatureRenderer) and \
isinstance(self.mVectorLayer, QgsVectorLayer):
self.mVectorLayer.setRendererV2(renderer)
def setVectorLayer(self, lyr:QgsVectorLayer):
"""
Sets a QgsVectorLayer that is shown on top of all other QgsRasterLayers
:param lyr:
:return:
"""
b = False
#remove last layer
if isinstance(self.mVectorLayer, QgsVectorLayer):
for mapCanvas in self.mapCanvases():
if self.mVectorLayer in mapCanvas.mLayerSources.remove(self.mVectorLayer):
mapCanvas.mLayerSources.remove(self.mVectorLayer)
b = True
# add new layer

Benjamin Jakimow
committed
if isinstance(lyr, QgsVectorLayer) and self.ui.gbVectorRendering.isChecked():
b = True
#add vector layer
self.mVectorLayer = lyr
self.mVectorLayer.rendererChanged.connect(self.sigVectorLayerChanged)

benjamin.jakimow@geo.hu-berlin.de
committed
for mapCanvas in self.mapCanvases():
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.mapLayerModel().addMapLayerSources([self.mVectorLayer])
if b:
self.sigVectorLayerChanged.emit()

benjamin.jakimow@geo.hu-berlin.de
committed
def applyStyles(self):
"""Applies all style changes to all sensor views."""
for sensorView in self.mSensorViews.values():

benjamin.jakimow@geo.hu-berlin.de
committed
sensorView.applyStyle()
def setTitle(self, title:str):
"""
Sets the widget title
:param title: str
"""
old = self.title()
if old != title:
self.ui.tbName.setText(title)

benjamin.jakimow@geo.hu-berlin.de
committed
def title(self)->str:
"""Returns the title."""

benjamin.jakimow@geo.hu-berlin.de
committed

Benjamin Jakimow
committed
def refreshMapView(self, sensor=None):
if isinstance(sensor, SensorInstrument):
sensorSettings = [self.mSensorViews[sensor]]
else:
#update all sensors
sensorSettings = self.mSensorViews.values()

Benjamin Jakimow
committed
for renderSetting in sensorSettings:
assert isinstance(renderSetting, MapViewRenderSettings)
renderSetting.applyStyle()
for mapCanvas in self.mapCanvases():
if isinstance(mapCanvas, MapCanvas):
mapCanvas.refresh()
def setCrosshairStyle(self, crosshairStyle:CrosshairStyle):
"""

Benjamin Jakimow
committed
Seths the CrosshairStyle of this MapView
:param crosshairStyle: CrosshairStyle
"""

Benjamin Jakimow
committed
self.onCrosshairChanged(crosshairStyle)

benjamin.jakimow@geo.hu-berlin.de
committed
def setHighlighted(self, b=True, timeout=1000):
"""
Activates or deactivates a red-line border of the MapCanvases
:param b: True | False to activate / deactivate the highlighted lines-
:param timeout: int, milliseconds how long the highlighted frame should appear
"""

benjamin.jakimow@geo.hu-berlin.de
committed
styleOn = """.MapCanvas {
border: 4px solid red;
border-radius: 4px;
}"""
styleOff = """"""
if b is True:
for mapCanvas in self.mapCanvases():
mapCanvas.setStyleSheet(styleOn)
if timeout > 0:
QTimer.singleShot(timeout, lambda : self.setHighlighted(False))
else:
for mapCanvas in self.mapCanvases():
mapCanvas.setStyleSheet(styleOff)
def rasterVisibility(self)->bool:
"""
Returns whether raster images should be visible.
:return: bool
"""
return self.ui.actionToggleRasterVisibility.isChecked()
def vectorVisibility(self)->bool:
"""
Returns whether vector images should be visible.
:return: bool
"""
return self.ui.actionToggleVectorVisibility.isChecked()
def setRasterVisibility(self, b:bool):
"""
Sets visibility of rasters.
:param b: bool
"""

Benjamin Jakimow
committed
self.ui.actionToggleRasterVisibility.setChecked(b)
for mapCanvas in self.mapCanvases():
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.setLayerVisibility(QgsRasterLayer, b)
def setVectorVisibility(self, b:bool):
"""
Sets the visibility of vector layers.
:param b:
:return:
"""
assert isinstance(b, bool)

Benjamin Jakimow
committed
self.mVectorsVisible = b
self.ui.actionToggleVectorVisibility.setChecked(b)

Benjamin Jakimow
committed
for mapCanvas in self.mapCanvases():
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.setLayerVisibility(QgsVectorLayer, b)
def removeSensor(self, sensor:SensorInstrument):
assert sensor in self.mSensorViews.keys()
self.mSensorViews.pop(sensor)

benjamin.jakimow@geo.hu-berlin.de
committed
def hasSensor(self, sensor):
assert type(sensor) is SensorInstrument
return sensor in self.mSensorViews.keys()

benjamin.jakimow@geo.hu-berlin.de
committed
def registerMapCanvas(self, sensor:SensorInstrument, mapCanvas:MapCanvas):
"""
Registers a new MapCanvas to this MapView
:param sensor:
:param mapCanvas:
:return:
"""
from timeseriesviewer.mapcanvas import MapCanvas
assert isinstance(mapCanvas, MapCanvas)
assert isinstance(sensor, SensorInstrument)
mapViewRenderSettings = self.mSensorViews[sensor]
assert isinstance(mapViewRenderSettings, MapViewRenderSettings)
mapViewRenderSettings.registerMapCanvas(mapCanvas)
#register signals sensor specific signals

Benjamin Jakimow
committed
mapCanvas.sigCrosshairVisibilityChanged.connect(self.onCrosshairChanged)
mapCanvas.sigCrosshairStyleChanged.connect(self.onCrosshairChanged)
#register non-sensor specific signals for this mpa view
self.sigMapViewVisibility.connect(mapCanvas.refresh)
self.sigVectorLayerChanged.connect(mapCanvas.refresh)
# self.sigVectorVisibility.connect(mapCanvas.refresh)

Benjamin Jakimow
committed
def crosshairStyle(self)->CrosshairStyle:
"""
Returns the CrosshairStyle
:return:
"""
for c in self.mapCanvases():
assert isinstance(c, MapCanvas)
style = c.crosshairStyle()
if isinstance(style, CrosshairStyle):
return style
return None

Benjamin Jakimow
committed
def onCrosshairChanged(self, obj):
"""
Synchronizes all crosshair positions. Takes care of CRS differences.
:param spatialPoint: SpatialPoint of the new Crosshair position
"""
from timeseriesviewer.crosshair import CrosshairStyle
srcCanvas = self.sender()
if isinstance(srcCanvas, MapCanvas):
dstCanvases = [c for c in self.mapCanvases() if c != srcCanvas]
else:
dstCanvases = [c for c in self.mapCanvases()]
if isinstance(obj, bool):
for mapCanvas in dstCanvases:
mapCanvas.setCrosshairVisibility(obj, emitSignal=False)

Benjamin Jakimow
committed
if isinstance(obj, CrosshairStyle):
for mapCanvas in dstCanvases:
mapCanvas.setCrosshairStyle(obj, emitSignal=False)

benjamin.jakimow@geo.hu-berlin.de
committed
def addSensor(self, sensor):
"""
:param sensor:
:return:
"""
if isinstance(sensor, SensorInstrument) and sensor not in self.mSensorViews.keys():

benjamin.jakimow@geo.hu-berlin.de
committed
#w.showSensorName(False)
w = self.ui.addSensor(sensor)

Benjamin Jakimow
committed
w.sigRendererChanged.connect(lambda s=sensor : self.refreshMapView(sensor=s))
#w.sigSensorRendererChanged.connect(self.onSensorRenderingChanged)
self.mSensorViews[sensor] = w
s =""
def onSensorRenderingChanged(self, renderer):
sensorSettings = self.sender()
assert isinstance(sensorSettings, MapViewSensorSettings)
for mapCanvas in sensorSettings.mapCanvases():
mapCanvas.setRenderer(renderer)
#mapCanvas.refresh()

benjamin.jakimow@geo.hu-berlin.de
committed
assert type(sensor) is SensorInstrument
return self.mSensorViews[sensor]

benjamin.jakimow@geo.hu-berlin.de
committed
1006
1007
1008
1009
1010
1011
1012
1013
1014
1015
1016
1017
1018
1019
1020
1021
1022
1023
1024
1025
1026
1027
1028
1029
1030
1031
1032
1033
1034
1035
1036
1037
1038
1039
1040
1041
1042
1043
1044
1045
1046
1047
def displayBandNames(provider_or_dataset, bands=None):
results = None
if isinstance(provider_or_dataset, QgsRasterLayer):
return displayBandNames(provider_or_dataset.dataProvider())
elif isinstance(provider_or_dataset, QgsRasterDataProvider):
if provider_or_dataset.name() == 'gdal':
ds = gdal.Open(provider_or_dataset.dataSourceUri())
results = displayBandNames(ds, bands=bands)
else:
# same as in QgsRasterRendererWidget::displayBandName
results = []
if bands is None:
bands = range(1, provider_or_dataset.bandCount() + 1)
for band in bands:
result = provider_or_dataset.generateBandName(band)
colorInterp ='{}'.format(provider_or_dataset.colorInterpretationName(band))
if colorInterp != 'Undefined':
result += '({})'.format(colorInterp)
results.append(result)
elif isinstance(provider_or_dataset, gdal.Dataset):
results = []
if bands is None:
bands = range(1, provider_or_dataset.RasterCount+1)
for band in bands:
b = provider_or_dataset.GetRasterBand(band)
descr = b.GetDescription()
if len(descr) == 0:
descr = 'Band {}'.format(band)
results.append(descr)
return results
class RasterDataProviderMockup(QgsRasterDataProvider):
def __init__(self):
super(RasterDataProviderMockup, self).__init__('')
class MapViewRenderSettings(QgsCollapsibleGroupBox, loadUI('mapviewrendersettings.ui')):
LUT_RENDERER = {QgsMultiBandColorRenderer:QgsMultiBandColorRendererWidget,
QgsSingleBandGrayRenderer:QgsSingleBandGrayRendererWidget,
QgsSingleBandPseudoColorRenderer:QgsSingleBandPseudoColorRendererWidget,
QgsPalettedRasterRenderer:QgsPalettedRendererWidget}
LUT_RENDERER[QgsMultiBandColorRenderer]=MultiBandColorRendererWidget
LUT_RENDERER[QgsSingleBandPseudoColorRenderer]=SingleBandPseudoColorRendererWidget
LUT_RENDERER[QgsSingleBandGrayRenderer]=SingleBandGrayRendererWidget
LUT_RENDERER[QgsPalettedRasterRenderer] = PalettedRendererWidget

Benjamin Jakimow
committed
sigRendererChanged = pyqtSignal()
def __init__(self, sensor, parent=None):
"""Constructor."""
1064
1065
1066
1067
1068
1069
1070
1071
1072
1073
1074
1075
1076
1077
1078
1079
1080
1081
1082
1083
1084
1085
# Set up the user interface from Designer.
# After setupUI you can access any designer object by doing
# self.<objectname>, and you can use autoconnect slots - see
# http://qt-project.org/doc/qt-4.8/designer-using-a-ui-file.html
# #widgets-and-dialogs-with-auto-connect
self.setupUi(self)
QApplication.clipboard().dataChanged.connect(self.onClipboardChanged)
assert isinstance(sensor, SensorInstrument)
self.mSensor = sensor
self.mSensor.sigNameChanged.connect(self.onSensorNameChanged)
self.setTitle(self.mSensor.name())
from timeseriesviewer.models import OptionListModel, Option
rasterRendererModel = OptionListModel()
#rasterRendererModel.addOption(Option(QgsMultiBandColorRendererWidget, name='multibandcolor (QGIS)', mRenderType = QgsMultiBandColorRenderer))
#rasterRendererModel.addOption(Option(QgsPalettedRendererWidget, name='paletted (QGIS)', mRenderType = QgsPalettedRasterRenderer))
#rasterRendererModel.addOption(Option(QgsSingleBandGrayRendererWidget, name='singlegray (QGIS)', mRenderType = QgsSingleBandGrayRenderer))
#rasterRendererModel.addOption(Option(QgsSingleBandPseudoColorRendererWidget, name='singlebandpseudocolor (QGIS)', mRenderType = QgsSingleBandPseudoColorRenderer))
rasterRendererModel.addOption(Option(MultiBandColorRendererWidget, name='Multibandcolor', mRenderType=QgsMultiBandColorRenderer))
rasterRendererModel.addOption(Option(SingleBandGrayRendererWidget, name='Singlegray', mRenderType=QgsSingleBandGrayRenderer))
rasterRendererModel.addOption(Option(SingleBandPseudoColorRendererWidget, name='Singleband Pseudocolor', mRenderType=QgsSingleBandPseudoColorRenderer))
rasterRendererModel.addOption(Option(PalettedRendererWidget, name='Paletted', mRenderType=QgsPalettedRasterRenderer))
self.mRasterRendererModel = rasterRendererModel
self.cbRenderType.setModel(self.mRasterRendererModel)
assert isinstance(self.stackedWidget, QStackedWidget)

Benjamin Jakimow
committed
self.mMockupCanvas = QgsMapCanvas(parent)
self.mMockupCanvas.setVisible(False)
self.mMockupRasterLayer = self.mSensor.mockupLayer()
self.mMockupCanvas.setLayers([self.mMockupRasterLayer])
for func in rasterRendererModel.optionValues():
#extent = self.canvas.extent()
#w = func.create(self.mMockupRasterLayer, self.mMockupRasterLayer.extent())
w = func(self.mMockupRasterLayer, self.mMockupRasterLayer.extent())

Benjamin Jakimow
committed
assert isinstance(w, QgsRasterRendererWidget)
w.setParent(self.stackedWidget)
w.setSizePolicy(QSizePolicy(QSizePolicy.Preferred, QSizePolicy.Preferred))
w.setMapCanvas(self.mMockupCanvas)
self.stackedWidget.addWidget(w)

Benjamin Jakimow
committed
s = ""
self.initActions()
def initActions(self):
self.btnPasteStyle.setDefaultAction(self.actionPasteStyle)
self.btnCopyStyle.setDefaultAction(self.actionCopyStyle)
self.btnApplyStyle.setDefaultAction(self.actionApplyStyle)
clipboardRenderer = rendererFromXml(QApplication.clipboard().mimeData())
self.actionPasteStyle.setEnabled(isinstance(clipboardRenderer, QgsRasterRenderer))
self.actionPasteStyle.triggered.connect(self.pasteStyleFromClipboard)
self.actionCopyStyle.triggered.connect(self.pasteStyleToClipboard)
self.actionApplyStyle.triggered.connect(self.applyStyle)
def mapCanvases(self):
return self.mMapCanvases[:]
def registerMapCanvas(self, mapCanvas):
assert isinstance(mapCanvas, MapCanvas)
self.mMapCanvases.append(mapCanvas)
mapCanvas.sigChangeSVRequest.connect(self.onMapCanvasRendererChangeRequest)
def onMapCanvasRendererChangeRequest(self, mapCanvas, renderer):
self.setRasterRenderer(renderer)
self.applyStyle()
s = ""
def onSensorNameChanged(self, newName):
self.setTitle(self.mSensor.name())
self.actionApplyStyle.setToolTip('Apply style to all map view images from "{}"'.format(self.mSensor.name()))
1148
1149
1150
1151
1152
1153
1154
1155
1156
1157
1158
1159
1160
1161
1162
1163
1164
1165
1166
1167
1168
1169
1170
1171
1172
1173
1174
1175
1176
1177
1178
1179
1180
1181
1182
def currentRenderWidget(self):
"""
Returns the current QgsRasterRendererWidget
:return: QgsRasterRendererWidget
"""
return self.stackedWidget.currentWidget()
def setRasterRenderer(self, renderer):
assert isinstance(renderer, QgsRasterRenderer)
assert isinstance(self.stackedWidget, QStackedWidget)
self.mMockupRasterLayer.setRenderer(renderer)
#find the widget class that fits
cls = None
for option in self.mRasterRendererModel:
if type(renderer) == option.mRenderType:
cls = option.value()
break
if cls == None:
return
widgets = []
for i in range(self.stackedWidget.count()):
w = self.stackedWidget.widget(i)
if isinstance(w, cls):
widgets.append(w)
if len(widgets) > 0:
for w in widgets:
assert isinstance(w, QgsRasterRendererWidget)
#w.setRasterLayer(self.mMockupRasterLayer)
#w.setFromRenderer(cloneRenderer(renderer))
w.setFromRenderer(renderer)
#w.doComputations()
w = widgets[0]
self.stackedWidget.setCurrentWidget(w)
self.cbRenderType.setCurrentIndex(self.stackedWidget.currentIndex())
def rasterRenderer(self):
return self.stackedWidget.currentWidget().renderer()
def apply(self):
mRendererWidget = self.currentRenderWidget()
mRendererWidget.doComputations()

benjamin.jakimow@geo.hu-berlin.de
committed
def onClipboardChanged(self):
mimeData = QApplication.clipboard().mimeData()
renderer = rendererFromXml(mimeData)
b = isinstance(renderer, QgsRasterRenderer)
#if b == False:
# print(mimeData.formats())
# s = ""
self.actionPasteStyle.setEnabled(b)

benjamin.jakimow@geo.hu-berlin.de
committed
def pasteStyleFromClipboard(self):
mimeData = QApplication.clipboard().mimeData()
renderer = rendererFromXml(mimeData)
if isinstance(renderer, QgsRasterRenderer):
self.setRasterRenderer(renderer)
def pasteStyleToClipboard(self):
xml = rendererToXml(self.rasterRenderer())
assert isinstance(xml, QDomDocument)
md = QMimeData()
#['application/qgis.style', 'text/plain']
md.setData('application/qgis.style', xml.toByteArray())
md.setData('text/plain', xml.toByteArray())
QApplication.clipboard().setMimeData(md)
if isinstance(r, QgsRasterRenderer):
for mapCanvas in self.mapCanvases():
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.addToRefreshPipeLine(MapCanvas.Command.RefreshRenderer)
RENDER_CLASSES = {}
RENDER_CLASSES['rasterrenderer'] = {
'singlebandpseudocolor': QgsSingleBandPseudoColorRenderer,
'singlebandgray': QgsSingleBandGrayRenderer,
'paletted': QgsPalettedRasterRenderer,
'multibandcolor': QgsMultiBandColorRenderer,
'hillshade': QgsHillshadeRenderer
}
RENDER_CLASSES['renderer-v2'] = {
'categorizedSymbol':QgsCategorizedSymbolRenderer,
'singleSymbol':QgsSingleSymbolRenderer
}
"""
Returns a renderer XML representation
:param renderer: QgsRasterRender | QgsFeatureRenderer
:return: QDomDocument
"""
if isinstance(renderer, QgsRasterRenderer):
1267
1268
1269
1270
1271
1272
1273
1274
1275
1276
1277
1278
1279
1280
1281
1282
1283
1284
1285
1286
1287
1288
1289
1290
1291
1292
#create a dummy raster layer
import uuid
from timeseriesviewer.virtualrasters import write_vsimem, read_vsimem
xml = """<VRTDataset rasterXSize="1" rasterYSize="1">
<GeoTransform> 0.0000000000000000e+00, 1.0000000000000000e+00, 0.0000000000000000e+00, 0.0000000000000000e+00, 0.0000000000000000e+00, -1.0000000000000000e+00</GeoTransform>
<VRTRasterBand dataType="Float32" band="1">
<Metadata>
<MDI key="STATISTICS_MAXIMUM">0</MDI>
<MDI key="STATISTICS_MEAN">0</MDI>
<MDI key="STATISTICS_MINIMUM">0</MDI>
<MDI key="STATISTICS_STDDEV">0</MDI>
</Metadata>
<Description>Band 1</Description>
<Histograms>
<HistItem>
<HistMin>0</HistMin>
<HistMax>0</HistMax>
<BucketCount>1</BucketCount>
<IncludeOutOfRange>0</IncludeOutOfRange>
<Approximate>0</Approximate>
<HistCounts>0</HistCounts>
</HistItem>
</Histograms>
</VRTRasterBand>
</VRTDataset>
"""

Benjamin Jakimow
committed
path = '/vsimem/{}.vrt'.format(uuid.uuid4())
drv = gdal.GetDriverByName('VRT')
assert isinstance(drv, gdal.Driver)
write_vsimem(path, xml)
lyr.setRenderer(renderer.clone())
lyr.exportNamedStyle(doc)
#remove dummy raster layer
lyr = None
drv.Delete(path)
elif isinstance(renderer, QgsFeatureRenderer):
#todo: distinguish vector type from requested renderer
lyr = QgsVectorLayer('Point?crs=epsg:4326&field=id:integer', 'dummy', 'memory')
lyr.setRenderer(renderer.clone())
lyr.exportNamedStyle(doc)
else:
raise NotImplementedError()
1315
1316
1317
1318
1319
1320
1321
1322
1323
1324
1325
1326
1327
1328
1329
1330
1331
1332
1333
1334
1335
1336
1337
1338
1339
1340
1341
1342
1343
1344
1345
1346
1347
1348
1349
1350
1351
1352
1353
1354
1355
1356
1357
1358
1359
1360
return doc
def rendererFromXml(xml):
"""
Reads a string `text` and returns the first QgsRasterRenderer or QgsFeatureRenderer (if defined).
:param text:
:return:
"""
if isinstance(xml, QMimeData):
for format in ['application/qgis.style', 'text/plain']:
if format in xml.formats():
dom = QDomDocument()
dom.setContent(xml.data(format))
return rendererFromXml(dom)
return None
elif isinstance(xml, str):
dom = QDomDocument()
dom.setContent(xml)
return rendererFromXml(dom)
assert isinstance(xml, QDomDocument)
root = xml.documentElement()
for baseClass, renderClasses in RENDER_CLASSES.items():
elements = root.elementsByTagName(baseClass)
if elements.count() > 0:
elem = elements.item(0).toElement()
typeName = elem.attributes().namedItem('type').nodeValue()
if typeName in renderClasses.keys():
rClass = renderClasses[typeName]
if baseClass == 'rasterrenderer':
return rClass.create(elem, DUMMY_RASTERINTERFACE)
elif baseClass == 'renderer-v2':
context = QgsReadWriteContext()
return rClass.load(elem, context)
#return rClass.create(elem)
else:
print(typeName)
s =""
return None

benjamin.jakimow@geo.hu-berlin.de
committed
class DatumViewUI(QFrame, loadUI('timeseriesdatumview.ui')):
1364
1365
1366
1367
1368
1369
1370
1371
1372
1373
1374
1375
1376
1377
1378
1379
1380
1381
1382
1383
1384
1385
1386
1387
1388
1389
1390
1391
"""
Widget to host the MapCanvases of all map views that relate to a single Datum-Sensor combinbation.
"""
def __init__(self, title='<#>', parent=None):
super(DatumViewUI, self).__init__(parent)
self.setupUi(self)
def sizeHint(self):
m = self.layout().contentsMargins()
s = QSize(0, 0)
map = None
widgets = [self.layout().itemAt(i).widget() for i in range(self.layout().count())]
widgets = [w for w in widgets if isinstance(w, QWidget)]
maps = [w for w in widgets if isinstance(w, MapCanvas)]
others = [w for w in widgets if not isinstance(w, MapCanvas)]
s = self.layout().spacing()
m = self.layout().contentsMargins()
def totalHeight(widgetList):
total = QSize(0,0)
for w in widgetList:
ws = w.size()
if ws.width() == 0:
ws = w.sizeHint()
if w.isVisible():
total.setWidth(max([total.width(), ws.width()]))
total.setHeight(total.height() + ws.height())
return total
baseSize = totalHeight(widgets)
if baseSize.width() == 0:
for o in others:
baseSize.setWidth(9999)
s = QSize(baseSize.width() + m.left() + m.right(),
baseSize.height() + m.top() + m.bottom())
1404
1405
1406
1407
1408
1409
1410
1411
1412
1413
1414
1415
1416
1417
1418
1419
1420
1421
1422
1423
1424
1425
1426
1427
1428
1429
1430
1431
1432
1433
return s
"""
def sizeHint(self):
if not self.ui.isVisible():
return QSize(0,0)
else:
#return self.ui.sizeHint()
size = self.ui.sizeHint()
s = self.ui.layout().spacing()
m = self.ui.layout().contentsMargins()
dx = m.left() + m.right() + s
dy = self.ui.layout().spacing()
n = len([m for m in self.mapCanvases.keys() if m.isVisible()])
if n > 0:
baseSize = self.mapCanvases.values()[0].size()
size = QSize(baseSize.width()+ dx, \
size.height()+ (n+1)*(dy+2*s))
else:
s = ""
return size
"""

benjamin.jakimow@geo.hu-berlin.de
committed
sigRenderProgress = pyqtSignal(int,int)
sigLoadingStarted = pyqtSignal(MapView, TimeSeriesDatum)
sigLoadingFinished = pyqtSignal(MapView, TimeSeriesDatum)
sigVisibilityChanged = pyqtSignal(bool)
def __init__(self, timeSeriesDatum:TimeSeriesDatum, stv, parent=None):
assert isinstance(timeSeriesDatum, TimeSeriesDatum)
assert isinstance(stv, SpatialTemporalVisualization)

benjamin.jakimow@geo.hu-berlin.de
committed

benjamin.jakimow@geo.hu-berlin.de
committed
self.ui.create()
self.showLoading(False)
self.wOffset = self.ui.layout().count()-1

benjamin.jakimow@geo.hu-berlin.de
committed
self.minHeight = self.ui.height()

benjamin.jakimow@geo.hu-berlin.de
committed
self.renderProgress = dict()
assert isinstance(stv, SpatialTemporalVisualization)
self.STV = stv
self.mSensor.sigNameChanged.connect(lambda :self.setColumnInfo())
self.TSD.sigVisibilityChanged.connect(self.setVisibility)

benjamin.jakimow@geo.hu-berlin.de
committed
self.setColumnInfo()
self.mRenderState = dict()

benjamin.jakimow@geo.hu-berlin.de
committed

benjamin.jakimow@geo.hu-berlin.de
committed
def setColumnInfo(self):
labelTxt = '{}\n{}'.format(str(self.TSD.mDate), self.TSD.mSensor.name())
tooltip = '\n'.join([tss.uri()for tss in self.TSD.sources()])

benjamin.jakimow@geo.hu-berlin.de
committed
self.ui.labelTitle.setText(labelTxt)
self.ui.labelTitle.setToolTip(tooltip)

benjamin.jakimow@geo.hu-berlin.de
committed
def setVisibility(self, b):
self.ui.setVisible(b)
self.sigVisibilityChanged.emit(b)

benjamin.jakimow@geo.hu-berlin.de
committed
def setHighlighted(self, b=True, timeout=1000):
styleOn = """.QFrame {
border: 4px solid red;
border-radius: 4px;
}"""
styleOff = """"""
if b is True:
self.ui.setStyleSheet(styleOn)
if timeout > 0:
QTimer.singleShot(timeout, lambda : self.setHighlighted(b=False))
else:
self.ui.setStyleSheet(styleOff)

benjamin.jakimow@geo.hu-berlin.de
committed
def removeMapView(self, mapView):
self.ui.layout().removeWidget(canvas)

benjamin.jakimow@geo.hu-berlin.de
committed
canvas.close()
def mapCanvases(self)->list:
"""
Retuns the MapCanvases of this DataView
:return: [list-of-MapCanvases]
"""
return self.mMapCanvases.values()
def refresh(self):
"""
Refreshes the MapCanvases in this DatumView, if they are not hidden behind a scroll area.
"""

benjamin.jakimow@geo.hu-berlin.de
committed
if self.ui.isVisible():
for c in self.mapCanvases():

benjamin.jakimow@geo.hu-berlin.de
committed
if c.isVisible():

benjamin.jakimow@geo.hu-berlin.de
committed
def insertMapView(self, mapView):
assert isinstance(mapView, MapView)
from timeseriesviewer.mapcanvas import MapCanvas

benjamin.jakimow@geo.hu-berlin.de
committed
mapCanvas.setObjectName('MapCanvas {} {}'.format(mapView.title(), self.TSD.mDate))

Benjamin Jakimow
committed
self.registerMapCanvas(mapView, mapCanvas)
mapCanvas.setMapView(mapView)
mapCanvas.setTSD(self.TSD)

benjamin.jakimow@geo.hu-berlin.de
committed
mapView.registerMapCanvas(self.mSensor, mapCanvas)
mapCanvas.renderComplete.connect(lambda : self.onRenderingChange(False))
mapCanvas.renderStarting.connect(lambda : self.onRenderingChange(True))
#mapCanvas.sigMapRefreshed[float, float].connect(
# lambda dt: self.STV.TSV.ui.dockSystemInfo.addTimeDelta('Map {}'.format(self.mSensor.name()), dt))
#mapCanvas.sigMapRefreshed.connect(
# lambda dt: self.STV.TSV.ui.dockSystemInfo.addTimeDelta('All Sensors', dt))

benjamin.jakimow@geo.hu-berlin.de
committed
def showLoading(self, b):
if b:
self.ui.progressBar.setRange(0, 0)
self.ui.progressBar.setValue(-1)
else:
self.ui.progressBar.setRange(0,1)
self.ui.progressBar.setValue(0)
def onRenderingChange(self, b):
mc = self.sender()
#assert isinstance(mc, QgsMapCanvas)
self.mRenderState[mc] = b
self.showLoading(any(self.mRenderState.values()))
def onRendering(self, *args):
renderFlags = [m.renderFlag() for m in self.mMapCanvases.values()]
drawFlags = [m.isDrawing() for m in self.mMapCanvases.values()]

benjamin.jakimow@geo.hu-berlin.de
committed
# print((renderFlags, drawFlags))
isLoading = any(renderFlags)
self.showLoading(isLoading)
s = ""

benjamin.jakimow@geo.hu-berlin.de
committed
def registerMapCanvas(self, mapView, mapCanvas):
from timeseriesviewer.mapcanvas import MapCanvas
assert isinstance(mapCanvas, MapCanvas)

Benjamin Jakimow
committed
mapCanvas.setVisible(mapView.isVisible())
#mapView.sigTitleChanged.connect(lambda title : mapCanvas.setSaveFileName('{}_{}'.format(self.TSD.date, title)))
mapCanvas.mapLayerModel().addMapLayerSources(self.TSD.qgsMimeDataUtilsUris())
#self.ui.layout().insertWidget(self.wOffset + len(self.mapCanvases), mapCanvas)
self.ui.layout().insertWidget(self.ui.layout().count() - 1, mapCanvas)
self.ui.update()
#register signals handled on (this) DV level
mapCanvas.renderStarting.connect(lambda: self.sigLoadingStarted.emit(mapView, self.TSD))
mapCanvas.mapCanvasRefreshed.connect(lambda: self.sigLoadingFinished.emit(mapView, self.TSD))
mapCanvas.sigShowProfiles.connect(lambda c, t : mapView.sigShowProfiles.emit(c,mapCanvas, t))
mapCanvas.sigChangeDVRequest.connect(self.onMapCanvasRequest)

benjamin.jakimow@geo.hu-berlin.de
committed
def onMapCanvasRequest(self, mapCanvas, key):
if key == 'hide_date':
self.TSD.setVisibility(False)
QApplication.clipboard().setText(self.TSD.mSensor.name())
QApplication.clipboard().setText('\n'.join(self.TSD.sourceUris()))

benjamin.jakimow@geo.hu-berlin.de
committed
def __lt__(self, other):
assert isinstance(other, DatumView)
return self.TSD < other.TSD

benjamin.jakimow@geo.hu-berlin.de
committed
def __eq__(self, other):
assert isinstance(other, DatumView)
return self.TSD == other.TSD

benjamin.jakimow@geo.hu-berlin.de
committed
class SpatialTemporalVisualization(QObject):
"""
"""
sigLoadingStarted = pyqtSignal(DatumView, MapView)
sigLoadingFinished = pyqtSignal(DatumView, MapView)
sigShowProfiles = pyqtSignal(SpatialPoint, MapCanvas, str)

benjamin.jakimow@geo.hu-berlin.de
committed
sigShowMapLayerInfo = pyqtSignal(dict)
sigSpatialExtentChanged = pyqtSignal(SpatialExtent)
sigCRSChanged = pyqtSignal(QgsCoordinateReferenceSystem)

benjamin.jakimow@geo.hu-berlin.de
committed
def __init__(self, timeSeriesViewer):
super(SpatialTemporalVisualization, self).__init__()

benjamin.jakimow@geo.hu-berlin.de
committed
#assert isinstance(timeSeriesViewer, TimeSeriesViewer), timeSeriesViewer
#default map settings
self.mSpatialExtent = SpatialExtent.world()
self.mCRS = self.mSpatialExtent.crs()

Benjamin Jakimow
committed
self.mSize = QSize(200, 200)
self.mColor = Qt.black

benjamin.jakimow@geo.hu-berlin.de
committed
self.ui = timeSeriesViewer.ui
#map-tool handling
self.mMapToolActivator = None
self.mMapTools = []

benjamin.jakimow@geo.hu-berlin.de
committed
self.scrollArea = self.ui.scrollAreaSubsets
assert isinstance(self.scrollArea, MapViewScrollArea)
#self.scrollArea.sigResized.connect(self.refresh())
#self.scrollArea.horizontalScrollBar().valueChanged.connect(self.mRefreshTimer.start)

benjamin.jakimow@geo.hu-berlin.de
committed
self.TSV = timeSeriesViewer

benjamin.jakimow@geo.hu-berlin.de
committed
self.targetLayout = self.ui.scrollAreaSubsetContent.layout()
#self.MVC = MapViewCollection(self)
#self.MVC.sigShowProfiles.connect(self.sigShowProfiles.emit)
self.MVC = self.ui.dockMapViews
assert isinstance(self.MVC, MapViewCollectionDock)

benjamin.jakimow@geo.hu-berlin.de
committed
self.MVC.sigShowProfiles.connect(self.sigShowProfiles.emit)
self.MVC.sigMapViewAdded.connect(self.onMapViewAdded)

benjamin.jakimow@geo.hu-berlin.de
committed
self.vectorOverlay = None
self.DVC.sigResizeRequired.connect(self.adjustScrollArea)

Benjamin Jakimow
committed
#self.DVC.sigLoadingStarted.connect(self.ui.dockRendering.addStartedWork)
#self.DVC.sigLoadingFinished.connect(self.ui.dockRendering.addFinishedWork)
#self.timeSeriesDateViewCollection.sigSpatialExtentChanged.connect(self.setSpatialExtent)
self.TS.sigTimeSeriesDatesAdded.connect(self.DVC.addDates)
self.TS.sigTimeSeriesDatesRemoved.connect(self.DVC.removeDates)

benjamin.jakimow@geo.hu-berlin.de
committed
#add dates, if already existing

benjamin.jakimow@geo.hu-berlin.de
committed
if len(self.TS) > 0:

benjamin.jakimow@geo.hu-berlin.de
committed
self.mMapRefreshTimer = QTimer(self)
self.mMapRefreshTimer.timeout.connect(self.timedCanvasRefresh)
self.mMapRefreshTimer.setInterval(500)
self.mNumberOfHiddenMapsToRefresh = 2
def timedCanvasRefresh(self, *args, force:bool=False):
#do refresh maps
assert isinstance(self.scrollArea, MapViewScrollArea)
visibleMaps = [m for m in self.mapCanvases() if m.isVisibleToViewport()]
hiddenMaps = sorted([m for m in self.mapCanvases() if not m.isVisibleToViewport()],
n = 0
#redraw all visible maps
for c in visibleMaps:
assert isinstance(c, MapCanvas)
c.timedRefresh()
n += 1
if n < 10:
#refresh up to mNumberOfHiddenMapsToRefresh maps which are not visible to the user
i = 0
for c in hiddenMaps:
assert isinstance(c, MapCanvas)
c.timedRefresh()
i += 1
if i >= self.mNumberOfHiddenMapsToRefresh and not force:
break

Benjamin Jakimow
committed
def mapViewFromCanvas(self, mapCanvas:MapCanvas)->MapView:
"""
Returns the MapView a mapCanvas belongs to
:param mapCanvas: MapCanvas
:return: MapView
"""
for mapView in self.MVC:
assert isinstance(mapView, MapView)
if mapCanvas in mapView.mapCanvases():
return mapView
return None
def onMapViewAdded(self, *args):
self.adjustScrollArea()
s = ""

benjamin.jakimow@geo.hu-berlin.de
committed
def createMapView(self):
self.MVC.createMapView()
def registerMapCanvas(self, mapCanvas:MapCanvas):
"""
Connects a MapCanvas and its signals
:param mapCanvas: MapCanvas
"""
from timeseriesviewer.mapcanvas import MapCanvas
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.setMapLayerStore(self.TSV.mMapLayerStore)
self.mMapCanvases.append(mapCanvas)
#set general canvas properties
mapCanvas.setFixedSize(self.mSize)
mapCanvas.setSpatialExtent(self.mSpatialExtent)

Benjamin Jakimow
committed
mapCanvas.setLabelingModel(self.TSV.ui.dockLabeling.mLabelAttributeModel)
#register on map canvas signals
def onChanged(e, mapCanvas0=None):
self.setSpatialExtent(e, mapCanvas0=mapCanvas0)
#mapCanvas.sigSpatialExtentChanged.connect(lambda e: self.setSpatialExtent(e, mapCanvas0=mapCanvas))
mapCanvas.sigSpatialExtentChanged.connect(lambda e: onChanged(e, mapCanvas0=mapCanvas))

Benjamin Jakimow
committed
mapCanvas.sigCrosshairPositionChanged.connect(self.onCrosshairChanged)
def onCrosshairChanged(self, spatialPoint:SpatialPoint):

Benjamin Jakimow
committed
"""
Synchronizes all crosshair positions. Takes care of CRS differences.
:param spatialPoint: SpatialPoint of new Crosshair position

Benjamin Jakimow
committed
"""
from timeseriesviewer.crosshair import CrosshairStyle

Benjamin Jakimow
committed
srcCanvas = self.sender()
if isinstance(srcCanvas, MapCanvas):
dstCanvases = [c for c in self.mapCanvases() if c != srcCanvas]
else:
dstCanvases = [c for c in self.mapCanvases()]
if isinstance(spatialPoint, SpatialPoint):

Benjamin Jakimow
committed
for mapCanvas in dstCanvases:
mapCanvas.setCrosshairPosition(spatialPoint, emitSignal=False)

Benjamin Jakimow
committed
def setCrosshairStyle(self, crosshairStyle:CrosshairStyle):
"""
Sets a crosshair style to all map canvas
:param crosshairStyle: CrosshairStyle
"""
for mapView in self.mapViews():
assert isinstance(mapView, MapView)
mapView.setCrosshairStyle(crosshairStyle)

benjamin.jakimow@geo.hu-berlin.de
committed
"""
Sets the Crosshair visiblity
:param b: bool
"""
assert isinstance(b, bool)
self.onCrosshairChanged(b)

benjamin.jakimow@geo.hu-berlin.de
committed
def setVectorLayer(self, lyr:QgsVectorLayer):
"""
Sets a QgsVectorLaye to be shown on top of raster images
:param lyr: QgsVectorLayer
"""

benjamin.jakimow@geo.hu-berlin.de
committed
self.MVC.setVectorLayer(lyr)

benjamin.jakimow@geo.hu-berlin.de
committed
assert isinstance(size, QSize)
from timeseriesviewer.mapcanvas import MapCanvas
for mapCanvas in self.mMapCanvases:
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.setFixedSize(size)
self.sigMapSizeChanged.emit(self.mSize)

benjamin.jakimow@geo.hu-berlin.de
committed
self.adjustScrollArea()
def mapSize(self)->QSize:
"""
Returns the MapCanvas size
:return: QSize
"""

benjamin.jakimow@geo.hu-berlin.de
committed
def refresh(self):
"""
Refreshes all visible MapCanvases
"""
for c in self.mapCanvases():
assert isinstance(c, MapCanvas)
c.refresh()
#self.mMapRefreshTimer.stop()

benjamin.jakimow@geo.hu-berlin.de
committed
def adjustScrollArea(self):
"""
Adjusts the scroll area widget to fit all visible widgets
"""

benjamin.jakimow@geo.hu-berlin.de
committed
m = self.targetLayout.contentsMargins()

benjamin.jakimow@geo.hu-berlin.de
committed
w = h = 0
s = QSize()
r = None
tsdViews = [v for v in self.DVC if v.ui.isVisible()]
mapViews = [v for v in self.MVC if v.isVisible()]
nX = len(tsdViews)
nY = len(mapViews)
spacing = self.targetLayout.spacing()
margins = self.targetLayout.contentsMargins()

benjamin.jakimow@geo.hu-berlin.de
committed
sizeX = 1
sizeY = 50
if nX > 0:
s = tsdViews[0].ui.sizeHint().width()
s = nX * (s + spacing) + margins.left() + margins.right()
sizeX = s
if nY > 0 and nX > 0:
s = tsdViews[0].ui.sizeHint().height()
s = s + margins.top() + margins.bottom()

benjamin.jakimow@geo.hu-berlin.de
committed
#s = tsdViews[0].ui.sizeHint()
#s = QSize(nX * (s.width() + spacing) + margins.left() + margins.right(),
# s.height() + margins.top() + margins.bottom())

benjamin.jakimow@geo.hu-berlin.de
committed
#print(sizeX, sizeY)
self.targetLayout.parentWidget().resize(QSize(sizeX, sizeY))

benjamin.jakimow@geo.hu-berlin.de
committed
def setMapTool(self, mapToolKey, *args, **kwds):
"""
Create a maptool instance to each MapCanvas
:param mapToolKey: str which MapTool is to create, or QgsMapTool instance
:param args: optional maptool arguments
:param kwds: optional maptool keywords
:return: [list-of-QgsMapTools]
"""
self.mMapToolActivator = self.sender()
del self.mMapTools[:]

Benjamin Jakimow
committed
from timeseriesviewer.maptools import MapTools, CursorLocationMapTool, SpectralProfileMapTool, TemporalProfileMapTool
for canvas in self.mMapCanvases:
mt = None
if mapToolKey in MapTools.mapToolKeys():
mt = MapTools.create(mapToolKey, canvas, *args, **kwds)
if isinstance(mapToolKey, QgsMapTool):
mt = MapTools.copy(mapToolKey, canvas, *args, **kwds)
if isinstance(mt, QgsMapTool):
canvas.setMapTool(mt)
self.mMapTools.append(mt)
#if required, link map-tool with specific slots
if isinstance(mt, CursorLocationMapTool):
mt.sigLocationRequest[SpatialPoint, QgsMapCanvas].connect(lambda c, m : self.sigShowProfiles.emit(c,m, mapToolKey))
return self.mMapTools

benjamin.jakimow@geo.hu-berlin.de
committed
def setSpatialCenter(self, center, mapCanvas0=None):
"""
Sets the spatial center of all MapCanvases
:param center: SpatialPoint
:param mapCanvas0:
"""
assert isinstance(center, SpatialPoint)
center = center.toCrs(self.mCRS)
if not isinstance(center, SpatialPoint):
return
self.mSpatialExtent.setCenter(center)
for mapCanvas in self.mMapCanvases:
if mapCanvas != mapCanvas0:
oldState = mapCanvas.blockSignals(True)
mapCanvas.setCenter(center)
mapCanvas.blockSignals(oldState)
self.sigSpatialExtentChanged.emit(self.mSpatialExtent)
def setSpatialCenter(self, center:SpatialPoint, mapCanvas0=None):
"""
Sets the MapCanvas center.
:param center: SpatialPoint
:param mapCanvas0: MapCanvas0 optional
"""
assert isinstance(center, SpatialPoint)
center = center.toCrs(self.mCRS)
if not isinstance(center, SpatialPoint):
return None
for mapCanvas in self.mapCanvases():
assert isinstance(mapCanvas, MapCanvas)
if mapCanvas != mapCanvas0:
center0 = mapCanvas.spatialCenter()
if center0 != center:
oldState = mapCanvas.blockSignals(True)
mapCanvas.setCenter(center)
mapCanvas.blockSignals(oldState)
self.mMapRefreshTimer.start()
def setSpatialExtent(self, extent, mapCanvas0=None):
"""
Sets the spatial extent of all MapCanvases
:param extent: SpatialExtent
:param mapCanvas0:
:return:
"""
assert isinstance(extent, SpatialExtent)
extent = extent.toCrs(self.mCRS)
if not isinstance(extent, SpatialExtent) \
or extent.isEmpty() or not extent.isFinite() \
or extent.width() <= 0 \
or extent.height() <= 0 \
or extent == self.mSpatialExtent:
if self.mSpatialExtent == extent:
return
for mapCanvas in self.mapCanvases():
assert isinstance(mapCanvas, MapCanvas)
extent0 = mapCanvas.spatialExtent()
if mapCanvas != mapCanvas0 and extent0 != extent:
self.sigSpatialExtentChanged.emit(extent)

benjamin.jakimow@geo.hu-berlin.de
committed
def setBackgroundColor(self, color:QColor):
"""
Sets the MapCanvas background color
:param color: QColor
"""
assert isinstance(color, QColor)
self.mColor = color
for mapCanvas in self.mMapCanvases:
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.setCanvasColor(color)
def backgroundColor(self)->QColor:
"""
Returns the MapCanvas background color
:return: QColor
"""
return self.mColor

Benjamin Jakimow
committed
def mapCanvases(self, mapView=None)->list:
"""
Returns MapCanvases
:param mapView: a MapView to return MapCanvases from only, defaults to None
:return: [list-of-MapCanvas]
"""
if isinstance(mapView, MapView):
s = ""

Benjamin Jakimow
committed
def mapViews(self)->list:
"""
Returns a list of all mapviews
:return [list-of-MapViews]:
"""
return self.MVC[:]
def setCrs(self, crs):
assert isinstance(crs, QgsCoordinateReferenceSystem)
transform = QgsCoordinateTransform()
transform.setSourceCrs(self.mCRS)
transform.setDestinationCrs(crs)
if transform.isValid() and not transform.isShortCircuited():

Benjamin Jakimow
committed
for mapCanvas in self.mapCanvases():
# print(('STV set CRS {} {}', str(mapCanvas), self.mCRS.description()))
mapCanvas.setDestinationCrs(QgsCoordinateReferenceSystem(crs))
"""
from timeseriesviewer.utils import saveTransform
if saveTransform(self.mSpatialExtent, self.mCRS, crs):
self.mCRS = crs
else:
pass
"""
self.sigCRSChanged.emit(self.crs())
def crs(self)->QgsCoordinateReferenceSystem:
"""
Returns the QgsCoordinateReferenceSystem
:return: QgsCoordinateReferenceSystem
"""
def spatialExtent(self)->SpatialExtent:
"""
Returns the SpatialExtent
:return: SpatialExtent
"""

benjamin.jakimow@geo.hu-berlin.de
committed
return self.mSpatialExtent

benjamin.jakimow@geo.hu-berlin.de
committed
def navigateToTSD(self, TSD:TimeSeriesDatum):
"""
Changes the viewport of the scroll window to show the requested TimeSeriesDatum
:param TSD: TimeSeriesDatum
"""

benjamin.jakimow@geo.hu-berlin.de
committed
assert isinstance(TSD, TimeSeriesDatum)
#get widget related to TSD

benjamin.jakimow@geo.hu-berlin.de
committed
assert isinstance(self.scrollArea, QScrollArea)
self.scrollArea.ensureWidgetVisible(tsdv.ui)

benjamin.jakimow@geo.hu-berlin.de
committed
sigResizeRequired = pyqtSignal()
sigLoadingStarted = pyqtSignal(MapView, TimeSeriesDatum)
sigLoadingFinished = pyqtSignal(MapView, TimeSeriesDatum)

benjamin.jakimow@geo.hu-berlin.de
committed
sigShowProfiles = pyqtSignal(SpatialPoint)
sigSpatialExtentChanged = pyqtSignal(SpatialExtent)

benjamin.jakimow@geo.hu-berlin.de
committed
def __init__(self, STViz):
assert isinstance(STViz, SpatialTemporalVisualization)
super(DateViewCollection, self).__init__()

benjamin.jakimow@geo.hu-berlin.de
committed
#self.tsv = tsv
#self.timeSeries = tsv.TS
self.STV = STViz
self.ui = self.STV.targetLayout.parentWidget()

benjamin.jakimow@geo.hu-berlin.de
committed
self.scrollArea = self.ui.parentWidget().parentWidget()
#potentially there are many more dates than views.
#therefore we implement the addinng/removing of mapviews here
#we reduce the number of layout refresh calls by
#suspending signals, adding the new map view canvases, and sending sigResizeRequired
self.STV.MVC.sigMapViewAdded.connect(self.addMapView)
self.STV.MVC.sigMapViewRemoved.connect(self.removeMapView)

benjamin.jakimow@geo.hu-berlin.de
committed
self.setFocusView(None)

benjamin.jakimow@geo.hu-berlin.de
committed
def tsdFromMapCanvas(self, mapCanvas):
assert isinstance(mapCanvas, MapCanvas)
assert isinstance(view, DatumView)
if mapCanvas in view.mMapCanvases.values():
return view.TSD
return None

benjamin.jakimow@geo.hu-berlin.de
committed
def tsdView(self, tsd):
r = [v for v in self.mViews if v.TSD == tsd]

benjamin.jakimow@geo.hu-berlin.de
committed
if len(r) == 1:
return r[0]
else:
raise Exception('TSD not in list')
def addMapView(self, mapView):
assert isinstance(mapView, MapView)
w = self.ui
#w.setUpdatesEnabled(False)
#for tsdv in self.mViews:
# tsdv.ui.setUpdatesEnabled(False)

benjamin.jakimow@geo.hu-berlin.de
committed

benjamin.jakimow@geo.hu-berlin.de
committed
tsdv.insertMapView(mapView)
#for tsdv in self.mViews:
# tsdv.ui.setUpdatesEnabled(True)

benjamin.jakimow@geo.hu-berlin.de
committed
#mapView.sigSensorRendererChanged.connect(lambda *args : self.setRasterRenderer(mapView, *args))
mapView.sigMapViewVisibility.connect(lambda: self.sigResizeRequired.emit())
mapView.sigShowProfiles.connect(self.sigShowProfiles.emit)

benjamin.jakimow@geo.hu-berlin.de
committed
self.sigResizeRequired.emit()
def removeMapView(self, mapView):
assert isinstance(mapView, MapView)

benjamin.jakimow@geo.hu-berlin.de
committed
tsdv.removeMapView(mapView)
self.sigResizeRequired.emit()

benjamin.jakimow@geo.hu-berlin.de
committed
def highlightDate(self, tsd):
"""
Highlights a time series data for a specific time our
:param tsd:
:return:
"""
tsdView = self.tsdView(tsd)
if isinstance(tsdView, DatumView):
tsdView.setHighlight(True)

benjamin.jakimow@geo.hu-berlin.de
committed
def setFocusView(self, tsd):
self.focusView = tsd
def orderedViews(self):
#returns the
if self.focusView is not None:
assert isinstance(self.focusView, DatumView)
return sorted(self.mViews, key=lambda v: np.abs(v.TSD.date - self.focusView.TSD.date))

benjamin.jakimow@geo.hu-berlin.de
committed
else:

benjamin.jakimow@geo.hu-berlin.de
committed

benjamin.jakimow@geo.hu-berlin.de
committed
def setSubsetSize(self, size):
assert isinstance(size, QSize)
self.subsetSize = size
for tsdView in self.orderedViews():
tsdView.blockSignals(True)
for tsdView in self.orderedViews():
tsdView.setSubsetSize(size)
for tsdView in self.orderedViews():
tsdView.blockSignals(False)

benjamin.jakimow@geo.hu-berlin.de
committed
def addDates(self, tsdList):
"""
Create a new TSDView
:param tsdList:
:return:
"""
for tsd in tsdList:
assert isinstance(tsd, TimeSeriesDatum)
#tsdView.setSubsetSize(self.subsetSize)
DV.sigLoadingStarted.connect(self.sigLoadingStarted.emit)
DV.sigLoadingFinished.connect(self.sigLoadingFinished.emit)
DV.sigVisibilityChanged.connect(lambda: self.STV.adjustScrollArea())

benjamin.jakimow@geo.hu-berlin.de
committed
for i, mapView in enumerate(self.STV.MVC):
DV.insertMapView(mapView)

benjamin.jakimow@geo.hu-berlin.de
committed
bisect.insort(self.mViews, DV)
i = self.mViews.index(DV)

benjamin.jakimow@geo.hu-berlin.de
committed
DV.ui.setParent(self.STV.targetLayout.parentWidget())
self.STV.targetLayout.insertWidget(i, DV.ui)
DV.ui.show()

benjamin.jakimow@geo.hu-berlin.de
committed
if len(tsdList) > 0:
self.sigResizeRequired.emit()
def removeDates(self, tsdList):
toRemove = [v for v in self.mViews if v.TSD in tsdList]

benjamin.jakimow@geo.hu-berlin.de
committed
removedDates = []
Loading
Loading full blame...