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

Implmented new quality control function checkZ4S.

Browse Source
This commit is contained in:
Ludger Küperkoch 2015-07-01 15:30:13 +02:00
parent 489199af10
commit fb956f381a
1 changed files with 24 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -78,6 +78,8 @@ def run_autopicking(wfstream, pickparam):
minsiglength = pickparam.getParam('minsiglength')
minpercent = pickparam.getParam('minpercent')
nfacsl = pickparam.getParam('noisefactor')
# parameter to check for spuriously picked S onset
zfac = pickparam.getParam('zfac')
# initialize output
Pweight = 4 # weight for P onset
@ -162,6 +164,27 @@ def run_autopicking(wfstream, pickparam):
z_copy[0].data = tr_filt.data
Pflag = checksignallength(z_copy, aicpick.getpick(), tsnrz, minsiglength, \
nfacsl, minpercent, iplot)
if Pflag == 1:
# check for spuriously picked S onset
# both horizontal traces needed
if len(ndat) == 0 or len(edat) == 0:
print 'One or more horizontal components missing!'
print 'Skip control function checkZ4S.'
else:
# filter and taper horizontal traces
trH1_filt = edat.copy()
trH2_filt = ndat.copy()
trH1_filt.filter('bandpass', freqmin=bph1[0], freqmax=bph1[1], \
zerophase=False)
trH2_filt.filter('bandpass', freqmin=bph1[0], freqmax=bph1[1], \
zerophase=False)
trH1_filt.taper(max_percentage=0.05, type='hann')
trH2_filt.taper(max_percentage=0.05, type='hann')
zne = z_copy
zne += trH1_filt
zne += trH2_filt
Pflag = checkZ4S(zne, aicpick.getpick(), zfac, \
tsnrz[3], iplot)
##############################################################
# go on with processing if AIC onset passes quality control
if (aicpick.getSlope() >= minAICPslope and
@ -272,7 +295,7 @@ def run_autopicking(wfstream, pickparam):
h_copy[1].data = trH2_filt.data
elif algoS == 'AR3':
print zdat, edat, ndat
# re-create stream object including both horizontal components
# re-create stream object including all components
hdat = zdat.copy()
hdat += edat
hdat += ndat