__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
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
        #!TODO: copy the printer settings from TSQ
Yan's avatar
Yan committed
66
67
68
69
        painter = QtGui.QPainter(printdevice)
        painter.drawImage(QtCore.QPoint(0,0), img)
        painter.end()
    #TODO: substitute the QPrintPreviewDialog with QPrintPreviewWidget
70
71
    printPreview = QtPrintSupport.QPrintPreviewDialog()
    printer = printPreview.printer()
Yan's avatar
Yan committed
72
73
74
    printer.setResolution(600)
    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
77
    printPreview.paintRequested.connect(lambda:
            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
84
85
86
87
88
        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
89
        gt.ann_spec(spect, ms)
Yan's avatar
Yan committed
90
91
92
        spect.get_figure().canvas.draw()


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


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


Yan's avatar
Yan committed
106
107
108
109
110
111
112
113
114
115
116
117
118
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
119
120
121
122
123
124
125
126
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
127
128
129
130
131
def about(parent):
    """constructs window with "about" info"""
    QtWidgets.QMessageBox.information(
            parent, "About Prasopes",
            "Prasopes Finnigan raw file viewer\n\n"
3Yan's avatar
3Yan committed
132
            "Version: 0.0.10 (alpha)")
Yan's avatar
Yan committed
133
134


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

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

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

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

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

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

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

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

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

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

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

Yan's avatar
Yan committed
304

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

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

3Yan's avatar
3Yan committed
347

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

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

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

367

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