Stabilized zero-crosings determination for source spectrum estimation from P pulse.

This commit is contained in:
Ludger Küperkoch 2016-09-20 10:40:21 +02:00
parent b3edeca81d
commit f46c618748

View File

@ -387,21 +387,11 @@ def calcsourcespec(wfstream, onset, inventory, vp, delta, azimuth, incidence, qp
ldat[0].stats.delta)) ldat[0].stats.delta))
# get window after P pulse for # get window after P pulse for
# calculating source spectrum # calculating source spectrum
if zdat[0].stats.sampling_rate <= 100:
winzc = zdat[0].stats.sampling_rate
elif zdat[0].stats.sampling_rate > 100 and \
zdat[0].stats.sampling_rate <= 200:
winzc = 0.5 * zdat[0].stats.sampling_rate
elif zdat[0].stats.sampling_rate > 200 and \
zdat[0].stats.sampling_rate <= 400:
winzc = 0.2 * zdat[0].stats.sampling_rate
elif zdat[0].stats.sampling_rate > 400:
winzc = zdat[0].stats.sampling_rate
tstart = UTCDateTime(zdat[0].stats.starttime) tstart = UTCDateTime(zdat[0].stats.starttime)
tonset = onset.timestamp - tstart.timestamp tonset = onset.timestamp - tstart.timestamp
impickP = tonset * zdat[0].stats.sampling_rate impickP = tonset * zdat[0].stats.sampling_rate
wfzc = Ldat[impickP: impickP + winzc] wfzc = Ldat[impickP: len(Ldat) - 1]
# get time array # get time array
t = np.arange(0, len(inttrz) * zdat[0].stats.delta, \ t = np.arange(0, len(inttrz) * zdat[0].stats.delta, \
zdat[0].stats.delta) zdat[0].stats.delta)