__main__.py 14.5 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
            d_set.clear()
Yan's avatar
Yan committed
47
            chrom_s['timesarg'].clear()
3Yan's avatar
3Yan committed
48
49
            [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
50
51
52
53
54
55
56
57
        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
58
        gt.populate(chrom_plot, spc, d_set, ms_s, chrom_s)
3Yan's avatar
3Yan committed
59
        fn[0] = filename
60
61
        parent.setWindowTitle("Prasopes - {}".format(
            os.path.basename(filename)))
62
        update.signal.emit()
3Yan's avatar
3Yan committed
63

Yan's avatar
Yan committed
64

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


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


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


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


Yan's avatar
Yan committed
108
109
110
111
112
113
114
115
116
117
118
119
120
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
121
122
123
124
125
126
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
127
128
129
130
131
132
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
133
134
    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
135
136


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


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

Yan's avatar
Yan committed
156
157
    config = cf.settings()

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

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

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

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

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

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

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

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

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

Yan's avatar
Yan committed
309

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

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

3Yan's avatar
3Yan committed
352

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

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

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

372

Yan's avatar
Yan committed
373
374
if __name__ == "__main__":
    main()