__main__.py 13.1 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
Yan's avatar
Yan committed
10
from rawprasslib import load_raw
Yan's avatar
Yan committed
11
from rawprasslib import rawprasslib
Yan's avatar
Yan committed
12
from prasopes.predictmz import predict as getmzpattern
Yan's avatar
Yan committed
13
14
import prasopes.graphtools as gt
import prasopes.filetools as ft
Yan's avatar
Yan committed
15
import prasopes.zcetools as zce
Yan's avatar
Yan committed
16
import prasopes.drltools as drl
Yan's avatar
Yan committed
17
import prasopes.config as cf
Yan's avatar
Yan committed
18
import prasopes.tangoicons
Yan's avatar
Yan committed
19
20
21
import sys
import matplotlib
import logging
3Yan's avatar
3Yan committed
22
import os.path
Yan's avatar
Yan committed
23
24
25
matplotlib.use("Qt5Agg")


26
27
28
29
class update_signal(QtCore.QObject):
    signal = QtCore.pyqtSignal()


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

Yan's avatar
Yan committed
59

60
def update_spectrum(condition, chromatogram, spect, ds,
Yan's avatar
Yan committed
61
                    ms, fn, chrom, config):
62
    if fn[0] is not None and condition:
Yan's avatar
Yan committed
63
64
65
66
67
68
69
70
71
        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()


72
def dropped(event, parent, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
73
            chrom, update, config):
74
    dropurl = event.mimeData().urls()[0].toLocalFile()
75
    load_file(parent, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
76
              chrom, update, config, filename=dropurl)
Yan's avatar
Yan committed
77
78
79
80


def drag_entered(event):
    if event.mimeData().hasUrls() and event.mimeData().urls()[0]\
81
        .toLocalFile().lower().endswith('.raw'):
Yan's avatar
Yan committed
82
83
84
        event.accept()


Yan's avatar
Yan committed
85
86
87
88
89
90
91
92
93
94
95
96
97
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
98
99
100
101
102
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
103
            "Version: 0.0.10 (alpha)")
Yan's avatar
Yan committed
104
105


Yan's avatar
Yan committed
106
def main():
Yan's avatar
Yan committed
107
108
    # ds for data_set, ms for mass_spec_dataset, chrom for chrom_dataset
    ds = dict(chrom_dat=None, masses=None, matrix=None)
109
    ms = dict(x=[0], y=[0], annotation=[],
Yan's avatar
Yan committed
110
              name="Spectrum", xlabel="m/z", ylabel="ion count",
Yan's avatar
Yan committed
111
              xtics=20, predict=None)
112
    chrom = dict(x=[0], y=[0], t_start=None, t_end=None,
Yan's avatar
Yan committed
113
                 name="Chromatogram", xlabel="time(min)",
114
                 ylabel="total ion count", timesarg=None)
Yan's avatar
Yan committed
115
    filename = [None]
3Yan's avatar
3Yan committed
116
    drlcache = [None, None]
117
    update = update_signal()
Yan's avatar
Yan committed
118

Yan's avatar
Yan committed
119
120
    config = cf.settings()

Yan's avatar
Yan committed
121
    p_logger = logging.getLogger('parseLogger')
Yan's avatar
Yan committed
122
    drl_logger = logging.getLogger('drlLogger')
123
    zce_logger = logging.getLogger('zceLogger')
Yan's avatar
Yan committed
124
125
    logging.basicConfig()
    p_logger.setLevel("WARN")
Yan's avatar
Yan committed
126
127
    #drl_logger.setLevel("INFO")
    drl_logger.setLevel("DEBUG")
128
    zce_logger.setLevel("DEBUG")
Yan's avatar
Yan committed
129

Yan's avatar
Yan committed
130
    graph = Figure(figsize=(5, 4), dpi=100, facecolor="None")
Yan's avatar
Yan committed
131
132
133
    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
134
135
136
137
138
139
140
141
    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)
142
    gt.arrow_factory(spectrum, chromatogram, ds, ms, chrom)
Yan's avatar
Yan committed
143
144
    gt.pan_factory(spectrum, ms)
    gt.zoom_factory(spectrum, 1.15, ms)
145
    mass_selector = SpanSelector(
Yan's avatar
Yan committed
146
        spectrum, lambda x_min, x_max: gt.pick_masses(
147
148
149
            x_min, x_max, spectrum, ms), 'horizontal', minspan=0.01,
        useblit=True, rectprops=dict(alpha=0.15, facecolor='purple'),
        button=3)
150
    time_selector = SpanSelector(
Yan's avatar
Yan committed
151
152
153
154
155
156
         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
157
    main_window = QtWidgets.QMainWindow(windowTitle="Prasopes")
Yan's avatar
Yan committed
158

Yan's avatar
Yan committed
159
160
161
    if QtGui.QIcon.themeName() is "":
        QtGui.QIcon.setThemeName("TangoMFK")

Yan's avatar
Yan committed
162
    openact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
163
        "document-open"), "&Open...", None)
Yan's avatar
Yan committed
164
165
    openact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_O)
    openact.triggered.connect(lambda: load_file(
166
        main_window, chromatogram, spectrum, ds, ms, filename, chrom,
167
        update, config))
Yan's avatar
Yan committed
168
    exportact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
169
        "document-save-as"), "&Export...", None)
Yan's avatar
Yan committed
170
171
172
173
    exportact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_E)
    exportact.triggered.connect(lambda: ft.export_dial(
        ms, chrom, filename, main_window))
    settingsact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
174
        "preferences-system"), "&Settings...", None)
Yan's avatar
Yan committed
175
176
    settingsact.triggered.connect(lambda: cf.dial(main_window))
    quitact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
177
        "application-exit"), "&Quit", None)
Yan's avatar
Yan committed
178
179
180
    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
181
        "applications-utilities"), "&TSQ zce...", None)
Yan's avatar
Yan committed
182
183
    zceact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_T)
    zceact.triggered.connect(lambda: zce.dialog(
184
        main_window, ds, filename, update))
Yan's avatar
Yan committed
185
    drlact = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
186
        "applications-utilities"), "&DRL...", None)
Yan's avatar
Yan committed
187
188
    drlact.setShortcut(QtCore.Qt.CTRL + QtCore.Qt.Key_D)
    drlact.triggered.connect(lambda: drl.main_window(
189
        main_window, ds, filename, drlcache, update))
Yan's avatar
Yan committed
190
191
    aboutact = QtWidgets.QAction("&About Prasopes", None)
    aboutact.triggered.connect(lambda: about(main_window))
Yan's avatar
Yan committed
192
    autozoomy = QtWidgets.QAction(QtGui.QIcon.fromTheme(
Yan's avatar
Yan committed
193
194
195
196
        "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()))
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
    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
212
213
214
215
216
    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()))
217
218
219
    filebrowseract.triggered.connect(
        lambda: treedock.show() if filebrowseract.isChecked() == True
                else treedock.hide())
Yan's avatar
Yan committed
220

Yan's avatar
Yan committed
221
222
223
224
    predictform = QtWidgets.QLineEdit(maximumWidth=150)
    predictform.editingFinished.connect(lambda: predictmz(
        predictform, chromatogram, spectrum, ds, ms, chrom))

Yan's avatar
Yan committed
225
226
    file_menu = QtWidgets.QMenu('&File', main_window)
    main_window.menuBar().addMenu(file_menu)
Yan's avatar
Yan committed
227
228
229
230
231
232
    file_menu.addAction(openact)
    file_menu.addAction(exportact)
    file_menu.addSeparator()
    file_menu.addAction(settingsact)
    file_menu.addSeparator()
    file_menu.addAction(quitact)
Yan's avatar
Yan committed
233
234
    tools_menu = QtWidgets.QMenu('&Tools', main_window)
    main_window.menuBar().addMenu(tools_menu)
Yan's avatar
Yan committed
235
236
    tools_menu.addAction(zceact)
    tools_menu.addAction(drlact)
237
    tools_menu.addSeparator()
Yan's avatar
Yan committed
238
239
240
    view_menu = QtWidgets.QMenu('&View', main_window)
    view_menu.addAction(filebrowseract)
    view_menu.addAction(autozoomy)
241
242
243
    view_menu.addAction(intensitiesact)
    view_menu.addSeparator()
    view_menu.addAction(oddevenact)
Yan's avatar
Yan committed
244
    main_window.menuBar().addMenu(view_menu)
Yan's avatar
Yan committed
245
246
247
    help_menu = QtWidgets.QMenu('&Help', main_window)
    main_window.menuBar().addMenu(help_menu)
    help_menu.addAction(aboutact)
Yan's avatar
Yan committed
248
249
250
251

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

Yan's avatar
Yan committed
252
253
254
255
256
    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
257

Yan's avatar
Yan committed
258
259
260
    toolBar.addAction(openact)
    toolBar.addAction(exportact)
    toolBar.addSeparator()
Yan's avatar
Yan committed
261
262
263
    toolBar.addWidget(QtWidgets.QLabel("Predict Formula:"))
    toolBar.addWidget(predictform)
    toolBar.addSeparator()
Yan's avatar
Yan committed
264
265
    toolBar.addAction(zceact)
    toolBar.addAction(drlact)
Yan's avatar
Yan committed
266
267
    toolBar.addSeparator()
    toolBar.addAction(autozoomy)
Yan's avatar
Yan committed
268

Yan's avatar
Yan committed
269

3Yan's avatar
3Yan committed
270
    fileModel = QtWidgets.QFileSystemModel()
271
272
    fileModel.setRootPath('')
    activeDir = fileModel.index(config.value("open_folder"))
3Yan's avatar
3Yan committed
273
274
275
276
    treeview = QtWidgets.QTreeView()
    treeview.setModel(fileModel)
    treeview.setCurrentIndex(activeDir)
    treeview.expand(activeDir)
277
278
    dirview = QtWidgets.QListView()
    dirview.setModel(fileModel)
3Yan's avatar
3Yan committed
279
    treedock = QtWidgets.QDockWidget()
Yan's avatar
Yan committed
280
281
    if not config.value("view/filebrowservisible",type=bool):
        treedock.hide()
282
283
284
285
    treedock.setWidget(QtWidgets.QWidget())
    treedock_layout = QtWidgets.QVBoxLayout(treedock.widget())
    treedock_layout.addWidget(treeview)
    treedock_layout.addWidget(dirview)
3Yan's avatar
3Yan committed
286
287
    update.signal.connect(
        lambda: treeview.setCurrentIndex(fileModel.index(filename[0])))
288
    update.signal.connect(
3Yan's avatar
3Yan committed
289
290
        lambda: dirview.setRootIndex(fileModel.index(
            os.path.dirname(os.path.realpath(filename[0])))))
3Yan's avatar
3Yan committed
291
292
293
    treeview.doubleClicked.connect(lambda index: load_file(
        main_window, chromatogram, spectrum, ds, ms, filename, chrom,
        update, config, filename=fileModel.filePath(index)))
294
295
296
297
298
299
    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)))
300
301
302
303
304
305
    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
306
307
    main_window.addDockWidget(QtCore.Qt.LeftDockWidgetArea, treedock)

Yan's avatar
Yan committed
308
    layout = QtWidgets.QVBoxLayout(main_widget)
Yan's avatar
Yan committed
309
    layout.addWidget(toolBar)
Yan's avatar
Yan committed
310
311
312
    layout.addWidget(mpl_canvas)
    mpl_canvas.setFocus()

3Yan's avatar
3Yan committed
313

314
    mpl_canvas.mpl_connect('key_press_event', lambda event:
315
        update_spectrum(event.key == "f5", chromatogram, spectrum, ds, ms,
Yan's avatar
Yan committed
316
                        filename, chrom, config))
Yan's avatar
Yan committed
317
318
    main_window.dragEnterEvent = lambda event: drag_entered(event)
    main_window.dropEvent = lambda event: dropped(
319
        event, main_window, chromatogram, spectrum, ds, ms, filename,
3Yan's avatar
3Yan committed
320
        chrom, update, config)
Yan's avatar
Yan committed
321
    main_window.setAcceptDrops(True)
Yan's avatar
Yan committed
322

Yan's avatar
Yan committed
323
    if len(sys.argv) == 2:
324
        load_file(main_window, chromatogram, spectrum, ds, ms,
3Yan's avatar
3Yan committed
325
                  filename, chrom, update, config, filename=sys.argv[1])
Yan's avatar
Yan committed
326
327
328
329
330
331
    else:
        gt.pop_plot(spectrum, ms)
        gt.pop_plot(chromatogram, chrom)

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

333

Yan's avatar
Yan committed
334
335
if __name__ == "__main__":
    main()