CN105232027A - 一种便携式心电信号处理方法和心电信号处理装置 - Google Patents
一种便携式心电信号处理方法和心电信号处理装置 Download PDFInfo
- Publication number
- CN105232027A CN105232027A CN201410253367.4A CN201410253367A CN105232027A CN 105232027 A CN105232027 A CN 105232027A CN 201410253367 A CN201410253367 A CN 201410253367A CN 105232027 A CN105232027 A CN 105232027A
- Authority
- CN
- China
- Prior art keywords
- ripple
- electrocardiosignal
- threshold value
- filtering
- slope
- 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)
Abstract
本发明公开了一种便携式心电信号处理方法和心电信号处理装置,其中,心电信号处理方法包括:通过IIR数字滤波器对输入的心电信号进行滤波,得到第一次滤波后的心电信号;根据需要滤除噪声的频率计算取样周期;若第一次滤波后的心电信号在取样周期内中的最大值和最小值所对应的两点间的直线的斜率小于或等于预设的滤波的斜率阈值,则将第一次滤波后的心电信号输出为心电图信号;若该斜率大于滤波的斜率阈值,则通过滤波器对第一次滤波后的心电信号进行至少一次滤波,直至得到第二心电信号并输出,其中,第二心电信号在取样周期内的最大值和最小值所对应的两点间的直线的斜率小于滤波的斜率阈值。本发明提供的技术方案能够有效消除心电信号中的噪声。
Description
技术领域
本发明涉及医学信号处理领域,具体涉及一种便携式心电信号处理方法和心电信号处理装置。
背景技术
心电信号是一种电压信号,它的伏值一般在0.5毫伏至5毫伏之间,相对于外界的工频干扰等信号,心电信号属于非常微弱的信号,因此,若不对心电信号进行滤波处理,则原始的心电信号(即不加处理的心电信号)容易淹没在噪声中,从而严重影响基于心电信号的医疗诊断的实行。
目前常采用平滑滤波技术对心电信号进行处理,平滑滤波是低频增强的空间域滤波技术,一般采用简单平均法进行,就是求邻近像元点的平均亮度值,以此消除心电信号中的噪声。
虽然平滑滤波技术能够较好的消除心电信号中的噪声,然而,平滑滤波技术也存在明显不足,例如,由于心电信号的噪声强度随环境的变化而变化,而平滑滤波技术并不能很好的适应环境的变化,这使得在某些环境(如噪声强度较大的环境)下,采用平滑滤波技术并不能有效滤除心电信号的噪声。
发明内容
本发明提供一种便携式心电信号处理方法和心电信号处理装置,用于消除心电信号中的噪声。
本发明第一方面提供一种心电信号处理方法,包括:
通过无限脉冲响应(IIR,InfiniteImpulseResponse)数字滤波器对输入的心电信号进行滤波,得到第一次滤波后的心电信号;
根据上述第一次滤波后的心电信号的波形频率计算取样周期,其中,上述取样周期为上述第一次滤波后的心电信号的波形频率的倒数;
根据上述第一次滤波后的心电信号在上述取样周期内中的最大值和最小值,计算上述最大值和上述最小值所对应的两点间的直线的斜率;
若上述斜率小于或等于预设的滤波的斜率阈值,则将上述第一次滤波后的心电信号输出为心电图信号;
若上述斜率大于上述滤波的斜率阈值,则通过低通滤波器对上述第一次滤波后的心电信号进行至少一次滤波,直至得到第二心电信号,其中,上述第二心电信号在上述取样周期内的最大值和最小值所对应的两点间的直线的斜率小于或等于上述滤波的斜率阈值;
将上述第二心电信号输出为心电图信号。
本发明第一方面提供一种心电信号处理装置,包括:
第一滤波单元,用于通过无限脉冲响应IIR数字滤波器对输入的心电信号进行滤波,得到第一次滤波后的心电信号;
第一计算单元,用于根据上述第一次滤波后的心电信号的波形频率计算取样周期,其中,上述取样周期为上述第一次滤波后的心电信号的波形频率的倒数;
第二计算单元,用于根据上述第一次滤波后的心电信号在上述取样周期内中的最大值和最小值,计算上述最大值和上述最小值所对应的两点间的直线的斜率;
输出单元,用于当上述第二计算单元计算得到的斜率小于或等于预设的滤波的斜率阈值时,将上述第一次滤波后的心电信号输出为心电图信号;
第二滤波单元,用于当上述第二计算单元计算得到的斜率大于上述滤波的斜率阈值时,通过低通滤波器对上述第一次滤波后的心电信号进行至少一次滤波,直至得到第二心电信号,其中,上述第二心电信号在上述取样周期内的最大值和最小值所对应的两点间的直线的斜率小于或等于上述滤波的斜率阈值;
上述输出单元还用于将上述第二滤波单元滤波得到的上述第二心电信号输出为心电图信号。
从本发明上述技术方案可知,本发明提供了一种复合滤波方案,首先通过IIR数字滤波器对心电信号进行滤波,根据经IIR滤波后得到的第一次滤波后的心电信号在取样周期内的最大值和最小值确定斜率,然后根据该斜率与滤波的斜率阈值的关系判定是否需要通过低通滤波器对第一次滤波后的心电信号进行进一步滤波,一方面,通过IIR滤波器和低通滤波器的复合滤波,能够有效消除心电信号中的噪声,另一方面,上述滤波的斜率阈值可以根据不同的环境进行设定,从而实现在不同的环境下对心电信号的噪声的有效滤除。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的一种便携式心电信号处理方法一个实施例流程示意图;
图2-a为本发明提供的一种心电图信号场景示意图;
图2-b为本发明提供的一种便携式心电信号处理方法另一个实施例流程示意图;
图2-c为本发明提供的一种R波定位方法实施例流程示意图;
图3为本发明提供的一种便携式心电信号处理装置一个实施例结构示意图。
具体实施方式
为使得本发明的发明目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而非全部实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明实施例提供一种便携式心电信号处理方法,上述方法包括:通过IIR数字滤波器对输入的心电信号进行滤波,得到第一次滤波后的心电信号;根据上述第一次滤波后的心电信号的波形频率计算取样周期,其中,上述取样周期为上述第一次滤波后的心电信号的波形频率的倒数;根据上述第一次滤波后的心电信号在上述取样周期内中的最大值和最小值,计算上述最大值和上述最小值所对应的两点间的直线的斜率;若上述斜率小于或等于预设的滤波的斜率阈值,则将上述第一次滤波后的心电信号输出为心电图信号;若上述斜率大于上述滤波的斜率阈值,则通过低通滤波器对上述第一次滤波后的心电信号进行至少一次滤波,直至得到第二心电信号,其中,上述第二心电信号在上述取样周期内的最大值和最小值所对应的两点间的直线的斜率小于或等于上述滤波的斜率阈值;将上述第二心电信号输出为心电图信号。本发明实施例还提供相应的心电信号处理装置。以下分别进行详细说明。
下面对本发明实施例中提供的一种便携式心电信号处理方法进行描述,请参阅图1,本发明实施例中的便携式心电信号处理方法包括:
101、通过IIR数字滤波器对输入的心电信号进行滤波,得到第一次滤波后的心电信号;
首先说明的是,本发明实施例中的心电信号是指心脏活动时所产生的生理电信号。
可选的,本发明实施例中的IIR数字滤波器的设计方法为:首先,在单位圆上均匀分布零点,其中,零点是指当输入幅度不为零且输入频率使输出幅值为零的点,即,零点表明与该零点所在位置相对应的频率处的幅值完全衰减;其次,选择位于单位圆上的极点,精确抵消单位圆上的一些零点,使得该零点所在位置相对应的频率处的幅值不会衰减,其中,极点是指当输入幅度不为零且输入频率使输出幅值为无穷大的点。由于单位圆上的每一点代表不同的频率,因此,零-极点位置决定了频率响应,IIR数字滤波器为递归滤波器的特定形式。当然,本发明实施例也可以采用IIR数字滤波器的常规设计方案对输入的心电信号进行滤波,此处不做限定。
102、根据上述第一次滤波后的心电信号的波形频率计算取样周期,其中,上述取样周期为上述第一次滤波后的心电信号的波形频率的倒数。
103、根据上述第一次滤波后的心电信号在上述取样周期内中的最大值和最小值,计算上述最大值和上述最小值所对应的两点间的直线的斜率;
本发明实施例中,心电信号处理装置根据上述第一次滤波后的心电信号在上述取样周期内中的最大值和最小值,根据直线斜率计算公式求出上述最大值对应的点与上述最小值对应的点连接形成的直线的斜率。
104、判断上述斜率是否小于或等于预设的滤波的斜率阈值;
本发明实施例中,上述滤波的斜率阈值可以根据实际环境进行设定,以适应不同环境的变化。例如,对于信号质量较高的原始心电信号(即为经过滤波处理的心电信号),上述滤波的斜率阈值可以设为较大值,对于信号质量较低的原始心电信号,上述滤波的斜率阈值可以设为较小值,即,上述滤波的斜率阈值的取值可以与原始心电信号的信号质量成正相关。
心电信号处理装置判断上述斜率是否小于或等于预置的滤波的斜率阈值,若是,则执行步骤105,若否,则执行步骤106。
105、将上述第一次滤波后的心电信号输出为心电图信号。
106、通过低通滤波器对上述第一次滤波后的心电信号进行至少一次滤波,直至得到第二心电信号;
其中,上述第二心电信号在上述取样周期内的最大值和最小值所对应的两点间的直线的斜率小于或等于上述滤波的斜率阈值。
本发明实施例中,心电信号处理装置通过低通滤波器(例如有限长单位冲激响应(FIR,FiniteImpulseResponse)滤波器或其它衰减较小的滤波器)对上述第一次滤波后的心电信号进行滤波,之后计算滤波后的第一次滤波后的心电信号(假设为心电信号A)在上述取样周期内的最大值和最小值所对应的两点间的连线是否小于等于上述滤波的斜率阈值,如果小于或等于上述滤波的斜率阈值,则将心电信号A确定为第二心电信号,如果大于上述滤波的斜率阈值,则通过上述低通滤波器对心电信号A进一步进行滤波,之后计算滤波后的心电信号A(假设为心电信号B)在上述取样周期内的最大值和最小值所对应的两点间的连线是否小于或等于上述滤波的斜率阈值,如果心电信号B小于或等于上述滤波的斜率阈值,则将心电信号B确定为第二心电信号,如果心电信号B大于上述滤波的斜率阈值,则通过上述低通滤波器对心电信号B进一步进行滤波,如此反复,直至得到第二心电信号。
107、将上述第二心电信号输出为心电图信号。
从本发明上述技术方案可知,本发明提供了一种复合滤波方案,首先通过IIR数字滤波器对心电信号进行滤波,根据经IIR滤波后得到的第一次滤波后的心电信号在取样周期内的最大值和最小值确定斜率,然后根据该斜率与滤波的斜率阈值的关系判定是否需要通过低通滤波器对第一次滤波后的心电信号进行进一步滤波,一方面,通过IIR滤波器和低通滤波器的复合滤波,能够有效消除心电信号中的噪声,另一方面,上述滤波的斜率阈值可以根据不同的环境进行设定,从而实现在不同的环境下对心电信号的噪声的有效滤除。
进一步,在图1所示实施例的基础上,本发明实施例还提供对心电图信号的QRS波群定位的方案,首先对QRS波群进行说明,QRS波群是心电图信号中电压幅值最高的一段波,由左、右心室除极而形成,主要用于反映电刺激自房室结传递到左、右心室,进而引起左、右心室收缩的过程。它可以细分为两个波峰朝下的Q波和S波,以及心电信号幅度最大的R波共三个小波,如图2-a所示,从发生时间看,首先出现的是低于等电位线波峰朝下的Q波,然后是整个心电图信号中幅度最大波峰朝上的R波,最后是和Q波同样方向的S波,通常,QRS波群的频率在7赫兹至27赫兹之间,电压变化幅度为5毫伏以下。下面对本发明实施例中的心电信号处理方法进一步描述,请参阅图2-b,本发明实施例中的心电信号处理方法包括:
步骤201~207,其中,步骤201~207与图1所示实施例中的步骤101~107相似,步骤201~207的具体实现过程可以参见步骤101~107中的描述,此处不再赘述。
208、对上述心电图信号进行QRS波群定位。
经过步骤201~207的滤波后,基本可以滤除50赫兹工频干扰和主要由呼吸引起的低频基线漂移,得到较干净和平滑的心电图信号,心电信号处理装置进一步对上述心电图信号进行QRS波群定位。
由于人的心跳每秒通常是40~120次,因此,正常情况下,1.5秒内应该存在一个R波,因此,本发明实施例的一种实现方式中,将一连串的心电图信号按照1.5秒为一个检测周期动态查找检测周期内的最高点,确定每个检测周期中的R波,进一步,由于人的QRS期间(即Q点和S点的时间间期)在150毫秒以内,而Q波、S波的方向与R波正好相反,因此,在每个上述检测周期中,将时间轴上位于R波波峰前75毫秒的位置所对应的信号点确定为S点,将时间轴上位于R波波峰后75毫秒的位置所对应的信号点确定为Q点。
可选地,确定每个检测周期中的R波的一种实现方式如图2-c所示,包括:
A1、根据第一公式初始化R波阈值;
其中,上述第一公式为:
thinitial=(max[0]+max[1]+...+max[N-1])/(2*N),式中thinitial表示初始化的R波阈值,max[N-1]表示在上述心电图信号的第N个检测周期内的最大值,其中,上述N为大于或等于2的自然数,可选地,N取4。
A2、在上述心电图信号的首个检测周期中:将大于初始化的R波阈值的信号序列确定为上述首个检测周期中的R波,并根据第二公式更新R波阈值;
其中,上述第二公式为:
thnew=0.8*tholder+0.2*y(max)last/2,式中thnew表示更新后的R波阈值,tholder表示在当前检测周期之前用于确定R波所使用的R波阈值(在步骤A2中,tholder为初始化的的R波阈值),y(max)last表示检测周期中的R波的波峰值。
A3、在上述心电图信号后续的每个检测周期中:将大于更新后的R波阈值的信号序列确定为当前检测周期中的R波,并根据上述第二公式更新R波阈值。
进一步,心电信号处理装置还可以根据以下原则剔除误判的R波和调整R波阈值:
情况一:心室收缩后有200毫秒的不应期,因此这段时间内心脏不可能再次收缩,如果在R波后的200毫秒内出现一个满足阈值的R波那应该将该R波剔除。
情况二:由于心律失常,QRS波群的幅度可能会很小,这时很有可能按照预设的R波阈值判定找不出符合条件的R波,因此,一旦出现前一个R波和后一个R波相隔大于2.5倍的平均RR间期就说明少检测一个R波,这时将R波阈值每次减小四分之一,以此找出漏检的R波,如果R波阈值降为原R波阈值的二分之一时还没有发现漏检的R波,则说明确实没有R波,这段RR过长可能是由于室颤造成,其中,RR间期是指相邻的两个R波之间的时间间期,平均RR间期是指心电图信号中的多个RR间期取平均值。
从本发明上述技术方案可知,本发明提供了一种复合滤波方案,首先通过IIR数字滤波器对心电信号进行滤波,根据经IIR滤波后得到的第一次滤波后的心电信号在取样周期内的最大值和最小值确定斜率,然后根据该斜率与滤波的斜率阈值的关系判定是否需要通过低通滤波器对第一次滤波后的心电信号进行进一步滤波,一方面,通过IIR滤波器和低通滤波器的复合滤波,能够有效消除心电信号中的噪声,另一方面,上述滤波的斜率阈值可以根据不同的环境进行设定,从而实现在不同的环境下对心电信号的噪声的有效滤除。进一步,本发明实施例还提供了对心电图信号的QRS波群定位的方案,以便根据心电图信号的QRS波群进行与心率有关的诊断。
本发明实施例还提供了根据QRS波群进行与心率有关的智能诊断方案。当通过图2-b所示实施例中的方法定位出QRS波群后,心电信号处理装置便可进一步计算出RR间期(即相邻的两个R波之间的时间间期)、RR间期差(即相邻两个RR间期之前的差)、RR间期的平均值以及QRS时限(即QRS波群的宽度),根据RR间期、RR间期差、RR间期的平均值以及QRS时限限定诊断标准,并结合心率,实现对心动过缓、心动过速、心律不齐、房颤、停搏、漏搏、早搏以及Ron-T现象。
可选地,假设诊断参数定义为:RRi表示RR间期,rr表示RR间期的平均值,RRdif表示RR间期差,W表示QRS时限,HR表示心率,则诊断依据设定如下:
若RR>1.5秒,则判定为心房颤动(房颤);
若RR>3秒,则判定为停搏;
若RR>1.5rr,则判定为漏搏;
若HR>100次/分,则判定为心动过速;
若HR<60次/分,则判定为心动过缓;
若RRdif>0.12秒,则判定为心律不齐;
若RR<0.75rr且W≤120ms,或RRi+RRi+1<2rr,则判定为房性早搏;
若RR<0.75rr且W>120ms且RRi+RRi+1≥2rr,则判定为室性早搏;
当室性早搏的QRS波群落在T波上时,称为室性早搏的Ron-T现象,一直被认为是心律失常领域内的危险信号。由于正常的R波和T波的间距差不多在0.2秒至0.4秒,而当R波叠加到T波上时,两个R波的间距将缩小到R波与T波的间距长度。因此,如果0.2秒<RRi<1/3rr时,认为出现了Ron-T现象。
进一步,在诊断出出现上述任意一种或两种以上状况时,心电信号处理装置根据不同的状况输出相应的报警信号。
本发明实施例还提供一种便携式心电信号处理装置,如图3所示,本发明实施例中的心电信号处理装置300包括:
第一滤波单元301,用于通过IIR数字滤波器对输入的心电信号进行滤波,得到第一次滤波后的心电信号;
第一计算单元302,用于根据上述第一次滤波后的心电信号的波形频率计算取样周期,其中,上述取样周期为上述第一次滤波后的心电信号的波形频率的倒数;
第二计算单元303,用于根据上述第一次滤波后的心电信号在上述取样周期内中的最大值和最小值,计算上述最大值和上述最小值所对应的两点间的直线的斜率;
输出单元304,用于当第二计算单元303计算得到的斜率小于或等于预设的滤波的斜率阈值时,将上述第一次滤波后的心电信号输出为心电图信号;
第二滤波单元305,用于当第二计算单元303计算得到的斜率大于上述滤波的斜率阈值时,通过低通滤波器(例如FIR滤波器或其它衰减较小的滤波器)对上述第一次滤波后的心电信号进行至少一次滤波,直至得到第二心电信号,其中,上述第二心电信号在上述取样周期内的最大值和最小值所对应的两点间的直线的斜率小于或等于上述滤波的斜率阈值;
输出单元304还用于将第二滤波单元305滤波得到的上述第二心电信号输出为心电图信号。
本发明实施例中,上述滤波的斜率阈值可以根据实际环境进行设定,以适应不同环境的变化。例如,对于信号质量较高的原始心电信号(即为经过滤波处理的心电信号),上述滤波的斜率阈值可以设为较大值,对于信号质量较低的原始心电信号,上述滤波的斜率阈值可以设为较小值,即,上述滤波的斜率阈值的取值可以与原始心电信号的信号质量成正相关。
可选地,在图3所示实施例的基础上,心电信号处理装置还包括:定位单元,用于对上述心电图信号进行QRS波群定位。
可选地,上述定位单元包括:R波确定单元,用于确定每个检测周期中的R波,其中,上述检测周期为1.5秒;S点确定单元,用于在每个上述检测周期中,将时间轴上位于R波波峰前75毫秒的位置所对应的信号点确定为S点;Q点确定单元,用于在每个上述检测周期中,将时间轴上位于R波波峰后75毫秒的位置所对应的信号点确定为Q点。
可选地,上述R波确定单元具体包括:
初始化单元,用于根据第一公式初始化R波阈值;
R波阈值更新单元,用于根据第二公式更新R波阈值;
R波确定子单元,用于在上述心电图信号的首个检测周期中:将大于初始化的R波阈值的信号序列确定为上述首个检测周期中的R波,并通过上述R波阈值更新单元更新R波阈值;在上述心电图信号后续的每个检测周期中:将大于更新后的R波阈值的信号序列确定为当前检测周期中的R波,并通过上述R波阈值更新单元更新R波阈值。
其中,上述第一公式为:
thinitial=(max[0]+max[1]+...+max[N-1])/(2*N),式中thinitial表示初始化的R波阈值,max[N-1]表示在上述心电图信号的第N个检测周期内的最大值,其中,上述N为大于或等于2的自然数,可选地,N取4。
上述第二公式为:
thnew=0.8*tholder+0.2*y(max)last/2,式中thnew表示更新后的R波阈值,thlast表示在当前检测周期之前用于定位R波所使用的R波阈值,y(max)last表示当前检测周期中的R波的波峰值。
进一步,心电信号处理装置还可以根据以下原则剔除误判的R波和调整R波阈值:
情况一:心室收缩后有200毫秒的不应期,因此这段时间内心脏不可能再次收缩,如果在R波后的200毫秒内出现一个满足阈值的R波那应该将该R波剔除。
情况二:由于心律失常,QRS波群的幅度可能会很小,这时很有可能按照预设的R波阈值判定找不出符合条件的R波,因此,一旦出现前一个R波和后一个R波相隔大于2.5倍的平均RR间期就说明少检测一个R波,这时将R波阈值每次减小四分之一,以此找出漏检的R波,如果R波阈值降为原R波阈值的二分之一时还没有发现漏检的R波,则说明确实没有R波,这段RR过长可能是由于室颤造成,其中,RR间期是指相邻的两个R波之间的时间间期,平均RR间期是指心电图信号中的多个RR间期取平均值。
需要说明的是,本发明实施例中的心电信号处理装置可以如上述方法实施例中的心电信号处理装置,可以用于实现上述方法实施例中的全部技术方案,其各个功能模块的功能可以根据上述方法实施例中的方法具体实现,其具体实现过程可参照上述实施例中的相关描述,此处不再赘述。
从本发明上述技术方案可知,本发明提供了一种复合滤波方案,首先通过IIR数字滤波器对心电信号进行滤波,根据经IIR滤波后得到的第一次滤波后的心电信号在取样周期内的最大值和最小值确定斜率,然后根据该斜率与滤波的斜率阈值的关系判定是否需要通过低通滤波器对第一次滤波后的心电信号进行进一步滤波,一方面,通过IIR滤波器和低通滤波器的复合滤波,能够有效消除心电信号中的噪声,另一方面,上述滤波的斜率阈值可以根据不同的环境进行设定,从而实现在不同的环境下对心电信号的噪声的有效滤除。进一步,本发明实施例还提供了对心电图信号的QRS波群定位的方案,以便根据心电图信号的QRS波群进行与心率有关的诊断。
在本申请所提供的几个实施例中,应该理解到,所揭露的装置和方法,可以通过其它的方式实现。例如,以上所描述的装置实施例仅仅是示意性的,例如,所述单元的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个单元或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。另一点,所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口,装置或单元的间接耦合或通信连接,可以是电性,机械或其它的形式。
需要说明的是,对于前述的各方法实施例,为了简便描述,故将其都表述为一系列的动作组合,但是本领域技术人员应该知悉,本发明并不受所描述的动作顺序的限制,因为依据本发明,某些步骤可以采用其它顺序或者同时进行。其次,本领域技术人员也应该知悉,说明书中所描述的实施例均属于优选实施例,所涉及的动作和模块并不一定都是本发明所必须的。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详述的部分,可以参见其它实施例的相关描述。
以上为对本发明所提供的一种便携式心电信号处理方法和心电信号处理装置的描述,对于本领域的一般技术人员,依据本发明实施例的思想,在具体实施方式及应用范围上均会有改变之处,综上,本说明书内容不应理解为对本发明的限制。
Claims (10)
1.一种便携式心电信号处理方法,其特征在于,包括:
通过无限脉冲响应IIR数字滤波器对输入的心电信号进行滤波,得到第一次滤波后的心电信号;
根据所述第一次滤波后的心电信号的波形频率计算取样周期,其中,所述取样周期为所述第一次滤波后的心电信号的波形频率的倒数;
根据所述第一次滤波后的心电信号在所述取样周期内中的最大值和最小值,计算所述最大值和所述最小值所对应的两点间的直线的斜率;
若所述斜率小于或等于预设的滤波的斜率阈值,则将所述第一次滤波后的心电信号输出为心电图信号;
若所述斜率大于所述滤波的斜率阈值,则通过低通滤波器对所述第一次滤波后的心电信号进行至少一次滤波,直至得到第二心电信号,其中,所述第二心电信号在所述取样周期内的最大值和最小值所对应的两点间的直线的斜率小于或等于所述滤波的斜率阈值;
将所述第二心电信号输出为心电图信号。
2.根据权利要求1所述的方法,其特征在于,所述方法还包括:
对所述心电图信号进行QRS波群定位。
3.根据权利要求2所述的方法,其特征在于,
所述对所述心电图信号进行QRS波群定位包括:
确定每个检测周期中的R波,其中,所述检测周期为1.5秒;
在每个所述检测周期中,将时间轴上位于R波波峰前75毫秒的位置所对应的信号点确定为S点,将时间轴上位于R波波峰后75毫秒的位置所对应的信号点确定为Q点。
4.根据权利要求3所述的方法,其特征在于,
所述确定每个检测周期中的R波,包括:
根据第一公式初始化R波阈值;
在所述心电图信号的首个检测周期中:将大于初始化的R波阈值的信号序列确定为所述首个检测周期中的R波,并根据第二公式更新R波阈值;
在所述心电图信号后续的每个检测周期中:将大于更新后的R波阈值的信号序列确定为当前检测周期中的R波,并根据所述第二公式更新R波阈值;
其中,所述第一公式为:
thinitial=(max[0]+max[1]+...+max[N-1])/(2*N),式中thinitial表示初始化的R波阈值,max[N-1]表示在所述心电图信号的第N个检测周期内的最大值,其中,所述N为大于或等于2的自然数;
所述第二公式为:
thnew=0.8*tholder+0.2*y(max)last/2,式中thnew表示更新后的R波阈值,tholder表示在当前检测周期之前用于确定R波所使用的R波阈值,y(max)last表示当前检测周期中的R波的波峰值。
5.根据权利要求4所述的方法,其特征在于,所述N等于4。
6.根据权利要求1至5任一项所述的方法,其特征在于,
所述低通滤波器具体为有限长单位冲激响应FIR滤波器。
7.一种便携式心电信号处理装置,其特征在于,包括:
第一滤波单元,用于通过无限脉冲响应IIR数字滤波器对输入的心电信号进行滤波,得到第一次滤波后的心电信号;
第一计算单元,用于根据所述第一次滤波后的心电信号的波形频率计算取样周期,其中,所述取样周期为所述第一次滤波后的心电信号的波形频率的倒数;
第二计算单元,用于根据所述第一次滤波后的心电信号在所述取样周期内中的最大值和最小值,计算所述最大值和所述最小值所对应的两点间的直线的斜率;
输出单元,用于当所述第二计算单元计算得到的斜率小于或等于预设的滤波的斜率阈值时,将所述第一次滤波后的心电信号输出为心电图信号;
第二滤波单元,用于当所述第二计算单元计算得到的斜率大于所述滤波的斜率阈值时,通过FIR滤波器对所述第一次滤波后的心电信号进行至少一次滤波,直至得到第二心电信号,其中,所述第二心电信号在所述取样周期内的最大值和最小值所对应的两点间的直线的斜率小于或等于所述滤波的斜率阈值;
所述输出单元还用于将所述第二滤波单元滤波得到的所述第二心电信号输出为心电图信号。
8.根据权利要求7所述的心电信号处理装置,其特征在于,所述心电信号处理装置还包括:
定位单元,用于对所述心电图信号进行QRS波群定位。
9.根据权利要求8所述的心电信号处理装置,其特征在于,
所述定位单元包括:
R波确定单元,用于确定每个检测周期中的R波,其中,所述检测周期为1.5秒;
S点确定单元,用于在每个所述检测周期中,将时间轴上位于R波波峰前75毫秒的位置所对应的信号点确定为S点;
Q点确定单元,用于在每个所述检测周期中,将时间轴上位于R波波峰后75毫秒的位置所对应的信号点确定为Q点。
10.根据权利要求9所述的心电信号处理装置,其特征在于,
所述R波确定单元具体包括:
初始化单元,用于根据第一公式初始化R波阈值;
R波阈值更新单元,用于根据第二公式更新R波阈值;
R波确定子单元,用于在所述心电图信号的首个检测周期中:将大于初始化的R波阈值的信号序列确定为所述首个检测周期中的R波,并通过所述R波阈值更新单元更新R波阈值;在所述心电图信号后续的每个检测周期中:将大于更新后的R波阈值的信号序列确定为当前检测周期中的R波,并通过所述R波阈值更新单元更新R波阈值;
其中,所述第一公式为:
thinitial=(max[0]+max[1]+...+max[N-1])/(2*N),式中thinitial表示初始化的R波阈值,max[N-1]表示在所述心电图信号的第N个检测周期内的最大值,其中,所述N为大于或等于2的自然数;
所述第二公式为:
thnew=0.8*tholder+0.2*y(max)last/2,式中thnew表示更新后的R波阈值,tholder表示在当前检测周期之前用于确定R波所使用的R波阈值,y(max)last表示当前检测周期中的R波的波峰值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410253367.4A CN105232027B (zh) | 2014-06-09 | 2014-06-09 | 一种便携式心电信号处理方法和心电信号处理装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410253367.4A CN105232027B (zh) | 2014-06-09 | 2014-06-09 | 一种便携式心电信号处理方法和心电信号处理装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105232027A true CN105232027A (zh) | 2016-01-13 |
CN105232027B CN105232027B (zh) | 2017-12-08 |
Family
ID=55030061
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410253367.4A Active CN105232027B (zh) | 2014-06-09 | 2014-06-09 | 一种便携式心电信号处理方法和心电信号处理装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105232027B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107080522A (zh) * | 2017-03-16 | 2017-08-22 | 深圳竹信科技有限公司 | 信号处理方法及装置 |
CN107647864A (zh) * | 2017-11-06 | 2018-02-02 | 上海联影医疗科技有限公司 | 一种心电信号分析方法及成像方法 |
CN109106397A (zh) * | 2017-06-25 | 2019-01-01 | 吴健康 | 一种胎心音监测和分析系统 |
CN109688916A (zh) * | 2016-09-07 | 2019-04-26 | 心脏起搏器股份公司 | 可植入心脏监视器的过缓暂停检测 |
CN109843165A (zh) * | 2016-10-17 | 2019-06-04 | 日本电信电话株式会社 | 心跳检测方法和心跳检测设备 |
CN112401906A (zh) * | 2020-11-10 | 2021-02-26 | 河北省科学院应用数学研究所 | 一种基于幅值的自动脑电去伪迹方法 |
US11138770B2 (en) | 2017-11-06 | 2021-10-05 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for medical imaging |
CN114533080A (zh) * | 2020-11-26 | 2022-05-27 | 深圳市科瑞康实业有限公司 | 一种对心电信号数据进行颤动波识别的方法和装置 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060155199A1 (en) * | 2002-05-31 | 2006-07-13 | Regis Logier | Frequency processing of an rr series in an analogue cardiac signal |
CN101049235A (zh) * | 2007-04-26 | 2007-10-10 | 中国科学技术大学 | 心室复极高频波体表检测方法与装置 |
CN101361650A (zh) * | 2007-08-07 | 2009-02-11 | 深圳市理邦精密仪器有限公司 | Iir滤波器的零相位实现方法及零相位iir滤波装置 |
CN101449971A (zh) * | 2008-12-30 | 2009-06-10 | 南京大学 | 基于节律模式的便携式心电诊断监测设备 |
CN201759552U (zh) * | 2010-09-04 | 2011-03-16 | 中南林业科技大学 | 便携式家用心电检测装置 |
US20110196246A1 (en) * | 2010-02-09 | 2011-08-11 | Angel Medical Systems Inc. | Qrs onset and offset detection with adaptive temporal windowing |
US20130237874A1 (en) * | 2012-03-12 | 2013-09-12 | Texas Instruments Incorporated | Real Time QRS Detection Using Adaptive Threshold |
WO2014058664A1 (en) * | 2012-10-11 | 2014-04-17 | Medtronic, Inc. | Determining onsets and offsets of cardiac depolarization and repolarization waves |
-
2014
- 2014-06-09 CN CN201410253367.4A patent/CN105232027B/zh active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20060155199A1 (en) * | 2002-05-31 | 2006-07-13 | Regis Logier | Frequency processing of an rr series in an analogue cardiac signal |
CN101049235A (zh) * | 2007-04-26 | 2007-10-10 | 中国科学技术大学 | 心室复极高频波体表检测方法与装置 |
CN101361650A (zh) * | 2007-08-07 | 2009-02-11 | 深圳市理邦精密仪器有限公司 | Iir滤波器的零相位实现方法及零相位iir滤波装置 |
CN101449971A (zh) * | 2008-12-30 | 2009-06-10 | 南京大学 | 基于节律模式的便携式心电诊断监测设备 |
US20110196246A1 (en) * | 2010-02-09 | 2011-08-11 | Angel Medical Systems Inc. | Qrs onset and offset detection with adaptive temporal windowing |
CN201759552U (zh) * | 2010-09-04 | 2011-03-16 | 中南林业科技大学 | 便携式家用心电检测装置 |
US20130237874A1 (en) * | 2012-03-12 | 2013-09-12 | Texas Instruments Incorporated | Real Time QRS Detection Using Adaptive Threshold |
WO2014058664A1 (en) * | 2012-10-11 | 2014-04-17 | Medtronic, Inc. | Determining onsets and offsets of cardiac depolarization and repolarization waves |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109688916A (zh) * | 2016-09-07 | 2019-04-26 | 心脏起搏器股份公司 | 可植入心脏监视器的过缓暂停检测 |
CN109688916B (zh) * | 2016-09-07 | 2022-01-11 | 心脏起搏器股份公司 | 可植入心脏监视器的过缓暂停检测 |
CN109843165B (zh) * | 2016-10-17 | 2021-10-01 | 日本电信电话株式会社 | 心跳检测方法和心跳检测设备 |
CN109843165A (zh) * | 2016-10-17 | 2019-06-04 | 日本电信电话株式会社 | 心跳检测方法和心跳检测设备 |
CN107080522A (zh) * | 2017-03-16 | 2017-08-22 | 深圳竹信科技有限公司 | 信号处理方法及装置 |
CN109106397A (zh) * | 2017-06-25 | 2019-01-01 | 吴健康 | 一种胎心音监测和分析系统 |
US11138770B2 (en) | 2017-11-06 | 2021-10-05 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for medical imaging |
CN107647864A (zh) * | 2017-11-06 | 2018-02-02 | 上海联影医疗科技有限公司 | 一种心电信号分析方法及成像方法 |
CN107647864B (zh) * | 2017-11-06 | 2024-04-26 | 上海联影医疗科技股份有限公司 | 一种心电信号分析方法及成像方法 |
US11989804B2 (en) | 2017-11-06 | 2024-05-21 | Shanghai United Imaging Healthcare Co., Ltd. | Systems and methods for medical imaging |
CN112401906A (zh) * | 2020-11-10 | 2021-02-26 | 河北省科学院应用数学研究所 | 一种基于幅值的自动脑电去伪迹方法 |
CN114533080A (zh) * | 2020-11-26 | 2022-05-27 | 深圳市科瑞康实业有限公司 | 一种对心电信号数据进行颤动波识别的方法和装置 |
CN114533080B (zh) * | 2020-11-26 | 2023-10-13 | 深圳市科瑞康实业有限公司 | 一种对心电信号数据进行颤动波识别的方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105232027B (zh) | 2017-12-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105232027A (zh) | 一种便携式心电信号处理方法和心电信号处理装置 | |
US8442624B2 (en) | System for cardiac medical condition detection | |
Thakor et al. | Applications of adaptive filtering to ECG analysis: noise cancellation and arrhythmia detection | |
JPS63226338A (ja) | Ecg波信号からペース・パルス信号を除去する方法とecg波信号からペース・パルス信号を除去する装置 | |
WO2002009582A2 (en) | Method and apparatus for a morphology-preserving smoothing | |
CN105899258B (zh) | 可变带宽ecg高通滤波器 | |
CN108113665A (zh) | 一种心电信号自动降噪方法 | |
US20120197088A1 (en) | Detection of waveform artifact | |
CN110226919B (zh) | 心电信号类型检测方法、装置、计算机设备及存储介质 | |
CN109171711A (zh) | 一种基于极值法的快速p波检测方法 | |
CN106691376A (zh) | 心电信号自适应滤波方法及装置 | |
CN110327032A (zh) | 一种单导心电信号pqrst波联合精准识别算法 | |
CN105078444A (zh) | 一种噪声检测方法、装置及医疗检测设备 | |
Prasad et al. | Heart rate detection using Hilbert transform | |
Rahman et al. | Baseline wandering removal from ECG signal by wandering path finding algorithm | |
CN109107041B (zh) | 起搏信号检测方法、装置及可穿戴健康设备 | |
JP7189210B2 (ja) | 植込み型医療装置用のハイブリッド信号処理回路及び方法 | |
US9020584B2 (en) | Method for the detection of subcutanous cardiac signals and a cardiac device for use in detecting subcutaneous cardiac signals | |
KR101048763B1 (ko) | 신호 검출 장치 및 방법 | |
Amin et al. | Separation of fetal electrocardiography (ECG) from composite ECG using adaptive linear neural network for fetal monitoring | |
KR20220042908A (ko) | 적응적 메디안 필터를 이용한 심전도 신호의 r 피크 검출 장치 및 방법 | |
CN110432885A (zh) | 一种光电容积脉搏波噪声去除方法 | |
CN111643053A (zh) | 一种减少脉搏波信号中运动伪影的方法及系统 | |
CN109770858B (zh) | 一种心电图的异常检测方法、装置、设备和存储介质 | |
EP3413793B1 (en) | Detecting stimulus pulses |
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 | ||
GR01 | Patent grant |