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

Debugging.

Browse Source
This commit is contained in:
Ludger Küperkoch 2015-06-24 15:43:59 +02:00
parent 0cc29802a4
commit c68597ce62
1 changed files with 8 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -13,7 +13,7 @@ import matplotlib.pyplot as plt
import numpy as np import numpy as np
from pylot.core.pick.Picker import * from pylot.core.pick.Picker import *
from pylot.core.pick.CharFuns import * from pylot.core.pick.CharFuns import *
import pdb
def run_autopicking(wfstream, pickparam): def run_autopicking(wfstream, pickparam):
""" """
:param: wfstream :param: wfstream
@ -157,10 +157,11 @@ def run_autopicking(wfstream, pickparam):
# of class AutoPicking # of class AutoPicking
aicpick = AICPicker(aiccf, tsnrz, pickwinP, iplot, None, tsmoothP) aicpick = AICPicker(aiccf, tsnrz, pickwinP, iplot, None, tsmoothP)
############################################################## ##############################################################
# check signal length to detect spuriously picked noise peaks if aicpick.getpick() is not None:
z_copy[0].data = tr_filt.data # check signal length to detect spuriously picked noise peaks
Pflag = checksignallength(z_copy, aicpick.getpick(), tsnrz, minsiglength, \ z_copy[0].data = tr_filt.data
nfacsl, minpercent, iplot) Pflag = checksignallength(z_copy, aicpick.getpick(), tsnrz, minsiglength, \
nfacsl, minpercent, iplot)
############################################################## ##############################################################
# 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
@ -320,7 +321,8 @@ def run_autopicking(wfstream, pickparam):
############################################################### ###############################################################
# go on with processing if AIC onset passes quality control # go on with processing if AIC onset passes quality control
if (aicarhpick.getSlope() >= minAICSslope and if (aicarhpick.getSlope() >= minAICSslope and
aicarhpick.getSNR() >= minAICSSNR): aicarhpick.getSNR() >= minAICSSNR and
aicarhpick.getpick() is not None):
aicSflag = 1 aicSflag = 1
print 'AIC S-pick passes quality control: Slope: %f, SNR: %f' \ print 'AIC S-pick passes quality control: Slope: %f, SNR: %f' \
% (aicarhpick.getSlope(), aicarhpick.getSNR()) % (aicarhpick.getSlope(), aicarhpick.getSNR())