微信公众号搜"智元新知"关注
微信扫一扫可直接关注哦!

实时心电图实时散点图显示但混乱,如何清理?

如何解决实时心电图实时散点图显示但混乱,如何清理?

我有一些数据字典,我已经成功地使python实时更新并绘制。取得线路是可行的,但是现在我需要绘制散点图。我不确定如何设置实时更新数据,以便基本上可以相互绘制x2和y5,我习惯于将每个数据集分开,所以目前我有scatter2(x2)和scatter5(y5)分别更新,但是如何将它们像ax.scatter(x2,y5)而不是ax.scatter(x2)和另一行ax.scatter(y5)一起绘制?同时绘制它们,但分别绘制,使我得到了一个填充过度的散点图(也许它也保存了太多的值),并且似乎没有以与成功绘制的信号相同的速度移动。

import heartpy as hp
import numpy as np
import matplotlib.pyplot as plt

import matplotlib.animation as animation 
from collections import deque
from matplotlib.ticker import FuncFormatter
#from . import config 





sample_rate = 40
def init(): # initialization function: plot the background of each frame
   line.set_ydata([np.nan]*len(x))
   line1.set_ydata([np.nan]*len(x1))
   scatter2.set_xdata([np.nan]*len(x2))
   # scatter3.set_ydata([np.nan]*len(x3))
   # scatter4.set_ydata([np.nan]*len(x4))
   scatter5.set_ydata([np.nan]*len(y5))

   return line,line1,scatter2,scatter5,#scatter4,scatter5

def animate(i):
    read_fname = 'temp.1D'
    datax = np.loadtxt('temp.1D')
    wd,m = hp.process(datax,40 ) #run analysis        
    hr = wd['hr'] 
    rm = wd['rolling_mean'] 
    yb = wd['ybeat']
    pl = wd['peaklist']
    rp = wd['removed_beats']
    rp_y = wd['removed_beats_y']
    #pl=ax.collections[pl].get_array()

    yhr = np.array(hr[:])
    yrm = np.array(rm[:])
    yyb = np.array(yb[:])
    ypl = np.array(pl[:])
    yrp = np.array(rp[:])
    yrp_y = np.array(rp_y[:])
    yhr = np.concatenate([yhr])
    yrm = np.concatenate([yrm])
    yyb = np.concatenate([yyb])
    ypl = np.concatenate([ypl])
    yrp = np.concatenate([rp])
    yrp_y = np.concatenate([rp_y])

    data.extend(yhr)
    data1.extend(yrm) 
    data2.extend(ypl)
  #  data3.extend(yrp)
   # data4.extend(yrp_y)
    data5.extend(yyb)
    

    line.set_ydata(data)
    line1.set_ydata(data1)
    scatter2.set_xdata(data2)
   # scatter3.set_ydata(data3)
    # scatter4.set_ydata(data4)
    scatter5.set_ydata(data5)


    print(yhr)
    print(yrm)
    print(ypl)

    print()
    return line,scatter5

max_x = 250
max_y = 1700
#max_x2 = 1000
#max_y5 = 1000
data = deque(np.zeros(max_x),maxlen=max_x)
data1 = deque(np.zeros(max_x),maxlen=max_x)
data2= deque(np.zeros(max_x),maxlen=max_x)
#data3 = deque(np.zeros(max_x),maxlen=max_x)
#data4 = deque(np.zeros(max_x),maxlen=max_x) #hold the last 10 values
data5 = deque(np.zeros(max_y),maxlen=max_y)

x = np.arange(0,max_x)
x1 = np.arange(0,max_x)   
x2 = np.arange(0,max_x)
#x3 = np.arange(0,max_x)
#x4 = np.arange(0,max_x)
y5 = np.arange(0,max_y)

fig,ax = plt.subplots()
#fig,ax1 = plt.subplots()
ax.set_ylim(-500,max_y)
ax.set_xlim(0,max_x-1)
line,= ax.plot(x,np.linspace(0,1000,250,endpoint = False))
line1,= ax.plot(x1,endpoint = False))
#ax = plt.gca()
scatter2,= ax.plot(x2,marker="o",ls="",color = "green")
#scatter3,= ax.scatter(x3,50,1000))
#scatter4 = ax.scatter(x4,1000))
scatter5,= ax.plot(y5,ls="")

#ax.grid()
#ax.relim()
#ax.autoscale_view(True,True,True)

ax.xaxis.set_major_formatter(FuncFormatter(lambda x,pos: '{:.0}s'.format(max_x -x -990)))

plt.xlabel('Seconds ago')



ani = animation.FuncAnimation(fig,animate,init_func=init,interval=25,blit=True,save_count = 10)




plt.show()

版权声明:本文内容由互联网用户自发贡献,该文观点与技术仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请发送邮件至 dio@foxmail.com 举报,一经查实,本站将立刻删除。