__main__.py 14.3 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
46
47
48
49
50
51
52
53
54
55
        try:
            d_set['chrom_dat'], d_set['masses'], d_set['matrix'] \
                = load_raw(filename, settings.value("tmp_location"))
        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
56
        gt.populate(chrom_plot, spc, d_set, ms_s, chrom_s)
3Yan's avatar
3Yan committed
57
        fn[0] = filename
58
59
        parent.setWindowTitle("Prasopes - {}".format(
            os.path.basename(filename)))
60
        update.signal.emit()
3Yan's avatar
3Yan committed
61

Yan's avatar
Yan committed
62

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


Yan's avatar
Yan committed
80
81
def update_spectrum(chromatogram, spect, ds, ms, fn, chrom, config):
    if fn[0] is not None:
Yan's avatar
Yan committed
82
83
84
85
86
87
        slims = [spect.get_xlim(), spect.get_ylim()]
        ds['chrom_dat'], ds['masses'], ds['matrix'] \
            = load_raw(fn[0], config.value("tmp_location"))
        gt.populate(chromatogram, spect, ds, ms, chrom)
        spect.set_xlim(slims[0])
        spect.set_ylim(slims[1])
Yan's avatar
Yan committed
88
        gt.ann_spec(spect, ms)
Yan's avatar
Yan committed
89
90
91
        spect.get_figure().canvas.draw()


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


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


Yan's avatar
Yan committed
105
106
107
108
109
110
111
112
113
114
115
116
117
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()


Yan's avatar
Yan committed
118
119
120
121
122
123
124
125
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
126
127
128
129
130
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
131
            "Version: 0.0.11 (alpha)")
Yan's avatar
Yan committed
132
133


Yan's avatar
Yan committed
134
def main():
Yan's avatar
Yan committed
135
136
    # ds for data_set, ms for mass_spec_dataset, chrom for chrom_dataset
    ds = dict(chrom_dat=None, masses=None, matrix=None)
137
    ms = dict(x=[0], y=[0], annotation=[],
Yan's avatar
Yan committed
138
              name="Spectrum", xlabel="m/z", ylabel="ion count",
Yan's avatar
Yan committed
139
              xtics=20, predict=None)
140
    chrom = dict(x=[0], y=[0], t_start=None, t_end=None,
Yan's avatar
Yan committed
141
                 name="Chromatogram", xlabel="time(min)",
142
                 ylabel="total ion count", timesarg=None)
Yan's avatar
Yan committed
143
    filename = [None]
3Yan's avatar
3Yan committed
144
    drlcache = [None, None]
145
    update = update_signal()
Yan's avatar
Yan committed
146

Yan's avatar
Yan committed
147
148
    config = cf.settings()

Yan's avatar
Yan committed
149
    p_logger = logging.getLogger('parseLogger')
Yan's avatar
Yan committed
150
    drl_logger = logging.getLogger('drlLogger')
151
    zce_logger = logging.getLogger('zceLogger')
Yan's avatar
Yan committed
152
153
    logging.basicConfig()
    p_logger.setLevel("WARN")
Yan's avatar
Yan committed
154
155
    #drl_logger.setLevel("INFO")
    drl_logger.setLevel("DEBUG")
156
    zce_logger.setLevel("DEBUG")
Yan's avatar
Yan committed
157

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

Yan's avatar
Yan committed
186
187
188
    if QtGui.QIcon.themeName() is "":
        QtGui.QIcon.setThemeName("TangoMFK")

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

Yan's avatar
Yan committed
253
254
255
256
    predictform = QtWidgets.QLineEdit(maximumWidth=150)
    predictform.editingFinished.connect(lambda: predictmz(
        predictform, chromatogram, spectrum, ds, ms, chrom))

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

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

Yan's avatar
Yan committed
286
287
288
289
290
    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
291

Yan's avatar
Yan committed
292
293
294
    toolBar.addAction(openact)
    toolBar.addAction(exportact)
    toolBar.addSeparator()
Yan's avatar
Yan committed
295
296
297
    toolBar.addWidget(QtWidgets.QLabel("Predict Formula:"))
    toolBar.addWidget(predictform)
    toolBar.addSeparator()
Yan's avatar
Yan committed
298
299
    toolBar.addAction(zceact)
    toolBar.addAction(drlact)
Yan's avatar
Yan committed
300
301
    toolBar.addSeparator()
    toolBar.addAction(autozoomy)
Yan's avatar
Yan committed
302

Yan's avatar
Yan committed
303

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

Yan's avatar
Yan committed
342
    layout = QtWidgets.QVBoxLayout(main_widget)
Yan's avatar
Yan committed
343
    layout.addWidget(toolBar)
Yan's avatar
Yan committed
344
345
    layout.addWidget(mpl_canvas)

3Yan's avatar
3Yan committed
346

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

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

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

366

Yan's avatar
Yan committed
367
368
if __name__ == "__main__":
    main()