CN103006209B - 心电信号检波方法及检波器 - Google Patents
心电信号检波方法及检波器 Download PDFInfo
- Publication number
- CN103006209B CN103006209B CN201210590401.8A CN201210590401A CN103006209B CN 103006209 B CN103006209 B CN 103006209B CN 201210590401 A CN201210590401 A CN 201210590401A CN 103006209 B CN103006209 B CN 103006209B
- Authority
- CN
- China
- Prior art keywords
- data
- ecgminl
- electrocardiogram
- ecg
- ecgmaxl
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 22
- 238000000718 qrs complex Methods 0.000 claims description 35
- 238000001514 detection method Methods 0.000 claims description 27
- 238000004458 analytical method Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 6
- 230000036651 mood Effects 0.000 claims description 2
- 230000002159 abnormal effect Effects 0.000 abstract description 3
- 239000002131 composite material Substances 0.000 abstract 2
- 238000002583 angiography Methods 0.000 abstract 1
- 238000007519 figuring Methods 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000001575 pathological effect Effects 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000000747 cardiac effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005284 excitation Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 238000011897 real-time detection Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000002861 ventricular Effects 0.000 description 1
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明涉及生物医学工程领域,提供一种心电信号检波方法,计算原始心电序列的差分绝对值序列,采用滑动窗对窗内差分绝对值求和,寻找大于特定值的位置区间,并在该区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL,根据利用EcgMaxL与EcgMinL大小关系、是否为avR导联,识别QRS复合波;还提供了一种心电信号检波器;本发明有效地避免了由于干扰等导致QRS复合波异常而不能准确识别,本发明简捷可行,准确性高,抗干扰能力强。
Description
技术领域
本发明涉及生物医学工程领域,特别涉及一种心电信号检波方法及检波器。
背景技术
理论上,如图1所示,一般一次心动周期就会在心电图上记录出一系列地高低宽窄不同的波形:包括P波、QRS复合波、T波和U波(U波有时没有)。
典型的QRS复合波包括三个相连的波动。第一个向下的波为Q波,继Q波后一个狭高向上的波为R波,与R波相连接的又一个向下的波为S波。由于这三个波紧密相连且总时间不超过0.10秒,故合称QRS波。QRS所占时间代表心室肌兴奋传播所需时间,正常人在0.06~0.10秒之间。
R检测是心电分析的基础与关键,现有关于R检测研究也比较热,如一种心电信号R波峰检测方法(ZL201010214626.4),通过对差向量进行对数极坐标转换来度量波形技术实现R波检测,又如心电图R波检测方法(CN101856225A),主要通过公式d(n)*d(n)*log(d(n)*d(n))对预处理后心电数据进行香农能量转换,在此基础上检测R波位置,还如基于波形特征匹配的心电信号R波峰检测方法:ZL 201010170490.1,通过模板匹配法来检测R波位置,这些方法对于一般的病理性不是很明显的心电信号检波,效果非常好,准确率可达99%,但是,一方面对于一些病理性强的QRS波或者QRS波群中R波可能不存在情况(如QS波),由于识别特征不再明显,现有检波算法误差较大,漏检和虚检比较严重,易误把QS波识别成R波,另一个方面,现有算法相对复杂,计算量大。并且,现在R波检波算法没有考虑aVR导联的特殊性,即与其他导联正常波形相比其波形是倒置的,适用范围受限。
发明内容
为解决以上问题,本发明提供了一种心电信号检波方法及检波器,根据QRS复合波识别特征,不需对心电数据进行基线漂移抑制和滤波等预处理,即能准确识别出QRS复合波,同时有效识别心脏停博信号,用于不同导联、不同采样率、不同数据长度的心电数据的QRS复合波实时和非实时检测。
本发明提供一种心电信号检波方法,包括:
获取原始心电序列EcgData;
计算原始心电序列的差分绝对值序列DiffAbs;
对滑动窗内差分绝对值求和得到窗口差分和序列DiffSumAbs;
在窗口差分和序列DiffSumAbs中,寻找大于DiffSumAbsMax×f1的第一个位置点,设为Is,寻找小于DiffSumAbsMax×f1的第一个位置点,设为Ie;其中f1为调节系数,在0.1–0.9之间取值;
在Is到Ie区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL;
利用EcgMaxL与EcgMinL大小关系、是否avR导联,识别QRS复合波。
本发明还提供一种心电信号检波器,包括分别与中央处理器相连接的存储器、差分绝对值计算单元、求和单元、比较判断单元、分析导联单元和复合波检测单元;
所述存储器用于存储数据,包括原始心电序列EcgData、原始心电序列的差分绝对值序列DiffAbs、窗口差分和序列DiffSumAbs;
所述差分绝对值计算单元用于计算原始心电序列的差分绝对值序列DiffAbs;
所述求和单元用于对序列值求和,包括对滑动窗内差分绝对值求和得到窗口差分和序列DiffSumAbs
所述比较判断单元用于对数据进行比较和判断,包括在窗口差分和序列DiffSumAbs中,寻找大于DiffSumAbsMax×f1的第一个位置点,设为Is,寻找小于DiffSumAbsMax×f1的第一个位置点,设为Ie;其中f1为调节系数,在0.1–0.9之间取值;在Is到Ie区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL;比较EcgMaxL与EcgMinL大小关系;
分析导联单元用于判断是否属于avR导联;
复合波检测单元,用于根据是否属于avR导联和EcgMaxL与EcgMinL大小关系识别QRS复合波
与现有技术相比,本发明计算原始心电序列的差分绝对值序列,采用滑动窗对窗内差分绝对值求和,寻找大于特定值的位置区间,并在该区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL,根据利用EcgMaxL与EcgMinL大小关系、avR导联,识别QRS复合波。本发明提供的心电信号检波方法和检波器简捷可行,准确性高,抗干扰能力强,具有广泛的应用场景。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简要介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域的普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为现有技术心电图示意图;
图2为本发明心电信号检波方法实施例流程示意图;
图3为本发明心电信号检波方法另一实施例流程示意图
图4为本发明心电信号检波器实施例结构示意图;
图5为本发明心电信号检波器另一实施例结构示意图。
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步地详细描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
作为本发明一个实施例,提供一种心电信号检波方法,如图2所示,包括:
301、获取原始心电序列EcgData;
所述获取原始心电序列数据EcgData可以通过心电采集装置采集,也可以通过来源于心电序列的离线数据,本发明对此不作限定。假设原始心电序列EcgData为一维数组,序列长度表示为SL,原始心电序列数据的采样率为F。
优选地,在301获取原始心电序列EcgData后,包括301A、无心博识别步骤,即求原始心电序列数据EcgData最大值,设为EcgMax,求原始心电序列数据EcgData最小值,设为EcgMin;如果EcgMax-EcgMin<0.1×Dv,所述Dv值在0.05-0.3毫伏(mv)之间,则认为无心博数据,结束检波分析;否则,进入步骤302。
302、计算原始心电序列的差分绝对值序列DiffAbs;
所述对原始心电序列数据求差分绝对值序列DiffAbs,可表示为:
DiffAbs[i]=abs(EcgData[i]-EcgData[i+1]),i=0,…,SL-2
303、采用滑动窗,对窗内差分序列绝对值求和,得到窗口差分和序列DiffSumAbs;
所述滑动窗口大小为WL,值为F×WLH秒,WLH为0.2-20秒;所述F为心电数据采样率,滑动步长为1;并把窗口中心数据下标作为差分和序列下标,其和值作为其值 j表示窗口起始数据下标;
本实施例WLH为0.2-20秒,主要是考虑分析的心电信号心率范围一般在20-300次范围之内,滑动窗口大小系数(60/300至60/30,即0.2-3),并经研究实验测试表明,当滑动窗口大小在(0.2-3)×F范围内时,检测效果最佳。
优选地,如图3所示,包括304、无心博重识别步骤,即计算特定区间内窗口差分和序列DiffSumAbs中的最大值DiffSumAbsMax与最小值DiffSumAbsMin,
所述特定区间为i到i+L区间,数据长度L值为F×LH秒,设起始点为i,开始分析时i=0,根据实验研究表明,LH值选取在2-10秒之间为佳;
如果DiffSumAbsMax-DiffSumAbsMin<Dv,则结束该特定区间的检波分析,i=i+L,若i<SL-2,回到该步骤之首,否则,结束检波分析;若DiffSumAbsMax-DiffSumAbsMin不小于Dv,继续步骤306;所述Dv值在0.05-0.3毫伏(mv)之间
306、在窗口差分和序列DiffSumAbs中,寻找大于DiffSumAbsMax×f1的第一个位置点,设为Is,寻找小于DiffSumAbsMax×f1的第一个位置点,设为Ie;所述f1为调节系数,大小在0.1–0.9之间;
如果DiffSumAbsMax-DiffSumAbsMin≧Dv,则从i开始到L之间,在窗口差分和序列DiffSumAbs中,依次寻找大于DiffSumAbsMax×f1的位置,当找到时记录第一个位置,设为Is,然后继续寻找小于DiffSumAbsMax×f1的位置,当找到后,记录第一个位置,设为Ie,或者在找到并记录大于DiffSumAbsMax×f1的第一个位置(设为Is)后,从Is往后跳一段数据,然后继续寻找小于DiffSumAbsMax×f1的位置,当找到后,记录第一个位置,设为Ie,后一种方式可以节省跳过的一段数据的检测,并且跳过的一段数据长度在0.05秒左右,因为一般QRS波群宽度不小于0.05秒,所以跳过这段数据不会造成检波错误,而且会加快了算法速度。
所述f1为调节系数,在0.1–0.9之间取值;
307、在Is到Ie区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL;
从Is开始,依次判断寻找到Ie区间内的最大值,并记录该最大值位置数据下标EcgMaxL,依次判断寻找到Ie区间内的最小值,并记录该最小值位置数据下标EcgMinL;
308、利用EcgMaxL与EcgMinL大小关系、avR导联,识别QRS复合波。
本发明计算原始心电序列的差分绝对值序列,采用滑动窗对窗内差分绝对值求和,寻找大于特定值的位置区间,并在该区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL,根据利用EcgMaxL与EcgMinL大小关系、是否是avR导联,识别QRS复合波,有效地避免了由于干扰等异常导致QRS复合波识别特征不明显导致的不能准确识别出QRS复合波的影响,经实验验证,本发明在干扰因素等异常情况情况下(如基线漂移很大),不需干扰预处理,也能够确保很高的识别准确率。
作为本发明一个实施例,所述利用EcgMaxL与EcgMinL大小关系、avR导联,识别QRS复合波包括,如果EcgMaxL<EcgMinL,若不是avR导联,则R波位置就是EcgMaxL到EcgMinL区间内的心电数据最大值点位置,若是avR导联,则R波位置就是EcgMaxL到EcgMinL区间内的心电数据最小值点位置,设为iR;
作为本发明一个实施例,所述利用EcgMaxL与EcgMinL大小关系、avR导联,识别QRS复合波包括,如果EcgMaxL>EcgMinL,则包括:
在EcgMinL到EcgMaxL区间内求心电数据最小值点位置,设为iQRS,然后以此点向前搜索到心电数据上升段顶点位置,设为iQRSs,以此点向后搜索到心电数据上升段顶点位置,设为iQRSe;根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测。
作为本发明一个实施例,所述根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测包括:
如果abs(EcgData[iQRSs]-EcgData[iQRSe])<tDv则R波消失,Q、S波重合,即QS波,若不是avR导联,则QS波顶点就是EcgMinL与EcgMaxL之间的心电数据最小值点位置,若是avR导联,则QS波顶点就是EcgMinL与EcgMaxL之间的心电数据最大值点位置,其中tDv值在0.1-0.9毫伏之间取值。
作为本发明一个实施例,所述根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测包括:
如果abs(EcgData[iQRSs]-EcgData[iQRSe])>=tDv,并且iQRSs点心电数据小于在iQRSe点心电数据,若不是avR导联,则Q波顶点就是EcgMinL与EcgMaxL之间的心电数据最小值点位置,R波顶点就是Is与EcgMaxL之间的心电数据离EcgMaxL最近的最大值点位置或者最大拐点位置,S波顶点就是EcgMinL与Ie之间的心电数据最小值点位置或者最大拐点位置,若是avR导联,则Q波顶点就是EcgMinL与EcgMaxL之间的心电数据最大值点位置,R波顶点就是Is与EcgMaxL之间的心电数据离EcgMaxL最近的最小值点位置或者最大拐点位置,S波顶点就是EcgMinL与Ie之间的心电数据最大值点位置或者最大拐点位置;
作为本发明一个实施例,所述根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测包括:
如果abs(EcgData[iQRSs]-EcgData[iQRSe])>=tDv并且在iQRSs点心电数据大于在iQRSe点心电数据,若是avR导联,则Q波顶点就是Is与EcgMaxL之间的心电数据最小值点位置,R波顶点就是Is与EcgMinL之间的心电数据离EcgMinL最近的最大值点位置或者最大拐点位置,S波顶点就是EcgMinL与Ie之间的心电数据最小值点位置或者最大拐点位置,若不是avR导联,则Q波顶点就是Is与EcgMaxL之间的心电数据最大值点位置,R波顶点就是Is与EcgMinL之间的心电数据离EcgMinL最近的最小值点位置或者最大拐点位置,S波顶点就是EcgMinL与Ie之间的心电数据最大值点位置或者最大拐点位置,其中tDv值在0.1-0.9mv之间取值。
本发明提供一种心电信号检波器,如图4所示,包括分别与中央处理器相连接的存储器、差分绝对值计算单元、求和单元、比较判断单元、分析导联单元和复合波检测单元;
所述中央处理器用于调度处理,具体调度过程见上述步骤301-308过程的描述,不再赘述。
所述存储器用于存储数据,包括原始心电序列EcgData、原始心电序列的差分绝对值序列DiffAbs、窗口差分和序列DiffSumAbs;
所述差分绝对值计算单元用于计算原始心电序列的差分绝对值序列DiffAbs;
所述求和单元用于对序列值求和,包括对滑动窗内差分绝对值求和得到窗口差分和序列DiffSumAbs
所述比较判断单元用于对数据进行比较和判断,包括在窗口差分和序列DiffSumAbs中,寻找大于DiffSumAbsMax×f1的第一个位置点,设为Is,寻找小于DiffSumAbsMax×f1的第一个位置点,设为Ie;其中f1为调节系数,在0.1-0.9之间取值;在Is到Ie区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL;比较EcgMaxL与EcgMinL大小关系;
所述分析导联单元用于判断是否属于avR导联,判断结果供复合波检测单元使用。
所述复合波检测单元,用于根据是否属于avR导联和EcgMaxL与EcgMinL大小关系识别QRS复合波。
本发明差分绝对值计算单元计算原始心电序列的差分绝对值序列,求和单元采用滑动窗对窗内差分绝对值求和,比较判断单元寻找大于特定值的位置区间,并在该区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL,复合波检测单元利用EcgMaxL与EcgMinL大小关系、avR导联,识别QRS复合波,有效地避免了由于干扰等异常导致QRS复合波识别特征不明显导致的不能准确识别出QRS复合波的影响,经实验验证,本发明在干扰因素等异常情况下,能够确保较成的识别成功率。
优选地,如图5所示,所述心电信号检波器还包括最值计算单元,与中央处理器相连接;
所述最值计算单元用于计算一段数据的最大值或/和最小值,包括计算特定区间内窗口差分和序列DiffSumAbs中的最大值DiffSumAbsMax与最小值DiffSumAbsMin。
作为本发明一个实施例,所述复合波检测单元根据是否属于avR导联和EcgMaxL与EcgMinL大小关系识别QRS复合波,包括:
如果EcgMaxL<EcgMinL,若不是avR导联,则R波位置就是EcgMaxL到EcgMinL区间内的心电数据最大值点位置,若是avR导联,则R波位置就是EcgMaxL到EcgMinL区间内的心电数据最小值点位置,设为iR;
作为本发明一个实施例,所述复合波检测单元根据是否属于avR导联和EcgMaxL与EcgMinL大小关系识别QRS复合波,包括:如果EcgMaxL>EcgMinL,则包括:
在EcgMinL到EcgMaxL区间内求心电数据最小值点位置,设为iQRS,然后以此点向前搜索到心电数据上升段顶点位置,设为iQRSs,以此点向后搜索到心电数据上升段顶点位置,设为iQRSe;根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测。
作为本发明一个实施例,所述根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测包括:
如果abs(EcgData[iQRSs]-EcgData[iQRSe])<tDv则R波消失,Q、S波重合,即QS波,若不是avR导联,则QS波顶点就是EcgMinL与EcgMaxL之间的心电数据最小值点位置,若是avR导联,则QS波顶点就是EcgMinL与EcgMaxL之间的心电数据最大值点位置,其中tDv值在0.1-0.9毫伏之间取值。
作为本发明一个实施例,所述根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测包括:
如果abs(EcgData[iQRSs]-EcgData[iQRSe])>=tDv,并且iQRSs点心电数据小于在iQRSe点心电数据,若不是avR导联,则Q波顶点就是EcgMinL与EcgMaxL之间的心电数据最小值点位置,R波顶点就是Is与EcgMaxL之间的心电数据离EcgMaxL最近的最大值点位置或者最大拐点位置,S波顶点就是EcgMinL与Ie之间的心电数据最小值点位置或者最大拐点位置,若分析导联是avR导联,则Q波顶点就是EcgMinL与EcgMaxL之间的心电数据最大值点位置,R波顶点就是Is与EcgMaxL之间的心电数据离EcgMaxL最近的最小值点位置或者最大拐点位置,S波顶点就是EcgMinL与Ie之间的心电数据最大值点位置或者最大拐点位置;
作为本发明一个实施例,所述根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测包括:
如果abs(EcgData[iQRSs]-EcgData[iQRSe])>=tDv并且在iQRSs点心电数据大于在iQRSe点心电数据,若分析导联不是avR导联,则Q波顶点就是Is与EcgMaxL之间的心电数据最小值点位置,R波顶点就是Is与EcgMinL之间的心电数据离EcgMinL最近的最大值点位置或者最大拐点位置,S波顶点就是EcgMinL与Ie之间的心电数据最小值点位置或者最大拐点位置,若分析导联不是avR导联,则Q波顶点就是Is与EcgMaxL之间的心电数据最大值点位置,R波顶点就是Is与EcgMinL之间的心电数据离EcgMinL最近的最小值点位置或者最大拐点位置,S波顶点就是EcgMinL与Ie之间的心电数据最大值点位置或者最大拐点位置,其中tDv值在0.1-0.9mv之间取值。
以上所举实施例,对本发明的目的、技术方案和优点进行了进一步的详细说明,所应理解的是,以上所举实施例仅为本发明的优选实施方式而已,并不用以限制本发明,凡在本发明的精神和原则之内对本发明所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种心电信号检波方法,其特征在于:包括:
获取原始心电序列EcgData;
计算原始心电序列的差分绝对值序列DiffAbs;
以某一窗宽为单位,对在该窗口内的序列差分绝对值序列求和,并不断滑动窗口,形成窗口差分和序列DiffSumAbs;
在窗口差分和序列DiffSumAbs中,寻找大于DiffSumAbsMax×f1的第一个位置点,设为Is,寻找小于DiffSumAbsMax×f1的第一个位置点,设为Ie;其中f1为调节系数,在0.1–0.9之间取值;所述DiffSumAbsMax表示窗口差分和序列DiffSumAbs中的最大值;
在Is到Ie区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL;
利用EcgMaxL与EcgMinL大小关系、是否avR导联,识别QRS复合波,包括:
如果EcgMaxL<EcgMinL,若不是avR导联,则R波位置就是EcgMaxL到EcgMinL区间内的心电数据最大值点位置,若是avR导联,则R波位置就是EcgMaxL到EcgMinL区间内的心电数据最小值点位置;
如果EcgMaxL>EcgMinL,在EcgMinL到EcgMaxL区间内求心电数据最小值点位置,设为iQRS,以此点向前搜索到心电数据上升段顶点位置,设为iQRSs,以此点向后搜索到心电数据上升段顶点位置,设为iQRSe;根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测。
2.根据权利要求1所述心电信号检波方法,其特征在于:在所述获取原始心电序列之后,在所述计算原始心电序列的差分绝对值序列DiffAbs之前,包括无心博识别步骤,包括:求原始心电序列数据EcgData最大值,设为EcgMax,求原始心电序列数据EcgData最小值,设为EcgMin;如果EcgMax-EcgMin<0.1×Dv,则认为无心博数据,结束检波分析;所述Dv值在0.05-0.3 mv之间。
3.根据权利要求1所述心电信号检波方法,其特征在于:在所述对滑动窗内差分绝对值求和得到窗口差分和序列DiffSumAbs之后,在所述在窗口差分和序列DiffSumAbs中,寻找大于DiffSumAbsMax×f1的第一个位置点之前,包括无心博重识别步骤,包括计算特定区间内窗口差分和序列DiffSumAbs中的最大值DiffSumAbsMax与最小值DiffSumAbsMin,判断若DiffSumAbsMax-DiffSumAbsMin<Dv,则结束该特定区间的检波分析,所述Dv值在0.05-0.3mv之间。
4.根据权利要求1所述心电信号检波方法,其特征在于:所述根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测,包括:
如果abs(EcgData[iQRSs]-EcgData[iQRSe])<tDv,则R波消失,Q、S波重合,即QS波,进一步地,若不是avR导联,则QS波顶点就是EcgMinL与EcgMaxL之间的心电数据最小值点位置,若是avR导联,则QS波顶点就是EcgMinL与EcgMaxL之间的心电数据最大值点位置,其中tDv值在0.1-0.9毫伏之间取值。
5.根据权利要求1所述心电信号检波方法,其特征在于:所述根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测,包括:
如果abs(EcgData[iQRSs]-EcgData[iQRSe])>=tDv并且iQRSs点心电数据小于在iQRSe点心电数据,进一步地,若不是avR导联,则Q波顶点就是EcgMinL与EcgMaxL之间的心电数据最小值点位置,R波顶点就是Is与EcgMaxL之间的心电数据离EcgMaxL最近的最大值点位置,S波顶点就是EcgMinL与Ie之间的心电数据最小值点位置,若分析导联是avR导联,则Q波顶点就是EcgMinL与EcgMaxL之间的心电数据最大值点位置,R波顶点就是Is与EcgMaxL之间的心电数据离EcgMaxL最近的最小值点位置,S波顶点就是EcgMinL与Ie之间的心电数据最大值点位置;其中tDv值在0.1-0.9mv之间取值。
6.根据权利要求1所述心电信号检波方法,其特征在于:所述根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测,包括:
如果abs(EcgData[iQRSs]-EcgData[iQRSe])>=tDv并且在iQRSs点心电数据大于在iQRSe点心电数据,进一步地,若分析导联不是avR导联,则Q波顶点就是Is与EcgMaxL之间的心电数据最小值点位置,R波顶点就是Is与EcgMinL之间的心电数据离EcgMinL最近的最大值点位置,S波顶点就是EcgMinL与Ie之间的心电数据最小值点位置,若分析导联不是avR导联,则Q波顶点就是Is与EcgMaxL之间的心电数据最大值点位置,R波顶点就是Is与EcgMinL之间的心电数据离EcgMinL最近的最小值点位置,S波顶点就是EcgMinL与Ie之间的心电数据最大值点位置,其中tDv值在0.1-0.9mv之间取值。
7.一种心电信号检波器,其特征在于:包括分别与中央处理器相连接的存储器、差分绝对值计算单元、求和单元、比较判断单元、分析导联单元和复合波检测单元;
所述存储器用于存储数据,包括原始心电序列EcgData、原始心电序列的差分绝对值序列DiffAbs、窗口差分和序列DiffSumAbs;
所述差分绝对值计算单元用于计算原始心电序列的差分绝对值序列DiffAbs;
所述求和单元用于对序列值求和,包括以某一窗宽为单位,对在该窗口内的序列差分绝对值序列求和,并不断滑动窗口,形成窗口差分和序列DiffSumAbs;;
所述比较判断单元用于对数据进行比较和判断,包括在窗口差分和序列DiffSumAbs中,寻找大于DiffSumAbsMax×f1的第一个位置点,设为Is,寻找小于DiffSumAbsMax×f1的第一个位置点,设为Ie;其中f1为调节系数,在0.1–0.9之间取值;在Is到Ie区间内寻找原始数据最大值点位置EcgMaxL和最小值点位置EcgMinL;比较EcgMaxL与EcgMinL大小关系;所述DiffSumAbsMax表示窗口差分和序列DiffSumAbs中的最大值;
分析导联单元用于判断是否属于avR导联;
复合波检测单元,用于根据是否属于avR导联、EcgMaxL与EcgMinL大小关系识别QRS复合波,包括:
如果EcgMaxL<EcgMinL,若不是avR导联,则R波位置就是EcgMaxL到EcgMinL区间内的心电数据最大值点位置,若是avR导联,则R波位置就是EcgMaxL到EcgMinL区间内的心电数据最小值点位置;
如果EcgMaxL>EcgMinL,在EcgMinL到EcgMaxL区间内求心电数据最小值点位置,设为iQRS,然后以此点向前搜索到心电数据上升段顶点位置,设为iQRSs,以此点向后搜索到心电数据上升段顶点位置,设为iQRSe;根据iQRSs与iQRSe处原始心电数据的大小关系进行QRS复合波的检测。
8.根据权利要求7所述心电信号检波器,其特征在于:还包括最值计算单元,与中央处理器相连接,用于计算一段数据的最大值或/和最小值,包括计算特定区间内窗口差分和序列DiffSumAbs中的最大值与最小值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210590401.8A CN103006209B (zh) | 2012-12-29 | 2012-12-29 | 心电信号检波方法及检波器 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210590401.8A CN103006209B (zh) | 2012-12-29 | 2012-12-29 | 心电信号检波方法及检波器 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103006209A CN103006209A (zh) | 2013-04-03 |
CN103006209B true CN103006209B (zh) | 2014-10-22 |
Family
ID=47955706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210590401.8A Expired - Fee Related CN103006209B (zh) | 2012-12-29 | 2012-12-29 | 心电信号检波方法及检波器 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103006209B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103705234B (zh) * | 2013-12-05 | 2015-08-26 | 深圳先进技术研究院 | 动态心电信号数据中的检波方法和装置 |
TWI539931B (zh) * | 2014-01-03 | 2016-07-01 | 緯創資通股份有限公司 | 心電圖量測方法與裝置 |
CN107303177B (zh) * | 2016-04-15 | 2020-04-14 | 中国科学院微电子研究所 | 一种心电图t波、p波的检测方法及系统 |
CN105877742B (zh) * | 2016-05-18 | 2019-03-22 | 四川长虹电器股份有限公司 | 一种针对心电信号ⅰ导联电极的错放检测方法 |
US11138770B2 (en) | 2017-11-06 | 2021-10-05 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for medical imaging |
CN110432894B (zh) * | 2019-08-09 | 2022-04-22 | 上海鹰瞳医疗科技有限公司 | 心电图关键点标注方法及电子设备 |
CN111513706B (zh) * | 2020-04-20 | 2022-09-27 | 重庆邮电大学 | 一种针对含有异常r波的心电信号的检测方法和装置 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101467879A (zh) * | 2007-12-26 | 2009-07-01 | 深圳迈瑞生物医疗电子股份有限公司 | 一种基于多导同步心电信号处理方法及装置 |
CN101522099A (zh) * | 2006-08-03 | 2009-09-02 | Bsp生物信号处理有限公司 | 利用高频qrs电位分析识别心肌缺血的装置和方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20090048528A1 (en) * | 2007-08-16 | 2009-02-19 | Bruce Hopenfeld | System and methods for detecting ischemia with a limited extracardiac lead set |
US8855752B2 (en) * | 2010-02-09 | 2014-10-07 | Angel Medical Systems, Inc. | QRS onset and offset detection with adaptive temporal windowing |
-
2012
- 2012-12-29 CN CN201210590401.8A patent/CN103006209B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101522099A (zh) * | 2006-08-03 | 2009-09-02 | Bsp生物信号处理有限公司 | 利用高频qrs电位分析识别心肌缺血的装置和方法 |
CN101467879A (zh) * | 2007-12-26 | 2009-07-01 | 深圳迈瑞生物医疗电子股份有限公司 | 一种基于多导同步心电信号处理方法及装置 |
Non-Patent Citations (3)
Title |
---|
一种新型心电信号采集及分析系统;王恩美;《仪器仪表学报》;20010831;第22卷(第4期);368-370 * |
王恩美.一种新型心电信号采集及分析系统.《仪器仪表学报》.2001,第22卷(第4期),368-370. |
田杰.基于Cortex的实时心电处理及检测技术.《《中国优秀硕士学位论文全文数据库,信息科技特辑》》.2011,(第8期),50-60. * |
Also Published As
Publication number | Publication date |
---|---|
CN103006209A (zh) | 2013-04-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103006209B (zh) | 心电信号检波方法及检波器 | |
CN101856225B (zh) | 一种心电信号r波峰检测方法 | |
WO1998009241A1 (en) | Ecg p-qrs-t onset, offset and peak detection method and apparatus | |
Chatterjee et al. | Real time P and T wave detection from ECG using FPGA | |
CN104173064A (zh) | 基于心率变异分析的测谎方法及测谎装置 | |
CN103810393A (zh) | 一种基于心电向量的心电波形特征点定位方法 | |
CN109171787A (zh) | 脉冲信号的采样方法、装置和计算机程序介质 | |
CN107622259B (zh) | 一种t波检测方法、心电数据分析方法及装置 | |
CN103705234A (zh) | 动态心电信号数据中的检波方法和装置 | |
Plesinger et al. | Robust multichannel QRS detection | |
Salih et al. | A novel approach for detecting QRS complex of ECG signal | |
Sadhukhan et al. | Detection of ECG characteristic features using slope thresholding and relative magnitude comparison | |
CN112952827A (zh) | 一种精确识别电动自行车充电的非侵入式全负荷识别技术 | |
CN111190049B (zh) | 主分量分析的混沌系统检测纳伏级微弱正弦信号的方法 | |
Wartak et al. | Computer program for pattern recognition of electrocardiograms | |
Arafat et al. | Automatic detection of ECG wave boundaries using empirical mode decomposition | |
Narsimha et al. | Denoising and QRS detection of ECG signals using Empirical Mode Decomposition | |
Espiritu-Santo-Rincon et al. | ECG feature extraction via waveform segmentation | |
Zhang et al. | Novel QRS detection by CWT for ECG sensor | |
Pal et al. | QRS Complex detection using Empirical Mode Decomposition based windowing technique | |
Li et al. | A novel abnormal ECG beats detection method | |
US8706203B2 (en) | Classification of a signal in a time domain | |
Rodríguez et al. | Adaptive threshold and principal component analysis for features extraction of electrocardiogram signals | |
CN114073534B (zh) | 一种心功能分析算法 | |
Huang et al. | Detecting QRS complexes of two-channel ECG signals by using combined wavelet entropy |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20141022 |