CN113907765B - 一种无创胎儿心电信号质量评估方法 - Google Patents

一种无创胎儿心电信号质量评估方法 Download PDF

Info

Publication number
CN113907765B
CN113907765B CN202111178668.1A CN202111178668A CN113907765B CN 113907765 B CN113907765 B CN 113907765B CN 202111178668 A CN202111178668 A CN 202111178668A CN 113907765 B CN113907765 B CN 113907765B
Authority
CN
China
Prior art keywords
electrocardiosignal
abdomen
signal
fetal
electrocardiosignals
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.)
Active
Application number
CN202111178668.1A
Other languages
English (en)
Other versions
CN113907765A (zh
Inventor
张松
谢佳兴
杨益民
郝冬梅
杨琳
李旭雯
张建楠
张雅俊
刘建利
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing University of Technology
Original Assignee
Beijing University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Beijing University of Technology filed Critical Beijing University of Technology
Priority to CN202111178668.1A priority Critical patent/CN113907765B/zh
Publication of CN113907765A publication Critical patent/CN113907765A/zh
Application granted granted Critical
Publication of CN113907765B publication Critical patent/CN113907765B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • A61B5/318Heart-related electrical modalities, e.g. electrocardiography [ECG]
    • A61B5/344Foetal cardiography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7203Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7221Determining signal validity, reliability or quality
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7225Details of analog processing, e.g. isolation amplifier, gain or sensitivity adjustment, filtering, baseline or drift compensation
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems

Abstract

本发明公开了一种无创胎儿心电信号质量评估方法,将胎儿心电明显的信号与受到强烈干扰或胎儿心电不清晰的信号区分。该方法对采集到的腹部电信号进行预处理和分离,通过计算预处理后的腹部心电信号和分离得到的胎儿心电信号的信息熵、Lempel‑Ziv复杂度、母胎能量比、R波个数等作为评估参数,运用支持向量机判断腹部心电信号和胎儿心电信号为“质量可接受”或“质量不可接受”。本发明实现了对腹部心电信号和胎儿心电信号的信号质量分类,能够避免因信号质量不高而导致提取错误胎心率的情况,在长时间胎心监护过程中,能够实时判断信号质量,对无创提取准确可靠的胎心率有一定的帮助。

Description

一种无创胎儿心电信号质量评估方法
技术领域
本发明属于医学信号处理技术领域,特别涉及一种无创胎儿心电信号质量评估方法。
背景技术
胎儿及新生儿的死亡率是衡量一个国家经济、文化及医疗等发展水平的重要依据。
胎儿在母亲体内出现胎儿窘迫等其他问题都会直接反映在胎心率的变化上。胎儿健康监护能提前检测到胎儿的某些心脏疾病、被脐带缠绕引起的窒息和宫内窘迫等异常情况,实现胎儿健康安全预防和优生,能大大降低胎儿和新生儿的死亡率。
目前,获取胎心率有两种主要方法:超声多普勒法和胎儿心电法。其中,超声多普勒法存在着许多的不足,首先胎心探头对于位置要求比较高,往往需要专业人员进行操作;在监测过程中孕妇不能移动,行动受限,无法做到长时间监护等等。并且有研究表明,长时间超声检测或许会对胎儿发育产生副作用,表现在对绒毛膜结构、角膜结构、神经发育的影响。而相比于传统的超声多普勒监护方式,胎儿心电法可以进行长时间监测,在监测过程中孕妇可以适当活动。综合来看,胎儿心电法是目前医护研究人员研究的重点。
胎儿心电法通过在母亲腹部体表放置采集电极来被动接收微弱的非入侵式生物电信号,但在采集过程中通常会受到许多因素的干扰,最明显的是母亲心电信号,另外非平稳噪声的存在会影响胎儿心电信号的信噪比,导致心率不易识别或错误识别的情况出现。因此,需要对预处理后的腹部心电信号以及分离得到的胎儿心电信号进行信号质量评估,将信号质量差的信号在后续分析过程中剔除,可以提高心率检测的准确性,及时避免错误获取母亲心率和胎心率。
现有的腹部心电信号和胎儿心电信号质量评估方法相对较少,大多是通过计算信号质量指数来对腹部心电信号或胎儿心电信号进行判断。同时由于采集到的腹部心电信号和分离得到的胎儿心电信号个体差异性较大,只计算某一个信号质量指数难以做到准确判断信号质量。因此,需要一种能够对腹部心电信号和胎儿心电信号进行质量评估的方法,在对腹部心电信号进行分离前,能够剔除掉信号质量差的信号;同时在计算胎心率前,能够评估胎儿心电信号质量,避免获取错误的胎心率。
发明内容
本发明的目的是提供一种无创胎儿心电信号质量评估方法,以解决上述背景技术中提到的问题。
本发明的无创胎儿心电信号质量评估方法,包括以下步骤:
步骤1、采集得到多组4路通道的腹部电信号;
步骤2、对采集到的信号进行预处理以及胎儿心电信号分离,将预处理后的腹部心电信号和分离得到的胎儿心电信号的质量评估结果分为可接受与不可接受两类;
步骤3、对预处理后的腹部心电信号的每一通道信号计算获取信息熵、Lempel-Ziv复杂度、母胎能量比以及相对功率比,作为对应通道信号的特征数据。并将得到的特征数据,通过基于支持向量机方法进行训练学习获得腹部心电信号分类模型;
步骤4、对分离得到的胎儿心电信号计算信息熵、Lempel-Ziv复杂度、R波个数作为特征数据,通过支持向量机方法进行训练学习获得胎儿心电信号分类模型;
步骤5、对新采集的腹部心电测试信号,截取预设的长度后,根据步骤3计算特征数据,放入训练好的腹部心电信号分类模型获得信号质量的评估结果;
步骤6、对分离得到的胎儿心电测试信号,根据步骤4计算特征数据,放入训练好的胎儿心电信号分类模型获得信号质量的评估结果。
本发明具有以下优点:
(1)本发明选取预处理后的腹部心电信号和胎儿心电信号的波形特征和非线性特征,将其作为特征数据对信号质量进行评估,质量评估结果可以很好地反映腹部心电信号和胎儿心电信号的生理特征。同时计算多个特征有利于消除个体差异性,提高分类准确率;
(2)本发明对腹部心电信号和胎儿心电信号两种信号的信号质量都进行了评估,不再只对一种信号进行质量评估。当腹部心电信号或胎儿心电信号质量很差时,能够实时获取模型评估结果并结束后续的模块处理,大大节省了时间和计算复杂度;
(3)本发明提出的质量评估方法同样适用于单通道的腹部心电信号,拥有较高的准确率;
附图说明
图1为本发明的总流程图
图2为腹部心电信号质量评估方法的流程图
图3为胎儿心电信号质量评估方法的流程图
图4为腹部心电信号质量评估实际操作流程图
图5为胎儿心电信号质量评估实际操作流程图
具体实施方式
下面将结合附图对本发明的具体实施方式进行详细描述。
本发明提供一种无创胎儿心电信号质量评估方法,其总流程如图1所示;其中腹部心电信号质量评估方法如图2所示,胎儿心电信号质量评估方法如图3所示,信号质量评估实际操作过程如图4和图5所示;具体实施步骤如下:
步骤1、采集得到多组4路通道的腹部电信号X;
步骤2、对采集得到的信号X进行预处理以及胎儿心电信号分离,并将预处理后的腹部心电信号和分离得到的胎儿心电信号的质量评估结果分为可接受与不可接受两类;
具体地,预处理算法主要包括:
5Hz高通滤波器去除腹部心电信号的基线漂移;
采用小波软阈值去噪的方式,去除腹部心电信号的高频噪声;
50Hz陷波器,去除腹部心电信号的工频干扰;
采用Maurizio Varanini在2013年提出的算法对胎儿心电信号进行分离,(该算法的文章发表在COMPUTING IN CARDIOLOGY CONFERENCE会议上,文章名称为“一种多步骤的胎儿心电分离方法”,是该会议上评分最高的开源算法)。从腹部心电信号分离得到胎儿心电信号;
步骤3、对预处理后的腹部心电信号的每一通道信号计算获取信息熵、Lempel-Ziv复杂度、母胎能量比以及相对功率比,作为对应通道信号的特征数据。并将得到的特征数据,通过基于支持向量机方法进行训练学习获得腹部心电信号分类模型;
具体地,特征数据计算公式如下包括:
特征1:各通道预处理后的腹部心电信号信息熵的计算,信息熵公式为:
H(X)=-∑χ∈xp(χ)logp(χ)
式中p(χ)代表随机事件X为χ的概率。
特征2:各通道预处理后的腹部心电信号Lempel-Ziv复杂度的计算,Lempel-Ziv复杂度公式为:
式中,r(i)为二进制符号序列,组成符号序列R,x(i)为当前信号幅值,Th为设定好的腹部心电信号阈值;
式中,n是腹部心电信号序列的长度,C(n)是标准化后的Lempel-Ziv复杂度,c(n)是包含于二进制符号序列R的不同子串个数。
特征3:各通道预处理后的腹部心电信号母胎能量比的计算,母胎能量比公式为:
其中,k是当前信号通道序号,l是拐点法确定的母亲QRS波区间长度,s是母亲前后R波中点确定的单波区间长度,E(i)是当前点位的幅值。
特征4:各通道预处理后的腹部心电信号相对功率比的计算,相对功率比公式为:
其中,PSD(k,f)为功率谱密度,k为当前信号通道序号,f为频率。
当有一个通道以上的信号预测为质量可接受,则此时段信号为有效信号;若所有通道均预测为质量不可接受,则此时段信号为异常信号,执行跳出算法操作。
步骤4、对分离得到的胎儿心电信号计算信息熵、Lempel-Ziv复杂度、R波个数等参数作为特征数据,通过支持向量机方法进行训练学习获得胎儿心电信号分类模型;
具体地,计算R波个数的步骤如下:
1)判断信号方向:将4秒数据等分为多份,分别计算各区间内的最大值均值meanMax和最小值均值meanMin,比较二者大小。若最大值均值大于最小值均值,则认定信号为正向;否则认为信号为反向,进行信号反转操作。
2)确定R波个数并去除异常值:确定信号方向后,设定一个信号幅值阈值,选取信号幅值大于阈值的极大值点作为识别到的R波。在识别过程中,可能会检测到两个距离相近的R波点,需要按照生理情况去除异常值,选定一个心率阈值,如果两个R波计算得到的心率高于心率阈值,则选择幅值更高的极大值点作为R波点,去除掉幅值较低的极大值点。
步骤5、对新采集的腹部心电测试信号,截取预设的长度后,根据步骤3计算特征数据,放入训练好的腹部心电信号分类模型获得信号质量的评估结果;
步骤6、对分离得到的胎儿心电测试信号,根据步骤4计算特征数据,放入训练好的胎儿心电信号分类模型获得信号质量的评估结果。
以上所述仅表达了本发明的优选实施方式,其描述较为具体和详细,但并不能因此而理解为对本发明专利范围的限制。应当指出的是,对于本领域的普通技术人员来说,在不脱离本发明构思的前提下,还可以做出若干变形、改进及替代,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。

Claims (1)

1.一种无创胎儿心电信号质量评估方法,其特征在于:包括以下步骤:
步骤1、采集得到多组4路通道的腹部电信号,并对采集到的信号进行预处理;
步骤2、将预处理后的腹部心电信号和分离得到的胎儿心电信号的质量评估结果分为可接受与不可接受两类;
步骤3、对预处理后的腹部心电信号的每一通道信号计算获取信息熵、Lempel-Ziv复杂度、母胎能量比以及相对功率比,作为对应通道信号的特征数据;并将得到的特征数据放入基于支持向量机方法进行训练学习获得腹部心电信号分类模型;
步骤4、对分离得到的胎儿心电信号计算信息熵、Lempel-Ziv复杂度、R波个数参数作为特征数据,放入支持向量机方法进行训练学习获得胎儿心电信号分类模型;
步骤5、对新采集的腹部心电测试信号,截取预设的长度后,根据步骤3计算特征数据,放入训练好的腹部心电信号分类模型获得信号质量的评估结果;
步骤6、对分离得到的胎儿心电测试信号,根据步骤4计算特征数据,放入训练好的胎儿心电信号分类模型获得信号质量的评估结果;
截取长度统一为1000个采样点,腹部心电信号采样频率为250赫兹;
步骤2所述的腹部心电信号预处理方法包括以下步骤:
1):设计滤波器对采集到的腹部心电信号进行高通滤波,得到去除基线漂移干扰的腹部心电信号;
2):采用小波阈值去噪的方式处理腹部心电信号,得到去除高频噪声干扰的腹部心电信号;
3):50Hz陷波器,去除腹部心电信号的工频干扰;
步骤3中所述的特征数据计算具体步骤如下:
特征1:各通道预处理后的腹部心电信号信息熵的计算
信息熵H(X)具体公式如下:
其中,p(χ)代表随机事件X为χ的概率;
特征2:各通道预处理后的腹部心电信号Lempel-Ziv复杂度的计算Lempel-Ziv复杂度是时间序列的非线性特征的测度,反映了时间序列的随机性;具体步骤如下:
其中,r(i)为二进制符号序列,组成符号序列R,x(i)为当前信号幅值,Th为设定好的腹部心电信号阈值;
其中,n是腹部心电信号序列的长度,C(n)是标准化后的Lempel-Ziv复杂度,c(n)是包含于二进制符号序列R的不同子串个数;
特征3:各通道预处理后的腹部心电信号母胎能量比的计算
其中,k是当前信号通道序号,l是拐点法确定的母亲QRS波区间长度,s是母亲前后R波中点确定的单波区间长度,E(i)是当前点位的幅值;
特征4:各通道预处理后的腹部心电信号相对功率比的计算
其中,PSD(k,f)为功率谱密度,k为当前信号通道序号,f为频率;当有一个通道以上的信号预测为质量可接受,则此时段信号为有效信号;若所有通道均预测为质量不可接受,则此时段信号为异常信号,执行跳出算法操作;
步骤4所述的胎儿R波个数计算具体步骤如下:
1)判断信号方向:将4秒数据等分为多份,分别计算各区间内的最大值均值meanMax和最小值均值meanMin,比较二者大小;若最大值均值大于最小值均值,则认定信号为正向;否则认为信号为反向,进行信号反转操作;
2)确定R波个数并去除异常值:确定信号方向后,设定一个信号幅值阈值,选取信号幅值大于阈值的极大值点作为识别到的R波;在识别过程中,可能会检测到两个距离相近的R波点,需要按照生理情况去除异常值,选定一个心率阈值,如果两个R波计算得到的心率高于心率阈值,则选择幅值更高的极大值点作为R波点,去除掉幅值较低的极大值点。
CN202111178668.1A 2021-10-10 2021-10-10 一种无创胎儿心电信号质量评估方法 Active CN113907765B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111178668.1A CN113907765B (zh) 2021-10-10 2021-10-10 一种无创胎儿心电信号质量评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111178668.1A CN113907765B (zh) 2021-10-10 2021-10-10 一种无创胎儿心电信号质量评估方法

Publications (2)

Publication Number Publication Date
CN113907765A CN113907765A (zh) 2022-01-11
CN113907765B true CN113907765B (zh) 2024-02-23

Family

ID=79238850

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111178668.1A Active CN113907765B (zh) 2021-10-10 2021-10-10 一种无创胎儿心电信号质量评估方法

Country Status (1)

Country Link
CN (1) CN113907765B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN115553787A (zh) * 2022-09-30 2023-01-03 哈尔滨理工大学 一种基于多尺度残差收缩U-Net的胎儿心电信号提取方法
CN116407133B (zh) * 2023-06-02 2023-08-29 中国医学科学院阜外医院 心电信号的质量评估方法、装置、存储介质及电子设备

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106874722A (zh) * 2016-12-23 2017-06-20 中国科学院深圳先进技术研究院 一种基于心电信号的身份识别方法及其装置
CN106889981A (zh) * 2017-01-26 2017-06-27 浙江铭众科技有限公司 一种用于提取胎儿心率的智能终端
CN107049259A (zh) * 2017-04-28 2017-08-18 深圳睿京科技有限公司 生理数据的监控方法及装置
CN108013872A (zh) * 2018-01-10 2018-05-11 北京大学第三医院 用于母体胎儿心率监测的系统
CN109117769A (zh) * 2018-07-31 2019-01-01 东南大学 一种针对穿戴式心电采集的实时质量评估反馈方法
CN110101378A (zh) * 2019-04-30 2019-08-09 索思(苏州)医疗科技有限公司 一种由母亲腹部混合ecg信号提取胎儿心率的算法
CN110507317A (zh) * 2019-09-03 2019-11-29 西安邮电大学 一种心电信号r波的自适应ca-cfar定位方法
CN110960211A (zh) * 2019-12-30 2020-04-07 江南大学 一种基于嵌入式的心电实时监测系统
CN111407315A (zh) * 2020-03-24 2020-07-14 浙江大学温州研究院 一种针对超声多普勒胎心信号的质量评估方法
CN112704503A (zh) * 2021-01-15 2021-04-27 苏州大学 心电信号噪声处理方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9492096B2 (en) * 2009-11-03 2016-11-15 Vivaquant Llc ECG sensing apparatuses, systems and methods

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106874722A (zh) * 2016-12-23 2017-06-20 中国科学院深圳先进技术研究院 一种基于心电信号的身份识别方法及其装置
CN106889981A (zh) * 2017-01-26 2017-06-27 浙江铭众科技有限公司 一种用于提取胎儿心率的智能终端
CN107049259A (zh) * 2017-04-28 2017-08-18 深圳睿京科技有限公司 生理数据的监控方法及装置
CN108013872A (zh) * 2018-01-10 2018-05-11 北京大学第三医院 用于母体胎儿心率监测的系统
CN109117769A (zh) * 2018-07-31 2019-01-01 东南大学 一种针对穿戴式心电采集的实时质量评估反馈方法
CN110101378A (zh) * 2019-04-30 2019-08-09 索思(苏州)医疗科技有限公司 一种由母亲腹部混合ecg信号提取胎儿心率的算法
CN110507317A (zh) * 2019-09-03 2019-11-29 西安邮电大学 一种心电信号r波的自适应ca-cfar定位方法
CN110960211A (zh) * 2019-12-30 2020-04-07 江南大学 一种基于嵌入式的心电实时监测系统
CN111407315A (zh) * 2020-03-24 2020-07-14 浙江大学温州研究院 一种针对超声多普勒胎心信号的质量评估方法
CN112704503A (zh) * 2021-01-15 2021-04-27 苏州大学 心电信号噪声处理方法

Also Published As

Publication number Publication date
CN113907765A (zh) 2022-01-11

Similar Documents

Publication Publication Date Title
US7949389B2 (en) Fetal ECG monitoring
US8892181B2 (en) Non-invasive fetal monitoring
Jezewski et al. Determination of fetal heart rate from abdominal signals: evaluation of beat-to-beat accuracy in relation to the direct fetal electrocardiogram
CN113907765B (zh) 一种无创胎儿心电信号质量评估方法
US20130102856A1 (en) Non-invasive detection of fetal or maternal illness
Chiu et al. Using correlation coefficient in ECG waveform for arrhythmia detection
CN108577834B (zh) 一种用于癫痫间期棘波自动检测的方法
Mischi et al. Dedicated entropy measures for early assessment of pregnancy progression from single-channel electrohysterography
CN103263262B (zh) 测量胎儿心率的系统及方法
CN112826513B (zh) 一种基于深度学习和特异性矫正在fecg上胎儿心率检测系统
CN110123304B (zh) 基于多模板匹配和相关系数矩阵的动态心电噪声滤除方法
Fanelli et al. Telefetalcare: A first prototype of a wearable fetal electrocardiograph
CN112971797A (zh) 连续生理信号质量评估方法
Liu et al. Systematic methods for fetal electrocardiographic analysis: Determining the fetal heart rate, RR interval and QT interval
Altini et al. Combining electrohysterography and heart rate data to detect labour
CN102274014A (zh) 一种基于小波去噪的柯氏音血压测量与信号分析方法
Fotiadou et al. Deep convolutional long short-term memory network for fetal heart rate extraction
Cheng et al. Novel multichannel entropy features and machine learning for early assessment of pregnancy progression using electrohysterography
Di Maria et al. An algorithm for the analysis of fetal ECGs from 4-channel non-invasive abdominal recordings
CN103876731B (zh) 一种胎儿心电信号提取装置及方法
Liang et al. Passive fetal movement signal detection system based on intelligent sensing technology
CN112971801A (zh) 连续生理信号质量评估装置
Pani et al. A real-time algorithm for tracking of foetal ECG sources obtained by block-on-line BSS techniques
TWM574469U (zh) Qrs波即時檢測裝置
Orjuela-Cañón et al. Onset and peak detection over pulse wave using supervised SOM network

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
GR01 Patent grant
GR01 Patent grant