CN107616782A - 一种心电信号质量检测方法与装置 - Google Patents
一种心电信号质量检测方法与装置 Download PDFInfo
- Publication number
- CN107616782A CN107616782A CN201710797429.1A CN201710797429A CN107616782A CN 107616782 A CN107616782 A CN 107616782A CN 201710797429 A CN201710797429 A CN 201710797429A CN 107616782 A CN107616782 A CN 107616782A
- Authority
- CN
- China
- Prior art keywords
- electrocardiosignal
- quality
- zero crossing
- window
- energy
- 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.)
- Pending
Links
Landscapes
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Abstract
本发明提供一种心电信号质量检测方法与装置,所述方法包括:对心电信号进行能量包络提取;对能量包络进行差分运算;分段统计窗口内过零点;根据阈值判定信号质量。本发明通过使用能量包络的差分方法,将过零点的数量与心电信号质量关联起来,即只需通过统计差分数据过零点的数量即可判定心电信号质量,从而实现准确判定心电信号的质量,相比于传统心电信号质量分析算法存在的复杂度高且必须依靠强大的硬件资源才能实现,本算法简单,易于实现且运算复杂度低。
Description
技术领域
本发明涉及医学信号处理技术领域,特别是一种心电信号质量检测方法与装置。
背景技术
当今医学上,心血管疾病属于发病率和死亡率最高的疾病之一,有效预防和诊断心脏病成为急需解决的重要问题之一。心电信号是诊断心血管疾病和评估心脏功能的主要依据,因此关于心电信号检测和处理的研究一直备受关注,心电图作为一种无创的检测技术,因诊断方法简单、无痛无害等特点,一直是诊断心脏病的主要技术之一,在心血管疾病的诊断中发挥了极其重要的作用。
但是心电信号的采集非常容易受到各种因素的影响,如采集设备、人为因素、环境因素以及高频杂音等造成的噪声,导致采集到的信号出现质量问题,严重影响心电信号的分析与诊断,因此在对心电信号分析诊断之前,使用有效的方式检测出心电信号的质量并提取出质量好的信号进行分析和诊断,具有极其迫切的现实需求。
当前心电信号质量检测方法有许多种,主要是通过信号波形的频域分析和时域分析来判断信号波形质量。其中频域分析主要有三种方法:
(1)通过提取原始信号频谱的高频分量来判断波形质量。
(2)对原始信号进行带通滤波处理,提取处理后的信号频谱的高频分量来判断波形质量。
(3)通过提取原始信号能量包络的频谱来判断波形质量。
然而正常波形和非正常波形在频谱上虽然存在一定差异,却没有很明显的界限,导致以上三种方法容易出现对信号的漏检和误检,无法对信号质量作出准确判断。
时域分析主要有以下两种方法:
(1)基于波形质量差的区域,香农能量包络在一定时间内无法回零的原理,通过统计香农能量包络中小于阈值点的个数来判断波形质量。但是这种方法容易出现误判。
(2)通过多阈值分布的方法来判断波形质量。缺点在于不易取得合适的阈值。
发明内容
本发明的目的是提供一种心电信号质量检测方法与装置,旨在解决当前心电信号质量检测算法复杂度较高,需要耗费大量的资源,而无法准确判定心电信号质量而影响分析和诊断的问题,实现准确度高、复杂性低的技术效果。
为达到上述技术目的,本发明提供了一种心电信号质量检测方法,包括以下步骤:
S1、对心电信号进行能量包络提取;
S2、对能量包络进行差分运算;
S3、分段统计窗口内过零点;
S4、根据阈值判定信号质量。
优选地,所述步骤S1具体包括以下步骤:
S101、对心电信号进行归一化处理;
S102、计算归一化心电信号的香农能量;
S103、对所述香农能量进行平滑处理;
S104、提取香农能量包络。
优选地,所述过零点根据以下公式确定:
y(t-1)<0且y(t)>0,或
y(t-1)>0且y(t)<0
y(t)为香农能量差分运算值。
优选地,所述步骤S4具体为:设定所述阈值为所有窗口内过零点个数的平均值,当窗口内过零点数量大于所述阈值,则判定信号质量不合格;当窗口内过零点数量小于所述阈值,则判定信号质量合格。
本发明还提供了一种心电信号质量检测装置,包括:
能量包络提取模块,用于对心电信号进行能量包络提取;
差分运算模块,用于对能量包络进行差分运算;
过零点统计模块,用于分段统计窗口内过零点;
信号质量判定模块,用于根据阈值判定信号质量。
优选地,所述能量包络提取模块包括:
归一化处理单元,用于对心电信号进行归一化处理;
香农能量计算单元,用于计算归一化心电信号的香农能量;
平滑处理单元,用于对所述香农能量进行平滑处理;
提取单元,用于提取香农能量包络。
优选地,所述过零点根据以下公式确定:
y(t-1)<0且y(t)>0,或
y(t-1)>0且y(t)<0
y(t)为香农能量差分运算值。
优选地,所述信号质量判定模块包括:
阈值设定单元,用于设定窗口过零点阈值;
判定单元,用于将窗口过零点数量与所述阈值进行对比,当窗口内过零点数量大于所述阈值,则判定信号质量不合格;当窗口内过零点数量小于所述阈值,则判定信号质量合格。
发明内容中提供的效果仅仅是实施例的效果,而不是发明所有的全部效果,上述技术方案中的一个技术方案具有如下优点或有益效果:
与现有技术相比,本发明通过使用能量包络的差分方法,将过零点的数量与心电信号质量关联起来,即只需通过统计差分数据过零点的数量即可判定心电信号质量,从而实现准确判定心电信号的质量,相比于传统心电信号质量分析算法存在的复杂度高且必须依靠强大的硬件资源才能实现,本算法简单,易于实现且运算复杂度低。
附图说明
图1为本发明实施例中所提供的一种心电信号质量检测方法流程图;
图2为本发明实施例中所提供的一种心电信号质量检测方法过程中能量包络提取方法流程图;
图3为本发明实施例中所提供的一种心电信号质量检测装置整体框图。
具体实施方式
为了能清楚说明本方案的技术特点,下面通过具体实施方式,并结合其附图,对本发明进行详细阐述。下文的公开提供了许多不同的实施例或例子用来实现本发明的不同结构。为了简化本发明的公开,下文中对特定例子的部件和设置进行描述。此外,本发明可以在不同例子中重复参考数字和/或字母。这种重复是为了简化和清楚的目的,其本身不指示所讨论各种实施例和/或设置之间的关系。应当注意,在附图中所图示的部件不一定按比例绘制。本发明省略了对公知组件和处理技术及工艺的描述以避免不必要地限制本发明。
下面结合附图对本发明实施例所提供的一种心电信号质量检测方法与装置进行详细说明。
如图1所示,本发明实施例公开了一种心电信号质量检测方法,包括以下步骤:
S1、对心电信号进行包络提取,具体步骤如图2所示:
S101、归一化处理;
从MIT-BIH(美国麻省理工学院提供的研究心律失常的数据库)中提取心电信号,对原始心电信号进行幅值归一化处理,计算原始心电信号绝对值的最大值,将原始心电信号除以得到的最大值,即完成信号幅值的归一化处理。公式如下:
x(t)为原始心电信号,xnorm(t)为归一化后的心电信号。归一化后的心电信号方便后续信号处理。
S102、计算归一化心电信号的香农能量;
计算归一化后心电信号的香农能量,计算公式如下:
E(t)为归一化后心电信号的香农能量,xnorm(t)为归一化后的心电信号。
S103、对心电信号的香农能量进行平滑处理;
设置平均区间N,对归一化心电信号的香农能量进行平滑处理,计算公式如下:
Es(t)为平滑后的香农能量。
S104、提取香农能量包络。
对平滑处理后的香农能量求平均值和方差,从而获得香农能量包络,计算公式如下:
P(t)为归一化香农能量包络,M(ES(t))为香农能量平均值,S(ES(t))为香农能量方差。
S2、对香农能量包络进行差分运算;
通过差分运算可以将波形的波峰和波谷转化为过零点,便于对波形的波峰及波谷也即极大值和极小值点进行统计,通过统计过零点的数量即可求得对应峰谷的数量,计算公式如下:
y(t)=P(t)-P(t-1)
采取前向差分的方式对香农能量包络进行差分运算,y(t)为香农能量差分运算值。
S3、分段统计过零点;
根据差分运算后的香农能量包络设置滑动窗口,每个滑动窗口的大小设置为两个分段点的大小,根据该滑动窗口,将差分运算后的数据按时间推移,将每一次窗口滑动后窗口内的过零点个数进行统计,窗口每次滑动一个数据点的位置。当存在以下情况时,则确定为过零点:
y(t-1)<0且y(t)>0,或
y(t-1)>0且y(t)<0
统计所有确定为过零点的数量。
S4、设定阈值,判定信号质量。
由于窗口内过零点的个数与窗口内心电信号的质量相关,窗口内的过零点的数量越多,说明心电信号噪声越大,反之,则心电信号噪声越小。因此,对窗口内过零点的数量设置一个阈值,阈值的大小可以根据经验来确定,当窗口内过零点的数量大于这个阈值,则判定该窗口内的心电信号数据质量不符合要求,并做相应的标记,当窗口内过零点的数量小于这个阈值,则判定该窗口内的心电信号数据符合要求。当需要对心电信号作进一步分析时,则直接跳过标记为不符合要求的心电信号数据段落,从而避免得到错误的分析或诊断结果。
对原始心电信号数据进行香农能量包络、差分运算、过零点分布以及质量不合格标记的效果仿真如图2所示。
本发明实施例通过使用能量包络的差分方法,将过零点的数量与心电信号质量关联起来,即只需通过统计差分数据过零点的数量即可判定心电信号质量,从而实现准确判定心电信号的质量,相比于传统心电信号质量分析算法存在的复杂度高且必须依靠强大的硬件资源才能实现,本算法简单,易于实现且运算复杂度低。
如图3所示,本发明实施例还公开了一种心电信号质量检测装置,包括能量包络提取模块、差分运算模块、过零点统计模块、信号质量判定模块。
所述能量包络提取模块,用于对心电信号进行能量包络提取,包括:
归一化处理单元,用于对心电信号进行归一化处理;
对原始心电信号进行幅值归一化处理,计算原始心电信号绝对值的最大值,将原始心电信号除以得到的最大值,即完成信号幅值的归一化处理。公式如下:
x(t)为原始心电信号,xnorm(t)为归一化后的心电信号。归一化后的心电信号方便后续信号处理。
香农能量计算单元,用于计算归一化心电信号的香农能量;
计算归一化后心电信号的香农能量,计算公式如下:
E(t)为归一化后心电信号的香农能量,xnorm(t)为归一化后的心电信号。
平滑处理单元,用于对所述香农能量进行平滑处理;
设置平均区间N,对归一化心电信号的香农能量进行平滑处理,计算公式如下:
Es(t)为平滑后的香农能量。
提取单元,用于提取香农能量包络。
对平滑处理后的香农能量求平均值和方差,从而获得香农能量包络,计算公式如下:
P(t)为归一化香农能量包络,M(Es(t))为香农能量平均值,S(Es(t))为香农能量方差。
所述差分运算模块,用于对能量包络进行差分运算;
通过差分运算可以将波形的波峰和波谷转化为过零点,便于对波形的波峰及波谷也即极大值和极小值点进行统计,通过统计过零点的数量即可求得对应峰谷的数量,计算公式如下:
y(t)=P(t)-P(t-1)
采取前向差分的方式对香农能量包络进行差分运算,y(t)为香农能量差分运算值。
所述过零点统计模块,用于分段统计过零点;
根据差分运算后的香农能量包络设置滑动窗口,每个滑动窗口的大小设置为两个分段点的大小,根据该滑动窗口,将差分运算后的数据按时间推移,将每一次窗口滑动后窗口内的过零点个数进行统计,窗口每次滑动一个数据点的位置。当存在以下情况时,则确定为过零点:
y(t-1)<0且y(t)>0,或
y(t-1)>0且y(t)<0
统计所有确定为过零点的数量。
所述信号质量判定模块,用于根据阈值判定信号质量,包括:
阈值设定单元,用于设定窗口过零点阈值,将阈值设定为所有窗口内过零点个数的平均值。
判定单元,用于将窗口过零点数量与所述阈值进行对比,当窗口内过零点数量大于所述阈值,则判定信号质量不合格;当窗口内过零点数量小于所述阈值,则判定信号质量合格。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种心电信号质量检测方法,其特征在于,包括以下步骤:
S1、对心电信号进行能量包络提取;
S2、对能量包络进行差分运算;
S3、分段统计窗口内过零点;
S4、根据阈值判定信号质量。
2.根据权利要求1所述的一种心电信号质量检测方法,其特征在于,所述步骤S1具体包括以下步骤:
S101、对心电信号进行归一化处理;
S102、计算归一化心电信号的香农能量;
S103、对所述香农能量进行平滑处理;
S104、提取香农能量包络。
3.根据权利要求1所述的一种心电信号质量检测方法,其特征在于,所述过零点根据以下公式确定:
y(t-1)<0且y(t)>0,或
y(t-1)>0且y(t)<0
y(t)为香农能量差分运算值。
4.根据权利要求1所述的一种心电信号质量检测方法,其特征在于,所述步骤S4具体为:设定所述阈值为所有窗口内过零点个数的平均值,当窗口内过零点数量大于所述阈值,则判定信号质量不合格;当窗口内过零点数量小于所述阈值,则判定信号质量合格。
5.一种心电信号质量检测装置,其特征在于,包括:
能量包络提取模块,用于对心电信号进行能量包络提取;
差分运算模块,用于对能量包络进行差分运算;
过零点统计模块,用于分段统计窗口内过零点;
信号质量判定模块,用于根据阈值判定信号质量。
6.根据权利要求5所述的一种心电信号质量检测装置,其特征在于,所述能量包络提取模块包括:
归一化处理单元,用于对心电信号进行归一化处理;
香农能量计算单元,用于计算归一化心电信号的香农能量;
平滑处理单元,用于对所述香农能量进行平滑处理;
提取单元,用于提取香农能量包络。
7.根据权利要求5所述的一种心电信号质量检测装置,其特征在于,所述过零点根据以下公式确定:
y(t-1)<0且y(t)>0,或
y(t-1)>0且y(t)<0
y(t)为香农能量差分运算值。
8.根据权利要求5所述的一种心电信号质量检测装置,其特征在于,所述信号质量判定模块包括:
阈值设定单元,用于设定窗口过零点阈值;
判定单元,用于将窗口过零点数量与所述阈值进行对比,当窗口内过零点数量大于所述阈值,则判定信号质量不合格;当窗口内过零点数量小于所述阈值,则判定信号质量合格。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710797429.1A CN107616782A (zh) | 2017-09-06 | 2017-09-06 | 一种心电信号质量检测方法与装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710797429.1A CN107616782A (zh) | 2017-09-06 | 2017-09-06 | 一种心电信号质量检测方法与装置 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107616782A true CN107616782A (zh) | 2018-01-23 |
Family
ID=61089662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710797429.1A Pending CN107616782A (zh) | 2017-09-06 | 2017-09-06 | 一种心电信号质量检测方法与装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107616782A (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109893121A (zh) * | 2019-03-26 | 2019-06-18 | 深圳理邦智慧健康发展有限公司 | 心电信号的采集方法、装置、终端和计算机可读存储介质 |
CN110179451A (zh) * | 2019-06-10 | 2019-08-30 | 深圳市是源医学科技有限公司 | 心电信号质量检测方法、装置、计算机设备和存储介质 |
CN112741623A (zh) * | 2020-12-09 | 2021-05-04 | 中国科学院深圳先进技术研究院 | 一种基于心电的糖化血红蛋白检测系统和方法 |
CN113349754A (zh) * | 2021-07-20 | 2021-09-07 | 成都芯跳医疗科技有限责任公司 | 一种基于便携式动态心电监护仪的心率计算方法 |
CN113499059A (zh) * | 2021-06-01 | 2021-10-15 | 武汉理工大学 | 基于光纤传感非接触式的bcg信号处理系统及方法 |
CN115343564A (zh) * | 2022-10-18 | 2022-11-15 | 青岛鼎信通讯股份有限公司 | 一种应用于电力领域的信号检测方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102028459A (zh) * | 2010-12-02 | 2011-04-27 | 广东宝莱特医用科技股份有限公司 | 一种心电图机通道起搏信号检测方法 |
CN104305991A (zh) * | 2014-11-18 | 2015-01-28 | 北京海思敏医疗技术有限公司 | 从心电信号中检测噪声的方法和设备 |
-
2017
- 2017-09-06 CN CN201710797429.1A patent/CN107616782A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102028459A (zh) * | 2010-12-02 | 2011-04-27 | 广东宝莱特医用科技股份有限公司 | 一种心电图机通道起搏信号检测方法 |
CN104305991A (zh) * | 2014-11-18 | 2015-01-28 | 北京海思敏医疗技术有限公司 | 从心电信号中检测噪声的方法和设备 |
Non-Patent Citations (1)
Title |
---|
M.SABARIMALAI MANIKANDAN,ET AL.: "A novel method for detecting R-peaks in electrocardiogram(ECG) signal", 《BIOMEDICAL SIGNAL PROCESSING AND CONTROL》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109893121A (zh) * | 2019-03-26 | 2019-06-18 | 深圳理邦智慧健康发展有限公司 | 心电信号的采集方法、装置、终端和计算机可读存储介质 |
CN110179451A (zh) * | 2019-06-10 | 2019-08-30 | 深圳市是源医学科技有限公司 | 心电信号质量检测方法、装置、计算机设备和存储介质 |
CN110179451B (zh) * | 2019-06-10 | 2021-10-29 | 深圳市是源医学科技有限公司 | 心电信号质量检测方法、装置、计算机设备和存储介质 |
CN112741623A (zh) * | 2020-12-09 | 2021-05-04 | 中国科学院深圳先进技术研究院 | 一种基于心电的糖化血红蛋白检测系统和方法 |
CN113499059A (zh) * | 2021-06-01 | 2021-10-15 | 武汉理工大学 | 基于光纤传感非接触式的bcg信号处理系统及方法 |
CN113349754A (zh) * | 2021-07-20 | 2021-09-07 | 成都芯跳医疗科技有限责任公司 | 一种基于便携式动态心电监护仪的心率计算方法 |
CN115343564A (zh) * | 2022-10-18 | 2022-11-15 | 青岛鼎信通讯股份有限公司 | 一种应用于电力领域的信号检测方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107616782A (zh) | 一种心电信号质量检测方法与装置 | |
Lozano et al. | Automatic differentiation of normal and continuous adventitious respiratory sounds using ensemble empirical mode decomposition and instantaneous frequency | |
Karimipour et al. | Real-time electrocardiogram P-QRS-T detection–delineation algorithm based on quality-supported analysis of characteristic templates | |
CN104173043B (zh) | 一种适合于移动平台的心电数据分析方法 | |
CN109077715A (zh) | 一种基于单导联的心电信号自动分类方法 | |
US20150105666A1 (en) | Narrow band feature extraction from cardiac signals | |
CN108992053B (zh) | 一种实时的无束缚检测心率和心跳间隔的方法 | |
CN104887220A (zh) | 一种由腹壁心电信号提取胎儿心电信号的方法及系统 | |
US10602944B2 (en) | Detecting artifacts in a signal | |
CN105550659A (zh) | 基于随机投影的实时心电分类方法 | |
Satija et al. | A simple method for detection and classification of ECG noises for wearable ECG monitoring devices | |
CN107811631A (zh) | 心电信号质量评估方法 | |
Chen et al. | Identification of Velcro rales based on Hilbert–Huang transform | |
Mahesh et al. | ECG arrhythmia classification based on logistic model tree | |
CN106361325A (zh) | 一种便携式心电仪所测单导联心电图的筛查识别系统 | |
Das et al. | Arduino-based noise robust online heart-rate detection | |
Li | A quality assessment method of single-lead ECG signal based on spectral analysis | |
Satija et al. | Low-complexity detection and classification of ECG noises for automated ECG analysis system | |
US7197358B2 (en) | Identifying infants at risk for sudden infant death syndrome | |
So-In et al. | Real-time ECG noise reduction with QRS complex detection for mobile health services | |
Pan et al. | Detection of ECG characteristic points using biorthogonal spline wavelet | |
US20080281216A1 (en) | Cardiac Arrhythmias Analysis of Electrophysiological Signals Based on Symbolic Dynamics | |
US10285613B2 (en) | Apparatus and method for determining the occurrence of a QRS complex in ECG data | |
US11944464B2 (en) | Methods and system for detecting inhalations and extracting measures of neural respiratory drive from an EMG signal | |
Chatterjee et al. | Real–time detection of electrocardiogram wave features using template matching and implementation in FPGA |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180123 |
|
RJ01 | Rejection of invention patent application after publication |