From f46c6187482a8a072df30fade0061373c27405ed Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Ludger=20K=C3=BCperkoch?= Date: Tue, 20 Sep 2016 10:40:21 +0200 Subject: [PATCH] Stabilized zero-crosings determination for source spectrum estimation from P pulse. --- pylot/core/analysis/magnitude.py | 14 ++------------ 1 file changed, 2 insertions(+), 12 deletions(-) diff --git a/pylot/core/analysis/magnitude.py b/pylot/core/analysis/magnitude.py index 8635fbc9..40a3643f 100644 --- a/pylot/core/analysis/magnitude.py +++ b/pylot/core/analysis/magnitude.py @@ -387,21 +387,11 @@ def calcsourcespec(wfstream, onset, inventory, vp, delta, azimuth, incidence, qp ldat[0].stats.delta)) # get window after P pulse for - # 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 + # calculating source spectrum tstart = UTCDateTime(zdat[0].stats.starttime) tonset = onset.timestamp - tstart.timestamp impickP = tonset * zdat[0].stats.sampling_rate - wfzc = Ldat[impickP: impickP + winzc] + wfzc = Ldat[impickP: len(Ldat) - 1] # get time array t = np.arange(0, len(inttrz) * zdat[0].stats.delta, \ zdat[0].stats.delta)