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

New key in pick dictionary: marker, indicating if pick has been rejected by checksignallength or checkZ4S.

Browse Source
This commit is contained in:
Ludger Küperkoch 2015-07-02 09:26:05 +02:00
parent 8282a7aa12
commit 0d8b6b4446
1 changed files with 10 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

9
pylot/core/pick/run_autopicking.py
View File

@ -102,6 +102,7 @@ def run_autopicking(wfstream, pickparam):
aicPflag = 0 aicPflag = 0
Pflag = 0 Pflag = 0
Sflag = 0 Sflag = 0
Pmarker = []
# split components # split components
zdat = wfstream.select(component="Z") zdat = wfstream.select(component="Z")
@ -185,6 +186,12 @@ def run_autopicking(wfstream, pickparam):
zne += trH2_filt zne += trH2_filt
Pflag = checkZ4S(zne, aicpick.getpick(), zfac, \ Pflag = checkZ4S(zne, aicpick.getpick(), zfac, \
tsnrz[3], iplot) tsnrz[3], iplot)
if Pflag == 0:
Pmarker = 'SinsteadP'
Pweight = 9
else:
Pmarker = 'shortsignallength'
Pweight = 9
############################################################## ##############################################################
# go on with processing if AIC onset passes quality control # go on with processing if AIC onset passes quality control
if (aicpick.getSlope() >= minAICPslope and if (aicpick.getSlope() >= minAICPslope and
@ -631,6 +638,8 @@ def run_autopicking(wfstream, pickparam):
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}
picks = {phase: phasepick} picks = {phase: phasepick}
# add P marker
picks[phase]['marked'] = Pmarker
# add S phase # add S phase
phase = 'S' phase = 'S'
phasepick = {'lpp': lpickS, 'epp': epickS, 'mpp': mpickS, 'spe': Serror, \ phasepick = {'lpp': lpickS, 'epp': epickS, 'mpp': mpickS, 'spe': Serror, \