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


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


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

Yan's avatar
Yan committed
61

62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
def print_graph(mass_spec, chrom_spec, spect, fn):
    printfig = Figure(figsize=(10, 4), dpi=300)
    printplot = printfig.add_subplot(111)
    if cf.settings().value("view/oddeven",type=bool):
        print("TODO")
    else:
        gt.pop_plot(printplot, mass_spec)
    printPreview = QtPrintSupport.QPrintPreviewDialog()
    printer = printPreview.printer()
    def ubercommand():
        cache_file = BytesIO()
        printfig.savefig(cache_file)
        cache_file.seek(0)
        image = QtGui.QImage.fromData(cache_file.read())
        printer.setResolution(600)
        painter = QtGui.QPainter(printPreview.printer())
        painter.drawImage(QtCore.QPoint(0,0), image)
        painter.end()
    printPreview.paintRequested.connect(lambda: ubercommand())
    printPreview.exec()


84
def update_spectrum(condition, chromatogram, spect, ds,
Yan's avatar
Yan committed
85
                    ms, fn, chrom, config):
86
    if fn[0] is not None and condition:
Yan's avatar
Yan committed
87
88
89
90
91
92
93
94
95
        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])
        spect.get_figure().canvas.draw()


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


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


Yan's avatar
Yan committed
109
110
111
112
113
114
115
116
117
118
119
120
121
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
122
123
124
125
126
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
127
            "Version: 0.0.10 (alpha)")
Yan's avatar
Yan committed
128
129


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

Yan's avatar
Yan committed
143
144
    config = cf.settings()

Yan's avatar
Yan committed
145
    p_logger = logging.getLogger('parseLogger')
Yan's avatar
Yan committed
146
    drl_logger = logging.getLogger('drlLogger')
147
    zce_logger = logging.getLogger('zceLogger')
Yan's avatar
Yan committed
148
149
    logging.basicConfig()
    p_logger.setLevel("WARN")
Yan's avatar
Yan committed
150
151
    #drl_logger.setLevel("INFO")
    drl_logger.setLevel("DEBUG")
152
    zce_logger.setLevel("DEBUG")
Yan's avatar
Yan committed
153

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

Yan's avatar
Yan committed
183
184
185
    if QtGui.QIcon.themeName() is "":
        QtGui.QIcon.setThemeName("TangoMFK")

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

Yan's avatar
Yan committed
249
250
251
252
    predictform = QtWidgets.QLineEdit(maximumWidth=150)
    predictform.editingFinished.connect(lambda: predictmz(
        predictform, chromatogram, spectrum, ds, ms, chrom))

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

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

Yan's avatar
Yan committed
282
283
284
285
286
    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
287

Yan's avatar
Yan committed
288
289
290
    toolBar.addAction(openact)
    toolBar.addAction(exportact)
    toolBar.addSeparator()
Yan's avatar
Yan committed
291
292
293
    toolBar.addWidget(QtWidgets.QLabel("Predict Formula:"))
    toolBar.addWidget(predictform)
    toolBar.addSeparator()
Yan's avatar
Yan committed
294
295
    toolBar.addAction(zceact)
    toolBar.addAction(drlact)
Yan's avatar
Yan committed
296
297
    toolBar.addSeparator()
    toolBar.addAction(autozoomy)
Yan's avatar
Yan committed
298

Yan's avatar
Yan committed
299

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

Yan's avatar
Yan committed
338
    layout = QtWidgets.QVBoxLayout(main_widget)
Yan's avatar
Yan committed
339
    layout.addWidget(toolBar)
Yan's avatar
Yan committed
340
341
342
    layout.addWidget(mpl_canvas)
    mpl_canvas.setFocus()

3Yan's avatar
3Yan committed
343

344
    mpl_canvas.mpl_connect('key_press_event', lambda event:
345
        update_spectrum(event.key == "f5", chromatogram, spectrum, ds, ms,
Yan's avatar
Yan committed
346
                        filename, chrom, config))
Yan's avatar
Yan committed
347
348
    main_window.dragEnterEvent = lambda event: drag_entered(event)
    main_window.dropEvent = lambda event: dropped(
349
        event, main_window, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
350
        chrom, update, config)
Yan's avatar
Yan committed
351
    main_window.setAcceptDrops(True)
Yan's avatar
Yan committed
352

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

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

363

Yan's avatar
Yan committed
364
365
if __name__ == "__main__":
    main()