# -*- coding: utf-8 -*-
# noinspection PyPep8Naming
"""
/***************************************************************************
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
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 *
import numpy as np
from .utils import *
from .import Option, OptionListModel
from .timeseries import SensorInstrument, TimeSeriesDate, TimeSeries, SensorProxyLayer
from .utils import loadUI
from .mapviewscrollarea import MapViewScrollArea
from .mapcanvas import MapCanvas, MapTools, MapCanvasInfoItem, MapCanvasMapTools
from .externals.qps.crosshair.crosshair import getCrosshairStyle, CrosshairStyle, CrosshairMapCanvasItem
from .externals.qps.layerproperties import showLayerPropertiesDialog
from .externals.qps.maptools import *
#assert os.path.isfile(dummyPath)
#lyr = QgsRasterLayer(dummyPath)
#assert lyr.isValid()
DUMMY_RASTERINTERFACE = QgsSingleBandGrayRenderer(None, 0)
KEY_LOCKED_LAYER = 'eotsv/locked'
KEY_SENSOR_GROUP = 'eotsv/sensorgroup'
KEY_SENSOR_LAYER = 'eotsv/sensorlayer'
def equalTextFormats(tf1:QgsTextFormat, tf2:QgsTextFormat)->True:
return tf1.toMimeData().text() == tf2.toMimeData().text()
class MapViewLayerTreeViewMenuProvider(QgsLayerTreeViewMenuProvider):
def __init__(self, mapView, view: QgsLayerTreeView, canvas: QgsMapCanvas):
super(MapViewLayerTreeViewMenuProvider, self).__init__()
assert isinstance(view, QgsLayerTreeView)
assert isinstance(canvas, QgsMapCanvas)
self.mLayerTreeView = view
self.mDummyCanvas = canvas
self.mDefActions = QgsLayerTreeViewDefaultActions(self.mLayerTreeView)
self.mMapView = mapView
self.actionAddGroup = self.mDefActions.actionAddGroup()
self.actionRename = self.mDefActions.actionRenameGroupOrLayer()
self.actionRemove = self.mDefActions.actionRemoveGroupOrLayer()
#self.actionZoomToLayer = self.mDefActions.actionZoomToGroup(self.mDummyCanvas)
self.actionCheckAndAllChildren = self.mDefActions.actionCheckAndAllChildren()
self.actionShowFeatureCount = self.mDefActions.actionShowFeatureCount()
#self.actionZoomToLayer = self.mDefActions.actionZoomToLayer(self.mDummyCanvas)
#self.actionZoomToSelected = self.mDefActions.actionZoomToSelection(self.mDummyCanvas)
#self.actionZoomToGroup = self.mDefActions.actionZoomToGroup(self.mDummyCanvas)
self.actionAddEOTSVSpectralProfiles = QAction('Add Spectral Profile Layer')
self.actionAddEOTSVTemporalProfiles = QAction('Add Temporal Profile Layer')
def mapView(self):
return self.mMapView
def layerTreeView(self)->QgsLayerTreeView:
return self.mLayerTreeView
def layerTree(self)->QgsLayerTree:
return self.layerTreeModel().rootGroup()
def layerTreeModel(self)->QgsLayerTreeModel:
return self.layerTreeView().model()
def createContextMenu(self)->QMenu:
model = self.layerTreeModel()
ltree = self.layerTree()
view = self.layerTreeView()
g = view.currentGroupNode()
l = view.currentLayer()
i = view.currentIndex()
#fixedNodes = len([l for l in view.selectedLayersRecursive() if l.property(KEY_LOCKED_LAYER) == True]) > 0 or \
# isinstance(g, QgsLayerTreeGroup) and g.property(KEY_LOCKED_LAYER) == True
# disable actions
#self.actionRemove.setEnabled(fixedNodes == False)
menu = QMenu(view)
isSensorGroup = isinstance(g, QgsLayerTreeGroup) and g.customProperty(KEY_SENSOR_GROUP) in [True, 'true']
isSensorLayer = isinstance(l, QgsRasterLayer) and l.customProperty(KEY_SENSOR_LAYER) in [True, 'true']
self.actionRemove.setEnabled(not (isSensorGroup or isSensorLayer))
self.actionAddGroup.setEnabled(not (isSensorGroup or isSensorLayer))
menu.addAction(self.actionAddGroup)
menu.addAction(self.actionRename)
menu.addAction(self.actionRemove)
#menu.addAction(self.actionZoomToGroup)
#menu.addAction(self.actionZoomToLayer)
#menu.addAction(self.actionZoomToSelected)
menu.addSeparator()
menu.addAction(self.actionAddEOTSVSpectralProfiles)
menu.addAction(self.actionAddEOTSVTemporalProfiles)
menu.addSeparator()
centerCanvas = None
if isinstance(self.mapView(), MapView):
visibleCanvases = self.mapView().visibleMapCanvases()
if len(visibleCanvases) > 0:
i = int(len(visibleCanvases) / 2)
centerCanvas = visibleCanvases[i]
a = menu.addAction('Set Properties')
a.triggered.connect(lambda *args,
canvas = centerCanvas,
lyr = l,
b = not isinstance(l, SensorProxyLayer):
showLayerPropertiesDialog(lyr, canvas, useQGISDialog=b))
a.setEnabled(isinstance(centerCanvas, QgsMapCanvas))
from .externals.qps.layerproperties import pasteStyleFromClipboard, pasteStyleToClipboard
a = menu.addAction('Copy Style')
a.setToolTip('Copy the current layer style to clipboard')
a.triggered.connect(lambda *args, lyr=l: pasteStyleToClipboard(lyr))
a = menu.addAction('Paste Style')
a.setEnabled('application/qgis.style' in QApplication.clipboard().mimeData().formats())
a.triggered.connect(lambda *args, lyr=l: pasteStyleFromClipboard(lyr))
#a = menu.addAction('Settings')
#from qps.layerproperties import showLayerPropertiesDialog
#a.triggered.connect(lambda *args, lyr=l:showLayerPropertiesDialog(lyr, self._canvas))
return menu
class MapViewLayerTreeModel(QgsLayerTreeModel):
"""
Layer Tree as shown in a MapView
"""
def __init__(self, rootNode, parent=None):
super(MapViewLayerTreeModel, self).__init__(rootNode, parent=parent)
def dataXXX(self, index:QModelIndex, role=Qt.DisplayRole):
node = self.index2node(index)
# if node.name() == 'testlayer':
# s = ""
if True:
if isinstance(node, QgsLayerTreeGroup) and node.customProperty(KEY_SENSOR_GROUP) in ['true', True]:
if role == Qt.FontRole:
f = super(MapViewLayerTreeModel, self).data(index, role=role)
f.setBold(True)
return f
if isinstance(node, QgsLayerTreeLayer) and node.customProperty(KEY_SENSOR_LAYER) in ['true', True]:
if role == Qt.FontRole:
f = super(MapViewLayerTreeModel, self).data(index, role=role)
assert isinstance(f, QFont)
f.setItalic(True)
return f
if role == Qt.DecorationRole:
return QIcon(':/timeseriesviewer/icons/icon.svg')
return super(MapViewLayerTreeModel, self).data(index, role=role)
def flagsXXX(self, index:QModelIndex):
f = super(MapViewLayerTreeModel, self).flags(index)
node = self.index2node(index)
if isinstance(node, QgsLayerTreeNode) and ( \
node.customProperty(KEY_SENSOR_LAYER) in ['true', True] or \
node.customProperty(KEY_SENSOR_GROUP) in ['true', True]):
f = f ^ Qt.ItemIsDragEnabled
f = f ^ Qt.ItemIsDropEnabled
return f
class MapView(QFrame, loadUIFormClass(jp(DIR_UI, 'mapview.ui'))):
"""
A MapView defines how a single map canvas visualizes sensor specific EOTS data plus additional vector overlays
"""
#sigVisibilityChanged = pyqtSignal(bool)
sigCanvasAppearanceChanged = pyqtSignal()
sigCrosshairChanged = pyqtSignal()
sigTitleChanged = pyqtSignal(str)
sigSensorRendererChanged = pyqtSignal(SensorInstrument, QgsRasterRenderer)
sigShowProfiles = pyqtSignal(SpatialPoint, MapCanvas, str)
def __init__(self, name='Map View', parent=None):
super(MapView, self).__init__(parent)
self.setupUi(self)
from eotimeseriesviewer.settings import DEFAULT_VALUES, Keys
self.mMapBackgroundColor = DEFAULT_VALUES[Keys.MapBackgroundColor]
self.mMapTextFormat = DEFAULT_VALUES[Keys.MapTextFormat]
self.mCurrentLayer = None
self.mMapWidget = None
self.mTimeSeries = None
self.mSensorLayerList = list()
self.mCrossHairStyle = CrosshairStyle()
m = QMenu(self.btnToggleCrosshair)
m.addAction(self.actionSetCrosshairStyle)
self.btnToggleCrosshair.setMenu(m)
self.btnToggleCrosshair.setDefaultAction(self.actionToggleCrosshairVisibility)
self.btnToggleCrosshair.setChecked(self.crosshairStyle().mShow)
self.btnToggleMapViewVisibility.setDefaultAction(self.actionToggleMapViewHidden)
self.tbName.textChanged.connect(self.onTitleChanged)
self.actionSetCrosshairStyle.triggered.connect(self.onChangeCrosshairStyle)
self.actionToggleMapViewHidden.toggled.connect(self.sigCanvasAppearanceChanged)
self.actionToggleCrosshairVisibility.toggled.connect(self.setCrosshairVisibility)
self.actionAddMapLayer.triggered.connect(lambda *args: self.onAddMapLayer())
self.actionAddVectorLayer.triggered.connect(lambda *args: self.onAddMapLayer(QgsMapLayerProxyModel.VectorLayer))
self.actionAddRasterLayer.triggered.connect(lambda *args: self.onAddMapLayer(QgsMapLayerProxyModel.RasterLayer))
self.btnAddLayer.setDefaultAction(self.actionAddMapLayer)
m = QMenu()
m.addAction(self.actionAddVectorLayer)
m.addAction(self.actionAddRasterLayer)
self.btnAddLayer.setMenu(m)
self.btnHighlightMapView.setDefaultAction(self.actionHighlightMapView)
self.actionHighlightMapView.triggered.connect(lambda: self.setHighlighted(True, timeout=500))
assert isinstance(self.mLayerTreeView, QgsLayerTreeView)
self.mDummyCanvas = QgsMapCanvas() # dummy map canvas for dummy layers
self.mDummyCanvas.setVisible(False)
self.mLayerTree = QgsLayerTree()
self.mLayerTreeMapCanvasBridget = QgsLayerTreeMapCanvasBridge(self.mLayerTree, self.mDummyCanvas)
# self.mLayerTreeModel = QgsLayerTreeModel(self.mLayerTree)
self.mLayerTreeModel = MapViewLayerTreeModel(self.mLayerTree)
self.mLayerTreeModel.setFlags(QgsLayerTreeModel.AllowNodeChangeVisibility |
QgsLayerTreeModel.AllowNodeRename |
QgsLayerTreeModel.AllowNodeReorder)
self._createSensorNode()
self.mLayerTreeView.setModel(self.mLayerTreeModel)
self.mMapLayerTreeViewMenuProvider = MapViewLayerTreeViewMenuProvider(self, self.mLayerTreeView, self.mDummyCanvas)
# register some actions that interact with other GUI elements
self.mMapLayerTreeViewMenuProvider.actionAddEOTSVSpectralProfiles.triggered.connect(self.addSpectralProfileLayer)
self.mMapLayerTreeViewMenuProvider.actionAddEOTSVTemporalProfiles.triggered.connect(self.addTemporalProfileLayer)
self.mLayerTreeView.setMenuProvider(self.mMapLayerTreeViewMenuProvider)
self.mLayerTreeView.currentLayerChanged.connect(self.setCurrentLayer)
self.mLayerTree.removedChildren.connect(self.onChildNodesRemoved)
self.mIsVisible = True
self.setTitle(name)
m = QMenu()
m.addAction(self.optionShowDate)
m.addAction(self.optionShowSensorName)
m.addAction(self.optionShowMapViewName)
self.btnInfoOptions.setMenu(m)
for action in m.actions():
action.toggled.connect(self.sigCanvasAppearanceChanged)
fixMenuButtons(self)
def setName(self, name:str):
self.setTitle(name)
def name(self)->str:
return self.title()
def setMapTextFormat(self, textformat:QgsTextFormat)->QgsTextFormat:
if not equalTextFormats(self.mapTextFormat(), textformat):
self.mMapTextFormat = textformat
self.sigCanvasAppearanceChanged.emit()
return self.mapTextFormat()
def mapTextFormat(self)->QgsTextFormat:
return self.mMapTextFormat
def mapBackgroundColor(self)->QColor:
"""
Returns the map background color
:return: QColor
"""
return self.mMapBackgroundColor
def setMapBackgroundColor(self, color:QColor)->QColor:
"""
Sets the map background color
:param color: QColor
:return: QColor
"""
if self.mMapBackgroundColor != color:
self.mMapBackgroundColor = color
self.sigCanvasAppearanceChanged.emit()
return self.mMapBackgroundColor
def visibleMapCanvases(self)->list:
"""
Returns the currently visible mapcanvases
:return: [list-of-MapCanvases]
"""
return [m for m in self.mapCanvases() if m.isVisibleToViewport()]
def onAddMapLayer(self, filter:QgsMapLayerProxyModel.Filter=QgsMapLayerProxyModel.All):
"""
Slot that opens a SelectMapLayersDialog for any kind of layer
"""
from .externals.qps.utils import SelectMapLayersDialog
d = SelectMapLayersDialog()
if filter == QgsMapLayerProxyModel.All:
title = 'Select Layer'
text = 'Layer'
elif filter == QgsMapLayerProxyModel.RasterLayer:
title = 'Select Raster Layer'
text = 'Raster'
elif filter == QgsMapLayerProxyModel.VectorLayer:
title = 'Select Vector Layer'
text = 'Vector'
d.setWindowTitle(title)
d.addLayerDescription(text, filter)
if d.exec() == QDialog.Accepted:
for l in d.mapLayers():
self.addLayer(l)
def setCurrentLayer(self, layer:QgsMapLayer):
"""
Sets the QgsMapCanvas.currentLayer() that is used by some QgsMapTools
:param layer: QgsMapLayer | None
:return:
"""
assert layer is None or isinstance(layer, QgsMapLayer)
self.mCurrentLayer = layer
if layer not in self.mSensorLayerList:
for c in self.mapCanvases():
c.setCurrentLayer(layer)
else:
s = ""
def addSpectralProfileLayer(self):
"""Adds the EOTSV Spectral Profile Layer"""
from eotimeseriesviewer.main import TimeSeriesViewer
tsv = TimeSeriesViewer.instance()
if isinstance(tsv, TimeSeriesViewer):
lyr = tsv.spectralLibrary()
if lyr not in self.layers():
self.addLayer(lyr)
def addTemporalProfileLayer(self):
"""Adds the EOTSV Temporal Profile Layer"""
from eotimeseriesviewer.main import TimeSeriesViewer
tsv = TimeSeriesViewer.instance()
if isinstance(tsv, TimeSeriesViewer):
lyr = tsv.temporalProfileLayer()
if lyr not in self.layers():
self.addLayer(lyr)
def addLayer(self, layer:QgsMapLayer):
"""
Add a QgsMapLayer to the MapView layer tree
:param layer: QgsMapLayer
"""
if isinstance(layer, QgsVectorLayer):
self.mLayerTree.insertLayer(0, layer)
else:
self.mLayerTree.addLayer(layer)
def _createSensorNode(self):
self.mLayerTreeSensorNode = QgsLayerTreeGroup(name='Raster Time Series', checked=True)
self.mLayerTreeSensorNode.setCustomProperty(KEY_LOCKED_LAYER, True)
self.mLayerTreeSensorNode.setCustomProperty(KEY_SENSOR_GROUP, True)
self.mLayerTree.addChildNode(self.mLayerTreeSensorNode)
def _containsSensorNode(self, root:QgsLayerTreeGroup)->bool:
assert isinstance(root, QgsLayerTreeGroup)
if root.customProperty(KEY_SENSOR_GROUP) in [True, 'true']:
return True
for grp in root.findGroups():
if self._containsSensorNode(grp):
return True
return False
def onChildNodesRemoved(self, node, idxFrom, idxTo):
if not self._containsSensorNode(self.mLayerTreeModel.rootGroup()):
self._createSensorNode()
def onChangeCrosshairStyle(self):
canvases = self.mapCanvases()
if len(canvases) > 0:
mapCanvas = canvases[0]
else:
mapCanvas = None
style = getCrosshairStyle(parent=self, crosshairStyle=self.crosshairStyle(), mapCanvas=mapCanvas)
if isinstance(style, CrosshairStyle):
self.setCrosshairStyle(style)
def setVisibility(self, b: bool):
"""
Sets the map view visibility
:param b: bool
"""
assert isinstance(b, bool)
changed = False
if self.actionToggleMapViewHidden.isChecked() == b:
self.actionToggleMapViewHidden.setChecked(not b)
if changed:
self.sigCanvasAppearanceChanged.emit()
def isVisible(self)->bool:
"""
Returns the map view visibility
:return: bool
"""
return not self.actionToggleMapViewHidden.isChecked()
def mapCanvases(self)->list:
"""
Returns the MapCanvases related to this map view. Requires that this mapview was added to a MapWidget
:return: [list-of-MapCanvases]
"""
if isinstance(self.mMapWidget, MapWidget):
return self.mMapWidget.mapViewCanvases(self)
else:
return []
def onTitleChanged(self, *args):
self.setWindowTitle('Map View "{}"'.format(self.title()))
self.sigTitleChanged.emit(self.title())
if self.optionShowMapViewName.isChecked():
self.sigCanvasAppearanceChanged.emit()
def setMapWidget(self, w):
if isinstance(w, MapWidget):
self.mMapWidget = w
else:
self.mMapWidget = None
def setTimeSeries(self, timeSeries:TimeSeries):
"""
Conntects the MapView with a TimeSeries.
:param timeSeries: TimeSeries
"""
assert isinstance(timeSeries, TimeSeries)
for s in self.sensors():
self.removeSensor(s)
self.mTimeSeries = timeSeries
self.mTimeSeries.sigSensorAdded.connect(self.addSensor)
self.mTimeSeries.sigSensorRemoved.connect(self.removeSensor)
for s in timeSeries.sensors():
self.addSensor(s)
def timeSeries(self)->TimeSeries:
"""
Returns the TimeSeries this mapview is connected with
:return: TimeSeries
"""
return self.mTimeSeries
def setTitle(self, title:str):
"""
Sets the widget title
:param title: str
"""
old = self.title()
if old != title:
self.tbName.setText(title)
def layers(self)->list:
"""
Returns the visible layers, including proxy layer for time-series data
:return: [list-of-QgsMapLayers]
"""
return [l for l in self.mLayerTree.checkedLayers() if isinstance(l, QgsMapLayer)]
def title(self, maskNewLines=True)->str:
"""
Returns the MapView title
:return: str
"""
if maskNewLines:
return self.tbName.text().replace('\\n', ' ').strip()
else:
return self.tbName.text().strip()
def setCrosshairStyle(self, crosshairStyle:CrosshairStyle)->CrosshairStyle:
"""
Seths the CrosshairStyle of this MapView
:param crosshairStyle: CrosshairStyle
"""
if self.mCrossHairStyle != crosshairStyle:
self.mCrossHairStyle = crosshairStyle
self.sigCrosshairChanged.emit()
return self.mCrossHairStyle
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
"""
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 crosshairStyle(self)->CrosshairStyle:
"""
Returns the CrosshairStyle
:return: CrosshairStyle
"""
return self.mCrossHairStyle
def setCrosshairVisibility(self, b:bool):
"""
Enables / diables the map canvas crosshair.
:param b: bool
"""
if b != self.actionToggleCrosshairVisibility.isChecked():
self.actionToggleCrosshairVisibility.setChecked(b)
else:
self.mCrossHairStyle.setVisibility(b)
self.sigCrosshairChanged.emit()
def sensorProxyLayers(self)->list:
layers = [n.layer() for n in self.mLayerTreeSensorNode.findLayers()]
return [l for l in layers if isinstance(l, SensorProxyLayer)]
def sensorProxyLayer(self, sensor:SensorInstrument)->SensorProxyLayer:
"""
Returns the proxy layer related to a SensorInstrument
:param sensor: SensorInstrument
:return: SensorLayer
"""
for l in self.sensorProxyLayers():
if l.sensor() == sensor:
return l
return None
def sensors(self)->list:
"""
Returns a list of SensorsInstruments
:return: [list-of-SensorInstruments]
"""
return [t[0] for t in self.mSensorLayerList]
def addSensor(self, sensor:SensorInstrument):
"""
Adds a SensorInstrument to be shown in this MapView. Each sensor will be represented as a Raster Layer in the
Tree Model.
:param sensor: SensorInstrument
"""
assert isinstance(sensor, SensorInstrument)
if sensor not in self.sensors():
sensor.sigNameChanged.connect(self.sigCanvasAppearanceChanged)
dummyLayer = sensor.proxyLayer()
assert isinstance(dummyLayer.renderer(), QgsRasterRenderer)
dummyLayer.rendererChanged.connect(lambda sensor=sensor: self.onSensorRendererChanged(sensor))
#QgsProject.instance().addMapLayer(dummyLayer)
layerTreeLayer = self.mLayerTreeSensorNode.addLayer(dummyLayer)
assert isinstance(layerTreeLayer, QgsLayerTreeLayer)
layerTreeLayer.setCustomProperty(KEY_LOCKED_LAYER, True)
layerTreeLayer.setCustomProperty(KEY_SENSOR_LAYER, True)
self.mSensorLayerList.append((sensor, dummyLayer))
def onSensorRendererChanged(self, sensor:SensorInstrument):
for c in self.sensorCanvases(sensor):
assert isinstance(c, MapCanvas)
c.addToRefreshPipeLine(MapCanvas.Command.RefreshRenderer)
def sensorCanvases(self, sensor:SensorInstrument)->list:
"""
Returns the MapCanvases that show a layer with data for the given ``sensor``
:param sensor: SensorInstrument
:return:
"""
assert isinstance(sensor, SensorInstrument)
return [c for c in self.mapCanvases() if isinstance(c, MapCanvas) and \
isinstance(c.tsd(), TimeSeriesDate) and c.tsd().sensor() == sensor]
def sensorLayer(self, sensor: SensorInstrument):
"""
Returns the QgsRasterLayer that is used a proxy to specify the QgsRasterRenderer for a sensor
:param sensor: SensorInstrument
:return: QgsRasterLayer
"""
assert isinstance(sensor, SensorInstrument)
for t in self.mSensorLayerList:
s, l = t
assert isinstance(s, SensorInstrument)
assert isinstance(l, QgsRasterLayer)
if s == sensor:
return l
raise Exception('Sensor "{}" not registered to MapView "{}"'.format(sensor.name(), self.title()))
def removeSensor(self, sensor:SensorInstrument):
"""
Removes a sensor from this map view
:param sensor:
:return:
"""
toRemove = []
for t in self.mSensorLayerList:
if t[0] == sensor:
toRemove.append(t)
for t in toRemove:
self.mLayerTreeSensorNode.removeLayer(t[1])
self.mSensorLayerList.remove(t)
def hasSensor(self, sensor)->bool:
"""
:param sensor:
:return:
"""
assert isinstance(sensor, SensorInstrument)
return sensor in self.sensors()
class MapViewListModel(QAbstractListModel):
"""
A model to store a list of map views.
"""
sigMapViewsAdded = pyqtSignal(list)
sigMapViewsRemoved = pyqtSignal(list)
def __init__(self, parent=None):
super(MapViewListModel, self).__init__(parent)
self.mMapViewList = []
def addMapView(self, mapView):
i = len(self.mMapViewList)
self.insertMapView(i, mapView)
def insertMapView(self, i, mapView):
self.insertMapViews(i, [mapView])
def insertMapViews(self, i, mapViews):
assert isinstance(mapViews, list)
assert i >= 0 and i <= len(self.mMapViewList)
self.beginInsertRows(QModelIndex(), i, i + len(mapViews) - 1)
for j in range(len(mapViews)):
mapView = mapViews[j]
assert isinstance(mapView, MapView)
mapView.sigTitleChanged.connect(
lambda : self.doRefresh([mapView])
)
self.mMapViewList.insert(i + j, mapView)
self.endInsertRows()
self.sigMapViewsAdded.emit(mapViews)
def doRefresh(self, mapViews):
for mapView in mapViews:
idx = self.mapView2idx(mapView)
self.dataChanged.emit(idx, idx)
def removeMapView(self, mapView):
self.removeMapViews([mapView])
def removeMapViews(self, mapViews):
assert isinstance(mapViews, list)
for mv in mapViews:
assert mv in self.mMapViewList
idx = self.mapView2idx(mv)
self.beginRemoveRows(idx.parent(), idx.row(), idx.row())
self.mMapViewList.remove(mv)
self.endRemoveRows()
self.sigMapViewsRemoved.emit(mapViews)
def rowCount(self, parent=None, *args, **kwargs):
return len(self.mMapViewList)
def columnCount(self, QModelIndex_parent=None, *args, **kwargs):
return 1
def idx2MapView(self, index):
if isinstance(index, QModelIndex):
if index.isValid():
index = index.row()
else:
return None
assert index >= 0 and index < len(self.mMapViewList)
return self.mMapViewList[index]
def mapView2idx(self, mapView):
assert isinstance(mapView, MapView)
row = self.mMapViewList.index(mapView)
return self.createIndex(row, 0, mapView)
def __len__(self):
return len(self.mMapViewList)
def __iter__(self):
return iter(self.mMapViewList)
def __getitem__(self, slice):
return self.mMapViewList[slice]
def data(self, index, role=Qt.DisplayRole):
if not index.isValid():
return None
if (index.row() >= len(self.mMapViewList)) or (index.row() < 0):
return None
mapView = self.idx2MapView(index)
assert isinstance(mapView, MapView)
value = None
if role == Qt.DisplayRole:
value = '{} {}'.format(index.row() +1 , mapView.title())
#if role == Qt.DecorationRole:
#value = classInfo.icon(QSize(20,20))
if role == Qt.UserRole:
value = mapView
return value
class MapWidget(QFrame, loadUIFormClass(jp(DIR_UI, 'mapwidget.ui'))):
"""
This widget contains all maps
"""
class ViewMode(enum.Enum):
MapViewByRows = 1,
MapViewByCols = 2
sigSpatialExtentChanged = pyqtSignal(SpatialExtent)
sigCrosshairPositionChanged = pyqtSignal([SpatialPoint], [SpatialPoint, MapCanvas])
sigCRSChanged = pyqtSignal(QgsCoordinateReferenceSystem)
sigMapSizeChanged = pyqtSignal(QSize)
sigMapBackgroundColorChanged = pyqtSignal(QColor)
sigMapTextColorChanged = pyqtSignal(QColor)
sigMapTextFormatChanged = pyqtSignal(QgsTextFormat)
sigMapsPerMapViewChanged = pyqtSignal(int)
sigMapViewsChanged = pyqtSignal()
sigMapViewAdded = pyqtSignal(MapView)
sigMapViewRemoved = pyqtSignal(MapView)
sigCurrentDateChanged = pyqtSignal(TimeSeriesDate)
sigCurrentLocationChanged = pyqtSignal(SpatialPoint, MapCanvas)
sigVisibleDatesChanged = pyqtSignal(list)
sigViewModeChanged = pyqtSignal(ViewMode)
def __init__(self, *args, **kwds):
super(MapWidget, self).__init__(*args, **kwds)
self.setupUi(self)
self.setContentsMargins(1, 1, 1, 1)
self.mGrid = self.gridFrame.layout()
assert isinstance(self.mGrid, QGridLayout)
self.mGrid.setSpacing(0)
self.mGrid.setContentsMargins(0, 0, 0, 0)
self.mSyncLock = False
self.mMaxNumberOfCachedLayers = 0
self.mMapLayerStore = QgsMapLayerStore()
self.mMapLayerCache = dict()
self.mCanvasCache = dict()
self.mMapViews = []
self.mCanvases = dict()
self.mCanvasSignals = dict()
self.mTimeSeries = None
self.mMapToolKey = MapTools.Pan
self.mViewMode = MapWidget.ViewMode.MapViewByRows
self.mMpMV = 3
self.mSpatialExtent = SpatialExtent.world()
self.mCrs = self.mSpatialExtent.crs()
self.mCurrentDate = None
self.mCrosshairPosition = None
self.mMapSize = QSize(200, 200)
from eotimeseriesviewer.settings import DEFAULT_VALUES, Keys
self.mMapTextFormat = DEFAULT_VALUES[Keys.MapTextFormat]
self.mMapRefreshTimer = QTimer(self)
self.mMapRefreshTimer.timeout.connect(self.timedRefresh)
self.mMapRefreshTimer.setInterval(500)
self.mMapRefreshTimer.start()
self.btnFirst.setDefaultAction(self.actionFirstDate)
self.btnLast.setDefaultAction(self.actionLastDate)
self.btnBackward.setDefaultAction(self.actionBackward)
self.btnForward.setDefaultAction(self.actionForward)
self.btnBackwardFast.setDefaultAction(self.actionBackwardFast)
self.btnForwardFast.setDefaultAction(self.actionForwardFast)
self.actionFirstDate.triggered.connect(self.moveToFirstTSD)
self.actionLastDate.triggered.connect(self.moveToLastTSD)
self.actionBackward.triggered.connect(self.moveToPreviousTSD)
self.actionForward.triggered.connect(self.moveToNextTSD)
self.actionBackwardFast.triggered.connect(self.moveToPreviousTSDFast)
self.actionForwardFast.triggered.connect(self.moveToNextTSDFast)
self.mTimeSlider.setTickInterval(0)
self.mTimeSlider.valueChanged.connect(self.onSliderReleased)
def messageBar(self)->QgsMessageBar:
"""
Returns the QgsMessageBar
:return: QgsMessageBar
"""
return self.mMessageBar
def refresh(self):
for c in self.mapCanvases():
assert isinstance(c, MapCanvas)
c.timedRefresh()
def setMapTextFormat(self, textFormat:QgsTextFormat)->QgsTextFormat:
if not equalTextFormats(textFormat, self.mMapTextFormat):
self.mMapTextFormat = textFormat
for mapView in self.mapViews():
assert isinstance(mapView, MapView)
mapView.setMapTextFormat(textFormat)
self.sigMapTextFormatChanged.emit(self.mapTextFormat())
return self.mapTextFormat()
def mapTextFormat(self)->QgsTextFormat:
return self.mMapTextFormat
def setMapTool(self, mapToolKey:MapTools):
if self.mMapToolKey != mapToolKey:
self.mMapToolKey = mapToolKey
for c in self.mapCanvases():
assert isinstance(c, MapCanvas)
mts = c.mapTools()
mts.activate(self.mMapToolKey)
def visibleTSDs(self)->list:
"""
Returns the list of currently shown TimeSeriesDates.
:return: [list-of-TimeSeriesDates]
"""
for mv in self.mMapViews:
tsds = []
for c in self.mCanvases[mv]:
if isinstance(c.tsd(), TimeSeriesDate):
tsds.append(c.tsd())
return sorted(tsds)
return []
def spatialExtent(self)->SpatialExtent:
"""
Returns the current SpatialExtent
:return: SpatialExtent
"""
return self.mSpatialExtent
def setSpatialExtent(self, extent:SpatialExtent)->SpatialExtent:
"""
Sets a SpatialExtent to all MapCanvases.
:param extent: SpatialExtent
:return: SpatialExtent the current SpatialExtent
"""
assert isinstance(extent, SpatialExtent)
if self.mSpatialExtent != extent:
self.mSpatialExtent = extent
for c in self.mapCanvases():
assert isinstance(c, MapCanvas)
c.addToRefreshPipeLine(extent)
self.sigSpatialExtentChanged.emit(self.mSpatialExtent.__copy__())
return self.spatialExtent()
def setSpatialCenter(self, centerNew:SpatialPoint):
"""
Sets the spatial center of all MapCanvases
:param centerNew: SpatialPoint
"""
assert isinstance(centerNew, SpatialPoint)
extent = self.spatialExtent()
if isinstance(extent, SpatialExtent):
centerOld = extent.center()
centerNew = centerNew.toCrs(extent.crs())
if centerNew != centerOld and isinstance(centerNew, SpatialPoint):
extent = extent.__copy__()
extent.setCenter(centerNew)
self.setSpatialExtent(extent)
def spatialCenter(self)->SpatialPoint:
"""
Return the center of all map canvas
:return: SpatialPoint
"""
return self.spatialExtent().spatialCenter()
def setCrs(self, crs:QgsCoordinateReferenceSystem)->QgsCoordinateReferenceSystem:
"""
Sets the MapCanvas CRS.
:param crs: QgsCoordinateReferenceSystem
:return: QgsCoordinateReferenceSystem
"""
self.mCrs = crs
if isinstance(crs, QgsCoordinateReferenceSystem):
for c in self.mapCanvases():
c.setCrs(crs)
return self.crs()
def timedRefresh(self):
"""
Calls the timedRefresh() routine for all MapCanvases
"""
for c in self.mapCanvases():
assert isinstance(c, MapCanvas)
c.timedRefresh()
def currentLayers(self):
layers = set()
for c in self.mapCanvases():
layers = layers.union(set(c.layers()))
return list(layers)
def crs(self)->QgsCoordinateReferenceSystem:
return self.mCrs
def setTimeSeries(self, ts:TimeSeries)->TimeSeries:
assert ts == None or isinstance(ts, TimeSeries)
self.mTimeSeries = ts
if isinstance(self.mTimeSeries, TimeSeries):
self.mTimeSeries.sigVisibilityChanged.connect(self._updateCanvasDates)
self.mTimeSeries.sigTimeSeriesDatesRemoved.connect(self._updateCanvasDates)
self.mTimeSeries.sigTimeSeriesDatesAdded.connect(self._updateSliderRange)
self.mTimeSeries.sigTimeSeriesDatesRemoved.connect(self._updateSliderRange)
if len(self.mTimeSeries) > 0:
self.mCurrentDate = self.mTimeSeries[0]
else:
self.mTimeSeries.sigTimeSeriesDatesAdded.connect(self.onSetInitialCurrentDate)
self._updateSliderRange()
return self.timeSeries()
def onSetInitialCurrentDate(self):
if len(self.timeSeries()) > 0:
self.setCurrentDate(self.timeSeries()[0])
self.mTimeSeries.sigTimeSeriesDatesAdded.disconnect(self.onSetInitialCurrentDate)
def _updateSliderRange(self):
slider = self.timeSlider()
assert isinstance(slider, QSlider)
n = len(self.timeSeries())
slider.setRange(0, n)
slider.setEnabled(n > 0)
if n > 10:
pageStep = int(n/100)*10
slider.setTickInterval(pageStep)
else:
pageStep = 5
slider.setTickInterval(0)
slider.setPageStep(pageStep)
if n > 0:
tsd = self.currentDate()
if isinstance(tsd, TimeSeriesDate) and tsd in self.timeSeries():
i = self.timeSeries()[:].index(tsd)
slider.setValue(i+1)
def onSliderReleased(self):
i = self.mTimeSlider.value() - 1
if isinstance(self.mTimeSeries, TimeSeries) and len(self.mTimeSeries) > 0:
i = min(i, len(self.mTimeSeries)-1)
i = max(i, 0)
tsd = self.mTimeSeries[i]
self.setCurrentDate(tsd)
def timeSeries(self)->TimeSeries:
return self.mTimeSeries
def setMode(self, mode:ViewMode):
if mode != self.mViewMode:
self.mViewMode = mode
self._updateGrid()
self.sigViewModeChanged.emit(self.mViewMode)
def setMapsPerMapView(self, n:int):
assert n >= 0
if n != self.mMpMV:
self.mMpMV = n
self._updateGrid()
self.timeSlider().setPageStep(max(1, n))
self.sigMapsPerMapViewChanged.emit(n)
def setMapSize(self, size:QSize)->QSize:
"""
Sets the MapCanvas size
:param size: QSite
:return: QSize
"""
if size != self.mMapSize:
for canvas in self.mapCanvases():
canvas.setFixedSize(size)
self.mMapSize = size
self._updateWidgetSize()
self.sigMapSizeChanged.emit(size)
return self.mMapSize
def mapSize(self)->QSize:
"""
Returns the MapCanvas size
:return: QSize
"""
return self.mMapSize
def mapCanvases(self)->list:
"""
Returns all MapCanvases
:return: [list-of-MapCanvases]
"""
return self.findChildren(MapCanvas)
def mapViewCanvases(self, mapView:MapView)->list:
"""
Returns the MapCanvases related to a MapView
:param mapView: MapView
:return: [list-of-MapCanvases]
"""
return self.mCanvases[mapView]
def moveToNextTSD(self):
for tsd in self.timeSeries()[:]:
assert isinstance(tsd, TimeSeriesDate)
if tsd > self.currentDate() and tsd.isVisible():
self.setCurrentDate(tsd)
return
s = ""
def moveToPreviousTSD(self):
for tsd in reversed(self.timeSeries()[:]):
if tsd < self.currentDate() and tsd.isVisible():
self.setCurrentDate(tsd)
return
s = ""
def moveToNextTSDFast(self):
visible = list([tsd for tsd in self.timeSeries() if tsd.isVisible() and tsd > self.currentDate()])
if len(visible) > 0 and self.mMpMV > 0:
i = min(self.mMpMV-1, len(visible)-1)
self.setCurrentDate(visible[i])
def moveToPreviousTSDFast(self):
visible = list(reversed([tsd for tsd in self.timeSeries() if tsd.isVisible() and tsd < self.currentDate()]))
if len(visible) > 0 and self.mMpMV > 0:
i = min(self.mMpMV - 1, len(visible)-1)
self.setCurrentDate(visible[i])
def moveToFirstTSD(self):
for tsd in self.timeSeries()[:]:
if tsd.isVisible():
self.setCurrentDate(tsd)
return
s = ""
def moveToLastTSD(self):
for tsd in reversed(self.timeSeries()[:]):
if tsd.isVisible():
self.setCurrentDate(tsd)
return
s = ""
def setCurrentDate(self, tsd:TimeSeriesDate)->TimeSeriesDate:
"""
Sets the current TimeSeriesDate, i.e. the "center" date of all dates to be shown
:param tsd: TimeSeriesDate
:return: TimeSeriesDate
"""
assert isinstance(tsd, TimeSeriesDate)
b = tsd != self.mCurrentDate or (len(self.mapCanvases()) > 0 and self.mapCanvases()[0].tsd() is None)
self.mCurrentDate = tsd
if b:
self._updateCanvasDates()
i = self.mTimeSeries[:].index(self.mCurrentDate) + 1
if self.mTimeSlider.value() != i:
self.mTimeSlider.setValue(i)
self.sigCurrentDateChanged.emit(self.mCurrentDate)
if isinstance(self.currentDate(), TimeSeriesDate):
i = self.timeSeries()[:].index(self.currentDate())
canForward = i < len(self.mTimeSeries) - 1
canBackward = i > 0
else:
canForward = canBackward = False
for a in [self.actionForward, self.actionForwardFast, self.actionLastDate]:
a.setEnabled(canForward)
for a in [self.actionBackward, self.actionBackwardFast, self.actionFirstDate]:
a.setEnabled(canBackward)
return self.mCurrentDate
def timeSlider(self)->QSlider:
return self.mTimeSlider
def currentDate(self)->TimeSeriesDate:
"""
Returns the current TimeSeriesDate
:return: TimeSeriesDate
"""
return self.mCurrentDate
def addMapView(self, mapView:MapView)->MapView:
"""
Adds a MapView
:param mapView: MapView
:return: MapView
"""
assert isinstance(mapView, MapView)
if mapView not in self.mMapViews:
self.mMapViews.append(mapView)
mapView.setMapWidget(self)
# connect signals
mapView.sigCanvasAppearanceChanged.connect(self._updateCanvasAppearance)
mapView.sigCrosshairChanged.connect(self._updateCrosshair)
self._updateGrid()
self._updateCrosshair(mapView=mapView)
self.sigMapViewsChanged.emit()
self.sigMapViewAdded.emit(mapView)
return mapView
def removeMapView(self, mapView:MapView)->MapView:
"""
Removes a MapView
:param mapView: Mapview
:return: MapView
"""
if mapView in self.mMapViews:
self.mMapViews.remove(mapView)
mapView.setMapWidget(None)
# disconnect signals
self._updateGrid()
self.sigMapViewsChanged.emit()
self.sigMapViewRemoved.emit(mapView)
return mapView
def mapViews(self)->list:
"""
Returns a list of all MapViews
:return: [list-of-MapViews]
"""
return self.mMapViews[:]
def syncQGISCanvasCenter(self, qgisChanged:bool):
if self.mSyncLock:
return
iface = qgis.utils.iface
assert isinstance(iface, QgisInterface)
c = iface.mapCanvas()
if not isinstance(c, QgsMapCanvas):
return
tsvCenter = self.spatialExtent().spatialCenter()
qgsCenter = SpatialExtent.fromMapCanvas(c).spatialCenter()
if qgisChanged:
# change EOTSV
if tsvCenter.crs().isValid():
self.mSyncLock = True
qgsCenter = qgsCenter.toCrs(tsvCenter.crs())
if isinstance(qgsCenter, SpatialPoint):
self.setSpatialCenter(qgsCenter)
QApplication.processEvents()
self.mSyncLock = False
else:
# change QGIS
if qgsCenter.crs().isValid():
self.mSyncLock = True
tsvCenter = tsvCenter.toCrs(qgsCenter.crs())
if isinstance(tsvCenter, SpatialPoint):
c.setCenter(tsvCenter)
QApplication.processEvents()
self.mSyncLock = False
def _createMapCanvas(self)->MapCanvas:
mapCanvas = MapCanvas()
mapCanvas.setMapLayerStore(self.mMapLayerStore)
mapCanvas.mInfoItem.setTextFormat(self.mapTextFormat())
# set general canvas properties
mapCanvas.setFixedSize(self.mMapSize)
mapCanvas.setDestinationCrs(self.mCrs)
mapCanvas.setSpatialExtent(self.mSpatialExtent)
# activate the current map tool
mapTools = mapCanvas.mapTools()
mapTools.activate(self.mMapToolKey)
# connect signals
self._connectCanvasSignals(mapCanvas)
return mapCanvas
def _connectCanvasSignals(self, mapCanvas:MapCanvas):
mapCanvas.sigSpatialExtentChanged.connect(self.setSpatialExtent)
mapCanvas.sigDestinationCrsChanged.connect(self.setCrs)
mapCanvas.sigCrosshairPositionChanged.connect(self.onCrosshairPositionChanged)
mapCanvas.mapTools().mtCursorLocation.sigLocationRequest[SpatialPoint, QgsMapCanvas].connect(self.sigCurrentLocationChanged)
def _disconnectCanvasSignals(self, mapCanvas:MapCanvas):
mapCanvas.sigSpatialExtentChanged.disconnect(self.setSpatialExtent)
mapCanvas.sigDestinationCrsChanged.disconnect(self.setCrs)
mapCanvas.sigCrosshairPositionChanged.disconnect(self.onCrosshairPositionChanged)
mapCanvas.mapTools().mtCursorLocation.sigLocationRequest[SpatialPoint, QgsMapCanvas].disconnect(
self.sigCurrentLocationChanged)
def onCrosshairPositionChanged(self, spatialPoint:SpatialPoint):
canvas = self.sender()
if self.mCrosshairPosition != spatialPoint:
self.setCrosshairPosition(spatialPoint)
self.sigCrosshairPositionChanged[SpatialPoint, MapCanvas].emit(self.mCrosshairPosition, canvas)
def setCurrentLayer(self, layer:QgsMapLayer):
for mapView in self.mapViews():
mapView.setCurrentLayer(layer)
def setCrosshairPosition(self, spatialPoint)->SpatialPoint:
spatialPoint = spatialPoint.toCrs(self.crs())
if self.mCrosshairPosition != spatialPoint:
self.mCrosshairPosition = spatialPoint
for canvas in self.mapCanvases():
assert isinstance(canvas, MapCanvas)
canvas.setCrosshairPosition(spatialPoint)
self.sigCrosshairPositionChanged[SpatialPoint].emit(self.mCrosshairPosition)
return self.crosshairPosition()
def crosshairPosition(self)->SpatialPoint:
return self.mCrosshairPosition
def _updateGrid(self):
import time
t0 = time.time()
oldCanvases = self._updateLayerCache()
# crop grid
if self.mViewMode == MapWidget.ViewMode.MapViewByRows:
nc = self.mMpMV
nr = len(self.mapViews())
else:
raise NotImplementedError()
toRemove = []
for row in range(nr, self.mGrid.rowCount()):
for col in range(self.mGrid.columnCount()):
item = self.mGrid.itemAtPosition(row, col)
if isinstance(item, QLayoutItem) and isinstance(item.widget(), QWidget):
toRemove.append(item.widget())
for col in range(nc, self.mGrid.columnCount()):
for row in range(self.mGrid.rowCount()):
item = self.mGrid.itemAtPosition(row, col)
if isinstance(item, QLayoutItem) and isinstance(item.widget(), QWidget):
toRemove.append(item.widget())
for w in toRemove:
self.mGrid.removeWidget(w)
w.setParent(None)
w.setVisible(False)
usedCanvases = []
self.mCanvases.clear()
if self.mViewMode == MapWidget.ViewMode.MapViewByRows:
for row, mv in enumerate(self.mMapViews):
assert isinstance(mv, MapView)
self.mCanvases[mv] = []
for col in range(self.mMpMV):
item = self.mGrid.itemAtPosition(row, col)
if isinstance(item, QLayoutItem) and isinstance(item.widget(), MapCanvas):
c = item.widget()
else:
c = self._createMapCanvas()
self.mGrid.addWidget(c, row, col)
assert isinstance(c, MapCanvas)
#c.setFixedSize(self.mMapSize)
c.setTSD(None)
c.setMapView(mv)
usedCanvases.append(c)
self.mCanvases[mv].append(c)
else:
raise NotImplementedError()
t1 = time.time()
self._updateCanvasDates()
t2 = time.time()
self._updateWidgetSize()
t3 = time.time()
s = ""
# remove old canvases
for c in oldCanvases:
if c not in usedCanvases:
try:
self._disconnectCanvasSignals(c)
except:
pass
t4 = time.time()
if False:
print('t1: {}'.format(t1 - t0))
print('t2: {}'.format(t2 - t1))
print('t3: {}'.format(t3 - t2))
print('t4: {}'.format(t4 - t3))
def _updateWidgetSize(self):
self.mGrid.update()
#self.resize(self.sizeHint())
# self.setMaximumSize(self.sizeHint())
# self.setFixedSize(self.sizeHint())
#if False and self.parentWidget():
if True:
w = self
assert isinstance(w, QWidget)
rect = QGuiApplication.primaryScreen().geometry()
maxw, maxh = 0.66*rect.width(), 0.66*rect.height()
hint = self.sizeHint()
minw, minh = min(hint.width(), maxw), min(hint.height(), maxh)
w.setMinimumSize(minw, minh)
#w.setFixedSize(self.sizeHint())
w.layout().update()
w.update()
def _updateLayerCache(self)->list:
canvases = self.findChildren(MapCanvas)
for c in canvases:
assert isinstance(c, MapCanvas)
self.mMapLayerCache[self._layerListKey(c)] = c.layers()
return canvases
def _layerListKey(self, canvas:MapCanvas):
return (canvas.mapView(), canvas.tsd())
def _updateCanvasDates(self, updateLayerCache=True):
visibleBefore = self.visibleTSDs()
bTSDChanged = False
if updateLayerCache:
self._updateLayerCache()
if not (isinstance(self.mCurrentDate, TimeSeriesDate) and isinstance(self.timeSeries(), TimeSeries)):
for c in self.findChildren(MapCanvas):
assert isinstance(c, MapCanvas)
c.setTSD(None)
bTSDChanged = True
else:
visible = [tsd for tsd in self.timeSeries() if tsd.isVisible()]
t = self.mCurrentDate.date()
visible = sorted(visible, key=lambda tsd: abs(tsd.date() - t))
visible = visible[0:min(len(visible), self.mMpMV)]
visible = sorted(visible)
# set TSD of remaining canvases to None
while len(visible) < self.mMpMV:
visible.append(None)
for mapView in self.mapViews():
for tsd, canvas in zip(visible, self.mCanvases[mapView]):
assert isinstance(tsd, TimeSeriesDate) or tsd is None
assert isinstance(canvas, MapCanvas)
if canvas.tsd() != tsd:
canvas.setTSD(tsd)
key = self._layerListKey(canvas)
if key in self.mMapLayerCache.keys():
canvas.setLayers(self.mMapLayerCache.pop(key))
bTSDChanged = True
# canvas.setLayers()
if bTSDChanged:
self._updateCanvasAppearance()
visible2 = self.visibleTSDs()
if visible2 != visibleBefore:
self.sigVisibleDatesChanged.emit(visible2)
self._freeUnusedMapLayers()
def _freeUnusedMapLayers(self):
layers = [l for l in self.mMapLayerStore.mapLayers().values() if isinstance(l, SensorProxyLayer)]
needed = self.currentLayers()
toRemove = [l for l in layers if isinstance(l, SensorProxyLayer) and l not in needed]
# todo: use a kind of caching
# remove layers from MapLayerCache and MapLayerStore
for mv in self.mMapLayerCache.keys():
layers = [l for l in self.mMapLayerCache[mv] if l not in toRemove]
self.mMapLayerCache[mv] = layers
self.mMapLayerStore.removeMapLayers(toRemove)
def _updateCrosshair(self, mapView=None):
if isinstance(mapView, MapView):
mapViews = [mapView]
else:
mapViews = self.mapViews()
for mapView in mapViews:
assert isinstance(mapView, MapView)
style = mapView.crosshairStyle()
assert isinstance(style, CrosshairStyle)
for canvas in self.mCanvases[mapView]:
assert isinstance(canvas, MapCanvas)
item = canvas.mCrosshairItem
item.setVisibility(style.mShow)
assert isinstance(item, CrosshairMapCanvasItem)
item.setCrosshairStyle(style)
canvas.addToRefreshPipeLine(MapCanvas.Command.UpdateMapItems)
def _updateCanvasAppearance(self, mapView=None):
if isinstance(mapView, MapView):
mapViews = [mapView]
else:
mapViews = self.mapViews()
for mapView in mapViews:
assert isinstance(mapView, MapView)
v = mapView.isVisible()
bg = mapView.mapBackgroundColor()
tf = mapView.mapTextFormat()
showDate = mapView.optionShowDate.isChecked()
showName = mapView.optionShowMapViewName.isChecked()
showSensor = mapView.optionShowSensorName.isChecked()
for canvas in self.mCanvases[mapView]:
assert isinstance(canvas, MapCanvas)
# set overall visibility
if canvas.isVisible() != v:
canvas.setVisible(v)
tsd = canvas.tsd()
if canvas.canvasColor() != bg:
canvas.addToRefreshPipeLine(mapView.mapBackgroundColor())
# set info text
info = canvas.infoItem()
assert isinstance(info, MapCanvasInfoItem)
info.setTextFormat(tf)
uc = []
lc = []
if isinstance(tsd, TimeSeriesDate):
if showDate:
uc += ['{}'.format(tsd.date())]
if showName:
lc += ['{}'.format(mapView.title(maskNewLines=False))]
if showSensor:
uc += ['{}'.format(tsd.sensor().name())]
uc = '\n'.join(uc)
lc = '\n'.join(lc)
info.setUpperCenter(uc)
info.setLowerCenter(lc)
canvas.addToRefreshPipeLine(MapCanvas.Command.UpdateMapItems)
class MapViewDock(QgsDockWidget, loadUI('mapviewdock.ui')):
sigMapViewAdded = pyqtSignal(MapView)
sigMapViewRemoved = pyqtSignal(MapView)
sigShowProfiles = pyqtSignal(SpatialPoint, MapCanvas, str)
sigMapCanvasColorChanged = pyqtSignal(QColor)
sigMapCanvasTextFormatChanged = pyqtSignal(QgsTextFormat)
sigSpatialExtentChanged = pyqtSignal(SpatialExtent)
sigCrsChanged = pyqtSignal(QgsCoordinateReferenceSystem)
sigMapSizeChanged = pyqtSignal(QSize)
sigMapsPerMapViewChanged = pyqtSignal(int)
sigMapTextFormatChanged = pyqtSignal(QgsTextFormat)
def setTimeSeries(self, timeSeries:TimeSeries):
assert isinstance(timeSeries, TimeSeries)
self.mTimeSeries = timeSeries
self.mTimeSeries.sigSensorAdded.connect(self.addSensor)
self.mTimeSeries.sigSensorRemoved.connect(self.removeSensor)
def __init__(self, parent=None):
super(MapViewDock, self).__init__(parent)
self.setupUi(self)
self.baseTitle = self.windowTitle()
self.btnAddMapView.setDefaultAction(self.actionAddMapView)
self.btnRemoveMapView.setDefaultAction(self.actionRemoveMapView)
self.btnCrs.setOptionVisible(QgsProjectionSelectionWidget.LayerCrs, True)
self.btnCrs.setOptionVisible(QgsProjectionSelectionWidget.ProjectCrs, True)
self.btnCrs.setOptionVisible(QgsProjectionSelectionWidget.CurrentCrs, True)
self.btnCrs.setOptionVisible(QgsProjectionSelectionWidget.DefaultCrs, True)
self.btnCrs.setOptionVisible(QgsProjectionSelectionWidget.RecentCrs, True)
#self.btnCrs.setOptionVisible(QgsProjectionSelectionWidget.CrsNotSet, True)
self.btnCrs.crsChanged.connect(self.sigCrsChanged)
self.btnMapCanvasColor.colorChanged.connect(self.sigMapCanvasColorChanged)
self.btnTextFormat.changed.connect(lambda *args: self.sigMapTextFormatChanged.emit(self.mapTextFormat()))
self.btnApplySizeChanges.clicked.connect(self.onApplyButtonClicked)
self.actionAddMapView.triggered.connect(self.createMapView)
self.actionRemoveMapView.triggered.connect(lambda: self.removeMapView(self.currentMapView()) if self.currentMapView() else None)
self.toolBox.currentChanged.connect(self.onToolboxIndexChanged)
self.spinBoxMapSizeX.valueChanged.connect(lambda: self.onMapSizeChanged('X'))
self.spinBoxMapSizeY.valueChanged.connect(lambda: self.onMapSizeChanged('Y'))
self.mLastMapSize = self.mapSize()
self.mLastNDatesPerMapView = self.sbMpMV.value()
self.mTimeSeries = None
self.mMapWidget = None
def onApplyButtonClicked(self):
self.sigMapSizeChanged.emit(QSize(self.spinBoxMapSizeX.value(), self.spinBoxMapSizeY.value()))
self.sigMapsPerMapViewChanged.emit(self.mapsPerMapView())
def setMapWidget(self, mw)->MapWidget:
"""
Connects this MapViewDock with a MapWidget
:param mw: MapWidget
:return:
"""
assert isinstance(mw, MapWidget)
assert mw.timeSeries() == self.mTimeSeries, 'Set the time series first!'
self.mMapWidget = mw
self.sigCrsChanged.connect(mw.setCrs)
mw.sigCRSChanged.connect(self.setCrs)
self.sigMapSizeChanged.connect(mw.setMapSize)
mw.sigMapSizeChanged.connect(self.setMapSize)
self.sigMapTextFormatChanged.connect(mw.setMapTextFormat)
mw.sigMapTextFormatChanged.connect(self.setMapTextFormat)
self.sigMapsPerMapViewChanged.connect(mw.setMapsPerMapView)
mw.sigMapsPerMapViewChanged.connect(self.setMapsPerMapView)
self.sigMapViewAdded.connect(mw.addMapView)
self.sigMapViewRemoved.connect(mw.removeMapView)
mw.sigMapViewAdded.connect(self.addMapView)
mw.sigMapViewRemoved.connect(self.removeMapView)
for mapView in mw.mapViews():
self.addMapView(mapView)
return self.mMapWidget
def mapWidget(self)->MapWidget:
"""
Returns the connected MapWidget
:return: MapWidget
"""
return self.mMapWidget
def mapViews(self)->list:
"""
Returns the defined MapViews
:return: [list-of-MapViews]
"""
assert isinstance(self.toolBox, QToolBox)
mapViews = []
for i in range(self.toolBox.count()):
item = self.toolBox.widget(i)
if isinstance(item, MapView):
mapViews.append(item)
return mapViews
def mapCanvases(self)->list:
"""
Returns all MapCanvases from all MapViews
:return: [list-of-MapCanvases]
"""
maps = []
for mapView in self.mapViews():
assert isinstance(mapView, MapView)
maps.extend(mapView.mapCanvases())
return maps
def setCrs(self, crs):
if isinstance(crs, QgsCoordinateReferenceSystem):
old = self.btnCrs.crs()
if old != crs:
self.btnCrs.setCrs(crs)
self.btnCrs.setLayerCrs(crs)
def mapsPerMapView(self)->int:
return self.sbMpMV.value()
def setMapsPerMapView(self, n:int):
assert n >= 0
if self.sbMpMV.value != n:
self.sbMpMV.setValue(n)
self.mLastNDatesPerMapView = n
def setMapTextFormat(self, textFormat:QgsTextFormat):
if not equalTextFormats(textFormat, self.mapTextFormat()):
self.btnTextFormat.setTextFormat(textFormat)
def mapTextFormat(self)->QgsTextFormat:
return self.btnTextFormat.textFormat()
def setMapSize(self, size):
assert isinstance(size, QSize)
ws = [self.spinBoxMapSizeX, self.spinBoxMapSizeY]
oldSize = self.mapSize()
b = oldSize != size
for w in ws:
w.blockSignals(True)
self.spinBoxMapSizeX.setValue(size.width()),
self.spinBoxMapSizeY.setValue(size.height())
self.mLastMapSize = QSize(size)
for w in ws:
w.blockSignals(False)
self.mLastMapSize = QSize(size)
if b:
self.sigMapSizeChanged.emit(size)
def onMapSizeChanged(self, dim):
newSize = self.mapSize()
#1. set size of other dimension accordingly
if dim is not None:
if self.checkBoxKeepSubsetAspectRatio.isChecked():
if dim == 'X':
vOld = self.mLastMapSize.width()
vNew = newSize.width()
targetSpinBox = self.spinBoxMapSizeY
elif dim == 'Y':
vOld = self.mLastMapSize.height()
vNew = newSize.height()
targetSpinBox = self.spinBoxMapSizeX
oldState = targetSpinBox.blockSignals(True)
targetSpinBox.setValue(int(round(float(vNew) / vOld * targetSpinBox.value())))
targetSpinBox.blockSignals(oldState)
newSize = self.mapSize()
if newSize != self.mLastMapSize:
self.btnApplySizeChanges.setEnabled(True)
else:
self.sigMapSizeChanged.emit(self.mapSize())
self.btnApplySizeChanges.setEnabled(False)
self.setMapSize(newSize)
def mapSize(self)->QSize:
return QSize(self.spinBoxMapSizeX.value(),
self.spinBoxMapSizeY.value())
def dummySlot(self):
s =""
def onMapViewsRemoved(self, mapViews):
for mapView in mapViews:
idx = self.stackedWidget.indexOf(mapView.ui)
if idx >= 0:
self.stackedWidget.removeWidget(mapView.ui)
mapView.ui.close()
else:
s = ""
self.actionRemoveMapView.setEnabled(len(self.mMapViews) > 0)
def mapBackgroundColor(self)->QColor:
"""
Returns the map canvas background color
:return: QColor
"""
return self.btnMapCanvasColor.color()
def setMapBackgroundColor(self, color:QColor):
"""
Sets the MapCanvas background color
:param color: QColor
"""
if color != self.mapBackgroundColor():
self.btnMapCanvasColor.setColor(color)
def setMapTextColor(self, color:QColor):
"""
Sets the map text color
:param color: QColor
:return: QColor
"""
if color != self.mapTextColor():
self.btnMapTextColor.setColor(color)
return self.mapTextColor()
def mapTextColor(self)->QColor:
"""
Returns the map text color.
:return: QColor
"""
return self.btnMapTextColor.color()
def onMapViewsAdded(self, mapViews):
nextShown = None
for mapView in mapViews:
mapView.sigTitleChanged.connect(self.updateTitle)
self.stackedWidget.addWidget(mapView.ui)
if nextShown is None:
nextShown = mapView
contents = mapView.ui.scrollAreaWidgetContents
size = contents.size()
hint = contents.sizeHint()
#mapView.ui.scrollArea.update()
s = ""
#setMinimumSize(mapView.ui.scrollAreaWidgetContents.sizeHint())
#hint = contents.sizeHint()
#contents.setMinimumSize(hint)
if isinstance(nextShown, MapView):
self.setCurrentMapView(nextShown)
for mapView in mapViews:
self.sigMapViewAdded.emit(mapView)
def updateButtons(self, *args):
b = len(self.mMapViews) > 0
self.actionRemoveMapView.setEnabled(b)
self.actionApplyStyles.setEnabled(b)
self.actionHighlightMapView.setEnabled(b)
def createMapView(self, name:str=None)->MapView:
"""
Create a new MapView
:return: MapView
"""
mapView = MapView()
n = len(self.mapViews()) + 1
if isinstance(name, str) and len(name) > 0:
title = name
else:
title = 'Map View {}'.format(n)
while title in [m.title() for m in self.mapViews()]:
n += 1
title = 'Map View {}'.format(n)
if n == 1:
mapView.optionShowDate.setChecked(True)
mapView.optionShowSensorName.setChecked(True)
mapView.setTitle(title)
mapView.sigShowProfiles.connect(self.sigShowProfiles)
mapView.setTimeSeries(self.mTimeSeries)
self.addMapView(mapView)
return mapView
def onInfoOptionToggled(self):
self.sigMapInfoChanged.emit()
s = ""
def addMapView(self, mapView:MapView):
"""
Adds a MapView
:param mapView: MapView
"""
assert isinstance(mapView, MapView)
if mapView not in self:
mapView.sigTitleChanged.connect(lambda *args, mv=mapView: self.onMapViewUpdated(mv))
#mapView.sigVisibilityChanged.connect(lambda *args, mv=mapView: self.onMapViewUpdated(mv))
mapView.sigCanvasAppearanceChanged.connect(lambda *args, mv=mapView: self.onMapViewUpdated(mv))
self.sigMapCanvasColorChanged.connect(mapView.setMapBackgroundColor)
self.sigMapCanvasTextFormatChanged.connect(mapView.setMapTextFormat)
i = self.toolBox.addItem(mapView, mapView.windowIcon(), mapView.title())
self.toolBox.setCurrentIndex(i)
self.onMapViewUpdated(mapView)
if len(self.mapViews()) == 1:
self.setMapTextFormat(mapView.mapTextFormat())
self.sigMapViewAdded.emit(mapView)
def onToolboxIndexChanged(self):
b = isinstance(self.toolBox.currentWidget(), MapView)
self.actionRemoveMapView.setEnabled(b)
def onMapViewUpdated(self, mapView:MapView):
"""
Handles updates that react on MapView changes
:param mapView: MapView to make the update for
"""
numMV = 0
for i in range(self.toolBox.count()):
item = self.toolBox.widget(i)
if isinstance(item, MapView):
numMV += 1
if item == mapView:
if mapView.isVisible():
icon = QIcon(":/timeseriesviewer/icons/mapview.svg")
else:
icon = QIcon(":/timeseriesviewer/icons/mapviewHidden.svg")
self.toolBox.setItemIcon(i, icon)
self.toolBox.setItemText(i, 'Map View {} "{}"'.format(numMV, mapView.title()))
break
def removeMapView(self, mapView:MapView)->MapView:
"""
Removes a MapView
:param mapView: MapView
:return: MapView
"""
if mapView in self.mapViews():
for i in range(self.toolBox.count()):
w = self.toolBox.widget(i)
if isinstance(w, MapView) and w == mapView:
self.toolBox.removeItem(i)
mapView.close()
if self.toolBox.count() >= i:
self.toolBox.setCurrentIndex(min(i, self.toolBox.count()-1))
self.sigMapViewRemoved.emit(mapView)
return mapView
def __len__(self)->int:
"""
Returns the number of MapViews
:return: int
"""
return len(self.mapViews())
def __iter__(self):
"""
Provides an iterator over all MapViews
:return:
"""
return iter(self.mapViews())
def __getitem__(self, slice):
return self.mapViews()[slice]
def __contains__(self, mapView):
return mapView in self.mapViews()
def index(self, mapView):
assert isinstance(mapView, MapView)
return self.mapViews().index(mapView)
def addSensor(self, sensor:SensorInstrument):
"""
Adds a new SensorInstrument
:param sensor: SensorInstrument
"""
for mapView in self.mapViews():
mapView.addSensor(sensor)
def removeSensor(self, sensor:SensorInstrument):
"""
Removes a Sensor
:param sensor: SensorInstrument
"""
for mapView in self.mapViews():
assert isinstance(mapView, MapView)
mapView.removeSensor(sensor)
def applyStyles(self):
for mapView in self.mMapViews:
mapView.applyStyles()
def setCrosshairStyle(self, crosshairStyle):
for mapView in self.mMapViews:
mapView.setCrosshairStyle(crosshairStyle)
def setShowCrosshair(self, b):
for mapView in self.mMapViews:
mapView.setCrosshairVisibility(b)
def index(self, mapView):
assert isinstance(mapView, MapView)
return self.mapViewsDefinitions.index(mapView)
def setCurrentMapView(self, mapView):
assert isinstance(mapView, MapView) and mapView in self.mapViews()
self.toolBox.setCurrentWidget(mapView)
self.updateTitle()
def updateTitle(self, *args):
# self.btnToggleMapViewVisibility.setChecked(mapView)
mapView = self.currentMapView()
if isinstance(mapView, MapView):
title = '{} | {}'.format(self.baseTitle, mapView.title())
else:
title = self.baseTitle
self.setWindowTitle(title)
def currentMapView(self):
w = self.toolBox.currentWidget()
if isinstance(w, MapView):
return w
return None