__main__.py 15.9 KB
Newer Older
Yan's avatar
Yan committed
1
2
3
4
5
#!/usr/bin/env python3

from matplotlib.backends.backend_qt5agg import FigureCanvasQTAgg as FigureCanvas
from matplotlib.figure import Figure
from matplotlib.widgets import SpanSelector
Yan's avatar
Yan committed
6
from matplotlib.gridspec import GridSpec
Yan's avatar
Yan committed
7
8
from PyQt5 import QtCore
from PyQt5 import QtWidgets
Yan's avatar
Yan committed
9
from PyQt5 import QtGui
10
from PyQt5 import QtPrintSupport
Yan's avatar
Yan committed
11
from rawprasslib import load_raw
Yan's avatar
Yan committed
12
from rawprasslib import rawprasslib
Yan's avatar
Yan committed
13
from prasopes.predictmz import predict as getmzpattern
Yan's avatar
Yan committed
14
15
16
17
18
try:
    from rawautoparams import load_params
    autoparams = True
except:
    autoparams = False
Yan's avatar
Yan committed
19
import numpy as np
Yan's avatar
Yan committed
20
import prasopes.config as cf
21
import prasopes.datatools as dt
Yan's avatar
Yan committed
22
import prasopes.drltools as drl
Yan's avatar
Yan committed
23
import prasopes.filetools as ft
Yan's avatar
Yan committed
24
25
import prasopes.graphtools as gt
import prasopes.imagetools as imgt
Yan's avatar
Yan committed
26
import prasopes.zcetools as zce
Yan's avatar
Yan committed
27
import prasopes.tangoicons
Yan's avatar
Yan committed
28
29
30
import sys
import matplotlib
import logging
3Yan's avatar
3Yan committed
31
import os.path
Yan's avatar
Yan committed
32
33
34
matplotlib.use("Qt5Agg")


35
36
37
38
class update_signal(QtCore.QObject):
    signal = QtCore.pyqtSignal()


39
def load_file(parent, chrom_plot, spc, d_set, ms_s, fn, chrom_s,
40
              update, settings, filename=None):
Yan's avatar
Yan committed
41
    """populates dataset and plots it"""
3Yan's avatar
3Yan committed
42
43
44
    directory=settings.value("open_folder")
    if fn[0] is not None:
        directory=fn[0]
Yan's avatar
Yan committed
45
46
    if filename == None:
        filename = QtWidgets.QFileDialog.getOpenFileName(
3Yan's avatar
3Yan committed
47
48
49
            caption="Open spectrum",
            directory=directory,
            filter="Finnigan RAW files (*.raw, *.RAW)")[0]
3Yan's avatar
3Yan committed
50
    if filename is not '' and os.path.isfile(filename):
Yan's avatar
Yan committed
51
        try:
3Yan's avatar
3Yan committed
52
            d_set.clear()
Yan's avatar
Yan committed
53
            chrom_s['timesarg'].clear()
3Yan's avatar
3Yan committed
54
55
            [d_set.append(dict(chrom_dat=i[0], masses=i[1], matrix=i[2])) for i
                    in load_raw(filename, settings.value("tmp_location"))]
Yan's avatar
Yan committed
56
57
58
59
60
61
62
        except rawprasslib.ParsingException as pex:
            QtWidgets.QMessageBox.critical(
                parent, "Opening of the file has failed!",
                "File is incompatible with the rawprasslib, "
                "canceling request!\n\n"
                "Error message:\n{}".format(pex.args[0]))
            return
Yan's avatar
Yan committed
63
64
65
66
67
68
69
70
71
72
73
74
75
        [i.clear() for i in (ms_s['params'],ms_s['headers'])]
        if autoparams == True:
            try:
                ms_s['params'], rawheaders = load_params(filename)
                segments = [len(subset['chrom_dat'][0]) for subset in d_set]
                indicies = [sum(segments[:i+1]) for i in range(len(segments))]
                ms_s['headers'] = np.split(rawheaders, indicies)[:-1]
            except Exception as pex:
                QtWidgets.QMessageBox.critical(
                    parent, "Parameters readout has failed!",
                    "File is incompatible with the rawautoparams, "
                    "no parameters loaded!\n\n"
                    "Error message:\n{}".format(pex.args[0]))
Yan's avatar
Yan committed
76
        gt.populate(chrom_plot, spc, d_set, ms_s, chrom_s)
3Yan's avatar
3Yan committed
77
        fn[0] = filename
78
79
        parent.setWindowTitle("Prasopes - {}".format(
            os.path.basename(filename)))
80
        update.signal.emit()
3Yan's avatar
3Yan committed
81

Yan's avatar
Yan committed
82

Yan's avatar
Yan committed
83
84
def print_graph(data_set, mass_spec, chrom_spec, spect, fn):
    def printimage(printdevice, img):
Yan's avatar
Yan committed
85
        printer.setResolution(600)
Yan's avatar
Yan committed
86
        painter = QtGui.QPainter(printdevice)
Yan's avatar
Yan committed
87
        painter.drawImage(0,0,img)
Yan's avatar
Yan committed
88
89
        painter.end()
    #TODO: substitute the QPrintPreviewDialog with QPrintPreviewWidget
90
91
    printPreview = QtPrintSupport.QPrintPreviewDialog()
    printer = printPreview.printer()
Yan's avatar
Yan committed
92
93
    printer.setPageSize(printer.A5)
    printer.setDuplex(printer.DuplexNone)
Yan's avatar
Yan committed
94
    image = imgt.paint_image(mass_spec, spect, printer)
Yan's avatar
Yan committed
95
    printPreview.paintRequested.connect(lambda:
Yan's avatar
Yan committed
96
                                        printimage(printer, image))
97
98
99
    printPreview.exec()


Yan's avatar
Yan committed
100
101
def update_spectrum(chromatogram, spect, ds, ms, fn, chrom, config):
    if fn[0] is not None:
Yan's avatar
Yan committed
102
        slims = [spect.get_xlim(), spect.get_ylim()]
3Yan's avatar
3Yan committed
103
104
105
        ds.clear()
        [ds.append(dict(chrom_dat=i[0], masses=i[1], matrix=i[2])) for i
                in load_raw(fn[0], config.value("tmp_location"))]
Yan's avatar
Yan committed
106
107
108
        gt.populate(chromatogram, spect, ds, ms, chrom)
        spect.set_xlim(slims[0])
        spect.set_ylim(slims[1])
Yan's avatar
Yan committed
109
        gt.ann_spec(spect, ms)
Yan's avatar
Yan committed
110
111
112
        spect.get_figure().canvas.draw()


113
def dropped(event, parent, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
114
            chrom, update, config):
115
    dropurl = event.mimeData().urls()[0].toLocalFile()
116
    load_file(parent, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
117
              chrom, update, config, filename=dropurl)
Yan's avatar
Yan committed
118
119
120
121


def drag_entered(event):
    if event.mimeData().hasUrls() and event.mimeData().urls()[0]\
122
        .toLocalFile().lower().endswith('.raw'):
Yan's avatar
Yan committed
123
124
125
        event.accept()


Yan's avatar
Yan committed
126
127
128
129
130
131
132
133
134
135
136
137
138
def predictmz(form, chromatogram, spect, ds, ms, chrom):
    text = form.text()
    if text == "":
        ms["predict"] = None
        return
    slims = [spect.get_xlim(), spect.get_ylim()]
    ms["predict"] = getmzpattern(text)
    gt.populate(chromatogram, spect, ds, ms, chrom)
    spect.set_xlim(slims[0])
    spect.set_ylim(slims[1])
    spect.get_figure().canvas.draw()


3Yan's avatar
3Yan committed
139
140
141
142
143
144
def oddeven_changed(chromatogram, spectrum, ds, ms, filename, chrom, config,
        oddevenact):
    config.setValue("view/oddeven", oddevenact.isChecked())
    update_spectrum(chromatogram, spectrum, ds, ms, filename, chrom, config)


Yan's avatar
Yan committed
145
146
147
def key_pressed(event, chrom, spect, ds, ms_ds, fn, chrom_ds, config):
    if event.key() == QtCore.Qt.Key_F5:
        update_spectrum(chrom, spect, ds, ms_ds, fn, chrom_ds, config)
148
149
150
151
152
153
    if event.key() == QtCore.Qt.Key_C:
        if event.modifiers().__int__() == QtCore.Qt.ControlModifier:
            imgt.clip_spect_img(ms_ds, spect)
        if event.modifiers().__int__() == QtCore.Qt.ControlModifier+\
                                          QtCore.Qt.ShiftModifier:
            dt.clip_spectstr(spect, chrom_ds, fn)
Yan's avatar
Yan committed
154
155
    if event.key() in (QtCore.Qt.Key_Left, QtCore.Qt.Key_Right):
        gt.shift_times(event, spect, chrom, ds, ms_ds, chrom_ds)
Yan's avatar
Yan committed
156
157


Yan's avatar
Yan committed
158
159
160
161
162
def about(parent):
    """constructs window with "about" info"""
    QtWidgets.QMessageBox.information(
            parent, "About Prasopes",
            "Prasopes Finnigan raw file viewer\n\n"
Yan's avatar
Yan committed
163
            "Version: 0.0.11 (alpha)")
Yan's avatar
Yan committed
164
165


Yan's avatar
Yan committed
166
def main():
Yan's avatar
Yan committed
167
    print(autoparams)
3Yan's avatar
3Yan committed
168
169
    ds = []
    ms = dict(annotation=[], name="Spectrum", xlabel="m/z",
Yan's avatar
Yan committed
170
171
              ylabel="ion count", xtics=20, predict=None,
              params=[], headers=[])
172
    chrom = dict(x=[0], y=[0], t_start=None, t_end=None,
Yan's avatar
Yan committed
173
                 name="Chromatogram", xlabel="time(min)",
3Yan's avatar
3Yan committed
174
                 ylabel="total ion count", timesarg=[])
Yan's avatar
Yan committed
175
    filename = [None]
3Yan's avatar
3Yan committed
176
    drlcache = [None, None]
177
    update = update_signal()
Yan's avatar
Yan committed
178

Yan's avatar
Yan committed
179
180
    config = cf.settings()

Yan's avatar
Yan committed
181
    p_logger = logging.getLogger('parseLogger')
Yan's avatar
Yan committed
182
    drl_logger = logging.getLogger('drlLogger')
183
    zce_logger = logging.getLogger('zceLogger')
Yan's avatar
Yan committed
184
185
    logging.basicConfig()
    p_logger.setLevel("WARN")
Yan's avatar
Yan committed
186
187
    #drl_logger.setLevel("INFO")
    drl_logger.setLevel("DEBUG")
188
    zce_logger.setLevel("DEBUG")
Yan's avatar
Yan committed
189

Yan's avatar
Yan committed
190
    graph = Figure(figsize=(5, 4), dpi=100, facecolor="None")
Yan's avatar
Yan committed
191
192
193
    grid = graph.add_gridspec(2, 1)
    chromatogram = graph.add_subplot(grid[0,0], facecolor=(1, 1, 1, 0.8))
    spectrum = graph.add_subplot(grid[1,0], facecolor=(1, 1, 1, 0.8))
Yan's avatar
Yan committed
194
195
196
197
198
199
200
201
202
    graph.tight_layout()
    mpl_canvas = FigureCanvas(graph)
    mpl_canvas.setStyleSheet("background-color:transparent;")
    mpl_canvas.setAutoFillBackground(False)
    
    gt.pan_factory(chromatogram)
    gt.zoom_factory(chromatogram, 1.15)
    gt.pan_factory(spectrum, ms)
    gt.zoom_factory(spectrum, 1.15, ms)
203
    mass_selector = SpanSelector(
Yan's avatar
Yan committed
204
        spectrum, lambda x_min, x_max: gt.pick_masses(
205
206
207
            x_min, x_max, spectrum, ms), 'horizontal', minspan=0.01,
        useblit=True, rectprops=dict(alpha=0.15, facecolor='purple'),
        button=3)
208
    time_selector = SpanSelector(
Yan's avatar
Yan committed
209
210
211
212
213
214
         chromatogram, lambda x_min, x_max: gt.pick_times(
             x_min, x_max, spectrum, ds, chromatogram, ms, chrom), 
         'horizontal', useblit=True, rectprops=dict(
             alpha=0.15, facecolor='purple'), button=3)

    app = QtWidgets.QApplication(sys.argv)
Yan's avatar
Yan committed
215
    main_window = QtWidgets.QMainWindow(windowTitle="Prasopes")
Yan's avatar
Yan committed
216

Yan's avatar
Yan committed
217
218
219
    if QtGui.QIcon.themeName() is "":
        QtGui.QIcon.setThemeName("TangoMFK")

Yan's avatar
Yan committed
220
    openact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
221
        "document-open"), "&Open...", None)
Yan's avatar
Yan committed
222
223
    openact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_O)
    openact.triggered.connect(lambda: load_file(
224
        main_window, chromatogram, spectrum, ds, ms, filename, chrom,
225
        update, config))
Yan's avatar
Yan committed
226
    exportact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
227
        "document-save-as"), "&Export...", None)
Yan's avatar
Yan committed
228
229
    exportact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_E)
    exportact.triggered.connect(lambda: ft.export_dial(
230
        spectrum, chrom, filename, main_window))
231
232
233
    printact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
        "document-print"), "&Print", None)
    printact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_P)
Yan's avatar
Yan committed
234
235
    printact.triggered.connect(lambda: print_graph(
        ds, ms, chrom, spectrum, filename))
Yan's avatar
Yan committed
236
    settingsact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
237
        "preferences-system"), "&Settings...", None)
Yan's avatar
Yan committed
238
239
    settingsact.triggered.connect(lambda: cf.dial(main_window))
    quitact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
240
        "application-exit"), "&Quit", None)
Yan's avatar
Yan committed
241
242
243
    quitact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_Q)
    quitact.triggered.connect(main_window.close)
    zceact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
244
        "applications-utilities"), "&TSQ zce...", None)
Yan's avatar
Yan committed
245
246
    zceact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_T)
    zceact.triggered.connect(lambda: zce.dialog(
247
        main_window, ds, filename, update))
Yan's avatar
Yan committed
248
    drlact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
249
        "applications-utilities"), "&DRL...", None)
Yan's avatar
Yan committed
250
251
    drlact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_D)
    drlact.triggered.connect(lambda: drl.main_window(
252
        main_window, ds, filename, drlcache, update))
Yan's avatar
Yan committed
253
254
    aboutact = QtWidgets.QAction("&About Prasopes", None)
    aboutact.triggered.connect(lambda: about(main_window))
Yan's avatar
Yan committed
255
    autozoomy = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
256
257
258
259
        "zoom-original"), "Auto Zoom Y", None, checkable=True,
        checked=config.value("view/autozoomy", type=bool))
    autozoomy.triggered.connect(lambda: config.setValue(
                                "view/autozoomy", autozoomy.isChecked()))
260
261
262
263
264
265
266
267
268
269
270
    autozoomy.triggered.connect(lambda: gt.autozoomy(spectrum))
    intensitiesact = QtWidgets.QAction("&Show intensities", None,
            checkable=True, checked=config.value("view/intensities",
                type=bool))
    intensitiesact.triggered.connect(lambda: config.setValue(
        "view/intensities", intensitiesact.isChecked()))
    intensitiesact.triggered.connect(lambda: gt.ann_spec(spectrum, ms))
    intensitiesact.triggered.connect(lambda: spectrum.figure.canvas.draw())
    oddevenact = QtWidgets.QAction("&Odd / even", None, checkable=True,
            checked=config.value("view/oddeven", type=bool))
    oddevenact.triggered.connect(lambda:
3Yan's avatar
3Yan committed
271
272
        oddeven_changed(chromatogram, spectrum, ds, ms, filename, chrom,
            config, oddevenact))
Yan's avatar
Yan committed
273
274
275
276
277
    filebrowseract = QtWidgets.QAction(
            "&File browser", None, checkable=True,
            checked=config.value("view/filebrowservisible", type=bool))
    filebrowseract.triggered.connect(lambda: config.setValue(
        "view/filebrowservisible", filebrowseract.isChecked()))
278
279
280
    filebrowseract.triggered.connect(
        lambda: treedock.show() if filebrowseract.isChecked() == True
                else treedock.hide())
Yan's avatar
Yan committed
281

Yan's avatar
Yan committed
282
283
284
285
    predictform = QtWidgets.QLineEdit(maximumWidth=150)
    predictform.editingFinished.connect(lambda: predictmz(
        predictform, chromatogram, spectrum, ds, ms, chrom))

Yan's avatar
Yan committed
286
287
    file_menu = QtWidgets.QMenu('&File', main_window)
    main_window.menuBar().addMenu(file_menu)
Yan's avatar
Yan committed
288
289
290
    file_menu.addAction(openact)
    file_menu.addAction(exportact)
    file_menu.addSeparator()
291
292
    file_menu.addAction(printact)
    file_menu.addSeparator()
Yan's avatar
Yan committed
293
294
295
    file_menu.addAction(settingsact)
    file_menu.addSeparator()
    file_menu.addAction(quitact)
Yan's avatar
Yan committed
296
297
    tools_menu = QtWidgets.QMenu('&Tools', main_window)
    main_window.menuBar().addMenu(tools_menu)
Yan's avatar
Yan committed
298
299
    tools_menu.addAction(zceact)
    tools_menu.addAction(drlact)
300
    tools_menu.addSeparator()
Yan's avatar
Yan committed
301
302
303
    view_menu = QtWidgets.QMenu('&View', main_window)
    view_menu.addAction(filebrowseract)
    view_menu.addAction(autozoomy)
304
305
306
    view_menu.addAction(intensitiesact)
    view_menu.addSeparator()
    view_menu.addAction(oddevenact)
Yan's avatar
Yan committed
307
    main_window.menuBar().addMenu(view_menu)
Yan's avatar
Yan committed
308
309
310
    help_menu = QtWidgets.QMenu('&Help', main_window)
    main_window.menuBar().addMenu(help_menu)
    help_menu.addAction(aboutact)
Yan's avatar
Yan committed
311
312
313
314

    main_widget = QtWidgets.QWidget(main_window)
    main_window.setCentralWidget(main_widget)

Yan's avatar
Yan committed
315
316
317
318
319
    toolBar = QtWidgets.QToolBar(main_window)
    toolBar.setAllowedAreas(QtCore.Qt.TopToolBarArea)
    toolBar.setFloatable(False)
    toolBar.setMovable(False)
    toolBar.setToolButtonStyle(QtCore.Qt.ToolButtonTextBesideIcon)
Yan's avatar
Yan committed
320

Yan's avatar
Yan committed
321
322
323
    toolBar.addAction(openact)
    toolBar.addAction(exportact)
    toolBar.addSeparator()
Yan's avatar
Yan committed
324
325
326
    toolBar.addWidget(QtWidgets.QLabel("Predict Formula:"))
    toolBar.addWidget(predictform)
    toolBar.addSeparator()
Yan's avatar
Yan committed
327
328
    toolBar.addAction(zceact)
    toolBar.addAction(drlact)
Yan's avatar
Yan committed
329
330
    toolBar.addSeparator()
    toolBar.addAction(autozoomy)
Yan's avatar
Yan committed
331

332
333
    #fileModel gots stuck sometimes, better give it a separate thread
    #it complains in separate thread but work as should, maybe fix later..
3Yan's avatar
3Yan committed
334
    fileModel = QtWidgets.QFileSystemModel()
335
336
    treedockThread = QtCore.QThread()
    treedockThread.start()
337
    fileModel.setRootPath('')
338
    fileModel.moveToThread(treedockThread)
339
    activeDir = fileModel.index(config.value("open_folder"))
3Yan's avatar
3Yan committed
340
341
342
343
    treeview = QtWidgets.QTreeView()
    treeview.setModel(fileModel)
    treeview.setCurrentIndex(activeDir)
    treeview.expand(activeDir)
344
345
    dirview = QtWidgets.QListView()
    dirview.setModel(fileModel)
3Yan's avatar
3Yan committed
346
    treedock = QtWidgets.QDockWidget()
Yan's avatar
Yan committed
347
348
    if not config.value("view/filebrowservisible",type=bool):
        treedock.hide()
349
350
351
352
    treedock.setWidget(QtWidgets.QWidget())
    treedock_layout = QtWidgets.QVBoxLayout(treedock.widget())
    treedock_layout.addWidget(treeview)
    treedock_layout.addWidget(dirview)
3Yan's avatar
3Yan committed
353
354
    update.signal.connect(
        lambda: treeview.setCurrentIndex(fileModel.index(filename[0])))
355
    update.signal.connect(
3Yan's avatar
3Yan committed
356
357
        lambda: dirview.setRootIndex(fileModel.index(
            os.path.dirname(os.path.realpath(filename[0])))))
3Yan's avatar
3Yan committed
358
359
360
    treeview.doubleClicked.connect(lambda index: load_file(
        main_window, chromatogram, spectrum, ds, ms, filename, chrom,
        update, config, filename=fileModel.filePath(index)))
361
362
363
364
365
366
    dirview.doubleClicked.connect(lambda index: load_file(
        main_window, chromatogram, spectrum, ds, ms, filename, chrom,
        update, config, filename=fileModel.filePath(index)))
    dirview.clicked.connect(lambda index: load_file(
        main_window, chromatogram, spectrum, ds, ms, filename, chrom,
        update, config, filename=fileModel.filePath(index)))
367
368
369
370
371
372
    def closeOverride(event, dock, action):
        action.setChecked(False)
        dock.hide()
        event.ignore()
    treedock.closeEvent = lambda event: closeOverride(
            event, treedock, filebrowseract)
3Yan's avatar
3Yan committed
373
374
    main_window.addDockWidget(QtCore.Qt.LeftDockWidgetArea, treedock)

Yan's avatar
Yan committed
375
    layout = QtWidgets.QVBoxLayout(main_widget)
Yan's avatar
Yan committed
376
    layout.addWidget(toolBar)
Yan's avatar
Yan committed
377
378
    layout.addWidget(mpl_canvas)

3Yan's avatar
3Yan committed
379

Yan's avatar
Yan committed
380
    main_window.setFocus()
Yan's avatar
Yan committed
381
382
    main_window.dragEnterEvent = lambda event: drag_entered(event)
    main_window.dropEvent = lambda event: dropped(
383
        event, main_window, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
384
        chrom, update, config)
Yan's avatar
Yan committed
385
    main_window.setAcceptDrops(True)
Yan's avatar
Yan committed
386
387
    main_window.keyPressEvent = lambda event: key_pressed(event, chromatogram,
            spectrum, ds, ms, filename, chrom, config)
Yan's avatar
Yan committed
388

Yan's avatar
Yan committed
389
    if len(sys.argv) == 2:
390
        load_file(main_window, chromatogram, spectrum, ds, ms,
3Yan's avatar
3Yan committed
391
                  filename, chrom, update, config, filename=sys.argv[1])
Yan's avatar
Yan committed
392
    else:
3Yan's avatar
3Yan committed
393
394
        gt.pop_plot(0, 0, spectrum, ms)
        gt.pop_plot(0, 0, chromatogram, chrom)
Yan's avatar
Yan committed
395
396
397

    main_window.show()
    sys.exit(app.exec_())
Yan's avatar
Yan committed
398

399

Yan's avatar
Yan committed
400
401
if __name__ == "__main__":
    main()