__main__.py 14.4 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
import numpy as np
Yan's avatar
Yan committed
15
16
import prasopes.config as cf
import prasopes.drltools as drl
Yan's avatar
Yan committed
17
import prasopes.filetools as ft
Yan's avatar
Yan committed
18
19
import prasopes.graphtools as gt
import prasopes.imagetools as imgt
Yan's avatar
Yan committed
20
import prasopes.zcetools as zce
Yan's avatar
Yan committed
21
import prasopes.tangoicons
Yan's avatar
Yan committed
22
23
24
import sys
import matplotlib
import logging
3Yan's avatar
3Yan committed
25
import os.path
Yan's avatar
Yan committed
26
27
28
matplotlib.use("Qt5Agg")


29
30
31
32
class update_signal(QtCore.QObject):
    signal = QtCore.pyqtSignal()


33
def load_file(parent, chrom_plot, spc, d_set, ms_s, fn, chrom_s,
34
              update, settings, filename=None):
Yan's avatar
Yan committed
35
    """populates dataset and plots it"""
3Yan's avatar
3Yan committed
36
37
38
    directory=settings.value("open_folder")
    if fn[0] is not None:
        directory=fn[0]
Yan's avatar
Yan committed
39
40
    if filename == None:
        filename = QtWidgets.QFileDialog.getOpenFileName(
3Yan's avatar
3Yan committed
41
42
43
            caption="Open spectrum",
            directory=directory,
            filter="Finnigan RAW files (*.raw, *.RAW)")[0]
3Yan's avatar
3Yan committed
44
    if filename is not '' and os.path.isfile(filename):
Yan's avatar
Yan committed
45
        try:
3Yan's avatar
3Yan committed
46
47
48
            d_set.clear()
            [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
49
50
51
52
53
54
55
56
        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
57
        gt.populate(chrom_plot, spc, d_set, ms_s, chrom_s)
3Yan's avatar
3Yan committed
58
        fn[0] = filename
59
60
        parent.setWindowTitle("Prasopes - {}".format(
            os.path.basename(filename)))
61
        update.signal.emit()
3Yan's avatar
3Yan committed
62

Yan's avatar
Yan committed
63

Yan's avatar
Yan committed
64
65
def print_graph(data_set, mass_spec, chrom_spec, spect, fn):
    def printimage(printdevice, img):
Yan's avatar
Yan committed
66
        printer.setResolution(600)
Yan's avatar
Yan committed
67
        painter = QtGui.QPainter(printdevice)
Yan's avatar
Yan committed
68
        painter.drawImage(0,0,img)
Yan's avatar
Yan committed
69
70
        painter.end()
    #TODO: substitute the QPrintPreviewDialog with QPrintPreviewWidget
71
72
    printPreview = QtPrintSupport.QPrintPreviewDialog()
    printer = printPreview.printer()
Yan's avatar
Yan committed
73
74
    printer.setPageSize(printer.A5)
    printer.setDuplex(printer.DuplexNone)
Yan's avatar
Yan committed
75
    image = imgt.paint_image(mass_spec, spect, printer)
Yan's avatar
Yan committed
76
    printPreview.paintRequested.connect(lambda:
Yan's avatar
Yan committed
77
                                        printimage(printer, image))
78
79
80
    printPreview.exec()


Yan's avatar
Yan committed
81
82
def update_spectrum(chromatogram, spect, ds, ms, fn, chrom, config):
    if fn[0] is not None:
Yan's avatar
Yan committed
83
        slims = [spect.get_xlim(), spect.get_ylim()]
3Yan's avatar
3Yan committed
84
85
86
        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
87
88
89
        gt.populate(chromatogram, spect, ds, ms, chrom)
        spect.set_xlim(slims[0])
        spect.set_ylim(slims[1])
Yan's avatar
Yan committed
90
        gt.ann_spec(spect, ms)
Yan's avatar
Yan committed
91
92
93
        spect.get_figure().canvas.draw()


94
def dropped(event, parent, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
95
            chrom, update, config):
96
    dropurl = event.mimeData().urls()[0].toLocalFile()
97
    load_file(parent, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
98
              chrom, update, config, filename=dropurl)
Yan's avatar
Yan committed
99
100
101
102


def drag_entered(event):
    if event.mimeData().hasUrls() and event.mimeData().urls()[0]\
103
        .toLocalFile().lower().endswith('.raw'):
Yan's avatar
Yan committed
104
105
106
        event.accept()


Yan's avatar
Yan committed
107
108
109
110
111
112
113
114
115
116
117
118
119
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
120
121
122
123
124
125
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
126
127
128
129
130
131
132
133
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)
    if event.key() == QtCore.Qt.Key_C\
        and event.modifiers().__int__() == QtCore.Qt.ControlModifier:
        imgt.clip_spect_img(ms_ds, spect)


Yan's avatar
Yan committed
134
135
136
137
138
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
139
            "Version: 0.0.11 (alpha)")
Yan's avatar
Yan committed
140
141


Yan's avatar
Yan committed
142
def main():
3Yan's avatar
3Yan committed
143
144
145
    ds = []
    ms = dict(annotation=[], name="Spectrum", xlabel="m/z",
              ylabel="ion count", xtics=20, predict=None)
146
    chrom = dict(x=[0], y=[0], t_start=None, t_end=None,
Yan's avatar
Yan committed
147
                 name="Chromatogram", xlabel="time(min)",
3Yan's avatar
3Yan committed
148
                 ylabel="total ion count", timesarg=[])
Yan's avatar
Yan committed
149
    filename = [None]
3Yan's avatar
3Yan committed
150
    drlcache = [None, None]
151
    update = update_signal()
Yan's avatar
Yan committed
152

Yan's avatar
Yan committed
153
154
    config = cf.settings()

Yan's avatar
Yan committed
155
    p_logger = logging.getLogger('parseLogger')
Yan's avatar
Yan committed
156
    drl_logger = logging.getLogger('drlLogger')
157
    zce_logger = logging.getLogger('zceLogger')
Yan's avatar
Yan committed
158
159
    logging.basicConfig()
    p_logger.setLevel("WARN")
Yan's avatar
Yan committed
160
161
    #drl_logger.setLevel("INFO")
    drl_logger.setLevel("DEBUG")
162
    zce_logger.setLevel("DEBUG")
Yan's avatar
Yan committed
163

Yan's avatar
Yan committed
164
    graph = Figure(figsize=(5, 4), dpi=100, facecolor="None")
Yan's avatar
Yan committed
165
166
167
    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
168
169
170
171
172
173
174
    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)
175
    gt.arrow_factory(spectrum, chromatogram, ds, ms, chrom)
Yan's avatar
Yan committed
176
177
    gt.pan_factory(spectrum, ms)
    gt.zoom_factory(spectrum, 1.15, ms)
178
    mass_selector = SpanSelector(
Yan's avatar
Yan committed
179
        spectrum, lambda x_min, x_max: gt.pick_masses(
180
181
182
            x_min, x_max, spectrum, ms), 'horizontal', minspan=0.01,
        useblit=True, rectprops=dict(alpha=0.15, facecolor='purple'),
        button=3)
183
    time_selector = SpanSelector(
Yan's avatar
Yan committed
184
185
186
187
188
189
         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
190
    main_window = QtWidgets.QMainWindow(windowTitle="Prasopes")
Yan's avatar
Yan committed
191

Yan's avatar
Yan committed
192
193
194
    if QtGui.QIcon.themeName() is "":
        QtGui.QIcon.setThemeName("TangoMFK")

Yan's avatar
Yan committed
195
    openact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
196
        "document-open"), "&Open...", None)
Yan's avatar
Yan committed
197
198
    openact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_O)
    openact.triggered.connect(lambda: load_file(
199
        main_window, chromatogram, spectrum, ds, ms, filename, chrom,
200
        update, config))
Yan's avatar
Yan committed
201
    exportact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
202
        "document-save-as"), "&Export...", None)
Yan's avatar
Yan committed
203
204
205
    exportact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_E)
    exportact.triggered.connect(lambda: ft.export_dial(
        ms, chrom, filename, main_window))
206
207
208
    printact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
        "document-print"), "&Print", None)
    printact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_P)
Yan's avatar
Yan committed
209
210
    printact.triggered.connect(lambda: print_graph(
        ds, ms, chrom, spectrum, filename))
Yan's avatar
Yan committed
211
    settingsact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
212
        "preferences-system"), "&Settings...", None)
Yan's avatar
Yan committed
213
214
    settingsact.triggered.connect(lambda: cf.dial(main_window))
    quitact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
215
        "application-exit"), "&Quit", None)
Yan's avatar
Yan committed
216
217
218
    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
219
        "applications-utilities"), "&TSQ zce...", None)
Yan's avatar
Yan committed
220
221
    zceact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_T)
    zceact.triggered.connect(lambda: zce.dialog(
222
        main_window, ds, filename, update))
Yan's avatar
Yan committed
223
    drlact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
224
        "applications-utilities"), "&DRL...", None)
Yan's avatar
Yan committed
225
226
    drlact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_D)
    drlact.triggered.connect(lambda: drl.main_window(
227
        main_window, ds, filename, drlcache, update))
Yan's avatar
Yan committed
228
229
    aboutact = QtWidgets.QAction("&About Prasopes", None)
    aboutact.triggered.connect(lambda: about(main_window))
Yan's avatar
Yan committed
230
    autozoomy = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
231
232
233
234
        "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()))
235
236
237
238
239
240
241
242
243
244
245
    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
246
247
        oddeven_changed(chromatogram, spectrum, ds, ms, filename, chrom,
            config, oddevenact))
Yan's avatar
Yan committed
248
249
250
251
252
    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()))
253
254
255
    filebrowseract.triggered.connect(
        lambda: treedock.show() if filebrowseract.isChecked() == True
                else treedock.hide())
Yan's avatar
Yan committed
256

Yan's avatar
Yan committed
257
258
259
260
    predictform = QtWidgets.QLineEdit(maximumWidth=150)
    predictform.editingFinished.connect(lambda: predictmz(
        predictform, chromatogram, spectrum, ds, ms, chrom))

Yan's avatar
Yan committed
261
262
    file_menu = QtWidgets.QMenu('&File', main_window)
    main_window.menuBar().addMenu(file_menu)
Yan's avatar
Yan committed
263
264
265
    file_menu.addAction(openact)
    file_menu.addAction(exportact)
    file_menu.addSeparator()
266
267
    file_menu.addAction(printact)
    file_menu.addSeparator()
Yan's avatar
Yan committed
268
269
270
    file_menu.addAction(settingsact)
    file_menu.addSeparator()
    file_menu.addAction(quitact)
Yan's avatar
Yan committed
271
272
    tools_menu = QtWidgets.QMenu('&Tools', main_window)
    main_window.menuBar().addMenu(tools_menu)
Yan's avatar
Yan committed
273
274
    tools_menu.addAction(zceact)
    tools_menu.addAction(drlact)
275
    tools_menu.addSeparator()
Yan's avatar
Yan committed
276
277
278
    view_menu = QtWidgets.QMenu('&View', main_window)
    view_menu.addAction(filebrowseract)
    view_menu.addAction(autozoomy)
279
280
281
    view_menu.addAction(intensitiesact)
    view_menu.addSeparator()
    view_menu.addAction(oddevenact)
Yan's avatar
Yan committed
282
    main_window.menuBar().addMenu(view_menu)
Yan's avatar
Yan committed
283
284
285
    help_menu = QtWidgets.QMenu('&Help', main_window)
    main_window.menuBar().addMenu(help_menu)
    help_menu.addAction(aboutact)
Yan's avatar
Yan committed
286
287
288
289

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

Yan's avatar
Yan committed
290
291
292
293
294
    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
295

Yan's avatar
Yan committed
296
297
298
    toolBar.addAction(openact)
    toolBar.addAction(exportact)
    toolBar.addSeparator()
Yan's avatar
Yan committed
299
300
301
    toolBar.addWidget(QtWidgets.QLabel("Predict Formula:"))
    toolBar.addWidget(predictform)
    toolBar.addSeparator()
Yan's avatar
Yan committed
302
303
    toolBar.addAction(zceact)
    toolBar.addAction(drlact)
Yan's avatar
Yan committed
304
305
    toolBar.addSeparator()
    toolBar.addAction(autozoomy)
Yan's avatar
Yan committed
306

Yan's avatar
Yan committed
307

3Yan's avatar
3Yan committed
308
    fileModel = QtWidgets.QFileSystemModel()
309
310
    fileModel.setRootPath('')
    activeDir = fileModel.index(config.value("open_folder"))
3Yan's avatar
3Yan committed
311
312
313
314
    treeview = QtWidgets.QTreeView()
    treeview.setModel(fileModel)
    treeview.setCurrentIndex(activeDir)
    treeview.expand(activeDir)
315
316
    dirview = QtWidgets.QListView()
    dirview.setModel(fileModel)
3Yan's avatar
3Yan committed
317
    treedock = QtWidgets.QDockWidget()
Yan's avatar
Yan committed
318
319
    if not config.value("view/filebrowservisible",type=bool):
        treedock.hide()
320
321
322
323
    treedock.setWidget(QtWidgets.QWidget())
    treedock_layout = QtWidgets.QVBoxLayout(treedock.widget())
    treedock_layout.addWidget(treeview)
    treedock_layout.addWidget(dirview)
3Yan's avatar
3Yan committed
324
325
    update.signal.connect(
        lambda: treeview.setCurrentIndex(fileModel.index(filename[0])))
326
    update.signal.connect(
3Yan's avatar
3Yan committed
327
328
        lambda: dirview.setRootIndex(fileModel.index(
            os.path.dirname(os.path.realpath(filename[0])))))
3Yan's avatar
3Yan committed
329
330
331
    treeview.doubleClicked.connect(lambda index: load_file(
        main_window, chromatogram, spectrum, ds, ms, filename, chrom,
        update, config, filename=fileModel.filePath(index)))
332
333
334
335
336
337
    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)))
338
339
340
341
342
343
    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
344
345
    main_window.addDockWidget(QtCore.Qt.LeftDockWidgetArea, treedock)

Yan's avatar
Yan committed
346
    layout = QtWidgets.QVBoxLayout(main_widget)
Yan's avatar
Yan committed
347
    layout.addWidget(toolBar)
Yan's avatar
Yan committed
348
349
    layout.addWidget(mpl_canvas)

3Yan's avatar
3Yan committed
350

Yan's avatar
Yan committed
351
    main_window.setFocus()
Yan's avatar
Yan committed
352
353
    main_window.dragEnterEvent = lambda event: drag_entered(event)
    main_window.dropEvent = lambda event: dropped(
354
        event, main_window, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
355
        chrom, update, config)
Yan's avatar
Yan committed
356
    main_window.setAcceptDrops(True)
Yan's avatar
Yan committed
357
358
    main_window.keyPressEvent = lambda event: key_pressed(event, chromatogram,
            spectrum, ds, ms, filename, chrom, config)
Yan's avatar
Yan committed
359

Yan's avatar
Yan committed
360
    if len(sys.argv) == 2:
361
        load_file(main_window, chromatogram, spectrum, ds, ms,
3Yan's avatar
3Yan committed
362
                  filename, chrom, update, config, filename=sys.argv[1])
Yan's avatar
Yan committed
363
    else:
3Yan's avatar
3Yan committed
364
365
        gt.pop_plot(0, 0, spectrum, ms)
        gt.pop_plot(0, 0, chromatogram, chrom)
Yan's avatar
Yan committed
366
367
368

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

370

Yan's avatar
Yan committed
371
372
if __name__ == "__main__":
    main()