marcel/pylot
3
0
Fork 0
Code Issues 17 Pull Requests 1 Releases Wiki Activity

Put additional parameters w0 and fc to picks dictionary.

Browse Source
This commit is contained in:
Ludger Küperkoch 2015-11-30 14:41:59 +01:00
parent 466e0020a6
commit 9f93c25aa8
1 changed files with 9 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

14
pylot/core/pick/autopick.py
View File

@ -18,7 +18,7 @@ from pylot.core.pick.utils import checksignallength, checkZ4S, earllatepicker,\
getSNR, fmpicker, checkPonsets, wadaticheck, crossings_nonzero_all getSNR, fmpicker, checkPonsets, wadaticheck, crossings_nonzero_all
from pylot.core.util.utils import getPatternLine from pylot.core.util.utils import getPatternLine
from pylot.core.read.data import Data from pylot.core.read.data import Data
from pylot.core.analysis.magnitude import WApp, DCfc from pylot.core.analysis.magnitude import WApp, w0fc
def autopickevent(data, param): def autopickevent(data, param):
stations = [] stations = []
@ -137,7 +137,9 @@ def autopickstation(wfstream, pickparam):
Pflag = 0 Pflag = 0
Sflag = 0 Sflag = 0
Pmarker = [] Pmarker = []
Ao = None Ao = None # Wood-Anderson peak-to-peak amplitude
w0 = None # plateau of source spectrum
fc = None # corner frequancy of source spectrum
# split components # split components
zdat = wfstream.select(component="Z") zdat = wfstream.select(component="Z")
@ -313,7 +315,7 @@ def autopickstation(wfstream, pickparam):
FM = 'N' FM = 'N'
############################################################## ##############################################################
# get DC value (w0) and corner frequency (fc) of source spectrum # get DC value and corner frequency (fc) of source spectrum
# from P pulse # from P pulse
# initialize Data object # initialize Data object
data = Data() data = Data()
@ -339,7 +341,7 @@ def autopickstation(wfstream, pickparam):
index = min([3, len(zc) - 1]) index = min([3, len(zc) - 1])
calcwin = (zc[index] - zc[0]) * z_copy[0].stats.delta calcwin = (zc[index] - zc[0]) * z_copy[0].stats.delta
# calculate source spectrum and get w0 and fc # calculate source spectrum and get w0 and fc
specpara = DCfc(z_copy, mpickP, calcwin, iplot) specpara = w0fc(z_copy, mpickP, calcwin, iplot)
w0 = specpara.getw0() w0 = specpara.getw0()
fc = specpara.getfc() fc = specpara.getfc()
@ -788,7 +790,8 @@ def autopickstation(wfstream, pickparam):
# for P phase # for P phase
phase = 'P' phase = 'P'
phasepick = {'lpp': lpickP, 'epp': epickP, 'mpp': mpickP, 'spe': Perror, phasepick = {'lpp': lpickP, 'epp': epickP, 'mpp': mpickP, 'spe': Perror,
'snr': SNRP, 'snrdb': SNRPdB, 'weight': Pweight, 'fm': FM} 'snr': SNRP, 'snrdb': SNRPdB, 'weight': Pweight, 'fm': FM,
'w0': w0, 'fc': fc}
picks = {phase: phasepick} picks = {phase: phasepick}
# add P marker # add P marker
picks[phase]['marked'] = Pmarker picks[phase]['marked'] = Pmarker
@ -800,6 +803,7 @@ def autopickstation(wfstream, pickparam):
# add Wood-Anderson amplitude # add Wood-Anderson amplitude
picks[phase]['Ao'] = Ao picks[phase]['Ao'] = Ao
return picks return picks