Newer
Older
self.actionSet453.triggered.connect(lambda: self.setBandSelection('nIR,swIR,R'))
def addBtnAction(action):
btn = QToolButton()
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)
self.ui.show()
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)
for sensor in self.mapViewCollection.TS.Sensors:
self.addSensor(sensor)
self.setTitle(name)
#forward actions with reference to this band view
def dummy(self, *args):
print(args)
def setIsVisible(self, b):
assert isinstance(b, bool)
changed = False
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)
def isVisible(self):
return not self.ui.actionToggleMapViewHidden.isChecked()
def mapCanvases(self):
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
def vectorLayerRenderer(self):
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)

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

Benjamin Jakimow
committed
if isinstance(lyr, QgsVectorLayer) and self.ui.gbVectorRendering.isChecked():
#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.layerModel().setVectorLayerSources([self.mVectorLayer])
#mapCanvas.setLayers([l for l in mapCanvas.layers() if isinstance(l, QgsRasterLayer)])
#mapCanvas.setLazyVectorSources([lyr])
mapCanvas.refresh()

benjamin.jakimow@geo.hu-berlin.de
committed
else:
#remove vector layers
self.mVectorLayer = None
for mapCanvas in self.mapCanvases():
mapCanvas.layerModel().setVectorLayerSources([])
#mapCanvas.setLayers([l for l in mapCanvas.mLayers if not isinstance(l, QgsVectorLayer)])
mapCanvas.refresh()

benjamin.jakimow@geo.hu-berlin.de
committed
def applyStyles(self):
for sensorView in self.mSensorViews.values():

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

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

benjamin.jakimow@geo.hu-berlin.de
committed
def title(self):

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()

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

Benjamin Jakimow
committed
self.onCrosshairChanged(crosshairStyle)

benjamin.jakimow@geo.hu-berlin.de
committed
def setHighlighted(self, b=True, timeout=1000):
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):
return self.ui.actionToggleRasterVisibility.isChecked()
def vectorVisibility(self):
return self.ui.actionToggleVectorVisibility.isChecked()
def setRasterVisibility(self, b):

Benjamin Jakimow
committed
self.mRastersVisible = b
self.ui.actionToggleRasterVisibility.setChecked(b)
for mapCanvas in self.mapCanvases():
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.layerModel().setRasterLayerVisibility(b)
mapCanvas.refresh()
def setVectorVisibility(self, b):
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.layerModel().setVectorLayerVisibility(self.mVectorsVisible)
mapCanvas.refresh()

benjamin.jakimow@geo.hu-berlin.de
committed
def removeSensor(self, sensor):
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, mapCanvas):
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
mapCanvas.setRenderer(mapViewRenderSettings.rasterRenderer())
mapCanvas.setRenderer(self.vectorLayerRenderer())

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
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
1361
1362
1363
1364
1365
1366
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."""
1383
1384
1385
1386
1387
1388
1389
1390
1391
1392
1393
1394
1395
1396
1397
1398
1399
1400
1401
1402
1403
1404
# 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)
self.mMockupCanvas = QgsMapCanvas(parent=parent)
self.mMockupRasterLayer = QgsRasterLayer(self.mSensor.pathImg)
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())
#w = func(QgsRasterLayer(), QgsRectangle())
w.setSizePolicy(QSizePolicy(QSizePolicy.Preferred, QSizePolicy.Preferred))
#w.sizeHint = lambda : QSize(300, 50)
w.setMapCanvas(self.mMockupCanvas)
self.stackedWidget.addWidget(w)
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()
def onSensorNameChanged(self, newName):
self.setTitle(self.mSensor.name())
self.actionApplyStyle.setToolTip('Apply style to all map view images from "{}"'.format(self.mSensor.name()))
1462
1463
1464
1465
1466
1467
1468
1469
1470
1471
1472
1473
1474
1475
1476
1477
1478
1479
1480
1481
1482
1483
1484
1485
1486
1487
1488
1489
1490
1491
1492
1493
1494
1495
1496
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)
1549
1550
1551
1552
1553
1554
1555
1556
1557
1558
1559
1560
1561
1562
1563
1564
1565
1566
1567
1568
1569
1570
for mapCanvas in self.mMapCanvases:
assert isinstance(mapCanvas, MapCanvas)
mapCanvas.layerModel().setRenderer(r)
mapCanvas.refresh()
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):
1582
1583
1584
1585
1586
1587
1588
1589
1590
1591
1592
1593
1594
1595
1596
1597
1598
1599
1600
1601
1602
1603
1604
1605
1606
1607
#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()
1630
1631
1632
1633
1634
1635
1636
1637
1638
1639
1640
1641
1642
1643
1644
1645
1646
1647
1648
1649
1650
1651
1652
1653
1654
1655
1656
1657
1658
1659
1660
1661
1662
1663
1664
1665
1666
1667
1668
1669
1670
1671
1672
1673
1674
1675
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')):
1679
1680
1681
1682
1683
1684
1685
1686
1687
1688
1689
1690
1691
1692
1693
1694
1695
1696
1697
1698
1699
1700
1701
1702
1703
1704
1705
1706
"""
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())
1719
1720
1721
1722
1723
1724
1725
1726
1727
1728
1729
1730
1731
1732
1733
1734
1735
1736
1737
1738
1739
1740
1741
1742
1743
1744
1745
1746
1747
1748
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)
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 = self.TSD.sensor
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.date), self.TSD.sensor.name())
tooltip = '{}'.format(self.TSD.pathImg)

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()

benjamin.jakimow@geo.hu-berlin.de
committed
def refresh(self):

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

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.date))

Benjamin Jakimow
committed
self.registerMapCanvas(mapView, mapCanvas)

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

Benjamin Jakimow
committed
mapView.registerMapCanvas(self.mSensor, mapCanvas)
# register MapCanvas on STV level
self.STV.registerMapCanvas(mapCanvas)
mapCanvas.renderComplete.connect(lambda : self.onRenderingChange(False))
mapCanvas.renderStarting.connect(lambda : self.onRenderingChange(True))

benjamin.jakimow@geo.hu-berlin.de
committed
mapCanvas.sigDataLoadingFinished.connect(
lambda dt: self.STV.TSV.ui.dockSystemInfo.addTimeDelta('Map {}'.format(self.mSensor.name()), dt))

benjamin.jakimow@geo.hu-berlin.de
committed
mapCanvas.sigDataLoadingFinished.connect(
lambda dt: self.STV.TSV.ui.dockSystemInfo.addTimeDelta('All Sensors', dt))
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.layerModel().setRasterLayerSources([self.TSD.pathImg])
#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)
if key == 'copy_sensor':
QApplication.clipboard().setText(self.TSD.sensor.name())
if key == 'copy_date':
QApplication.clipboard().setText(str(self.TSD.date))
if key == 'copy_path':
QApplication.clipboard().setText(str(self.TSD.pathImg))

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()
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.mRefreshTimer = QTimer(self)
self.mRefreshTimer.setInterval(1000)
self.mRefreshTimer.timeout.connect(self.refresh)
self.scrollArea.sigResized.connect(self.mRefreshTimer.start)
self.scrollArea.horizontalScrollBar().valueChanged.connect(self.mRefreshTimer.start)

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

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:
self.setSpatialExtent(self.TS.getMaxSpatialExtent())

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

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