From fe8128a61a0eeb4bfd1c65ace042b06379089dae Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Ludger=20K=C3=BCperkoch?= Date: Tue, 15 Aug 2017 16:57:52 +0200 Subject: [PATCH] Use mean for noise level determination instead of max. --- pylot/core/pick/utils.py | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/pylot/core/pick/utils.py b/pylot/core/pick/utils.py index c3d8c103..f622cfa3 100644 --- a/pylot/core/pick/utils.py +++ b/pylot/core/pick/utils.py @@ -476,7 +476,7 @@ def getsignalwin(t, t1, tsignal): ''' # get signal window - isignal, = np.where((t <= min([t1 + tsignal, t[-1]])) \ + isignal, = np.where((t <= min([t1 + tsignal, len(t)])) \ & (t >= t1)) if np.size(isignal) < 1: print("getsignalwin: Empty array isignal, check signal window!") @@ -754,7 +754,7 @@ def checksignallength(X, pick, TSNR, minsiglength, nfac, minpercent, iplot=0, fi # get signal window isignal = getsignalwin(t, pick, minsiglength) # calculate minimum adjusted signal level - minsiglevel = max(rms[inoise]) * nfac + minsiglevel = np.mean(rms[inoise]) * nfac # minimum adjusted number of samples over minimum signal level minnum = len(isignal) * minpercent / 100 # get number of samples above minimum adjusted signal level