Source code for taurus.qt.qtgui.plot.arrayedit

#!/usr/bin/env python

#############################################################################
##
## This file is part of Taurus
## 
## http://taurus-scada.org
##
## Copyright 2011 CELLS / ALBA Synchrotron, Bellaterra, Spain
## 
## Taurus is free software: you can redistribute it and/or modify
## it under the terms of the GNU Lesser General Public License as published by
## the Free Software Foundation, either version 3 of the License, or
## (at your option) any later version.
## 
## Taurus is distributed in the hope that it will be useful,
## but WITHOUT ANY WARRANTY; without even the implied warranty of
## MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
## GNU Lesser General Public License for more details.
## 
## You should have received a copy of the GNU Lesser General Public License
## along with Taurus.  If not, see <http://www.gnu.org/licenses/>.
##
#############################################################################

"""
arrayedit.py: Widget for editing a spectrum/array via control points 
"""


import numpy
from taurus.external.qt import Qt, Qwt5
from taurus.qt.qtgui.util.ui import UILoadable
from curvesAppearanceChooserDlg import CurveAppearanceProperties 


@UILoadable
class ControllerBox(Qt.QWidget):
    
    def __init__(self, parent=None, x=0, y=0, corr=0):
        Qt.QWidget.__init__(self, parent)
        self.loadUi()
        self._x = x
        self.setY(y)
        self.box.setTitle('x=%6g'%self._x)
        self.corrSB.setValue(corr)
        self.ctrlObj = self.corrSB.ctrlObj = self.lCopyBT.ctrlObj= self.rCopyBT.ctrlObj = self.lScaleBT.ctrlObj = self.rScaleBT.ctrlObj = self
        self.corrSB.focusInEvent = self.corrSB_focusInEvent #reimplementing the focusInEvent method for the spinbox
        self.box.mousePressEvent = self.mousePressEvent
        #self.connect(self.corrSB, Qt.SIGNAL('valueChanged(double)'), self.enableScale)
        
    def mousePressEvent(self, event):
        self.emit(Qt.SIGNAL('selected'),self._x)
        #print 'SELECTED', self
        #Qt.QDoubleSpinBox.focusInEvent(self.corrSB, event)
    
    def corrSB_focusInEvent(self, event):
        self.emit(Qt.SIGNAL('selected'),self._x)
        #print 'GOT FOCUS', self
        Qt.QDoubleSpinBox.focusInEvent(self.corrSB, event)
        
    def setY(self, y):
        self._y = y
        self.enableScale()
    
    def enableScale(self, *args):
        enable = (self._y + self.corrSB.value()) != 0
        self.lScaleBT.setEnabled(enable)
        self.rScaleBT.setEnabled(enable)


@UILoadable
class EditCPointsDialog(Qt.QDialog):
    
    def __init__(self, parent=None, x=0):
        Qt.QDialog.__init__(self, parent)
        self.setupUi(self)

@UILoadable
class AddCPointsDialog(Qt.QDialog):
    
    def __init__(self, parent=None, x=0):
        Qt.QDialog.__init__(self, parent)
        self.loadUi()
        

@UILoadable
[docs]class ArrayEditor(Qt.QWidget): def __init__(self, parent=None): Qt.QWidget.__init__(self, parent) self.loadUi() self._controllers = [] #construct the layout for controllers container self.ctrlLayout = Qt.QHBoxLayout(self.controllersContainer) self.ctrlLayout.setContentsMargins ( 5, 0, 5, 0 ) self.ctrlLayout.setSpacing(1) #implement scroll bars for the controllers container self.scrollArea = Qt.QScrollArea(self) self.scrollArea.setWidget(self.controllersContainer) self.scrollArea.setVerticalScrollBarPolicy(Qt.Qt.ScrollBarAlwaysOff) self.scrollArea.setWidgetResizable(True) self.cpointsGroupBox.layout().insertWidget(0, self.scrollArea) #initialize data cpoints = 2 self.x = numpy.arange(256, dtype='double') self.y = numpy.zeros(256, dtype='double') self.xp = numpy.linspace(self.x[0],self.x[-1], cpoints) self.corrp = numpy.zeros(cpoints) self.yp = numpy.interp(self.xp, self.x, self.y) self.corr = numpy.zeros(self.x.size) #markers self.markerPos=self.xp[0] self.marker1 = Qwt5.QwtPlotMarker() self.marker1.setSymbol(Qwt5.QwtSymbol(Qwt5.QwtSymbol.Rect, Qt.QBrush(Qt.Qt.NoBrush), Qt.QPen(Qt.Qt.green), Qt.QSize(8,8))) self.marker1.attach(self.plot1) self.marker2 = Qwt5.QwtPlotMarker() self.marker2.setSymbol(Qwt5.QwtSymbol(Qwt5.QwtSymbol.Rect, Qt.QBrush(Qt.Qt.NoBrush), Qt.QPen(Qt.Qt.green), Qt.QSize(8,8))) self.marker2.attach(self.plot2) #cpointsPickers self._cpointMovingIndex = None self._cpointsPicker1 = Qwt5.QwtPicker(self.plot1.canvas()) self._cpointsPicker1.setSelectionFlags(Qwt5.QwtPicker.PointSelection) self._cpointsPicker2 = Qwt5.QwtPicker(self.plot2.canvas()) self._cpointsPicker2.setSelectionFlags(Qwt5.QwtPicker.PointSelection ) self._cpointsPicker1.widgetMousePressEvent = self.plot1MousePressEvent self._cpointsPicker1.widgetMouseReleaseEvent = self.plot1MouseReleaseEvent self._cpointsPicker2.widgetMousePressEvent = self.plot2MousePressEvent self._cpointsPicker2.widgetMouseReleaseEvent = self.plot2MouseReleaseEvent self._cpointsPicker1.widgetMouseDoubleClickEvent = self.plot1MouseDoubleClickEvent self._cpointsPicker2.widgetMouseDoubleClickEvent = self.plot2MouseDoubleClickEvent self._populatePlots() self.resetCorrection() self._selectedController =self._controllers[0] self._addCPointsDialog = AddCPointsDialog(self) #Launch low-priority initializations (to speed up load time) #Qt.QTimer.singleShot(0, <method>) #connections self.connect(self.addCPointsBT, Qt.SIGNAL('clicked (bool)'), self._addCPointsDialog.show) self.connect(self._addCPointsDialog.editBT, Qt.SIGNAL('clicked (bool)'), self.showEditCPointsDialog) self.connect(self._addCPointsDialog.cleanBT, Qt.SIGNAL('clicked (bool)'), self.resetCorrection) self.connect(self._addCPointsDialog.addSingleCPointBT, Qt.SIGNAL('clicked (bool)'), self.onAddSingleCPointBT) self.connect(self._addCPointsDialog.addRegEspCPointsBT, Qt.SIGNAL('clicked (bool)'), self.onAddRegEspCPointsBT)
[docs] def plot1MousePressEvent(self, event): self.plotMousePressEvent(event, self.plot1)
[docs] def plot2MousePressEvent(self, event): self.plotMousePressEvent(event, self.plot2)
[docs] def plotMousePressEvent(self, event, taurusplot): picked, pickedCurveName, pickedIndex = taurusplot.pickDataPoint(event.pos(), scope=20, showMarker=False, targetCurveNames=['Control Points']) if picked is not None: self.changeCPointSelection(picked.x()) self.makeControllerVisible(self._controllers[pickedIndex]) self._cpointMovingIndex = pickedIndex self._movingStartYPos = event.y() taurusplot.canvas().setCursor(Qt.Qt.SizeVerCursor)
[docs] def plot1MouseReleaseEvent(self, event): self.plotMouseReleaseEvent(event, self.plot1)
[docs] def plot2MouseReleaseEvent(self, event): self.plotMouseReleaseEvent(event, self.plot2)
[docs] def plotMouseReleaseEvent(self, event, taurusplot): if self._cpointMovingIndex is None: return #if no cpoint was picked, do nothing on release #no motion s performed if the y position is unchanged or if the mouse release is out of the canvas validMotion = (self._movingStartYPos != event.pos().y()) and taurusplot.canvas().rect().contains(event.pos()) if validMotion: #calculate the new correction newCorr = taurusplot.invTransform(taurusplot.getCurve('Control Points').yAxis(), event.y()) if taurusplot is self.plot1: newCorr -= self.yp[self._cpointMovingIndex] #apply new correction self._controllers[self._cpointMovingIndex].corrSB.setValue(newCorr) #reset the moving state self._cpointMovingIndex = None taurusplot.canvas().setCursor(Qt.Qt.CrossCursor)
[docs] def plot1MouseDoubleClickEvent(self, event): self.plotMouseDoubleClickEvent(event, self.plot1)
[docs] def plot2MouseDoubleClickEvent(self, event): self.plotMouseDoubleClickEvent(event, self.plot2)
[docs] def plotMouseDoubleClickEvent(self, event, taurusplot): picked, pickedCurveName, pickedIndex = taurusplot.pickDataPoint(event.pos(), scope=20, showMarker=False, targetCurveNames=['Control Points']) if picked is not None: return #we dont want to create a control point too close of an existing one xp = taurusplot.invTransform(taurusplot.getCurve('Control Points').xAxis(), event.x()) if xp<self.xp[0] or xp>self.xp[-1]: return #we dont want to create control points out of the curve range if Qt.QMessageBox.question(self, 'Create Control Point?', 'Insert a new control point at x=%g?'%xp, 'Yes', 'No') == 0: self.insertController(xp) self.changeCPointSelection(xp) Qt.QTimer.singleShot(1, self.makeControllerVisible) #singleshot is used as a hack to get out of the eventhandler
[docs] def makeControllerVisible(self, ctrl = None): if ctrl is None: ctrl = self._selectedController self.scrollArea.ensureWidgetVisible(ctrl)
[docs] def connectToController(self, ctrl): self.connect(ctrl, Qt.SIGNAL('selected'), self.changeCPointSelection) self.connect(ctrl.corrSB, Qt.SIGNAL('valueChanged (double)'), self.onCorrSBChanged) self.connect(ctrl.lCopyBT, Qt.SIGNAL('clicked (bool)'), self.onLCopy) self.connect(ctrl.rCopyBT, Qt.SIGNAL('clicked (bool)'), self.onRCopy) self.connect(ctrl.lScaleBT, Qt.SIGNAL('clicked (bool)'), self.onLScale) self.connect(ctrl.rScaleBT, Qt.SIGNAL('clicked (bool)'), self.onRScale)
[docs] def onAddSingleCPointBT(self): x=self._addCPointsDialog.singleCPointXSB.value() self.insertController(x)
[docs] def onAddRegEspCPointsBT(self): cpoints = self._addCPointsDialog.regEspCPointsSB.value() positions = numpy.linspace(self.x[0],self.x[-1], cpoints+2)[1:-1] for xp in positions: self.insertController(xp)
[docs] def showEditCPointsDialog(self): dialog = EditCPointsDialog(self) table = dialog.tableTW table.setRowCount(self.xp.size) for i,(xp,corrp) in enumerate(zip(self.xp, self.corrp)): table.setItem(i,0, Qt.QTableWidgetItem(str(xp))) table.setItem(i,1, Qt.QTableWidgetItem(str(corrp))) #show dialog and update values if it is accepted if dialog.exec_(): #delete previous controllers for c in self._controllers: c.setParent(None) self._controllers = [] #and create them anew new_xp = numpy.zeros(table.rowCount()) new_corrp = numpy.zeros(table.rowCount()) try: for i in xrange(table.rowCount()): new_xp[i] = float(table.item(i,0).text()) new_corrp[i] = float(table.item(i,1).text()) self.setCorrection(new_xp, new_corrp) except: Qt.QMessageBox.warning(self, 'Invalid data', 'Some values were not valid. Edition is ignored.')
def _getInsertionIndex(self, xp): i=0 while (self.xp[i]<xp):i+=1 return i
[docs] def insertControllers(self, xplist): for xp in xplist: self.insertController(xp)
[docs] def insertController(self, xp, index = None): if xp in self.xp: return None if index is None: index = self._getInsertionIndex(xp) #updating data (not in the most efficient way, but at least is clean) old_xp = self.xp self.xp = numpy.concatenate((self.xp[:index],(xp,), self.xp[index:])) self.yp = numpy.interp(self.xp, self.x, self.y) self.corrp = numpy.interp(self.xp, old_xp, self.corrp) #the new correction is obtained by interpolation from the neighbouring ones #creating the controller ctrl = ControllerBox(parent=None, x=xp, y=self.yp[index], corr=self.corrp[index]) self.ctrlLayout.insertWidget(index, ctrl) self._controllers.insert(index, ctrl) self.connectToController(ctrl) self.updatePlots() return index
[docs] def delController(self, index): self._controllers.pop(index).setParent(None) self.xp = numpy.concatenate((self.xp[:index], self.xp[index+1:])) self.yp = numpy.interp(self.xp, self.x, self.y) self.corrp = numpy.concatenate((self.corrp[:index], self.corrp[index+1:]))
def _populatePlots(self): #Curves appearance self._yAppearance= CurveAppearanceProperties( sStyle=Qwt5.QwtSymbol.NoSymbol, lStyle=Qt.Qt.SolidLine, lWidth=2, lColor='black', cStyle=Qwt5.QwtPlotCurve.Lines, yAxis=Qwt5.QwtPlot.yLeft) self._correctedAppearance= CurveAppearanceProperties( sStyle=Qwt5.QwtSymbol.NoSymbol, lStyle=Qt.Qt.DashLine, lWidth=1, lColor='blue', cStyle=Qwt5.QwtPlotCurve.Lines, yAxis=Qwt5.QwtPlot.yLeft) self._cpointsAppearance= CurveAppearanceProperties( sStyle=Qwt5.QwtSymbol.Rect, sSize=5, sColor='blue', sFill=True, lStyle=Qt.Qt.NoPen, yAxis=Qwt5.QwtPlot.yLeft) self._corrAppearance= CurveAppearanceProperties( sStyle=Qwt5.QwtSymbol.NoSymbol, lStyle=Qt.Qt.SolidLine, lWidth=1, lColor='blue', cStyle=Qwt5.QwtPlotCurve.Lines, yAxis=Qwt5.QwtPlot.yLeft) self.plot1.attachRawData({'x':self.x, 'y':self.y, 'title':'Master'}) self.plot1.setCurveAppearanceProperties({'Master':self._yAppearance}) self.plot1.attachRawData({'x':self.xp, 'y':self.yp+self.corrp, 'title':'Control Points'}) self.plot1.setCurveAppearanceProperties({'Control Points':self._cpointsAppearance}) self.plot1.attachRawData({'x':self.x, 'y':self.y+self.corr, 'title':'Corrected'}) self.plot1.setCurveAppearanceProperties({'Corrected':self._correctedAppearance}) self.plot2.attachRawData({'x':self.x, 'y':self.corr, 'title':'Correction'}) self.plot2.setCurveAppearanceProperties({'Correction':self._corrAppearance}) self.plot2.attachRawData({'x':self.xp, 'y':self.corrp, 'title':'Control Points'}) self.plot2.setCurveAppearanceProperties({'Control Points':self._cpointsAppearance})
[docs] def updatePlots(self): self.plot1.getCurve('Control Points').setData(self.xp, self.yp + self.corrp) self.plot1.getCurve('Corrected').setData(self.x, self.y + self.corr) self.plot2.getCurve('Correction').setData(self.x, self.corr) self.plot2.getCurve('Control Points').setData(self.xp, self.corrp) index = self._getInsertionIndex(self.markerPos) self.marker1.setValue(self.xp[index],self.yp[index]+self.corrp[index]) self.marker2.setValue(self.xp[index],self.corrp[index]) self.plot1.replot() self.plot2.replot()
[docs] def onLCopy(self, checked): sender = self.sender().ctrlObj index = self._getInsertionIndex(sender._x) v=sender.corrSB.value() for ctrl in self._controllers[:index]: ctrl.corrSB.setValue(v)
[docs] def onRCopy(self, checked): sender=self.sender().ctrlObj index = self._getInsertionIndex(sender._x) v=sender.corrSB.value() for ctrl in self._controllers[index+1:]: ctrl.corrSB.setValue(v)
[docs] def onLScale(self, checked): sender=self.sender().ctrlObj index = self._getInsertionIndex(sender._x) #y=numpy.interp(self.xp, self.x, self.y) #values of the master at the control points if self.yp[index] == 0: Qt.QMessageBox.warning(self, 'Scaling Error', 'The master at this control point is zero-valued. This point cannot be used as reference for scaling') return v=sender.corrSB.value()/(self.yp[index]) for i in range(0,index): self._controllers[i].corrSB.setValue(v*self.yp[i])
[docs] def onRScale(self, checked): sender=self.sender().ctrlObj index = self._getInsertionIndex(sender._x) #y=numpy.interp(self.xp, self.x, self.y) #values of the master at the control points if self.yp[index] == 0: Qt.QMessageBox.warning(self, 'Scaling Error', 'The master at this control point is zero-valued. This point cannot be used as reference for scaling') return v=sender.corrSB.value()/(self.yp[index]) for i in range(index+1,self.xp.size): self._controllers[i].corrSB.setValue(v*self.yp[i])
[docs] def changeCPointSelection(self, newpos): index = self._getInsertionIndex(newpos) old_index = self._getInsertionIndex(self.markerPos) self.markerPos = self.xp[index] self.marker1.setValue(self.xp[index],self.yp[index]+self.corrp[index]) self.marker2.setValue(self.xp[index],self.corrp[index]) self.plot1.replot() self.plot2.replot() self._controllers[old_index].corrSB.setStyleSheet('') self._controllers[index].corrSB.setStyleSheet('background:lightgreen') self._selectedController = self._controllers[index]
[docs] def onCorrSBChanged(self, value=None): '''recalculates the position and value of the control points (self.xp and self.corrp) as well as the correction curve (self.corr)''' ctrl = self.sender().ctrlObj self.corrp[self._getInsertionIndex(ctrl._x)] = value self.corr = numpy.interp(self.x, self.xp, self.corrp) # recalculate the correction curve self.updatePlots()
[docs] def setMaster(self, x, y, keepCP=False, keepCorr=False): #make sure that x,y are numpy arrays and that the values are sorted for x x,y = numpy.array(x), numpy.array(y) if x.shape != y.shape or x.size == 0 or y.size == 0: raise ValueError('The master curve is not valid' ) sortedindexes = numpy.argsort(x) self.x, self.y = x[sortedindexes], y[sortedindexes] self.plot1.getCurve('Master').setData(self.x, self.y) xp = None corrp = None if self.x[0] == self.xp[0] and self.x[-1] == self.x[-1]: if keepCP: xp = self.xp if keepCorr: corrp = self.corrp self.setCorrection(xp, corrp) self._addCPointsDialog.singleCPointXSB.setRange(self.x[0],self.x[-1])
[docs] def getMaster(self): '''returns x,m where x and m are numpy arrays representing the abcissas and ordinates for the master, respectively''' return self.x.copy(),self.y.copy()
[docs] def resetMaster(self): x = numpy.arange(256, dtype='double') y = numpy.zeros(256, dtype='double') self.setMaster(x,y)
[docs] def getCorrected(self): '''returns x,c where x and c are numpy arrays representing the abcissas and ordinates for the corrected curve, respectively''' return self.x.copy(), self.y+self.corr
[docs] def getCorrection(self): '''returns xp,cp where xp and cp are numpy arrays representing the abcissas and ordinates for the correction points, respectively''' return self.xp.copy(), self.corrp.copy()
[docs] def setCorrection(self, xp=None, corrp=None): '''sets control points at the points specified by xp and with the values specified by corrp. Example:: setCorrection([1,2,8,9], [0,0,0,0]) would set 4 control points with initial value 0 at x=1, 2, 8 and 9s ''' for c in self._controllers: c.setParent(None) #destroy previous controllers self._controllers = [] if xp is None: xp = numpy.array((self.x[0], self.x[-1])) corrp = numpy.zeros(2) if corrp is None: corrp = numpy.zeros(xp.size) #make sure that the extremes are there if xp[0] > self.x[0]: xp = numpy.concatenate(((self.x[0],), xp)) corrp = numpy.concatenate(((self.corrp[0],), corrp)) if xp[-1] < self.x[-1]: xp = numpy.concatenate((xp, (self.x[-1],))) corrp = numpy.concatenate((corrp, (self.corrp[-1],))) #now create everything self.xp = numpy.unique(xp) #make sure there are no repetitions and that it is sorted self.corrp = numpy.interp(self.xp, xp, corrp) #in case of repeated xp, take only one corrp self.yp = numpy.interp(self.xp, self.x, self.y) for i,(x,c) in enumerate(zip(xp,corrp)): ctrl = ControllerBox(parent=None, x=xp[i], y=self.yp[i]) self.ctrlLayout.insertWidget(i, ctrl) self._controllers.insert(i, ctrl) self.connectToController(ctrl) self.corr = numpy.interp(self.x, self.xp, self.corrp) # recalculate the correction curve self.updatePlots() self.changeCPointSelection(self.xp[0])
[docs] def resetCorrection(self): self.setCorrection()
#self.xp = numpy.linspace(self.x[0],self.x[-1], self.cpoints) #self.corrp = numpy.zeros(self.cpoints) #self._populateControllers() if __name__ == "__main__": import sys app = Qt.QApplication(sys.argv) form = ArrayEditor() #x = numpy.arange(100)-20 #y = -(x-50)**2+50**2 x = numpy.linspace(0.1,0.9,100) y = (x)**2-5*x form.setMaster(x,y) #form.setCorrection([1,30,70,90,100],[0,44,88,22,-100]) form.show() #sys.exit(app.exec_()) status=app.exec_() #x,y=form.getCorrected() #print "CORRECTED:",x,y sys.exit(status)