Closedown
This commit is contained in:
@@ -22,30 +22,21 @@ def blm_remove_spikes(x):
|
||||
:param x: input array
|
||||
:return: output array
|
||||
"""
|
||||
print("0")
|
||||
x = copy.copy(x)
|
||||
if x.size > 5: # Must have enough sample points
|
||||
print ("A")
|
||||
d_x = x[1:] - x[:-1]
|
||||
maximum = x.max() - x.min()
|
||||
print ("B")
|
||||
for i in range(x.size-3):
|
||||
print ("I=", i)
|
||||
if d_x[i+1] > 0.5 * maximum and d_x[i+2] < -0.5 * maximum:
|
||||
print ("E")
|
||||
# 1 point spike
|
||||
x[i+2] = (x[i:i+5].sum() - x[i+2])/4
|
||||
print ("F")
|
||||
print ("p=", i)
|
||||
if d_x[i+1] > 0.5 * maximum and d_x[i+2] >= -0.5 * maximum:
|
||||
print ("c")
|
||||
# 2 point spikes
|
||||
x[i+2] = (x[i:i+2].sum() + x[i+4])/3
|
||||
print ("C")
|
||||
x[i+3] = (x[i+1] + x[i+4:i+6].sum())/3
|
||||
print ("D")
|
||||
if i < x(.size-4): #TODO: FIX BY AG, CHECK
|
||||
x[i+2] = (x[i:i+2].sum() + x[i+4])/3
|
||||
if i < x(.size-6): #TODO: FIX BY AG, CHECK
|
||||
x[i+3] = (x[i+1] + x[i+4:i+6].sum())/3
|
||||
# Handle edge points
|
||||
print ("G")
|
||||
if d_x[-1] > 0.5 * maximum:
|
||||
x[-1] = x[-3:-1].sum() / 2
|
||||
|
||||
|
||||
Reference in New Issue
Block a user