问题描述
我正在寻找这个紫色圆圈的位置,以及在传感器开始收集数据之前存在(主要是)水平滞后的不同数据集的类似情况。
有什么想法吗?
解决方法
numpy.diff 为您提供数组中每个连续项目之间的差异。它可能非常有用,但您可能需要添加一些平滑功能以防止噪声产生误报。这是一个易于使用的平滑函数:
def moving_average(x,w):
return np.convolve(x,np.ones(w),'valid') / w
(来自here)
我建议你绘制 np.diff(data) 并看看它是什么样的。在该图上找到紫色点会容易得多,如果您执行两次 np.diff 甚至可能更容易找到它。
对于这种情况,我会这样做:
purple_index = np.argmax(np.diff(data) < -1)