如何使用Python在这个实时图中进行峰值检测

时间:2016-11-01 03:36:48

标签: python python-3.x matplotlib plot real-time

这是简单的代码matplotlib。 (用于实时情节)(python 3.x)

在这个例子中,我如何应用峰值检测?

我想获得MAX峰值点指数和绘制点。实时

下面这样的图片

最初的实时情节图

original

我希望得到峰值点数

enter image description here

我怎么能?

这是示例代码 但是peakutils方法中的错误。

http://matplotlib.org/examples/animation/animate_decay.html

import numpy as np
import matplotlib.pyplot as plt
import matplotlib.animation as animation
import peakutils 
from peakutils.plot import plot as pplot


def data_gen(t=0):
    cnt = 0
    while cnt < 1000:
        cnt += 1
        t += 0.1
        yield t, np.sin(2*np.pi*t) * np.exp(-t/10.)


def init():
    ax.set_ylim(-1.1, 1.1)
    ax.set_xlim(0, 10)
    del xdata[:]
    del ydata[:]
    line.set_data(xdata, ydata)
    return line,

fig, ax = plt.subplots()
line, = ax.plot([], [], lw=2)
ax.grid()
xdata, ydata = [], []


def run(data):
    # update the data
    t, y = data
    xdata.append(t)
    ydata.append(y)
    indexes = peakutils.indexes(ydata, thres=0.5, min_dist=5)  #error line
    print(indexes) 
    xmin, xmax = ax.get_xlim()

    if t >= xmax:
        ax.set_xlim(xmin, 2*xmax)
        ax.figure.canvas.draw()
    line.set_data(xdata, ydata)

    return line,

ani = animation.FuncAnimation(fig, run, data_gen, blit=False, interval=10,
                              repeat=False, init_func=init)
plt.show()

0 个答案:

没有答案