From 0d8b6b444663ce8b7747dc62670e83bed14fe972 Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Ludger=20K=C3=BCperkoch?= Date: Thu, 2 Jul 2015 09:26:05 +0200 Subject: [PATCH] New key in pick dictionary: marker, indicating if pick has been rejected by checksignallength or checkZ4S. --- pylot/core/pick/run_autopicking.py | 11 ++++++++++- 1 file changed, 10 insertions(+), 1 deletion(-) diff --git a/pylot/core/pick/run_autopicking.py b/pylot/core/pick/run_autopicking.py index 94ea1637..3295785d 100755 --- a/pylot/core/pick/run_autopicking.py +++ b/pylot/core/pick/run_autopicking.py @@ -102,6 +102,7 @@ def run_autopicking(wfstream, pickparam): aicPflag = 0 Pflag = 0 Sflag = 0 + Pmarker = [] # split components zdat = wfstream.select(component="Z") @@ -185,6 +186,12 @@ def run_autopicking(wfstream, pickparam): zne += trH2_filt Pflag = checkZ4S(zne, aicpick.getpick(), zfac, \ 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 if (aicpick.getSlope() >= minAICPslope and @@ -631,10 +638,12 @@ def run_autopicking(wfstream, pickparam): phasepick = {'lpp': lpickP, 'epp': epickP, 'mpp': mpickP, 'spe': Perror, \ 'snr': SNRP, 'snrdb': SNRPdB, 'weight': Pweight, 'fm': FM} picks = {phase: phasepick} + # add P marker + picks[phase]['marked'] = Pmarker # add S phase phase = 'S' phasepick = {'lpp': lpickS, 'epp': epickS, 'mpp': mpickS, 'spe': Serror, \ 'snr': SNRS, 'snrdb': SNRSdB, 'weight': Sweight, 'fm': None} picks[phase] = phasepick - + return picks