[bugfix] determine the SNR according to the demeaned waveform between the beginning of the noise window and the end of the signal window

2015-06-23 08:24:21 +02:00 · 2015-06-23 08:24:21 +02:00 · dc281c6c0e
commit dc281c6c0e
parent bb6ebc7ee1
1 changed files with 3 additions and 0 deletions
--- a/pylot/core/pick/utils.py
+++ b/pylot/core/pick/utils.py
@ -322,6 +322,9 @@ def getSNR(X, TSNR, t1):
        print 'getSNR: Empty array isignal, check signal window!'
        return

+    # demean over entire snr window
+    x -= x[inoise[0][0]:isignal[0][-1]].mean()
+
    # calculate ratios
    noiselevel = np.sqrt(np.mean(np.square(x[inoise])))
    signallevel = np.sqrt(np.mean(np.square(x[isignal])))