本发明涉及一种心电检测系统,尤其涉及一种基于通用平台的心电采集系统及其采集方法。
背景技术
心血管病患人数众多,而且每10秒就有1人死于心血管病,但心血管病的诊治情况却不容乐观,由于很多心脏疾病非常需要心电图数据来指导诊断,并且在疾病发作时才可以收集到相关性高的心电图数据,如果没有随时随地的收集处理设备,将会为诊断带来巨大麻烦。
现有技术主要有如下类别:从平台分,有专用平台设备(如嵌入式系统)和与通用平台(android或苹果);从测量时间分,有24小时不间断监护和按需测量;从采集模块与数据处理模块关联分,有分离模式与一体模式;从市场定位分有医院专业设备与家庭日常保健设备,以及非便携与便携的区分。一般来说,医院专业设备采用专用平台,采用一体模式,开发成本高,升级难度大;而家用日常保健设备更多选择通用平台,采集模块独立设备,采用蓝牙进行通信。从检测时间看,家用设备不间断监护与按需测量都有其市场需求,而按需测量应用相比更广泛。
在通用平台(android或苹果)上开发的家用便携保健设备,因为设备的非专用特点,用户打开检测软件,并点击选取心电测量功能,往往显得操作复杂,对于行动不便者,尤其老人,是其严重的缺陷。
同时心电模块工作时,需要一定的电能,即便是待机,电能消耗也不容忽视。传统的心电独立模块,通常采用手动电源开关来关闭电源,以达到节电目的。此外设备分离,在保管上也有其不便之处。
因此,亟待解决上述问题。
技术实现要素:
发明目的:本发明的目的是提供一种可实现心电图的实时快速采集、解析、绘制、传输、显示,且同时可确保待检测心电信号的可用性的基于通用平台的心电采集系统。
本发明的另一个的目的是提供该基于通用平台的心电采集系统的采集方法。
技术方案:为实现以上目的,本发明公开了一种基于通用平台的心电采集系统,包括集成于android平台内的android主控板、心电采集模块和传感器,所述心电采集模块与多个传感器相连接,所述传感器分别置于人体不同部位以实现心电信号测量,传感器的输出经心电采集模块处理后形成心电数据完成采集,并经串口通信将其发送至android主控板;
所述android主控板用于运行android应用程序,并对心电数据进行处理和储存,android主控板通过串口采集到心电数据后,多线程同步处理:
1.1、进行心电信号质量辨识;
1.2、根据辨识后的心电信号进行心电分析,计算得到心电参数、心率、qrs时间和st高度参数;
1.3、本地数据保存:包括原始心电数据保存,计算得到的心电参数和心电结果保存;
1.4、网络服务器数据保存:本地数据同步上传服务器;
android主控板连接有显示屏,并通过显示屏即时将采集的心电图波形及心电参数数据显示出。
优选的,所述步骤1.1中,所述心电信号质量辨识方法包括如下步骤:
s1、获取采样率为fs的原始的心电波形幅值数据,并将心电波形幅值数据分成以n秒为时间单位的处理段x,n为8~120间的任意一个值;
s2、在处理段x内,将心电信号以1秒为单位进行片段分割,分割段为xi,其中i=1…n,xi的数据长度为m=1×fs=fs;
s3、以处理段x为单位,取出每一心电波形片段xi的幅值最大值mmax和最小值mmin构成包络点,通过比较各个包络点得到分割段xi的包络差ei=mmax–mmin,再求取处理段x内波形包络差的平均值
s4、在处理段x内,获取每一分割段xi的幅值方差值
s5、在处理段x内,利用快速傅里叶变换,将心电片段信号xi的时域信号转化为频域信号,即为功率谱信号;对1~5hz的幅值进行积分得到功率
s6、在处理段x内,得到步骤s5中的心电片段信号xi的功率谱信号后,对5~40hz的幅值进行积分得到功率
s7、根据参数包络差、方差值、信号在1~5hz之间的功率和信噪比的合格情况对心电波形进行质量等级划分,心电波形的质量分成4个等级,分别用0~3表示质量等级系数;质量等级系数0表示心电波形良好,质量等级系数1表示心电波形较差,质量等级系数2表示心电波形差,质量等级系数3表示心电波形存在严重噪声。即从波形质量的程度“0”表示良好,“1”表示较差,“2”表示差,“3”表示严重噪声;所述心电波形的质量等级划分方法为:当信噪比参数不合格时,该处理段心电波形的质量等级系数直接评定为3,即波形存在严重噪声;在信噪比参数合格的前提下,根据包络差、方差值、信号在1~5hz之间的功率的合格情况来评定心电波形质量,当三个参数全部合格时,该处理段心电波形质量等级系数评定为0,即波形良好;当三个参数中有两个或一个合格时,心电波形质量等级系数评定为1,即波形较差;当三个参数全部不合格时,心电波形质量等级系数评定为2,即波形差;
进一步,所述步骤1.2中心电分析方法,包括如下步骤:
s11、获取n秒的待分析的心电波形,n为8~120间的任意一个值;
s12、依据心电信号质量辨识方法对待分析的心电波形进行质量评估,得到n秒的心电波形质量等级,其中心电波形质量评估以1s为单位,共得到n个数据分别代表第1秒到第n秒的质量等级系数;
s13、当心电信号质量等级系数为3且持续时间大于0.3n秒时,提示“信号质量噪声过大,无法分析”返回至步骤a重新获取待分析的心电波形;否则进行心电参数的计算和心电结果判断;心电参数的计算和心电结果判断的步骤为:
s1301、依据心电qrs波群宽度的检测方法找出qrs波群宽度、qrs波群中各个子波形的起止时间点、幅度参数和代表qrs波群形态的bn;然后再找出心电特征包括p、t波的起止时间点及其幅度参数,构成了一次心跳;
s1302、取心电波形质量等级系数小于3的心电特征用于计算心电参数rr间期和心率,用心电波形质量等级为0的心电特征计算心电参数pr间期、qrs宽度、qt间期、st段高度、p波平均电轴、qrs平均电轴、t波平均电轴、q波宽度和r波高度;
s1303、利用心电特征和上述心电参数进行心电异常判断得出心电图异常结果
s14、当心电信号质量等级系数为1或2且持续时间大于0.3n秒时,提示“信号质量较差”并输出计算得到的心电参数和心电结果;否则认为心电波形质量较好,直接输出计算得到的心电参数和心电结果。
优选的,所述步骤s1301中心电qrs波群宽度的检测方法,包括如下步骤:
a、获取采样率为fs的原始心电波形幅值数据s,并将心电波形s分成以n秒为时间单位的处理段信号x,n为8~120间的任意一值;
b、将信号x经过包络处理得到包络信号xe,然后在包络信号xe上利用阈值法定位qrs波群的位置得到集合bxe;
其中信号x的处理方法和利用阀值法的定位方法包括如下步骤:
b1、将信号x经过一个带通滤波器,其频率范围为5hz~25hz,得到凸显出qrs波群特征且滤除干扰波的信号xfiltered;对信号xfiltered做微分后再平方得到信号xdiff,对xdiff进行积分得到包络信号xe;
b2、将包络信号xe按0.2n秒的长度分段得到5个包络片段xe(i),其中i=1…5,在每个包络片段内找到最大值xemax(i),然后求出这5个最大值的中值
c、在集合bxe中选取第n个qrs波群位置bxe(n),求出该qrs波群的rr间期rrxe(n)=bxe(n)-bxe(n-1);从信号x中截取一片段xios,截取的起始位置为bxe(n-1)+rrxe(n)×k3,终止位置为bxe(n)-rrxe(n)×(0.5-k3),其中0.3<k3<0.5,即第n-1个心跳的t波之后到第n个心跳的p波之前的一段数据;作出片段xios中幅值范围的直方图gp,得到xios数据取值的统计分布;在直方图gp中找出xios幅值分布个数最多的子区间b,该子区间b代表的电位即为第n次心跳之前的等电位段电位vios;
d、从信号x中第n个qrs波群位置bxe(n)周围截取一数据段qn,在数据段qn中找到所有极大值和极小值点的集合为极值点p和qrs波群中绝对幅值最大的波qm;
其中获取极值点p和qm的方法包括如下步骤:
d1、信号x在bxe(n)处往前ω-=(0.07~0.15)fs、往后ω+=(0.1~0.23)fs个点,截取信号x中的qrs波群的数据段qn=(q1,…,qj,…,qω),其中ω=ω-+ω+表示qn的数据长度,qj表示qn中第j个数据;
d2、在数据段qn中找到所有极大值pk和极小值点lo,极大值pk和极小值点lo合称为极值点p;然后再找出极大值pk中的最大值pmax=(vmax,imax)和极小值lo中的最小值lmin=(vmin,imin),i、v分别表示点在qn中的位置和数据幅值大小;
d3、定义qrs波群中最明显的波顶点为qrs波群中绝对幅值最大的波qm,如果|vmax-vios|>k4×|vmin-vios|则qm点为lmin点,否则qm点为pmax点,其中2<k4<10;
e、极值点p中满足最小qrs可识别波条件的为主极点dn,再从主极点dn中筛选出特征点
其中筛选特征点
e1、从极值点p中选择一个极值点pj,用以下公式在处理段qn中搜索极值点pj的左右支撑区间
δqj,x=|qj-qx|
其中j和qj为极值点pj在qn中的位置和幅值,τ为qrs最大的有意义的生理时间宽度,取80~160ms;i、a、b和k均是为了求解支撑区间的辅助变量,没有特别含义;qx是qn中第x个数据,与qj一个意义;
e2、若极值点pj的左右支撑区间
e3、重复步骤e1、e2直到筛选出极值点p中所有可识别波,记为主极值点集合dn={pj};
e4、极值点pj的左斜率
e5、如
f、在数据段qn中,以步骤d中得到的qm为顶点,利用三角形法确定qrs的起止点qrson和qrsoff;所述三角形法包括如下步骤:
f1、三角形法的作用是寻找波形上的转折点,转折点包括极值点和非极值拐点;一心电波形线段,从顶点x开始向前取一个点z做一条辅助线l,求该线段上所有数据点到直线l的距离d,距离最大的点即为转折点y;使用三角形法找qrson和qrsoff的前提是:qrson和qrsoff是波形的非极值拐点,而不是峰值点;
f2、在数据段数据qn中首先以qm点作为顶点,向前50~200ms之间寻找一点作辅助线s1,寻找数据段数据qn中到辅助线s1的距离最大的点,即为拐点1,判断拐点1是否在
f3、在数据段数据qn中首先以qm点作为顶点,向后50~200ms之间寻找一点作辅助线l1,寻找数据段数据qn中到辅助线l1的距离最大的点,即为拐点1',判断拐点1'是否在
g、将特征点
h、采用新的rn或者
h1、当0<rn中特征点的个数≤6时,用rn作为分配qrs波群的集合on;当
h2、若集合on中的特征点oj是极大值点且大于电势值vios,则oj是正向波;如果特征点oj是极大值点且小于电势值vios,则抛弃oj;如果集合on中的特征点oj是极小值点且小于电势值vios,则oj是负向波;如果特征点oj是极小值点且大于电势值vios,则抛弃oj;直至得到包含qrs波群中的所有正向波和负向波的集合
h3、若集合
h4、集合
h5、用
rr间期与心率:当前r波所在位置与上一个qrs波群中r波的时间差即为rr间期,一般以ms(毫秒)为单位,心率=60000/rr间期;
q波的幅度与宽度:如果有q波,在数据qn(qrson:qrsoff)中,从位置q波顶点向后寻找第一个大于幅度qn(qrson)的点,此点即为q波终末点qoff,也是r波起始点ron;q波的幅度为q波顶点幅度减去幅度qn(qrson),q波宽度即是qoff减去ron;
r波的幅度与宽度:如果由q波,r波起点ron就是qoff,如果没有q波,ron就是qrson。在数据qn(qrson:qrsoff)中,从位置r波顶点向后寻找第一个小于幅度qn(ron)的点,即为r波终末点roff,也是s波起始点son;r波的幅度为r波顶点幅度减去幅度ron,r波宽度即是roff减去ron;
s波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置s波顶点向后寻找第一个小于幅度qn(son)的点,即为s波终末点soff,也是r'波起始点r'on;s波的幅度为s波顶点幅度减去幅度son,s波宽度即是soff减去son;
r'波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置r'波顶点向后寻找第一个小于幅度qn(r'on)的点,即为r'波终末点r'off,也是s'波起始点s'on;r'波的幅度为r'波顶点幅度减去幅度r'on,r'波宽度即是r'off减去r'on;
s'波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置s'波顶点向后寻找第一个小于幅度qn(s'on)的点,即为s'波终末点s'off,也是r〞波起始点r〞on;s'波的幅度为s'波顶点幅度减去幅度s'on,s'波宽度即是s'off减去s'on;
r〞波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置r〞波顶点向后寻找第一个小于幅度qn(r〞on)的点,即为r〞波终末点r〞off,也是s〞波起始点s〞on;r〞波的幅度为r〞波顶点幅度减去幅度r〞on,r〞波宽度即是r〞off减去r〞on;
s〞波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置s〞波顶点向后寻找第一个小于幅度qn(s〞on)的点,即为s〞波终末点r〞off;s〞波的幅度为s〞波顶点幅度减去幅度s〞on,s〞波宽度即是s〞off减去s〞on;
qsr波群的第一个子波的起始点即为qrson,qsr波群的最后一个子波的终末点即为qrsoff;
i、重复c到i步,直至检出集合bxe中所有qrs的特征点,用bn={bn(j)}表示;
j、重复a到j步,直到识别完全部原始心电波形幅值数据s。
优选的,所述android主控板采用极值抽样显示法将心电图波形显示在显示屏上,该极值抽样显示法的步骤包括:
a、首先按照显示屏的分辨率计算抽样点步长
b、将心电图波形的原始采样数据按照抽点步长片段分成若干组;
c、每一组分别选取最大值和最小值作为抽取点,且保持两点的先后顺序关系;
d、显示当前一组的抽取点,重复步骤c和d直至显示完所有分组数据。
优选的,所述android主控板中预设有快速采集指令,该快速采集指令中包括默认用户信息、默认测量导联方式和默认测量时长。
进一步,所述传感器的数量为3个。
再者,还包括供电模块,所述供电模块与心电采集模块之间设有电子开关。
本发明一种基于通用平台的心电采集系统的采集方法,包括如下步骤:
a、参数设置:启动android主控板,设定检测的用户、测量导联方式和测量时长;
b、传感器的固定:用户根据测量导联方式将传感器固定于身体相应位置进行测量;
c、心电数据采集:传感器将心电信号输出至心电采集模块,经由心电采集模块处理形成心电数据发送至android主控板;
d、心电数据的接收:android主控板接收心电采集模块发送的心电数据,并存储;
e、心电数据的显示:android主控板将心电数据发送至显示屏,并显示。
本发明另一种基于通用平台的心电采集系统的采集方法,包括如下步骤:
a、参数设置:启动android主控板,调用快速采集指令,获得默认用户、默认测量导联方式和默认测量时长;
b、传感器的固定:用户根据测量导联方式将传感器固定于手指相应位置进行测量;
c、心电数据采集:传感器将心电信号输出至心电采集模块,经由心电采集模块处理形成心电数据发送至android主控板;
d、心电数据的接收:android主控板接收心电采集模块发送的心电数据,并存储;
e、心电数据的显示:android主控板将心电数据发送至显示屏,并显示。
有益效果:与现有技术相比,本发明具有以下显著优点:
(1)本发明提供一种基于通用平台的集成心电采集模块的心电采集系统,所采用平台保留了android平板的通用性,在作为健康保健产品同时,还允许播放视频,游戏娱乐等功能;
(2)本发明设有快速采集指令,可快速激发测量,使得部分身体状况不好,或年龄大不擅操作平板人士,轻易使用设备进入测量状态;再者,内置电子开关可自动启停心电采集模块的电源,减少电能消耗;
(3)本发明集成于android平台内,其整体体积小,方便携带,可应用于医院,保险,家庭保健等诸多领域,具有广泛的应用前景;
(4)本发明利用心电波形幅值数据计算心电波形的包络差、方差值、信号在1~5hz之间的功率和信噪比四个参数,分析四个参数值是否合格,根据参数包络差、方差值、信号在1~5hz之间的功率和信噪比的合格情况对心电波形进行质量等级划分,为心电信号自动分析筛选合格的心电数据,确保待检测心电信号的可用性;从而进一步保障后续精准检测到qrs波群的位置;
(5)本发明利用心电信号质量辨识方法对待分析的心电波形进行质量评估,对心电波形质量等级好的心电特征进行心电参数的计算和异常结果判断,对于心电波形质量等级差的心电信号发出提示,提高心电信号分析的可信度以及准确度;本发明最后结合波形质量分析法对心电波质量的甄别,有效地分析心电波形。
(6)本发明中qrs波群宽度的检测方法中通过对信号进行包络处理,凸显qrs波群位置信息,削弱杂波对其影响,有效提高qrs波群的精准位置的准确度;
(7)本发明中qrs波群宽度的检测方法中统计法得到等电位段电位,再利用幅度、斜率和极值点信息识别出qrs波群中个各关键点,各关键点、等电位段电位和预估的qrs波群的起始点和终末点最终确定qrs波群中各子波的位置,并准确勾勒出qrs波群的形态特征;同时本发明可求出各个子波的幅度与宽度,进而可以为心电qrs波群形态特征分析提供详细的特征参数;
(8)本发明中qrs波群宽度的检测方法中使用极值法寻找qrs波中的波峰和波谷,再通过找到的极值点及其结合斜率、三角形法、等电位段电位值、幅值信息来确定qrs波群的起止点,该qrs波群的起止点时段即是qrs波群的宽度;其中三角形法中通过多次重复使用三角形法将待判断拐点与特征点进行对比确认,可提高起止点的精准度,保障心电图qrs波群的宽度的准确检测;心电图qrs波群的宽度的准确检测对于心脏功能检查、心血管疾病的诊断和预防有着重要的意义。
附图说明
图1为本发明的系统示意图;
图2为本发明中快速测量的流程图;
图3为本发明的流程图;
图4为本发明中极值抽样显示法的示意图;
图5为本发明中心电信号质量辨识方法的流程框图;
图6为本发明中心电波形中各波的示意图;
图7为本发明中心电分析方法的流程框图;
图8为本发明实施例中的实验结果示意图一;
图9为本发明实施例中的实验结果示意图二;
图10为本发明的流程框图;
图11为本发明步骤b中信号x的处理方法的流程框图;
图12为本发明步骤b中信号x处理后的示意图;
图13为本发明步骤c中的直方图gp的示意图;
图14为本发明步骤e中极值点的左右区间示意图;
图15为本发明步骤f中三角形法的原理示意图;
图16为本发明步骤f中采用三角形法确定qrson和qrsoff点的示意图;
图17为本发明步骤h中正向波和负向波的示意图;
图18为本发明中qrs命名规则示意图一;
图19为本发明中qrs命名规则示意图二;
图20为本发明中qrs命名规则示意图三。
具体实施方式
下面结合附图对本发明的技术方案作进一步说明。
本发明一种基于通用平台的心电采集系统,包括android主控板、心电采集模块和传感器,所述心电采集模块与多个传感器相连接,所述传感器的数量为3个,所述传感器分别置于人体不同部位以实现心电信号测量,传感器的输出经心电采集模块处理后形成心电数据完成采集,并经串口通信将其发送至android主控板;
所述android主控板用于运行android应用程序,并对心电数据进行处理和储存,android主控板连接有显示屏,并通过显示屏即时将采集的心电图波形及心电参数数据显示出;android主控板通过串口采集到心电数据后,多线程同步实现如下不同处理:
1.1、进行心电信号质量辨识;所述心电信号质量辨识方法,包括如下步骤:
s1、获取采样率为fs的原始的心电波形幅值数据,并将心电波形幅值数据分成以n秒为时间单位的处理段x,n为8~120间的任意一个值;
s2、在处理段x内,将心电信号以1秒为单位进行片段分割,分割段为xi,其中i=1…n,xi的数据长度为m=1×fs=fs;
s3、以处理段x为单位,取出每一心电波形片段xi的幅值最大值mmax和最小值mmin构成包络点,通过比较各个包络点得到分割段xi的包络差ei=mmax–mmin,再求取处理段x内波形包络差的平均值
s4、在处理段x内,获取每一分割段xi的幅值方差值
s5、在处理段x内,利用快速傅里叶变换,将心电片段信号xi的时域信号转化为频域信号,即为功率谱信号;对1~5hz的幅值进行积分得到功率
s6、在处理段x内,得到步骤s5中的心电片段信号xi的功率谱信号后,对5~40hz的幅值进行积分得到功率
s7、根据参数包络差、方差值、信号在1~5hz之间的功率和信噪比的合格情况对心电波形进行质量等级划分,心电波形的质量分成4个等级,分别用0~3表示质量等级系数;其中心电波形的质量等级划分方法为:当信噪比参数不合格时,该处理段心电波形的质量等级系数直接评定为3,即波形存在严重噪声;在信噪比参数合格的前提下,根据包络差、方差值、信号在1~5hz之间的功率的合格情况来评定心电波形质量,当三个参数全部合格时,该处理段心电波形质量等级系数评定为0,即波形良好;当三个参数中有两个或一个合格时,心电波形质量等级系数评定为1,即波形较差;当三个参数全部不合格时,心电波形质量等级系数评定为2,即波形差;质量等级系数0表示心电波形良好,质量等级系数1表示心电波形较差,质量等级系数2表示心电波形差,质量等级系数3表示心电波形存在严重噪声。即从波形质量的程度“0”表示良好,“1”表示较差,“2”表示差,“3”表示严重噪声。
1.2、根据辨识后的心电信号进行心电分析,计算得到心电参数、心率、qrs时间和st高度参数;所述心电分析方法,包括如下步骤:
s11、获取n秒的待分析的心电波形,n为8~120间的任意一个值;
s12、依据心电信号质量辨识方法对待分析的心电波形进行质量评估,得到n秒的心电波形质量等级,其中心电波形质量评估以1s为单位,共得到n个数据分别代表第1秒到第n秒的质量等级系数;
s13、当心电信号质量等级系数为3且持续时间大于0.3n秒时,提示“信号质量噪声过大,无法分析”返回至步骤a重新获取待分析的心电波形;否则进行心电参数的计算和异常结果判断;其中心电参数的计算和异常结果判断的步骤为:
s1301、依据心电qrs波群宽度的检测方法找出qrs波群宽度、qrs波群中各个子波形的起止时间点、幅度参数和代表qrs波群形态的bn;然后再找出心电特征包括p、t波的起止时间点及其幅度参数,构成了一次心跳;
如图10所示,其中心电qrs波群宽度的检测方法,包括如下步骤:
a、获取采样率为fs的原始心电波形幅值数据s,并将心电波形s分成以n秒为时间单位的处理段信号x,n为8~120间的任意一值;
b、将信号x经过包络处理得到包络信号xe,然后在包络信号xe上利用阈值法定位qrs波群的位置得到集合bxe;
如图11所示,其中信号x的处理方法和利用阀值法的定位方法包括如下步骤:
b1、将信号x经过一个带通滤波器,其频率范围为5hz~25hz,得到凸显出qrs波群特征且滤除干扰波的信号xfiltered;对信号xfiltered做微分后再平方得到信号xdiff,对xdiff进行积分得到包络信号xe,如图12所示;
b2、将包络信号xe按0.2n秒的长度分段得到5个包络片段xe(i),其中i=1…5,在每个包络片段内找到最大值xemax(i),然后求出这5个最大值的中值
c、在集合bxe中选取第n个qrs波群位置bxe(n),求出该qrs波群的rr间期rrxe(n)=bxe(n)-bxe(n-1);从信号x中截取一片段xios,截取的起始位置为bxe(n-1)+rrxe(n)×k3,终止位置为bxe(n)-rrxe(n)×(0.5-k3),其中0.3<k3<0.5,即第n-1个心跳的t波之后到第n个心跳的p波之前的一段数据;作出片段xios中幅值范围的直方图gp,得到xios数据取值的统计分布,如图13所示;在直方图gp中找出xios幅值分布个数最多的子区间b,该子区间b代表的电位即为第n次心跳之前的等电位段电位vios;
d、从信号x中第n个qrs波群位置bxe(n)周围截取一数据段qn,在数据段qn中找到所有极大值和极小值点的集合为极值点p和qrs波群中绝对幅值最大的波qm;
其中获取极值点p和qm的方法包括如下步骤:
d1、信号x在bxe(n)处往前ω-=(0.07~0.15)fs、往后ω+=(0.1~0.23)fs个点,截取信号x中的qrs波群的数据段qn=(q1,…,qj,…,qω),其中ω=ω-+ω+表示qn的数据长度,qj表示qn中第j个数据;
d2、在数据段qn中找到所有极大值pk和极小值点lo,极大值pk和极小值点lo合称为极值点p;然后再找出极大值pk中的最大值pmax=(vmax,imax)和极小值lo中的最小值lmin=(vmin,imin),i、v分别表示点在qn中的位置和数据幅值大小;
d3、定义qrs波群中最明显的波顶点为qrs波群中绝对幅值最大的波qm,如果|vmax-vios|>k4×|vmin-vios|则qm点为lmin点,否则qm点为pmax点,其中2<k4<10;
e、极值点p中满足最小qrs可识别波条件的为主极点dn,再从主极点dn中筛选出特征点
其中筛选特征点
e1、从极值点p中选择一个极值点pj,用以下公式在处理段qn中搜索极值点pj的左右支撑区间
δqj,x=|qj-qx|
其中j和qj为极值点pj在qn中的位置和幅值,τ为qrs最大的有意义的生理时间宽度,取80~160ms;i、a、b和k均是为了求解支撑区间的辅助变量,没有特别含义;qx是qn中第x个数据,与qj一个意义,如图14所示;
e2、若极值点pj的左右支撑区间
e3、重复步骤e1、e2直到筛选出极值点p中所有可识别波,记为主极值点集合dn={pj};
e4、极值点pj的左斜率
e5、如
f、在数据段qn中,以步骤d中得到的qm为顶点,利用三角形法确定qrs的起止点qrson和qrsoff;所述三角形法包括如下步骤:
f1、如图15所示,三角形法的作用是寻找波形上的转折点,转折点包括极值点和非极值拐点;一心电波形线段,从顶点x开始向前取一个点z做一条辅助线l,求该线段上所有数据点到直线l的距离d,距离最大的点即为转折点y;使用三角形法找qrson和qrsoff的前提是:qrson和qrsoff是波形的非极值拐点,而不是峰值点;
f2、在数据段数据qn中首先以qm点作为顶点,向前50~200ms之间寻找一点作辅助线s1,寻找数据段数据qn中到辅助线s1的距离最大的点,即为拐点1,判断拐点1是否在
f3、在数据段数据qn中首先以qm点作为顶点,向后50~200ms之间寻找一点作辅助线l1,寻找数据段数据qn中到辅助线l1的距离最大的点,即为拐点1',判断拐点1'是否在
g、将特征点
h、采用新的rn或者
h1、当0<rn中特征点的个数≤6时,用rn作为分配qrs波群的集合on;当
h2、若集合on中的特征点oj是极大值点且大于电势值vios,则oj是正向波;如果特征点oj是极大值点且小于电势值vios,则抛弃oj;如果集合on中的特征点oj是极小值点且小于电势值vios,则oj是负向波;如果特征点oj是极小值点且大于电势值vios,则抛弃oj;直至得到包含qrs波群中的所有正向波和负向波的集合
h3、若集合
h4、如图18所示,集合
h5、用
rr间期与心率:当前r波所在位置与上一个qrs波群中r波的时间差即为rr间期,一般以ms(毫秒)为单位,心率=60000/rr间期;
q波的幅度与宽度:如果有q波,在数据qn(qrson:qrsoff)中,从位置q波顶点向后寻找第一个大于幅度qn(qrson)的点,此点即为q波终末点qoff,也是r波起始点ron;q波的幅度为q波顶点幅度减去幅度qn(qrson),q波宽度即是qoff减去ron;
r波的幅度与宽度:如果由q波,r波起点ron就是qoff,如果没有q波,ron就是qrson。在数据qn(qrson:qrsoff)中,从位置r波顶点向后寻找第一个小于幅度qn(ron)的点,即为r波终末点roff,也是s波起始点son;r波的幅度为r波顶点幅度减去幅度ron,r波宽度即是roff减去ron;
s波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置s波顶点向后寻找第一个小于幅度qn(son)的点,即为s波终末点soff,也是r'波起始点r'on;s波的幅度为s波顶点幅度减去幅度son,s波宽度即是soff减去son;
r'波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置r'波顶点向后寻找第一个小于幅度qn(r'on)的点,即为r'波终末点r'off,也是s'波起始点s'on;r'波的幅度为r'波顶点幅度减去幅度r'on,r'波宽度即是r'off减去r'on;
s'波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置s'波顶点向后寻找第一个小于幅度qn(s'on)的点,即为s'波终末点s'off,也是r〞波起始点r〞on;s'波的幅度为s'波顶点幅度减去幅度s'on,s'波宽度即是s'off减去s'on;
r〞波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置r〞波顶点向后寻找第一个小于幅度qn(r〞on)的点,即为r〞波终末点r〞off,也是s〞波起始点s〞on;r〞波的幅度为r〞波顶点幅度减去幅度r〞on,r〞波宽度即是r〞off减去r〞on;
s〞波的幅度与宽度:在数据qn(qrson:qrsoff)中,从位置s〞波顶点向后寻找第一个小于幅度qn(s〞on)的点,即为s〞波终末点r〞off;s〞波的幅度为s〞波顶点幅度减去幅度s〞on,s〞波宽度即是s〞off减去s〞on;
qsr波群的第一个子波的起始点即为qrson,qsr波群的最后一个子波的终末点即为qrsoff;
i、重复c到i步,直至检出集合bxe中所有qrs的特征点,用bn={bn(j)}表示;
j、重复a到j步,直到识别完全部原始心电波形幅值数据s;
s1302、取心电波形质量等级系数小于3的心电特征用于计算心电参数rr间期和心率,用心电波形质量等级为0的心电特征计算心电参数pr间期、qrs宽度、qt间期、st段高度、p波平均电轴、qrs平均电轴、t波平均电轴、q波宽度和r波高度;
s1303、利用心电特征和上述心电参数进行心电异常判断得出心电图异常结果。
s14、当心电信号质量等级系数为1或2且持续时间大于0.3n秒时,提示“信号质量较差”并输出计算得到的心电参数和心电图异常结果;否则认为心电波形质量较好,直接输出计算得到的心电参数和心电图异常结果;
1.3、本地数据保存:包括原始心电数据保存,计算得到的心电参数和心电结果保存;
1.4、网络服务器数据保存:本地数据同步上传服务器;
所述android主控板、心电采集模块和传感器均集成于android平台内。本发明提供一种基于通用平台的集成心电采集模块的心电采集系统,所采用平台保留了android平板的通用性,在作为健康保健产品同时,还允许播放视频,游戏娱乐等功能。
本发明的android主控板采用极值抽样显示法将心电图波形显示在显示屏上,该极值抽样显示法的步骤包括:
a、首先按照显示屏的分辨率计算抽样点步长
b、将心电图波形的原始采样数据按照抽点步长片段分成若干组;
c、每一组分别选取最大值和最小值作为抽取点,且保持两点的先后顺序关系;
d、显示当前一组的抽取点,重复步骤c和d直至显示完所有分组数据。
本发明android主控板中预设有快速采集指令,该快速采集指令中包括默认用户信息、默认测量导联方式和默认测量时长。本发明设有快速采集指令,可快速激发测量,使得部分身体状况不好,或年龄大不擅操作平板人士,轻易使用设备进入测量状态。
还包括供电模块,所述供电模块与心电采集模块之间设有电子开关。内置电子开关可自动启停心电采集模块的电源,减少电能消耗。
第一种基于通用平台的心电采集系统的采集方法,包括如下步骤:
a、参数设置:启动android主控板,设定检测的用户、测量导联方式和测量时长;
b、传感器的固定:用户根据测量导联方式将传感器固定于身体相应位置进行测量;
c、心电数据采集:传感器将心电信号输出至心电采集模块,经由心电采集模块处理形成心电数据发送至android主控板;
d、心电数据的接收:android主控板接收心电采集模块发送的心电数据,并存储;
e、心电数据的显示:android主控板将心电数据发送至显示屏,并显示。
第二种基于通用平台的心电采集系统的采集方法,包括如下步骤:
a、参数设置:启动android主控板,调用快速采集指令,获得默认用户、默认测量导联方式和默认测量时长;
b、传感器的固定:用户根据测量导联方式将传感器固定于手指相应位置进行测量;
c、心电数据采集:传感器将心电信号输出至心电采集模块,经由心电采集模块处理形成心电数据发送至android主控板;
d、心电数据的接收:android主控板接收心电采集模块发送的心电数据,并存储;
e、心电数据的显示:android主控板将心电数据发送至显示屏,并显示。
实施例
如图5所示,本发明一种心电信号质量辨识方法,包括如下步骤:
s1、获取采样率为fs的原始的心电波形幅值数据,并将心电波形幅值数据分成以10秒为时间单位的处理段x;
s2、在处理段x内,将心电信号以1秒为单位进行片段分割,分割段为xi,其中i=1…10,xi的数据长度为m=1×fs=fs;
s3、以处理段x为单位,取出每一心电波形片段xi的幅值最大值mmax和最小值mmin构成包络点,通过比较各个包络点得到分割段xi的包络差ei=mmax–mmin,再求取处理段x内波形包络差的平均值
s4、在处理段x内,获取每一分割段xi的幅值方差值
s5、在处理段x内,利用快速傅里叶变换,将心电片段信号xi的时域信号转化为频域信号,即为功率谱信号;对1~5hz的幅值进行积分得到功率
s6、在处理段x内,得到步骤s5中的心电片段信号xi的功率谱信号后,对5~40hz的幅值进行积分得到功率
s7、根据参数包络差、方差值、信号在1~5hz之间的功率和信噪比的合格情况对心电波形进行质量等级划分,心电波形的质量分成4个等级,分别用0~3表示质量等级系数;其中心电波形的质量等级划分方法为:当信噪比参数不合格时,该处理段心电波形的质量等级系数直接评定为3,即波形存在严重噪声;在信噪比参数合格的前提下,根据包络差、方差值、信号在1~5hz之间的功率的合格情况来评定心电波形质量,当三个参数全部合格时,该处理段心电波形质量等级系数评定为0,即波形良好;当三个参数中有两个或一个合格时,心电波形质量等级系数评定为1,即波形较差;当三个参数全部不合格时,心电波形质量等级系数评定为2,即波形差;质量等级系数0表示心电波形良好,质量等级系数1表示心电波形较差,质量等级系数2表示心电波形差,质量等级系数3表示心电波形存在严重噪声。即从波形质量的程度“0”表示良好,“1”表示较差,“2”表示差,“3”表示严重噪声。
如图6所示,本发明中所提及的心电波形p波为心房激动波,也称为心房的除极波;qrs波群为心室激动的综合波群,也称为心室的除极波,其第一个向下的波为q波、任何向上的波都叫r波、在r波以后的任何向下的波都叫s波;t波称为心室的复极波,而心房的复极波很小,且埋于心室的除极波中,不易辨认,未给于专门命名。
如图7所示,本发明一种心电信号质量辨识方法的心电分析方法,包括如下步骤:
s11、获取10秒的待分析的心电波形;
s12、依据心电信号质量辨识方法对待分析的心电波形进行质量评估,得到10秒的心电波形质量等级,其中心电波形质量评估以1s为单位,共得到10个数据分别代表第1秒到第10秒的质量等级系数;
s13、当心电信号质量等级系数为3且持续时间大于3秒时,提示“信号质量噪声过大,无法分析”返回至步骤a重新获取待分析的心电波形;否则进行心电参数的计算和异常结果判断;其中心电参数的计算和异常结果判断的步骤为:
s1301、根据心电波形提取心电特征,心电特征包括p、q、r、s、t波的起止时间点及其幅度参数,构成了一次心跳;
s1302、取心电波形质量等级系数小于3的心电特征用于计算心电参数rr间期和心率,用心电波形质量等级为0的心电特征计算心电参数pr间期、qrs宽度、qt间期、st段高度、p波平均电轴、qrs平均电轴、t波平均电轴、q波宽度和r波高度;
s1303、利用心电特征和上述心电参数进行心电异常判断得出心电图异常结果。
s14、当心电信号质量等级系数为1或2且持续时间大于3秒时,提示“信号质量较差”并输出计算得到的心电参数和心电图异常结果;否则认为心电波形质量较好,直接输出计算得到的心电参数和心电图异常结果。
如图8和图9所示,本发明以10为时间单位进行心电分析的实验数据图和结果,图中方波表示质量辩识结果的标识,标识方波中:0表示波形质量好和良好;1、2表示波形质量较差和差;3表示波形质量严重差。
本发明中还预设有用户管理模块,该用户管理模块中设置多个用户,且在每一用户下设置对应的的信息属性。心电采集系统还包括血压采集模块、血糖采集模块、血氧采集模块、血脂采集模块、筛查采集模块和体脂采集模块,其采集方式为usb自动采集、蓝牙自动采集或人工输入采集中的一种或多种组合。采集得到的数据的显示方式有单次测量变动趋势图、散点图和30天趋势图,其中单次测量变动趋势图是以心跳次数为横轴,心电参数为纵轴显示的参数变化趋势图,散点图是由rr间期和差分rr间期绘制出的有关心率变异性的图表,30天趋势图是利用一天内测量得到的心电参数平均值和天数,绘制的30天内的参数变化趋势图。本发明提供包括微信在内的多种即时通讯发布方式提醒关联人员。具体方式如下:将心电参数、心电波形、自动判断的结果组成一定格式文件以利于通信工具传输;选择目标接收对象;调用系统api接口或是通信工具发布的接口,将心电文件发送出去;系统将通过因特网将数据上传至通信工具服务器,而后由服务器将数据转发给选定的接收对象。本发明可提供多人定时服药提醒,在各种状态下,定时闹铃提醒服药信息,包括人员名称,药品名,服药剂量,服药方法;在包括关机,锁屏,待机,普通操作等状态,弹出服药提醒对话框,并语音播报服药信息;在具有功能排他性的心电测量时,不弹出服药提醒对话框,只设置通知,并语音播报服药信息;各项数据通过wifi或数据流量上传至数据服务器,实现备份,网络诊断等功能。