getSNR: demean using only mean determined from noise window.

2015-06-25 10:07:35 +02:00 · 2015-06-25 10:07:35 +02:00 · c9f07b6540
commit c9f07b6540
parent 2cd7572ee4
1 changed files with 2 additions and 2 deletions
--- a/pylot/core/pick/utils.py
+++ b/pylot/core/pick/utils.py
@ -334,8 +334,8 @@ def getSNR(X, TSNR, t1):
        print 'getSNR: Empty array isignal, check signal window!'
        return

-    # demean over entire snr window
-    x = x - np.mean(x[np.hstack([inoise, isignal])])
+    # demean over entire waveform
+    x = x - np.mean(x[inoise])

    # calculate ratios
    noiselevel = np.sqrt(np.mean(np.square(x[inoise])))