2014-03-13 13:27:34 +01:00
|
|
|
#!/usr/bin/env python
|
2014-03-31 12:57:08 +02:00
|
|
|
# -*- coding: utf-8 -*-
|
|
|
|
"""
|
|
|
|
PyLoT: Main program
|
|
|
|
===================
|
|
|
|
PyLoT is a seismic data processing software capable of picking seismic
|
|
|
|
phases (symmetric and asymmetric error assignment), exporting these to
|
|
|
|
several common phase data formats and post process the data, e.g. locating
|
|
|
|
events, via external localization software.
|
|
|
|
Additionally PyLoT is meant as an interface to autoPyLoT which can
|
|
|
|
automatically pick seismic phases, if the parameters have properly been
|
|
|
|
chosen for the particular data set.
|
|
|
|
|
|
|
|
:author:
|
|
|
|
Sebastian Wehling-Benatelli
|
|
|
|
:copyright:
|
|
|
|
The PyLoT Development Team (https://ariadne.geophysik.rub.de/trac/PyLoT)
|
|
|
|
:license:
|
|
|
|
GNU Lesser General Public License, Version 3
|
|
|
|
(http://www.gnu.org/copyleft/lesser.html)
|
|
|
|
"""
|
2014-01-09 10:43:40 +01:00
|
|
|
|
|
|
|
import os
|
|
|
|
import platform
|
|
|
|
import sys
|
2014-03-13 13:27:34 +01:00
|
|
|
from PySide.QtCore import *
|
|
|
|
from PySide.QtGui import *
|
2014-04-01 16:46:10 +02:00
|
|
|
from obspy.core import (read, UTCDateTime)
|
2014-11-06 15:07:05 +01:00
|
|
|
from pylot import *
|
2014-01-09 10:43:40 +01:00
|
|
|
from pylot.core.util import _getVersionString
|
2014-07-16 12:06:05 +02:00
|
|
|
from pylot.core.read import (Data,
|
|
|
|
FilterOptions)
|
2014-03-13 13:27:34 +01:00
|
|
|
from pylot.core.util import FILTERDEFAULTS
|
2014-03-31 12:57:08 +02:00
|
|
|
from pylot.core.util import checkurl
|
2014-06-10 16:37:07 +02:00
|
|
|
from pylot.core.util import (PickDlg,
|
2014-10-27 11:56:44 +01:00
|
|
|
FilterOptionsDialog,
|
2014-06-10 16:37:07 +02:00
|
|
|
PropertiesDlg,
|
|
|
|
MPLWidget,
|
|
|
|
HelpForm)
|
2014-11-13 11:27:52 +01:00
|
|
|
from pylot.core.util import layoutStationButtons
|
2014-11-26 08:46:16 +01:00
|
|
|
import qrc_resources
|
2014-01-09 10:43:40 +01:00
|
|
|
|
|
|
|
# Version information
|
|
|
|
__version__ = _getVersionString()
|
|
|
|
|
2014-03-13 13:27:34 +01:00
|
|
|
|
2014-01-10 05:45:03 +01:00
|
|
|
class MainWindow(QMainWindow):
|
2014-03-13 13:27:34 +01:00
|
|
|
|
2014-01-10 05:45:03 +01:00
|
|
|
def __init__(self, parent=None):
|
|
|
|
super(MainWindow, self).__init__(parent)
|
2014-03-13 13:27:34 +01:00
|
|
|
|
2014-11-06 15:07:05 +01:00
|
|
|
# initialize filter parameter
|
2014-03-13 13:27:34 +01:00
|
|
|
filterOptionsP = FILTERDEFAULTS['P']
|
|
|
|
filterOptionsS = FILTERDEFAULTS['S']
|
|
|
|
self.filterOptionsP = FilterOptions(**filterOptionsP)
|
|
|
|
self.filterOptionsS = FilterOptions(**filterOptionsS)
|
|
|
|
|
2014-11-06 15:07:05 +01:00
|
|
|
# initialize data
|
|
|
|
self.data = None
|
2014-03-28 05:25:46 +01:00
|
|
|
self.loadData()
|
|
|
|
self.updateFilterOptions()
|
2014-11-26 08:46:16 +01:00
|
|
|
try:
|
|
|
|
self.startTime = min([tr.stats.starttime for tr in self.data.wfdata])
|
|
|
|
except:
|
|
|
|
self.startTime = UTCDateTime()
|
2014-03-28 05:25:46 +01:00
|
|
|
|
|
|
|
self.setupUi()
|
|
|
|
|
2014-10-27 11:57:34 +01:00
|
|
|
def _getCurrentPlotType(self):
|
|
|
|
return 'TestType'
|
|
|
|
|
2014-03-28 22:26:15 +01:00
|
|
|
def loadData(self):
|
2014-11-26 08:46:16 +01:00
|
|
|
self.data = None
|
|
|
|
|
|
|
|
def getComponent(self):
|
|
|
|
return self
|
2014-03-28 22:26:15 +01:00
|
|
|
|
2014-11-13 11:27:52 +01:00
|
|
|
def getData(self):
|
|
|
|
return self.data
|
|
|
|
|
|
|
|
def getDataWidget(self):
|
|
|
|
return self.DataPlot
|
|
|
|
|
2014-03-28 05:25:46 +01:00
|
|
|
def setupUi(self):
|
2014-11-06 15:07:05 +01:00
|
|
|
self.setWindowIcon(QIcon(":/pylot.ico"))
|
2014-03-28 05:25:46 +01:00
|
|
|
|
2014-11-06 15:07:05 +01:00
|
|
|
xlab = self.startTime.strftime('seconds since %d %b %Y %H:%M:%S (%Z)')
|
2014-10-27 11:55:06 +01:00
|
|
|
plottitle = self._getCurrentPlotType()
|
2014-07-16 12:06:05 +02:00
|
|
|
|
2014-03-28 05:25:46 +01:00
|
|
|
# create central matplotlib figure widget
|
2014-07-16 12:06:05 +02:00
|
|
|
self.DataPlot = MPLWidget(parent=self,
|
|
|
|
xlabel=xlab,
|
2014-10-27 11:55:06 +01:00
|
|
|
ylabel=None,
|
|
|
|
title=plottitle)
|
2014-03-28 05:25:46 +01:00
|
|
|
|
2014-04-01 16:46:10 +02:00
|
|
|
self.eventLabel = QLabel()
|
2014-06-11 05:48:21 +02:00
|
|
|
self.eventLabel.setFrameStyle(QFrame.StyledPanel | QFrame.Sunken)
|
2014-04-01 16:46:10 +02:00
|
|
|
status = self.statusBar()
|
|
|
|
status.setSizeGripEnabled(False)
|
|
|
|
status.addPermanentWidget(self.eventLabel)
|
|
|
|
status.showMessage("Ready", 5000)
|
|
|
|
|
2014-11-13 11:27:52 +01:00
|
|
|
statLayout = layoutStationButtons(self.getData(), self.getComponent())
|
|
|
|
dataLayout = self.getDataWidget()
|
2014-03-28 05:25:46 +01:00
|
|
|
|
|
|
|
maingrid = QGridLayout()
|
|
|
|
maingrid.setSpacing(10)
|
|
|
|
maingrid.addLayout(statLayout, 0, 0)
|
2014-11-26 08:46:16 +01:00
|
|
|
maingrid.addWidget(dataLayout, 1, 0)
|
|
|
|
self.setLayout(maingrid)
|
|
|
|
#self.setCentralWidget(dataLayout)
|
2014-03-28 05:25:46 +01:00
|
|
|
|
2014-07-16 12:06:05 +02:00
|
|
|
def plotData(self):
|
2014-11-26 08:46:16 +01:00
|
|
|
pass #self.data.plotData(self.DataPlot)
|
2014-03-19 12:24:41 +01:00
|
|
|
|
2014-11-06 15:07:05 +01:00
|
|
|
def adjustFilterOptions(self):
|
2014-11-13 11:27:52 +01:00
|
|
|
fstring = "Filter Options ({0})".format(self.getSeismicPhase())
|
|
|
|
filterDlg = FilterOptionsDialog(titleString=fstring,
|
2014-11-06 15:07:05 +01:00
|
|
|
parent=self,
|
2014-11-13 11:27:52 +01:00
|
|
|
filterOptions=self.getFilterOptions())
|
|
|
|
filterDlg.accepted.connect(filterDlg.getFilterOptions)
|
|
|
|
|
|
|
|
def getFilterOptions(self):
|
|
|
|
return self.filteroptions
|
|
|
|
|
|
|
|
def setFilterOptions(self, filterOptions):
|
|
|
|
cases = {'P':self.filterOptionsP,
|
|
|
|
'S':self.filterOptionsS}
|
|
|
|
cases[self.getSeismicPhase()] = filterOptions
|
|
|
|
self.updateFilterOptions()
|
2014-11-06 15:07:05 +01:00
|
|
|
|
2014-03-28 05:25:46 +01:00
|
|
|
def updateFilterOptions(self):
|
2014-10-27 12:04:01 +01:00
|
|
|
try:
|
|
|
|
self.filteroptions = [self.filterOptionsP
|
|
|
|
if not self.seismicPhase == 'S'
|
|
|
|
else self.filterOptionsS]
|
2014-11-06 15:07:05 +01:00
|
|
|
except Exception, e:
|
2014-10-27 12:04:01 +01:00
|
|
|
self.updateStatus('Error: %s' % e + ' ... no filteroptions loaded')
|
|
|
|
else:
|
|
|
|
self.updateStatus('Filteroptions succesfully loaded ...')
|
2014-03-28 22:26:15 +01:00
|
|
|
|
2014-11-13 11:27:52 +01:00
|
|
|
def getSeismicPhase(self):
|
|
|
|
return self.seismicPhase
|
|
|
|
|
|
|
|
def setSeismicPhase(self, phase):
|
|
|
|
self.seismicPhase = self.seismicPhaseButton.getValue()
|
|
|
|
|
2014-04-01 16:46:10 +02:00
|
|
|
def updateStatus(self, message):
|
|
|
|
self.statusBar().showMessage(message, 5000)
|
|
|
|
|
2014-03-28 22:26:15 +01:00
|
|
|
def helpHelp(self):
|
2014-06-11 15:19:37 +02:00
|
|
|
if checkurl():
|
2014-07-03 10:18:42 +02:00
|
|
|
form = HelpForm('https://ariadne.geophysik.ruhr-uni-bochum.de/trac/PyLoT/wiki')
|
2014-06-11 15:19:37 +02:00
|
|
|
else:
|
|
|
|
form = HelpForm(':/help.html')
|
|
|
|
form.show()
|
2014-03-28 22:26:15 +01:00
|
|
|
|
2014-01-10 05:45:03 +01:00
|
|
|
|
2014-03-28 22:26:15 +01:00
|
|
|
def main():
|
2014-11-06 15:07:05 +01:00
|
|
|
# create the Qt application
|
2014-01-30 13:11:44 +01:00
|
|
|
pylot_app = QApplication(sys.argv)
|
2014-01-09 10:43:40 +01:00
|
|
|
|
2014-03-28 22:26:15 +01:00
|
|
|
# set Application Information
|
|
|
|
pylot_app.setOrganizationName("Ruhr-University Bochum / MAGS2")
|
|
|
|
pylot_app.setOrganizationDomain("rub.de")
|
|
|
|
pylot_app.setApplicationName("PyLoT")
|
|
|
|
pylot_app.setWindowIcon(QIcon(":/pylot.ico"))
|
|
|
|
|
|
|
|
# create the main window
|
|
|
|
pylot_form = MainWindow()
|
2014-01-09 10:43:40 +01:00
|
|
|
|
2014-03-13 13:27:34 +01:00
|
|
|
# Show main window and run the app
|
2014-03-28 22:26:15 +01:00
|
|
|
pylot_form.show()
|
2014-02-03 12:51:23 +01:00
|
|
|
pylot_app.exec_()
|
2014-03-28 22:26:15 +01:00
|
|
|
|
2014-04-11 19:39:29 +02:00
|
|
|
if __name__ == "__main__":
|
|
|
|
main()
|