CN105125199A - 一种心率检测方法与装置 - Google Patents
一种心率检测方法与装置 Download PDFInfo
- Publication number
- CN105125199A CN105125199A CN201510455291.8A CN201510455291A CN105125199A CN 105125199 A CN105125199 A CN 105125199A CN 201510455291 A CN201510455291 A CN 201510455291A CN 105125199 A CN105125199 A CN 105125199A
- Authority
- CN
- China
- Prior art keywords
- ripple
- maximum
- threshold value
- heart rate
- sample values
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
Abstract
本发明实施例公开了一种心率检测方法与装置,属于心电数据处理技术领域,基于传统方法实现心率的准确计算,降低心率计算的复杂性。本发明实施例提供的方法包括:循环执行第一流程直至获取若干个R波;所述第一流程包括:采样心电数据以获取新的采样值,记所述新采样值的前一采样值为当前采样值;对当前采样值执行R波判定流程;根据所述获取的若干个R波确定当前心率。所述R波判定流程主要是根据R波的特点,排除与R波类似波形的影响,找出正确的R波,为心率的计算提供正确的数据。
Description
技术领域
本发明涉及心电数据处理技术领域,尤其涉及一种心率检测方法与装置。
背景技术
心率(HeartRate)是指心脏每分钟跳动的次数,可以表示心脏跳动快慢的意思。正常成年人安静时与运动状态的心率会有差异,不同成年人在同一状态下也会有显著的个体差异。而心率的数值在一定程度上可以体现人的身体健康状况。
心电数据具有周期性,在一个周期内,主要分为Q波、R波和S波,参考图1所示的心电图数据,目前最为常用的是检测R波的方式来确定心率。
比如句法的方法,具体为:将ECG(electrocardiogram,心电图)中不同的波形轮廓和线段(即模式)用一系列符号代表,然后检测这些符号所构成的序列。当某一序列符合QRS波所具有的符号序列时,则判定该序列所对应的ECG段为一QRS波。但句法的分析方法的处理速度较慢,分析规则的直观性不强,且在进一步的心率失常分析处理中不便沿用医生传统的分析方法与步骤,因而应用不够广泛。
除了句法的方法外,目前以检测R波的方式来确定心率的方法中还有非句法方法,比如在神经网络的方法、小波的方法以及软硬件结合的方法等等。小波的方法主要应用小波变换的多尺度特性,但多处于离线分析处理阶段,尽管目前微处理器和各种DSP芯片发展迅速,但实现起来仍显得复杂。寻求一种简单有效的方法确定心率具有重要的意义。
发明内容
为了解决上述技术问题,本发明实施例提供一种心率检测方法与装置,该方法基于传统的心电数据处理方法,实现简单。
有鉴于此,本发明实施例提供一种心率检测方法,该方法可以包括:
循环执行第一流程直至获取若干个R波;所述第一流程包括:采样心电数据以获取新的采样值,记所述新采样值的前一采样值为当前采样值;对当前采样值执行R波判定流程;
根据所述获取的若干个R波确定当前心率;
所述R波判定流程包括:若当前采样值为极大值,则标记极大值并返回执行第一流程;若当前采样值为极小值且极大值标记存在,则判断当前采样值和相邻的极大值采样间隔是否小于第一预设阈值;若小于,判定当前采样值相邻的极大值对应的波为R波;删除极大值标记;
其中,极大值的条件为大于左右相邻的采样值且大于波峰阈值;极小值的条件为小于左右相邻的采样值且小于波谷阈值。
优选的,根据若干个R波确定当前心率,计算公式为:
H=Sample*60/Raverage;
其中,H表示心率,Sample为心电数据采样率,Raverage为由若干个R波确定的两个R波之间平均间隔。
优选的,所述方法还包括:
记与当前采样值相邻的极大值为第一极大值,与第一极大值相邻的极大值为第二极大值;
当前采样值和第一极大值采样间隔小于第一预设阈值时,判断第一极大值和第二极大值的采样间隔是否小于第二预设阈值,若大于,则判定第一极大值对应的波为R波;若小于,则对R波进行修正;
所述对R波进行修正包括:将第一极大值和第二极大值中较大者对应的波确定为R波。
优选的,所述方法还包括:
在确定一个R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/4;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(3/8);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
优选的,所述方法还包括:
在确定一个R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/2;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(7/16);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
本发明还提供一种心率检测装置,所述装置包括:
执行模块,用于循环执行第一流程直至获取若干个R波;所述第一流程包括:采样心电数据以获取新的采样值,记所述新采样值的前一采样值为当前采样值;对当前采样值执行R波判定流程;
心率计算模块,用于根据所述获取的若干个R波确定当前时刻的心率;
R波判定模块,用于执行所述R波判定流程,所述R波判定流程包括:若当前采样值为极大值,则标记极大值并返回执行第一流程;若当前采样值为极小值且极大值标记存在,则判断当前采样值和相邻的极大值采样间隔是否小于第一预设阈值;若小于,判定当前采样值相邻的极大值对应的波为R波;删除极大值标记;
其中,极大值的条件为大于左右相邻的采样值且大于波峰阈值;极小值的条件为小于左右相邻的采样值且小于波谷阈值。
优选的,心率计算模块确定心率的计算公式为:
H=Sample*60/Raverage;
其中,H表示心率,Sample为心电数据采样率,Raverage为由若干个R波确定的R波平均间隔。
优选的,记与当前采样值相邻的极大值为第一极大值,与第一极大值相邻的极大值为第二极大值;所述装置还包括:
判定模块,用于当前采样值和第一极大值采样间隔小于第一预设阈值时,判断第一极大值和第二极大值的采样间隔是否小于第二预设阈值;
相应的,第一极大值和第二极大值的采样间隔大于第二预设阈值时,R波判定模块判定第一极大值对应的波为R波;
R波修正模块,用于第一极大值和第二极大值的采样间隔小于第二预设阈值,对R波进行修正,包括:将第一极大值和第二极大值中较大者对应的波确定为R波。
优选的,所述装置还包括:
第一阈值更新模块,用于在确定R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/4;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(3/8);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
优选的,所述装置还包括:
第二阈值更新模块,用于在确定一个R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/2;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(7/16);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
与现有技术相比,本发明实施例提供的技术方案具有以下优点:
本发明提供的技术方案,根据R波的特点,对心电图数据进行采样,并对采样数据进行极大值或极小值判断,确定出间隔小于第一预设阈值的极大值和极小值,初步确定该极大值对应的波为R波。本方案不仅可以实现R波的判定,且算法简单,同时通过合理设定第一预设阈值,可以排除心电数据中的T波(参考图1以及对图6的说明)对R波的影响,造成R波误判,提高R波判定的准确性,进而可以提高心率计算的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为一个周期内心电图数据示意图;
图2为本发明提供的一种心率检测方法实施例1的流程图;
图3为本发明提供的R波判定流程图;
图4为本发明提供的R波判定优化流程图;
图5为一种心电数据示意图;
图6为一种心电图中P波骤变示意图;
图7为本发明提供的一种心率检测装置实施例1的结构示意图。
具体实施方式
本发明实施例提供了一种心率检测方法与装置,基于传统方法实现R波的准确判定,实现简单,计算复杂度低,且可以保证心率计算结果的可靠性和准确性。
为了使本技术领域的人员更好地理解本发明方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分的实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
参考图2所示,为本发明提供的一种心率检测方法实施例1的流程图,本实施例主要可以包括如下步骤:
S201:循环执行第一流程直至获取若干个R波。
本步骤中要获取的R波的个数,可以提前设定,具体数值可以根据需要设定,这里不作具体限定。
其中,所述第一流程包括如下步骤:
步骤A1:采样心电数据以获取新的采样值。
实际操作时,心电数据采样频率可以自定义设置,比如200HZ。
步骤A2:记所述新采样值的前一采样值为当前采样值。
步骤A3:对当前采样值执行R波判定流程。
执行完步骤A1-步骤A3,则第一流程执行完成。
当获取若干个R波之后,进入步骤S202。
S202:根据所述获取的若干个R波确定当前心率。
下面对R波判定流程进行详细说明,参考图3所示的流程图,所述R波判定流程具体可以包括:
步骤B1:判断当前采样值是否为一个极大值。
其中,极大值为大于左右相邻的采样值和波峰阈值的采样值,即当前采样值要大于前一采样值和后一采样值(新的采样值),同时要大于波峰阈值(R波的波峰阈值),如果当前采样值是一个极大值,则对极大值进行标记,然后返回执行第一流程。
如果当前采样值不是一个极大值,则进入步骤B2。
步骤B2:若极大值标记存在,判断当前采样值是否是一个极小值。
极小值为小于左右相邻的采样值和波谷阈值(R波的波谷阈值)的采样值。若当前采样值是一个极小值且极大值标记存在,则执行步骤B3,反之则返回执行第一流程。需要说明的是,每确定一个极大值,则更新极大值标记,所以本步骤中的极大值是指最新确定的极大值。
步骤B3:判断当前采样值和相邻极大值的采样间隔是否小于第一预设阈值。
如果当前采样值和相邻的极大值的采样间隔小于第一预设阈值,则可以认为当前采样值相邻的极大值对应的波为R波。
步骤B4:确定R波后,删除极大值标记。
由于心电数据采样具有一定的频率,所以不同采样值之间的时间间隔除了可以直接采用时间表示,也可以采用采样点数进行表示,通过采样频率可以确定不同采样点数对应的时间。
本发明提供的技术方案,根据R波的特点,对心电图数据进行采样,并对采样数据进行极大值或极小值判断,确定出间隔小于第一预设阈值的极大值和极小值,初步确定该极大值对应的波为R波。本方案不仅可以实现R波的判定,同时通过合理设定第一预设阈值,可以排除心电数据中的T波(参考图1以及对图6的说明)对R波的影响,造成R波误判,提高R波判定的准确性,进而可以提高心率计算的准确性。
在一些实施例中,根据若干个R波确定当前心率,优选采用以下计算公式:H=Sample*60/Raverage;其中,H表示心率,Sample为心电数据采样率,Raverage为由若干个R波确定的R波平均间隔。
比如:假设若干个R波为Rwave[0]、Rwave[1]、Rwave[2],则Raverage=((Rwave[1]-Rwave[0])+(Rwave[2]-Rwave[1]))/2。
考虑到实际情况中存在人的心电信号骤变的现象,比如人受到某种刺激时心电信号出现骤然变小或变大(波形较为尖锐),这种现象会对R波的判定造成一定的影响,比如误判等。所以为了在实施例1的基础上,进一步解决上述问题,本发明还提供一种心率检测方法实施例2,本实施例主要对R波判定流程进行优化,参考图4所示的流程图,为了便于描述,本实施例中记与当前采样值相邻的极大值为第一极大值,与第一极大值相邻的极大值为第二极大值。
本实施例除了包括实施例1中的步骤外,进一步的R波判定流程还可以包括如下步骤:
步骤B5:当前采样值和第一极大值采样间隔小于第一预设阈值时,判断第一极大值和第二极大值的采样间隔是否小于第二预设阈值。
优选的,第二预设阈值可以定义为连个R波之间的时间间隔。
第一极大值和第二极大值的采样间隔小于第二预设阈值时,进入步骤B6。
步骤B6:对R波进行修正;
所述对R波进行修正包括:将第一极大值和第二极大值中较大者对应的波确定为R波。
步骤B7:在确定一个R波之后,删除极大值标记,更新波峰阈值和波谷阈值。
在确定一个R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值,更新公式如下:
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/4;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(3/8);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
为了更好的适应心电数据的骤变,需要阈值的设定更加宽泛,优选的,阈值还可以采用以下公式进行更新:
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/2;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(7/16);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
需要说明的是,步骤B4的执行与步骤B7-B8的执行没有先后顺序的限定。
本实施例不仅可以带来方法实施例1的有益效果,进一步的,对R波的判定进行了流程优化,并且合理设定了阈值的更新方式,可以避免由于心电数据骤变对R判定的影响,提高心率计算结果的准确性。
对于R波判定以及阈值更新的步骤,下面举例说明,记第i个采样值为N(i),每获取一个新的采样值,则采样个数计数加1(采用采样点数表示时间),采用极大值标志位Maxflag标记极大值。
步骤1:读取新采样值N(i+1),采样点数加1。
步骤2:判断Ni是否满足(1)式:Ni>N(i-1)且Ni>N(i+1)。
若N(i)满足(1)式,且N(i)大于波峰阈值,则N(i)对应的波段可能为R波,记N(i)的值和采样位置分别为Max和Rtime,极大值标志位maxflag置1。
i=i+1,返回执行步骤1。
如果N(i)不满足(1)式,或N(i)不小于波峰阈值,则进入步骤3。
步骤3:判断极大值标志是否为1。
若极大值标志为1,判断N(i)是否满足式(2):N(i)<N(i-1)且N(i)<N(i+1)。
若N(i)满足式(2)且N(i)小于波谷阈值,记N(i)的值和采样位置分别为Min和Stime。进入步骤4。
若N(i)不满足式(2)或N(i)不小于波谷阈值。
i=i+1,返回步骤1,继续采样心电数据。
若极大值标志位不为1,则i=i+1,返回步骤1继续采样。
步骤4:判断相邻两次极值Stime和Rtime对应时间差是否小于第一预设阈值。
若小于第一预设阈值,则进入步骤5,反之则i=i+1,返回步骤1继续采样。
步骤5:判断(Rwavenew-Rwaveold)是否小于两次R波间的时间长度。
其中,Rwaveold为上一R波对应的采样点数,Rwavenew为新的R波的采样点数,上一R波在Rwaveold采样时刻的采样值为Maxold,新的R波在Rwavenew采样时刻上的采样值为Maxnew。
若(Rwavenew-Rwaveold)不小于两次R波间的时间长度,则认为Maxnew是一个R波,根据R波对应的Max和Min更新波峰阈值和波谷阈值。
若(Rwavenew-Rwaveold)小于两次R波间的时间长度,则将Maxold和Maxnew中较大者对应的波认为是R波,更新波峰阈值和波谷阈值,极大值标志清零。
确定R波位置后,根据R波的位置来确定心率。
比如,记录第一个R波为Rwave[0],第二个R波为Rwave[1],则心率计算如下:H=Sample*60/(Rwave[1]-Rwave[0])。
在开始检测的2秒内,要确定最初的一个波峰阈值和波谷阈值。因为人的心率最慢一分钟30下,即在2秒内至少有一次心电图信号,这是第一次设置阈值的依据。找到一个极大值和一个极小值,并根据这两个值分别算出两个阈值。
阈值的计算公式分别是:
波峰阈值:Maxthreshold=Min+(Max-Min)/4;
波谷阈值:Minthreshold=Min+(Max-Min)*(3/8)。
其中,3/8(8分之3)是1/3的近似值,这样便于在C语言中直接用右移指令完成,以确保其实时性。
确定R波后,计算相邻R波间期:RRinterval=Rwave[1]-Rwave[0],在实际操作中可以计算多个相邻R波间期,并取平均值Raverage。此时,计算心率的公式是:H=sample*60/Raverage。其中,sample是采样率,根据设置的不同,可以有不同的采样率,优选的,采样率可以为200Hz。
基于硬件实现操作时,每计算一次心率后,Raverage和HRcountflag(心率计算标志位)都清零,以备下次使用。在开始采样时,由于数据不足,还不能计算心率,或是心率计算已经完成,需要为下一次计算做准备,则调用找最(极)大值函数,即采样值要大于左右相邻的两个采样值,且要大于波峰阈值,如图5中的D点,就不会被误认为是一个R波。采样点重新回到峰阈值以下,如F点,立即寻找一个极小值。这也是为了保证先找到可能的R波,再分析后续的采样是否为S波。只有满足RS复波的条件,才可以确定R波。如果极小值有效,判断极大和极小值的间隔。满足条件的,则可以认为这个波有可能是一个R波。否则,如图5中所示,当H和I点也分别大于波峰阈值和小于波谷阈值,但是它们之间的间隔,远大于窗的宽度(第一预设阈值),由此可以排除。
也有一些心电信号如图1中的P波较为尖锐时,和R波非常接近,参考图6所示。这时,有可能将X点也认为是一个R波。所以,当检测到第二个R波时,实际上是真正的R波(Y点),和上一次认为的R波X点比较一下之间的距离,如果小于特定值,则这个R波(Y点)和上次的R波(X点)是同一个R波。只记一次R波的位置,Y点的位置。
另外,考虑到基线漂移等因素,阈值设为可变的,根据新的R波重新设定阈值。由于,有的病人心电信号可能会出现骤然变小的情况,阈值的设定就变得更加宽泛。
R波峰阈值:Maxthreshold=Min+(Max-Min)/2,
R波谷阈值:Minthrshold=Min+(Max-Min)*(7/16)。
对应于上述方法实施例,下面本发明提供的装置实施例,参考图7所示,为本发明提供的一种心率检测装置实施例1的结构示意图,本实施例中的装置包括:执行模块701、心率计算模块702、R波判定模块703。
其中,所述执行模块701,用于循环执行第一流程直至获取若干个R波;所述第一流程包括:采样心电数据以获取新的采样值,记所述新采样值的前一采样值为当前采样值;对当前采样值执行R波判定流程;
心率计算模块702,用于根据所述获取的若干个R波确定当前时刻的心率;
R波判定模块703,用于执行所述R波判定流程,所述R波判定流程包括:若当前采样值为极大值,则标记极大值并返回执行第一流程;若当前采样值为极小值且极大值标记存在,则判断当前采样值和相邻的极大值采样间隔是否小于第一预设阈值;若小于,判定当前采样值相邻的极大值对应的波为R波;删除极大值标记。
其中,极大值的条件为大于左右相邻的采样值且大于波峰阈值;极小值的条件为小于左右相邻的采样值且小于波谷阈值。
本实施例中,根据R波的特点,对心电图数据进行采样,并对采样数据进行极大值或极小值判断,确定出间隔小于第一预设阈值的极大值和极小值,初步确定该极大值对应的波为R波。本方案不仅可以实现R波的判定,同时通过合理设定第一预设阈值,可以排除心电数据中的T波(参考图1以及对图6的说明)对R波的影响,造成R波误判,提高R波判定的准确性,进而可以提高心率计算的准确性。
在一些实施例中,优选的心率计算公式为:H=Sample*60/Raverage;其中,H表示心率,Sample为心电数据采样率,Raverage为由若干个R波确定的R波平均间隔。
在上述装置实施例1的基础上,本发明提供的一种心率检测装置实施例2,本实施例除了包括装置实施例1的组成模块之外,进一步的,为了便于描述,本实施例中记与当前采样值相邻的极大值为第一极大值,与第一极大值相邻的极大值为第二极大值。本实施例中的装置还可以包括:
判定模块,用于当前采样值和第一极大值采样间隔小于第一预设阈值时,判断第一极大值和第二极大值的采样间隔是否小于第二预设阈值;
R波修正模块,用于第一极大值和第二极大值的采样间隔小于第二预设阈值,对R波进行修正,包括:将第一极大值和第二极大值中较大者对应的波确定为R波。
为了适应心电数据骤变的情况,在装置实施例1或实施例2的基础上,心率检测装置还可以包括:
第一阈值更新模块,用于在确定R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/4;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(3/8);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
第二阈值更新模块,用于在确定一个R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/2;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(7/16);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
本实施例不仅可以带来装置实施例1的有益效果,进一步的,对R波的判定进行了优化,并且合理设定了阈值的更新方式,可以避免由于心电数据骤变对R波判定的影响,提高心率计算结果的准确性。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (10)
1.一种心率检测方法,其特征在于,所述方法包括:
循环执行第一流程直至获取若干个R波;所述第一流程包括:采样心电数据以获取新的采样值,记所述新采样值的前一采样值为当前采样值;对当前采样值执行R波判定流程;
根据所述获取的若干个R波确定当前心率;
所述R波判定流程包括:若当前采样值为极大值,则标记极大值并返回执行第一流程;若当前采样值为极小值且极大值标记存在,则判断当前采样值和相邻的极大值采样间隔是否小于第一预设阈值;若小于,判定当前采样值相邻的极大值对应的波为R波;删除极大值标记;
其中,极大值的条件为大于左右相邻的采样值且大于波峰阈值;极小值的条件为小于左右相邻的采样值且小于波谷阈值。
2.根据权利要求1所述的方法,其特征在于,根据若干个R波确定当前心率,计算公式为:
H=Sample*60/Raverage;
其中,H表示心率,Sample为心电数据采样率,Raverage为由若干个R波确定的R波平均间隔。
3.根据权利要求1所述的方法,其特征在于,所述方法还包括:
记与当前采样值相邻的极大值为第一极大值,与第一极大值相邻的极大值为第二极大值;
当前采样值和第一极大值采样间隔小于第一预设阈值时,判断第一极大值和第二极大值的采样间隔是否小于第二预设阈值,若大于,则判定第一极大值对应的波为R波;若小于,则对R波进行修正;
所述对R波进行修正包括:将第一极大值和第二极大值中较大者对应的波确定为R波。
4.根据权利要求1-3任一项所述的方法,其特征在于,所述方法还包括:
在确定一个R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/4;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(3/8);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
5.根据权利要求1-3任一项所述的方法,其特征在于,所述方法还包括:
在确定一个R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/2;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(7/16);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
6.一种心率检测装置,其特征在于,所述装置包括:
执行模块,用于循环执行第一流程直至获取若干个R波;所述第一流程包括:采样心电数据以获取新的采样值,记所述新采样值的前一采样值为当前采样值;对当前采样值执行R波判定流程;
心率计算模块,用于根据所述获取的若干个R波确定当前时刻的心率;
R波判定模块,用于执行所述R波判定流程,所述R波判定流程包括:若当前采样值为极大值,则标记极大值并返回执行第一流程;若当前采样值为极小值且极大值标记存在,则判断当前采样值和相邻的极大值采样间隔是否小于第一预设阈值;若小于,判定当前采样值相邻的极大值对应的波为R波;删除极大值标记;
其中,极大值的条件为大于左右相邻的采样值且大于波峰阈值;极小值的条件为小于左右相邻的采样值且小于波谷阈值。
7.根据权利要求6所述的装置,其特征在于,心率计算模块确定心率的计算公式为:
H=Sample*60/Raverage;
其中,H表示心率,Sample为心电数据采样率,Raverage为由若干个R波确定的R波平均间隔。
8.根据权利要求6所述的装置,其特征在于,记与当前采样值相邻的极大值为第一极大值,与第一极大值相邻的极大值为第二极大值;所述装置还包括:
判定模块,用于当前采样值和第一极大值采样间隔小于第一预设阈值时,判断第一极大值和第二极大值的采样间隔是否小于第二预设阈值;
R波修正模块,用于第一极大值和第二极大值的采样间隔小于第二预设阈值,对R波进行修正,包括:将第一极大值和第二极大值中较大者对应的波确定为R波。
9.根据权利要求6-8任一项所述的装置,其特征在于,所述装置还包括:
第一阈值更新模块,用于在确定R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/4;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(3/8);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
10.根据权利要求6-8任一项所述的装置,其特征在于,所述装置还包括:
第二阈值更新模块,用于在确定一个R波之后,根据该R波对应的极大值和极小值更新波峰阈值和波谷阈值;
波峰阈值更新公式为:Maxthreshold=Min+(Max-Min)/2;
波谷阈值更新公式为:Minthreshold=Min+(Max-Min)*(7/16);
其中,Maxthreshold为波峰阈值,Minthreshold为波谷阈值,Min为R波对应的极小值,Max为R波对应的极大值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510455291.8A CN105125199B (zh) | 2015-07-29 | 2015-07-29 | 一种心率检测方法与装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510455291.8A CN105125199B (zh) | 2015-07-29 | 2015-07-29 | 一种心率检测方法与装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105125199A true CN105125199A (zh) | 2015-12-09 |
CN105125199B CN105125199B (zh) | 2017-12-26 |
Family
ID=54711050
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510455291.8A Active CN105125199B (zh) | 2015-07-29 | 2015-07-29 | 一种心率检测方法与装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105125199B (zh) |
Cited By (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106491118A (zh) * | 2016-09-30 | 2017-03-15 | 武汉吉星医疗科技有限公司 | 基于安卓系统的心电图机的实时心率计算系统及其方法 |
CN106725450A (zh) * | 2017-02-20 | 2017-05-31 | 武汉中旗生物医疗电子有限公司 | Qrs波群检测方法 |
CN107348957A (zh) * | 2016-05-10 | 2017-11-17 | 中国科学院微电子研究所 | 一种心电图r波检测方法及装置、便携式智能移动终端 |
WO2018082190A1 (zh) * | 2016-11-03 | 2018-05-11 | 华为技术有限公司 | 一种ecg信号处理方法及装置 |
TWI629049B (zh) * | 2017-04-26 | 2018-07-11 | 啟德電子股份有限公司 | A method for analyzing a heart shock signal for calculating a short-term heart rate value |
CN108433715A (zh) * | 2018-03-08 | 2018-08-24 | 芜湖泰领信息科技有限公司 | 智能心率预警方法 |
CN108451521A (zh) * | 2018-03-08 | 2018-08-28 | 芜湖泰领信息科技有限公司 | 智能心率预警系统 |
CN108597336A (zh) * | 2018-02-28 | 2018-09-28 | 天津天堰科技股份有限公司 | 心电波形仿真方法 |
CN108685561A (zh) * | 2017-04-10 | 2018-10-23 | 中国科学院微电子研究所 | 一种信号分析方法及装置 |
CN110931127A (zh) * | 2019-11-26 | 2020-03-27 | 佳禾智能科技股份有限公司 | 有效rri数值采集方法、智能穿戴电子设备、计算机可读存储介质 |
US11989804B2 (en) | 2017-11-06 | 2024-05-21 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for medical imaging |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003000561A (ja) * | 2001-06-18 | 2003-01-07 | Canon Inc | R波認識方法及びr−r間隔測定方法及び心拍数測定方法及びr−r間隔測定装置及び心拍数測定装置 |
CN1525395A (zh) * | 2003-02-24 | 2004-09-01 | 深圳迈瑞生物医疗电子股份有限公司 | 信号时域波形极值点和周期的检测方法 |
JP2009112625A (ja) * | 2007-11-08 | 2009-05-28 | Almedio Inc | 心拍測定装置 |
CN103549950A (zh) * | 2013-11-19 | 2014-02-05 | 上海理工大学 | 移动心电监测差分阈值检测改进算法 |
CN103690156A (zh) * | 2013-11-22 | 2014-04-02 | 东软熙康健康科技有限公司 | 一种心率获取方法及心电信号的处理方法 |
CN104586384A (zh) * | 2013-10-30 | 2015-05-06 | 深圳华清心仪医疗电子有限公司 | 一种心电信号实时心率检测方法及系统 |
KR20150081763A (ko) * | 2014-01-06 | 2015-07-15 | 계명대학교 산학협력단 | 심전도 신호의 저전력 고효율 r파 검출 방법 및 시스템 |
-
2015
- 2015-07-29 CN CN201510455291.8A patent/CN105125199B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2003000561A (ja) * | 2001-06-18 | 2003-01-07 | Canon Inc | R波認識方法及びr−r間隔測定方法及び心拍数測定方法及びr−r間隔測定装置及び心拍数測定装置 |
CN1525395A (zh) * | 2003-02-24 | 2004-09-01 | 深圳迈瑞生物医疗电子股份有限公司 | 信号时域波形极值点和周期的检测方法 |
JP2009112625A (ja) * | 2007-11-08 | 2009-05-28 | Almedio Inc | 心拍測定装置 |
CN104586384A (zh) * | 2013-10-30 | 2015-05-06 | 深圳华清心仪医疗电子有限公司 | 一种心电信号实时心率检测方法及系统 |
CN103549950A (zh) * | 2013-11-19 | 2014-02-05 | 上海理工大学 | 移动心电监测差分阈值检测改进算法 |
CN103690156A (zh) * | 2013-11-22 | 2014-04-02 | 东软熙康健康科技有限公司 | 一种心率获取方法及心电信号的处理方法 |
KR20150081763A (ko) * | 2014-01-06 | 2015-07-15 | 계명대학교 산학협력단 | 심전도 신호의 저전력 고효율 r파 검출 방법 및 시스템 |
Non-Patent Citations (1)
Title |
---|
田絮资,杨建,黄力宇: "基于数学形态学的心电信号R波检测算法", 《现代电子技术》 * |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107348957A (zh) * | 2016-05-10 | 2017-11-17 | 中国科学院微电子研究所 | 一种心电图r波检测方法及装置、便携式智能移动终端 |
CN106491118A (zh) * | 2016-09-30 | 2017-03-15 | 武汉吉星医疗科技有限公司 | 基于安卓系统的心电图机的实时心率计算系统及其方法 |
WO2018082190A1 (zh) * | 2016-11-03 | 2018-05-11 | 华为技术有限公司 | 一种ecg信号处理方法及装置 |
US11337638B2 (en) | 2016-11-03 | 2022-05-24 | Huawei Technologies Co., Ltd. | ECG signal processing method and apparatus |
CN106725450A (zh) * | 2017-02-20 | 2017-05-31 | 武汉中旗生物医疗电子有限公司 | Qrs波群检测方法 |
CN108685561A (zh) * | 2017-04-10 | 2018-10-23 | 中国科学院微电子研究所 | 一种信号分析方法及装置 |
TWI629049B (zh) * | 2017-04-26 | 2018-07-11 | 啟德電子股份有限公司 | A method for analyzing a heart shock signal for calculating a short-term heart rate value |
US11989804B2 (en) | 2017-11-06 | 2024-05-21 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for medical imaging |
CN108597336B (zh) * | 2018-02-28 | 2021-11-05 | 天津天堰科技股份有限公司 | 心电波形仿真方法 |
CN108597336A (zh) * | 2018-02-28 | 2018-09-28 | 天津天堰科技股份有限公司 | 心电波形仿真方法 |
CN108451521A (zh) * | 2018-03-08 | 2018-08-28 | 芜湖泰领信息科技有限公司 | 智能心率预警系统 |
CN108433715A (zh) * | 2018-03-08 | 2018-08-24 | 芜湖泰领信息科技有限公司 | 智能心率预警方法 |
CN110931127A (zh) * | 2019-11-26 | 2020-03-27 | 佳禾智能科技股份有限公司 | 有效rri数值采集方法、智能穿戴电子设备、计算机可读存储介质 |
CN110931127B (zh) * | 2019-11-26 | 2023-09-12 | 佳禾智能科技股份有限公司 | 有效rri数值采集方法、智能穿戴电子设备、计算机可读存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN105125199B (zh) | 2017-12-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105125199A (zh) | 一种心率检测方法与装置 | |
CN102835954B (zh) | 一种心拍波形模板生成方法及模块 | |
CN103690156B (zh) | 一种心率获取方法及心电信号的处理方法 | |
CN101828918B (zh) | 基于波形特征匹配的心电信号r波峰检测方法 | |
Lin et al. | P-and T-wave delineation in ECG signals using a Bayesian approach and a partially collapsed Gibbs sampler | |
Kathirvel et al. | An efficient R-peak detection based on new nonlinear transformation and first-order Gaussian differentiator | |
CN103584854B (zh) | 心电信号r波的提取方法 | |
CN105030228A (zh) | 在心电信号中确定其p波位置的方法及装置 | |
CN105212922A (zh) | 面向fpga实现心电信号r波自动检测的方法及系统 | |
US10172531B2 (en) | Heartbeat detection method and heartbeat detection device | |
CN103810393A (zh) | 一种基于心电向量的心电波形特征点定位方法 | |
CN101933811A (zh) | 一种血氧饱和度检测方法及系统 | |
CN102579039B (zh) | 一种心电图中的t波交替检测的方法 | |
CN104905785A (zh) | 室性/室上性早搏检测方法和系统 | |
CN104921722A (zh) | 一种双导联融合心电qrs波检测方法 | |
CN101940476B (zh) | 一种血氧饱和度检测方法及系统 | |
CN110236538B (zh) | 一种基于肌电信号的运动起始点自动实时检测方法 | |
CN116350238A (zh) | 一种自适应阈值r波检测方法、装置、设备及存储介质 | |
CN113440145A (zh) | 一种心电信号t波及其终点自动检测方法 | |
CN109009087A (zh) | 一种心电信号r波的快速检测方法 | |
CN111248895A (zh) | 心电信号特征检测方法及其系统 | |
CN107361763B (zh) | 一种心电图数据r波检测方法及装置 | |
CN105411579A (zh) | 一种心电图r波检测方法及装置 | |
CN109875550A (zh) | 一种心室除极关键点检测方法 | |
CN111053552A (zh) | 一种基于深度学习的qrs波检测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
CP03 | Change of name, title or address |
Address after: Room 1310, 13 / F, building 7, No. 37, Chaoqian Road, science and Technology Park, Changping District, Beijing 102200 Patentee after: Beijing Lepu Mingshi Technology Co.,Ltd. Address before: No. 37, Chaoqian Road, Changping District, Beijing 102200 Patentee before: BEIJING E-CARE CENTURY TECHNOLOGY Co.,Ltd. |
|
CP03 | Change of name, title or address |