diff --git a/saqc/funcs/spikes_detection.py b/saqc/funcs/spikes_detection.py index ad55bebfa8933a61ac2798eefe8ce80ec5889d21..601d4dfbf78da2f63010c05f8c833ccde71ddad2 100644 --- a/saqc/funcs/spikes_detection.py +++ b/saqc/funcs/spikes_detection.py @@ -942,6 +942,7 @@ def spikes_flagBasic(data, field, flagger, thresh, tolerance, window, numba_kick # define spike testing function to roll with: def spikeTester(chunk, thresh=thresh, tol=tolerance): # signum change!!! + # chunk_stair=(np.sign(chunk[-2] - chunk[-1])*(chunk - chunk[-1]) < thresh)[::-1].cumsum() chunk_stair = (np.abs(chunk - chunk[-1]) < thresh)[::-1].cumsum() initial = np.searchsorted(chunk_stair, 2) if initial == len(chunk): @@ -960,7 +961,7 @@ def spikes_flagBasic(data, field, flagger, thresh, tolerance, window, numba_kick result = customRolling(to_roll, window, spikeTester, roll_mask, closed='both', engine=engine) # correct the result: only those values define plateaus, that do not have - # values at theire left starting point, that belong to other plateaus themselfs: + # values at their left starting point, that belong to other plateaus themself: def correctResult(result): for k in range(len(result)): if result[k] > 0: