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

renaming run_autopicking to autopickstation

Browse Source
This commit is contained in:
Sebastian Wehling-Benatelli 2015-07-13 06:40:51 +02:00
parent d79165a9f2
commit da22dcce90
1 changed files with 10 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

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

@ -35,7 +35,7 @@ def autopickevent(data, param):
for station in stations: for station in stations:
topick = data.select(station=station) topick = data.select(station=station)
all_onsets[station] = run_autopicking(topick, param) all_onsets[station] = autopickstation(topick, param)
# quality control # quality control
# median check and jackknife on P-onset times # median check and jackknife on P-onset times
@ -43,7 +43,7 @@ def autopickevent(data, param):
# check S-P times (Wadati) # check S-P times (Wadati)
return wadaticheck(jk_checked_onsets, wdttolerance, iplot) return wadaticheck(jk_checked_onsets, wdttolerance, iplot)
def run_autopicking(wfstream, pickparam): def autopickstation(wfstream, pickparam):
""" """
:param: wfstream :param: wfstream
:type: `~obspy.core.stream.Stream` :type: `~obspy.core.stream.Stream`
@ -144,7 +144,7 @@ def run_autopicking(wfstream, pickparam):
if algoP == 'HOS' or algoP == 'ARZ' and zdat is not None: if algoP == 'HOS' or algoP == 'ARZ' and zdat is not None:
print '##########################################' print '##########################################'
print 'run_autopicking: Working on P onset of station %s' % zdat[ print 'autopickstation: Working on P onset of station %s' % zdat[
0].stats.station 0].stats.station
print 'Filtering vertical trace ...' print 'Filtering vertical trace ...'
print zdat print zdat
@ -161,7 +161,7 @@ def run_autopicking(wfstream, pickparam):
Lwf = zdat[0].stats.endtime - zdat[0].stats.starttime Lwf = zdat[0].stats.endtime - zdat[0].stats.starttime
Ldiff = Lwf - Lc Ldiff = Lwf - Lc
if Ldiff < 0: if Ldiff < 0:
print 'run_autopicking: Cutting times are too large for actual ' \ print 'autopickstation: Cutting times are too large for actual ' \
'waveform!' 'waveform!'
print 'Use entire waveform instead!' print 'Use entire waveform instead!'
pstart = 0 pstart = 0
@ -234,7 +234,7 @@ def run_autopicking(wfstream, pickparam):
(aicpick.getSlope(), aicpick.getSNR()) (aicpick.getSlope(), aicpick.getSNR())
print 'Go on with refined picking ...' print 'Go on with refined picking ...'
# re-filter waveform with larger bandpass # re-filter waveform with larger bandpass
print 'run_autopicking: re-filtering vertical trace ...' print 'autopickstation: re-filtering vertical trace ...'
z_copy = zdat.copy() z_copy = zdat.copy()
tr_filt = zdat[0].copy() tr_filt = zdat[0].copy()
tr_filt.filter('bandpass', freqmin=bpz2[0], freqmax=bpz2[1], tr_filt.filter('bandpass', freqmin=bpz2[0], freqmax=bpz2[1],
@ -292,7 +292,7 @@ def run_autopicking(wfstream, pickparam):
else: else:
FM = 'N' FM = 'N'
print 'run_autopicking: P-weight: %d, SNR: %f, SNR[dB]: %f, ' \ print 'autopickstation: P-weight: %d, SNR: %f, SNR[dB]: %f, ' \
'Polarity: %s' % (Pweight, SNRP, SNRPdB, FM) 'Polarity: %s' % (Pweight, SNRP, SNRPdB, FM)
Sflag = 1 Sflag = 1
@ -302,7 +302,7 @@ def run_autopicking(wfstream, pickparam):
Sflag = 0 Sflag = 0
else: else:
print 'run_autopicking: No vertical component data availabler!, ' \ print 'autopickstation: No vertical component data availabler!, ' \
'Skip station!' 'Skip station!'
if edat is not None and ndat is not None and len(edat) > 0 and len( if edat is not None and ndat is not None and len(edat) > 0 and len(
@ -395,7 +395,7 @@ def run_autopicking(wfstream, pickparam):
cuttimesh2 = [round(aicarhpick.getpick() - Srecalcwin), cuttimesh2 = [round(aicarhpick.getpick() - Srecalcwin),
round(aicarhpick.getpick() + Srecalcwin)] round(aicarhpick.getpick() + Srecalcwin)]
# re-filter waveform with larger bandpass # re-filter waveform with larger bandpass
print 'run_autopicking: re-filtering horizontal traces...' print 'autopickstation: re-filtering horizontal traces...'
h_copy = hdat.copy() h_copy = hdat.copy()
# filter and taper data # filter and taper data
if algoS == 'ARH': if algoS == 'ARH':
@ -498,7 +498,7 @@ def run_autopicking(wfstream, pickparam):
elif Serror > timeerrorsS[3]: elif Serror > timeerrorsS[3]:
Sweight = 4 Sweight = 4
print 'run_autopicking: S-weight: %d, SNR: %f, SNR[dB]: %f' % ( print 'autopickstation: S-weight: %d, SNR: %f, SNR[dB]: %f' % (
Sweight, SNRS, SNRSdB) Sweight, SNRS, SNRSdB)
else: else:
@ -507,7 +507,7 @@ def run_autopicking(wfstream, pickparam):
'AIC-Slope=', aicarhpick.getSlope() 'AIC-Slope=', aicarhpick.getSlope()
else: else:
print 'run_autopicking: No horizontal component data available or ' \ print 'autopickstation: No horizontal component data available or ' \
'bad P onset, skipping S picking!' 'bad P onset, skipping S picking!'
############################################################## ##############################################################