[go: up one dir, main page]

CN107684423A - A kind of Fetal ECG separation method and device - Google Patents

A kind of Fetal ECG separation method and device Download PDF

Info

Publication number
CN107684423A
CN107684423A CN201610639848.8A CN201610639848A CN107684423A CN 107684423 A CN107684423 A CN 107684423A CN 201610639848 A CN201610639848 A CN 201610639848A CN 107684423 A CN107684423 A CN 107684423A
Authority
CN
China
Prior art keywords
qrs complex
signal
maternal
fetal
ecg
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
Application number
CN201610639848.8A
Other languages
Chinese (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.)
Shenzhen Institute of Advanced Technology of CAS
Original Assignee
Shenzhen Institute of Advanced Technology of CAS
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 Shenzhen Institute of Advanced Technology of CAS filed Critical Shenzhen Institute of Advanced Technology of CAS
Priority to CN201610639848.8A priority Critical patent/CN107684423A/en
Publication of CN107684423A publication Critical patent/CN107684423A/en
Pending legal-status Critical Current

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/24Detecting, measuring or recording bioelectric or biomagnetic signals of the body or parts thereof
    • A61B5/316Modalities, i.e. specific diagnostic methods
    • 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/346Analysis of electrocardiograms
    • A61B5/349Detecting specific parameters of the electrocardiograph cycle
    • A61B5/366Detecting abnormal QRS complex, e.g. widening
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/43Detecting, measuring or recording for evaluating the reproductive systems
    • A61B5/4306Detecting, measuring or recording for evaluating the reproductive systems for evaluating the female reproductive systems, e.g. gynaecological evaluations
    • A61B5/4343Pregnancy and labour monitoring, e.g. for labour onset detection
    • A61B5/4362Assessing foetal parameters
    • 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/7235Details of waveform analysis
    • A61B5/725Details of waveform analysis using specific filters therefor, e.g. Kalman or adaptive filters

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Public Health (AREA)
  • Molecular Biology (AREA)
  • Veterinary Medicine (AREA)
  • General Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Animal Behavior & Ethology (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Cardiology (AREA)
  • Signal Processing (AREA)
  • Physiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Psychiatry (AREA)
  • Pediatric Medicine (AREA)
  • Pregnancy & Childbirth (AREA)
  • Gynecology & Obstetrics (AREA)
  • Reproductive Health (AREA)
  • Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)

Abstract

本发明提供了一种胎儿心电分离方法及装置。所述方法包括:对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;利用平滑窗对所述QRS波群间的信号段进行平滑处理;对所述QRS波群进行线性平均,生成平均模板;根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。本发明的胎儿心电分离算法简单,且只需要一路母体腹壁心电信号,故可以完成单导联胎儿心电信号分离。

The invention provides a method and device for fetal electrocardiography separation. The method includes: performing QRS wave group detection on the received maternal abdominal wall electrocardiographic signal, determining the starting position of the QRS wave group and the starting position between the QRS wave groups; The segment is smoothed; the QRS wave group is linearly averaged to generate an average template; according to the correlation between the mother QRS wave group in the average template and the mother QRS wave group corresponding to the current ECG cycle, the The maternal QRS complex in the average template is dynamically amended to generate an estimated maternal abdominal wall ECG signal; the maternal abdominal wall ECG estimated signal is used as a reference signal, and the maternal abdominal wall ECG signal is used as an expected signal. An algorithm is used to separate the fetal ECG signal from the maternal abdominal wall ECG signal. The fetal electrocardiographic separation algorithm of the present invention is simple, and only one maternal abdominal wall electrocardiographic signal is needed, so the single-lead fetal electrocardiographic signal separation can be completed.

Description

一种胎儿心电分离方法及装置A method and device for separating fetal heart electricity

技术领域technical field

本发明涉及生物医学工程技术领域,尤其涉及一种胎儿心电分离方法及装置,具体的讲是一种基于平滑窗的胎儿心电分离方法及装置。The present invention relates to the technical field of biomedical engineering, in particular to a method and a device for separating fetal electrocardiograms, in particular to a method and a device for separating electrocardiograms of fetuses based on a smoothing window.

背景技术Background technique

本部分旨在为权利要求书中陈述的本发明的实施方式提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。This section is intended to provide a background or context for implementations of the invention that are recited in the claims. The descriptions herein are not admitted to be prior art by inclusion in this section.

胎儿心电信号(Fetal electrocardiogram,FECG)是胎儿宫内生理活动最主要的客观指标之一,通过分析胎儿心电信号,能够早期诊断胎儿宫内缺氧、宫内窘迫以及先天性心脏病等,有效地降低新生儿各类疾病的发病率和死亡率。因此,围产期监护仪的临床应用将是保证围产期各个阶段的胎儿安全,保障胎儿优生的有力措施之一。Fetal electrocardiogram (FECG) is one of the most important objective indicators of fetal intrauterine physiological activities. By analyzing fetal electrocardiogram signals, it can be used for early diagnosis of fetal hypoxia, intrauterine distress and congenital heart disease. Effectively reduce the morbidity and mortality of various diseases in newborns. Therefore, the clinical application of perinatal monitors will be one of the powerful measures to ensure the safety of the fetus in all stages of the perinatal period and ensure the eugenics of the fetus.

在围产期,胎儿在母体腹部由几个不同的电导率的包裹层包围,如胎盘、羊水、母体子宫和母体腹部皮下组织等,其中羊水的导电能力最高,而胎脂的导电能力最低,此外孕妇的皮肤和皮下脂肪也有微弱的导电性,但是导电能力比肌肉组织约小十倍。在采集胎儿心电时,由于胎儿心电信号需要经过上述导电层,本来就很微弱的胎儿心电信号更容易受外界的干扰。这些不同的组织,形成了所谓的“人体容积导体”,这个容积导体不是一个稳定的导体,其导电性和几何形状随着整个妊娠期不断改变。因此,这将导致母体腹部心电混合信号具有如下特点:1)信噪比低,胎儿心电分量太弱,仅为10~50uV左右,而母亲心电信号则可达750uV,两者可差十倍以上。肌电噪声与其它干扰也很强,信号完全淹没于噪声之中。2)时域中胎儿QRS波与母体QRS波有10~30%的重叠。3)频域中胎儿频谱与母体的频谱大部分重叠。4)信号随机性很强,是非平稳过程,特别是由于胎儿的自由移动造成胎儿心电传输通道的时变性。During the perinatal period, the fetus is surrounded by several layers of different electrical conductivity in the abdomen, such as placenta, amniotic fluid, maternal uterus and subcutaneous tissue of the maternal abdomen, among which amniotic fluid has the highest conductivity, while vernix has the lowest conductivity. In addition, the skin and subcutaneous fat of pregnant women also have weak conductivity, but the conductivity is about ten times smaller than that of muscle tissue. When collecting fetal ECG, since the fetal ECG signal needs to pass through the conductive layer, the already weak fetal ECG signal is more susceptible to external interference. These different tissues form the so-called "body volume conductor", which is not a stable conductor whose conductivity and geometry change continuously throughout gestation. Therefore, this will lead to the following characteristics of the maternal abdominal ECG mixed signal: 1) The signal-to-noise ratio is low, the fetal ECG component is too weak, only about 10-50uV, while the mother’s ECG signal can reach 750uV, the difference between the two can be different. More than ten times. Myoelectric noise and other interference are also very strong, and the signal is completely submerged in the noise. 2) In the time domain, the fetal QRS wave overlaps with the maternal QRS wave by 10-30%. 3) In the frequency domain, the fetal spectrum mostly overlaps with the maternal spectrum. 4) The signal is very random and is a non-stationary process, especially the time-varying nature of the fetal ECG transmission channel due to the free movement of the fetus.

母体腹部心电混合信号的上述特点,为获得高精度的胎儿心电信号制造了比较大的障碍,因此从母体腹壁混合信号中提取出胎儿心电信号,以实现围产期胎儿生命信息的监护,是目前广泛研究和关注的问题。The above-mentioned characteristics of the maternal abdominal ECG mixed signal create relatively large obstacles for obtaining high-precision fetal ECG signals. Therefore, the fetal ECG signal is extracted from the maternal abdominal wall mixed signal to realize the monitoring of fetal life information during the perinatal period , is currently a widely researched and concerned issue.

目前,从母体腹壁混合信号中提取出胎儿心电信号的方法主要有以下几种:At present, there are mainly the following methods for extracting fetal ECG signals from the mixed signals of the maternal abdominal wall:

1)匹配滤波:匹配滤波主要利用的是信号的自相关和互相关等统计特性,但是在实际环境中所要面临的情况远比假设情况要复杂的多,采用匹配滤波获得的胎儿心电信号不是很理想。1) Matched filtering: Matched filtering mainly uses statistical characteristics such as autocorrelation and cross-correlation of the signal, but the situation to be faced in the actual environment is far more complicated than the hypothetical situation. The fetal ECG signal obtained by using matched filtering is not Ideal.

2)独立成分分析法(ICA):此方法需要多路腹壁信号,无法实现单导联胎儿心电信号的提取。2) Independent component analysis (ICA): This method requires multiple abdominal wall signals and cannot realize the extraction of single-lead fetal ECG signals.

3)小波阈值去噪的胎儿心电信号提取方法滤波:利用小波变换结合先验信号知识能够获得效果不错的胎儿心电信号,但是提取过程中容易受到幅度和胎儿心电信号类似信号的干扰,往往导致胎儿心电信号中R波漏检和误检等情况产生,在计算胎儿心率时会引起误差,不利于于临床诊断。3) Fetal ECG signal extraction method filtering with wavelet threshold denoising: Using wavelet transform combined with prior signal knowledge can obtain a good fetal ECG signal, but the extraction process is easily disturbed by similar signals of amplitude and fetal ECG signal. It often leads to missed detection and false detection of the R wave in the fetal ECG signal, which will cause errors when calculating the fetal heart rate, which is not conducive to clinical diagnosis.

4)神经网络法:此方法可以根据神经网络的预测误差检测心电信号的QRS复合波。PT段信号平缓,可以被较好的预测,而对于QRS复合波,相邻点差异较大,会产生预测误差,利用这个预测误差训练网络,对QRS复合波进行预测,预测时的人工神经网络具有多层结构,利用神经网络学习训练过程需要花费较长时间,而且要求训练样本具有一定的代表性,这在实际腹壁心电信号降噪应用中很难实现。4) Neural network method: This method can detect the QRS complex wave of the ECG signal according to the prediction error of the neural network. The PT segment signal is flat and can be predicted better. For the QRS complex wave, the difference between adjacent points is large, which will cause a prediction error. Use this prediction error to train the network to predict the QRS complex wave. The artificial neural network used in the prediction With a multi-layer structure, it takes a long time to use the neural network to learn the training process, and the training samples are required to be representative, which is difficult to achieve in the actual application of abdominal wall ECG signal noise reduction.

5)自适应滤波法:自适应滤波方法包括使用一个或几个参考通道训练自适应滤波消除母体心电图。现有的自适应滤波方法去除母体ECG成分,需要形态上类似于被淹没FECG信号的母体ECG参考通道,但是现有的自适应对消方法要求额外提供参考信号,也就是说,除了腹壁电极外,还需要配备额外的母体胸导连电极,自适应方法对消掉胸导联和腹壁导联中相关的母体分量,这意味着母体腹壁心电信号和母体胸导联心电信号需要相关,才能够被较好的对消掉,并且,母体心电信号的传导路径非常复杂,某些情况下,腹壁中的母体心电信号和胸导联中的母体心电信号相关性差。因此,利用现有的自适应滤波算法对消提取胎儿心电信号会有大量的母体心电泄露,对正确识别胎儿心电分量造成了很大干扰。并且,现有的自适应滤波对消无法实现单导联的胎儿心电信号提取。5) Adaptive filtering method: The adaptive filtering method includes using one or several reference channels to train adaptive filtering to eliminate maternal ECG. Existing adaptive filtering methods to remove maternal ECG components require a maternal ECG reference channel morphologically similar to the submerged FECG signal, but existing adaptive cancellation methods require an additional reference signal, that is, in addition to abdominal wall electrodes , also needs to be equipped with additional electrodes connected to the maternal chest lead, and the adaptive method can eliminate the relevant maternal component in the chest lead and abdominal wall lead, which means that the maternal abdominal wall ECG signal and the maternal chest lead ECG signal need to be correlated, In addition, the conduction path of the maternal ECG signal is very complicated. In some cases, the correlation between the maternal ECG signal in the abdominal wall and the maternal ECG signal in the chest leads is poor. Therefore, using the existing adaptive filtering algorithm to cancel and extract the fetal ECG signal will cause a large amount of maternal ECG leakage, which will cause great interference to the correct identification of fetal ECG components. Moreover, the existing adaptive filter cancellation cannot realize single-lead fetal ECG signal extraction.

发明内容Contents of the invention

本发明提出一种胎儿心电分离方法及装置,用以解决现有的胎儿心电分离算法计算复杂且无法实现单导联的问题。The present invention proposes a fetal electrocardiographic separation method and device, which are used to solve the problem that the existing fetal electrocardiographic separation algorithm is complex in calculation and unable to realize a single lead.

为了达到上述目的,本发明实施例提供一种胎儿心电分离方法,包括:对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;利用平滑窗对所述QRS波群间的信号段进行平滑处理;对所述QRS波群进行线性平均,生成平均模板;根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。In order to achieve the above object, an embodiment of the present invention provides a fetal electrocardiographic separation method, including: performing QRS wave group detection on the received maternal abdominal wall electrocardiogram signal, determining the starting position of the QRS wave group and the starting point between the QRS wave groups start position; use the smoothing window to smooth the signal segments between the QRS complexes; carry out linear averaging to the QRS complexes to generate an average template; Correlation of the maternal QRS wave group corresponding to the electrical cycle, modifying the dynamic amplitude of the mother's QRS wave group in the average template to generate an estimated maternal abdominal wall electrocardiogram signal; using the maternal abdominal wall electrocardiographic estimated signal as a reference signal, Taking the maternal abdominal wall electrocardiographic signal as an expected signal, and separating the fetal electrocardiographic signal from the maternal abdominal wall electrocardiographic signal through an adaptive algorithm.

为了达到上述目的,本发明实施例还提供一种胎儿心电分离装置,包括:波群检测模块,用于对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;平滑处理模块,用于利用平滑窗对所述QRS波群间的信号段进行平滑处理;平均处理模块,用于对所述QRS波群进行线性平均,生成平均模板;波群修正模块,用于根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;分离模块,用于将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。In order to achieve the above object, the embodiment of the present invention also provides a fetal electrocardiogram separation device, including: a wave group detection module, which is used to detect the QRS wave group of the received maternal abdominal wall ECG signal, and determine the start of the QRS wave group The starting position between the position and the QRS complex; the smoothing processing module, which is used to smooth the signal segment between the QRS complexes by using a smoothing window; the average processing module, which is used to linearly average the QRS complexes, Generating an average template; a wave group correction module, used to modify the maternal QRS wave group in the average template according to the correlation between the mother QRS wave group in the average template and the mother QRS wave group corresponding to the current ECG cycle Perform dynamic amplitude correction to generate an estimated maternal abdominal wall ECG signal; a separation module is configured to use the maternal abdominal wall ECG estimated signal as a reference signal, and use the maternal abdominal wall ECG signal as an expected signal, through an adaptive algorithm, from The fetal electrocardiographic signal is separated from the maternal abdominal wall electrocardiographic signal.

本发明提出的胎儿心电分离方法及装置,通过平滑母体QRS波群间信号,消除母体QRS波群间信号段的胎儿分量,以及对母体QRS波群进行幅值修正,从而估计母体腹壁心电信号,此方法避免了胸导联和腹壁导联中母体分量非相关所带来的母体心电抑制不彻底的影响,具有更好的分离效果;本发明的实现只需要一路母体腹壁心电信号,故可以完成单导联胎儿心电信号分离;并且,本发明的胎儿分离算法设计思路简单直观,针对于异常的临床数据,易于查找原因和添加修正模块。The fetal electrocardiogram separation method and device proposed by the present invention can estimate the maternal abdominal wall ECG by smoothing the maternal QRS inter-group signal, eliminating the fetal component of the maternal QRS inter-group signal segment, and correcting the amplitude of the maternal QRS wave group. signal, this method avoids the influence of incomplete maternal ECG suppression caused by the non-correlation of the maternal component in the chest lead and the abdominal wall lead, and has a better separation effect; the realization of the present invention only needs one maternal abdominal wall ECG signal , so single-lead fetal ECG signal separation can be completed; moreover, the design idea of the fetal separation algorithm of the present invention is simple and intuitive, and it is easy to find the cause and add a correction module for abnormal clinical data.

附图说明Description of drawings

为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings that need to be used in the description of the embodiments or the prior art. Obviously, the accompanying drawings in the following description are only These are some embodiments of the present invention. For those skilled in the art, other drawings can also be obtained according to these drawings on the premise of not paying creative efforts.

图1为本发明实施例的胎儿心电分离方法的系统结构示意图;1 is a schematic diagram of the system structure of a method for fetal electrocardiography separation according to an embodiment of the present invention;

图2为本发明实施例的基于平滑窗的腹壁母体心电估计机制的示意框图;Fig. 2 is the schematic block diagram of the abdominal wall maternal ECG estimation mechanism based on the smoothing window of the embodiment of the present invention;

图3为本发明实施例的胎儿心电分离方法的处理流程图;Fig. 3 is the processing flow diagram of the method for fetal electrocardiography separation according to the embodiment of the present invention;

图4为本发明实施例的胎儿心电分离装置的结构示意图;4 is a schematic structural view of a fetal electrocardiography separation device according to an embodiment of the present invention;

图5为图4所示实施例的波群修正模块104的结构示意图;FIG. 5 is a schematic structural diagram of the wave group correction module 104 of the embodiment shown in FIG. 4;

图6为本发明的具体实施例采用的5路母体腹壁混合心电信号的波形图;Fig. 6 is the oscillogram of the mixed electrocardiogram of 5 road maternal abdominal wall that the specific embodiment of the present invention adopts;

图7为本发明的具体实施例的实验结果分析图。Fig. 7 is an analysis diagram of experimental results of a specific embodiment of the present invention.

具体实施方式detailed description

下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。The following will clearly and completely describe the technical solutions in the embodiments of the present invention with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the described embodiments are only some, not all, embodiments of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by persons of ordinary skill in the art without making creative efforts belong to the protection scope of the present invention.

本领域技术技术人员知道,本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此,本公开可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。Those skilled in the art know that the embodiments of the present invention can be implemented as a system, device, device, method or computer program product. Therefore, the present disclosure may be embodied in the form of complete hardware, complete software (including firmware, resident software, microcode, etc.), or a combination of hardware and software.

下面参考本发明的若干代表性实施方式,详细阐释本发明的原理和精神。The principle and spirit of the present invention will be explained in detail below with reference to several representative embodiments of the present invention.

本发明提出了一种基于平滑窗的自适应单导联胎儿心电分离方法及装置,对母体腹壁信号中的QRS波群间信号采用平滑窗进行平滑处理,以平滑掉QRS波群间信号段的胎儿心电节拍点,对腹壁心电信号中的母体QRS波群分割提取,对提取后的N个连续母体QRS波群以母体R波波峰位置为定位点,进行线性平均,生成平均QRS波群模板,截取当前输入的腹壁信号中的母体QRS波群和平均模板中的母体QRS波群,利用二者的相关系数,对平均模板中的母体QRS波群进行复制修正,以得到动态的母体腹壁心电估计信号。并且,通过自适应自主学习,对消掉腹壁信号中的母体分量,该方案避免了胸导联中母体分量和腹壁导联中母体分量非相关性带来的母体分量泄露影响,同时可以实现单导联胎儿心电分离。The present invention proposes an adaptive single-lead fetal electrocardiographic separation method and device based on a smoothing window, which uses a smoothing window to smooth the QRS inter-group signal in the maternal abdominal wall signal, so as to smooth out the QRS inter-group signal segment The maternal QRS wave group in the abdominal wall ECG signal is segmented and extracted, and the extracted N consecutive maternal QRS wave groups are linearly averaged with the peak position of the maternal R wave as the positioning point to generate an average QRS wave Group template, which intercepts the maternal QRS complex in the currently input abdominal wall signal and the maternal QRS complex in the average template, and uses the correlation coefficient between the two to copy and correct the maternal QRS complex in the average template to obtain a dynamic maternal Abdominal wall ECG estimation signal. Moreover, through self-adaptive learning, the maternal component in the abdominal wall signal is eliminated. This scheme avoids the influence of the maternal component leakage caused by the non-correlation between the maternal component in the chest lead and the maternal component in the abdominal wall lead. At the same time, it can achieve a single Lead fetal electrocardiogram separation.

图1为本发明实施例的胎儿心电分离方法的系统结构示意图。如图1所示,y(n)表示采集到的母体腹壁心电信号,f(n)为腹壁信号中的胎儿心电信号,x(n)为腹壁信号中的母体心电信号,为自适应滤波器的期望信号,x'(n)为通过腹壁母体心电估计机制估计的腹壁母体心电,自适应滤波器对输入信号不断的进行“学习”和“跟踪”,使输出信号最大限度的逼近腹壁心电信号,从而使误差信号最小,即为分离出的胎儿心电信号。FIG. 1 is a schematic diagram of the system structure of the fetal electrocardiography separation method according to the embodiment of the present invention. As shown in Figure 1, y(n) represents the collected maternal abdominal wall ECG signal, f(n) represents the fetal ECG signal in the abdominal wall signal, x(n) represents the maternal ECG signal in the abdominal wall signal, and is the The expected signal of the adaptive filter, x'(n) is the abdominal wall maternal ECG estimated by the abdominal wall maternal ECG estimation mechanism, and the adaptive filter continuously "learns" and "tracks" the input signal, so that the output signal Maximize the approximation of the abdominal wall ECG signal, so that the error signal minimum, That is, the isolated fetal ECG signal.

也就是说,基于自适应的胎儿心电分离算法,需要两路输入信号,即:一路母体腹壁心电信号(主输入信号),一路估计出的腹壁母体心电信号(参考输入信号),这相对于现有技术中的需要接入一路胸导联心电信号作为参考输入信号,实现了单导联。That is to say, the adaptive fetal ECG separation algorithm requires two input signals, namely: one maternal abdominal ECG signal (main input signal), and one estimated abdominal wall maternal ECG signal (reference input signal). Compared with the prior art that requires access to one ECG signal of a chest lead as a reference input signal, a single lead is realized.

本发明的核心是如何计算出腹壁母体心电估计信号,以作为自适应算法的参考输入信号。图2为本发明的腹壁母体心电估计机制的示意框图,如图2所示,将QRS波群信号和QRS波群间信号分开处理,采用基于平滑窗的方法进行母体QRS波群间信号进行平滑处理,消除母体QRS波群间信号段的胎儿心电节拍点,利用平均相关法对母体QRS波群进行幅值修正,从而估计出腹壁母体信号波形。The core of the present invention is how to calculate the estimated abdominal wall maternal electrocardiogram signal as the reference input signal of the adaptive algorithm. Fig. 2 is a schematic block diagram of the abdominal wall maternal ECG estimation mechanism of the present invention. As shown in Fig. 2, the QRS wave group signal and the inter-QRS inter-group signal are processed separately, and the method based on the smoothing window is used to carry out the maternal QRS inter-group signal. The smoothing process eliminates the fetal ECG beat points in the signal segment between the maternal QRS complexes, and uses the average correlation method to correct the amplitude of the maternal QRS complexes, thereby estimating the abdominal wall maternal signal waveform.

图3为本发明实施例的胎儿心电分离方法的处理流程图。如图3所示,包括:Fig. 3 is a processing flow chart of the fetal electrocardiographic separation method according to the embodiment of the present invention. As shown in Figure 3, including:

步骤S101,对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;Step S101, performing QRS wave group detection on the received maternal abdominal wall ECG signal, and determining the starting position of the QRS wave group and the starting position between the QRS wave groups;

步骤S102,利用平滑窗对所述QRS波群间的信号段进行平滑处理;Step S102, using a smoothing window to smooth the signal segments between the QRS complexes;

步骤S103,对所述QRS波群进行线性平均,生成平均模板;Step S103, linearly averaging the QRS complexes to generate an average template;

步骤S104,根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;Step S104, according to the correlation between the maternal QRS complex in the average template and the maternal QRS complex corresponding to the current ECG cycle, perform dynamic amplitude correction on the maternal QRS complex in the average template to generate a maternal Abdominal wall ECG estimation signal;

步骤S105,将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。Step S105, using the estimated maternal abdominal wall ECG signal as a reference signal and the maternal abdominal wall ECG signal as an expected signal, and using an adaptive algorithm to separate the fetal ECG signal from the maternal abdominal wall ECG signal.

在步骤S102中,采用固定长度窗,对所述QRS波群间的信号段进行平滑处理,所述固定长度窗的取值根据信号采样率和胎儿心电波群长度来确定。In step S102, the signal segment between the QRS complexes is smoothed by using a fixed-length window, and the value of the fixed-length window is determined according to the signal sampling rate and the length of the fetal electrocardiogram group.

如图2所示,用固定长度窗M,对QRS波群间信号段进行平滑处理,固定窗长度约为胎儿心电QRS波长度,胎儿心电波群长度一般为47ms至85ms,故M取值可为(0.047*d,0.085*d),其中参数d为信号采样率。As shown in Figure 2, a fixed length window M is used to smooth the signal segment between QRS wave groups. The length of the fixed window is about the length of the fetal ECG QRS wave. The length of the fetal ECG wave group is generally 47ms to 85ms, so the value of M It can be (0.047*d, 0.085*d), where the parameter d is the signal sampling rate.

QRS波群间信号可表示为:母体心电幅值一般是胎儿心电幅值的十几倍,经过平滑后,母体QRS波群间几乎不再含有胎儿QRS波群。The inter-QRS signal can be expressed as: The amplitude of the maternal ECG is generally more than ten times that of the fetal ECG. After smoothing, there is almost no fetal QRS complex between the maternal QRS complexes.

在步骤S103中,对所述QRS波群进行线性平均,生成平均片段,包括:从所述母体腹壁心电信号中,提取出母体心电R波波峰位置;以及以所述R波波峰位置为中心定位点,对所述QRS波群进行线性平均,生成所述平均模板。In step S103, the QRS wave group is linearly averaged to generate an average segment, including: extracting the maternal ECG R wave peak position from the maternal abdominal wall ECG signal; and taking the R wave peak position as The central location point is used to linearly average the QRS complex to generate the average template.

具体实施时,需要确定母体腹壁心电信号中的R波波峰位置,该方法为:对接收到的母体腹壁心电信号进行差分运算,以获取腹壁心电信号的斜率曲线,对斜率曲线做绝对滑动积分运算,记录每秒钟内的最大值,对前8秒内的最大值做平均,则R波判定阈值=系数比值*前8秒内的最大值平均,对差分值大于R波判定阈值的段,进行R波波峰位置最大值查找,以确定最终的R波波峰位置点。During specific implementation, it is necessary to determine the peak position of the R wave in the maternal abdominal wall ECG signal. The method is: perform differential calculation on the received maternal abdominal wall ECG signal to obtain the slope curve of the abdominal wall ECG signal, and perform an absolute calculation of the slope curve. Sliding integral operation, record the maximum value in each second, and average the maximum value in the first 8 seconds, then the R wave judgment threshold = coefficient ratio * average of the maximum value in the previous 8 seconds, if the difference value is greater than the R wave judgment threshold In the segment, the maximum value of the R wave peak position is searched to determine the final R wave peak position point.

然后以腹壁心电信号中R波波峰位置为中心定位点,截取母体QRS波群心电波形y1,y1=[y1(1),y1(2),...,y1(P)],其中P为母体QRS波群长度。Then take the peak position of the R wave in the abdominal wall ECG signal as the central positioning point, and intercept the maternal QRS complex ECG waveform y 1 , y 1 =[y 1 (1),y 1 (2),...,y 1 ( P)], where P is the length of the maternal QRS complex.

对当前心电周期对应的心电信号片段以及前面M-1个连续的心电信号片段进行线性平均,得到平均模板为: Linearly average the ECG signal segment corresponding to the current ECG cycle and the previous M-1 consecutive ECG signal segments, and the average template is obtained as:

在本发明实施例中,能否有效的抑制掉母体QRS波群分量,对正确识别分离结果中的胎儿心电分量起着决定性作用。平均模板中的母体心电波形,是对较早的心电波形和最新波形的等权累加和线性平均,估计出的母体心电片段不具有动态性。因为,本发明的步骤S104根据最新母体QRS波群和平均片段中母体QRS波群的相关性,对平均模板中的母体QRS波群进行动态幅值修正。In the embodiment of the present invention, whether the maternal QRS complex component can be effectively suppressed plays a decisive role in correctly identifying the fetal ECG component in the separation result. The maternal ECG waveform in the average template is the equal-weighted accumulation and linear average of the earlier ECG waveform and the latest waveform, and the estimated maternal ECG segment is not dynamic. Because, step S104 of the present invention performs dynamic amplitude correction on the maternal QRS complex in the average template according to the correlation between the latest maternal QRS complex and the maternal QRS complex in the average segment.

具体实施时,方法如下:In specific implementation, the method is as follows:

首先,计算所述当前心电周期对应的QRS波群片段与所述平均模板中的QRS波群片段的相关性,生成相关系数,具体步骤为:First, calculate the correlation between the QRS complex segment corresponding to the current electrocardiographic cycle and the QRS complex segment in the average template, and generate a correlation coefficient. The specific steps are:

(1)提取当前心电周期中的母体QRS波群片段:y3=[y3(1),y3(2),...,y3(P)],其中P为母体QRS波群片段的长度;(1) Extract the maternal QRS complex segment in the current ECG cycle: y 3 =[y 3 (1), y 3 (2),...,y 3 (P)], where P is the maternal QRS complex the length of the segment;

(2)提取出平均模板中母体QRS波群片段:y2=[y2(1),y2(2),...,y2(P)];(2) Extract the maternal QRS complex segment in the average template: y 2 =[y 2 (1),y 2 (2),...,y 2 (P)];

(3)计算量片段的相关系数r23(3) Calculate the correlation coefficient r 23 of the volume segment:

其次,根据所述相关系数,用所述当前心电周期对应的QRS波群片段对所述平均模板中的母体QRS波群进行动态幅值修正,所采用的表达式如下:Secondly, according to the correlation coefficient, use the QRS wave group segment corresponding to the current electrocardiographic cycle to carry out dynamic amplitude correction to the maternal QRS wave group in the average template, and the adopted expression is as follows:

由于平均段生成的平均模板具有当前腹壁心电周期内所对应母体心电波形的信号,并且基本消除掉了其中的胎儿心电波形,通过相关系数,对平均后的母体心电QRS波群段进行动态幅值修正,能最大限度的估计母体腹壁心电信号。Since the average template generated by the average section has the signal of the corresponding maternal ECG waveform in the current abdominal wall ECG cycle, and basically eliminated the fetal ECG waveform, the averaged maternal ECG QRS wave group segment The dynamic amplitude correction can estimate the maternal abdominal wall ECG signal to the greatest extent.

在步骤S105中,腹部电极取得主输入信号,将估计出的母体腹壁心电信号作为输入信号,此信号中基本不含有胎儿心电信号,只包含欲消除的噪声和母体心电信号;腹壁心电信号为期望信号,通过自适应算法,消除了输入信号与期望信号中相关的部分,保留不相关的部分,因此腹壁心电信号中的胎儿心电信号分离可以被最大化的提取出来。In step S105, the abdominal electrode obtains the main input signal, and takes the estimated maternal abdominal wall ECG signal as the input signal, which basically does not contain the fetal ECG signal, but only includes the noise to be eliminated and the maternal ECG signal; The electrical signal is an expected signal. Through an adaptive algorithm, the relevant part of the input signal and the expected signal is eliminated, and the irrelevant part is retained. Therefore, the separation of the fetal ECG signal in the abdominal wall ECG signal can be extracted to the maximum extent.

本发明实施例中,自适应算法可以为自适应的最小均方(LMS)算法、自适应的递归最小二乘方(RLS)算法、最小二乘格形(LSL)自适应算法等。In the embodiment of the present invention, the adaptive algorithm may be an adaptive least mean square (LMS) algorithm, an adaptive recursive least squares (RLS) algorithm, a least squares lattice (LSL) adaptive algorithm, and the like.

应当注意,尽管在附图中以特定顺序描述了本发明方法的操作,但是,这并非要求或者暗示必须按照该特定顺序来执行这些操作,或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地,可以省略某些步骤,将多个步骤合并为一个步骤执行,和/或将一个步骤分解为多个步骤执行。It should be noted that, although operations of the methods of the present invention are depicted in the drawings in a particular order, this does not require or imply that the operations must be performed in that particular order, or that all illustrated operations must be performed to achieve the desired results. . Additionally or alternatively, certain steps may be omitted, multiple steps may be combined into one step for execution, and/or one step may be decomposed into multiple steps for execution.

在介绍了本发明示例性实施方式的方法之后,接下来,参考图6对本发明示例性实施方式的胎儿心电分离装置进行介绍。该装置的实施可以参见上述方法的实施,重复之处不再赘述。以下所使用的术语“模块”和“单元”,可以是实现预定功能的软件和/或硬件。尽管以下实施例所描述的模块较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。After introducing the method of the exemplary embodiment of the present invention, next, the fetal electrocardiographic separation device of the exemplary embodiment of the present invention will be introduced with reference to FIG. 6 . For the implementation of the device, reference may be made to the implementation of the above method, and repeated descriptions will not be repeated. The terms "module" and "unit" used below may be software and/or hardware that realize predetermined functions. Although the modules described in the following embodiments are preferably implemented in software, implementations in hardware, or a combination of software and hardware are also possible and contemplated.

图4为本发明实施例的胎儿心电分离装置的结构示意图。如图所示,包括:Fig. 4 is a schematic structural diagram of a fetal electrocardiography separation device according to an embodiment of the present invention. As shown, including:

波群检测模块101,用于对接收到的母体腹壁心电信号进行QRS波群检测,确定QRS波群的起始位置和QRS波群间的起始位置;平滑处理模块102,用于利用平滑窗对所述QRS波群间的信号段进行平滑处理;平均处理模块103,用于对所述QRS波群进行线性平均,生成平均模板;波群修正模块104,用于根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号;分离模块105,用于将所述母体腹壁心电估计信号作为参考信号,将所述母体腹壁心电信号作为期望信号,通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号。The wave group detection module 101 is used for performing QRS wave group detection on the received maternal abdominal wall electrocardiographic signal to determine the starting position of the QRS wave group and the starting position between the QRS wave groups; the smoothing processing module 102 is used to use smoothing The window smoothes the signal segments between the QRS complexes; the average processing module 103 is used to linearly average the QRS complexes to generate an average template; Correlation between the maternal QRS wave group and the maternal QRS wave group corresponding to the current ECG cycle, the dynamic amplitude correction is performed on the maternal QRS wave group in the average template, and the maternal abdominal wall ECG estimation signal is generated; the separation module 105 , for using the estimated maternal abdominal wall ECG signal as a reference signal and the maternal abdominal wall ECG signal as an expected signal, and using an adaptive algorithm to separate the fetal ECG signal from the maternal abdominal wall ECG signal.

进一步地,在本发明实施例中,所述平滑处理模块102利用平滑窗对所述QRS波群间的信号段进行平滑处理,具体包括:采用固定长度窗,对所述QRS波群间的信号段进行平滑处理,所述固定长度窗的取值根据信号采样率和胎儿心电波群长度来确定。Further, in the embodiment of the present invention, the smoothing processing module 102 performs smoothing processing on the signal segments between the QRS complexes by using a smoothing window, which specifically includes: using a fixed length window to smooth the signal segments between the QRS complexes The segment is smoothed, and the value of the fixed-length window is determined according to the signal sampling rate and the length of the fetal electrocardiogram group.

进一步地,在本发明实施例中,所述平均处理模块103对所述QRS波群进行线性平均,生成平均模板,包括:从所述母体腹壁心电信号中,提取出母体心电R波波峰位置;以所述R波波峰位置为中心定位点,对所述QRS波群进行线性平均,生成所述平均模板。Further, in the embodiment of the present invention, the average processing module 103 performs linear average on the QRS complex to generate an average template, including: extracting the maternal ECG R wave peak from the maternal abdominal wall ECG signal Position: taking the position of the R wave peak as the central positioning point, performing linear average on the QRS complex to generate the average template.

进一步地,在本发明实施例中,所述波群修正模块104根据所述平均模板中的母体QRS波群与所述当前心电周期对应的母体QRS波群的相关性,对所述平均模板中的母体QRS波群行动态幅值修正,生成母体腹壁心电估计信号,如图5所示,具体包括:Further, in the embodiment of the present invention, the complex correction module 104 adjusts the average template according to the correlation between the maternal QRS complex in the average template and the maternal QRS complex corresponding to the current ECG cycle. The dynamic amplitude correction of the maternal QRS complex is performed to generate the maternal abdominal wall ECG estimation signal, as shown in Figure 5, which specifically includes:

相关系数计算单元1041,用于计算所述当前心电周期对应的QRS波群片段与所述平均模板中的QRS波群片段的相关性,生成相关系数;A correlation coefficient calculation unit 1041, configured to calculate the correlation between the QRS complex segment corresponding to the current ECG cycle and the QRS complex segment in the average template, and generate a correlation coefficient;

动态修正单元1042,用于根据所述相关系数,用所述当前心电周期对应的QRS波群片段对所述平均模板中的母体QRS波群进行动态幅值修正,生成母体腹壁心电估计信号。The dynamic correction unit 1042 is configured to use the QRS wave group segment corresponding to the current ECG cycle to perform dynamic amplitude correction on the maternal QRS wave group in the average template according to the correlation coefficient, so as to generate a maternal abdominal wall ECG estimation signal .

进一步地,在本发明实施例中,所述分离模块105通过自适应算法,从所述母体腹壁心电信号中分离出胎儿心电信号,所述自适应算法包括自适应的最小均方算法、自适应的递归最小二乘方算法、最小二乘格形自适应算法。Further, in the embodiment of the present invention, the separation module 105 separates the fetal ECG signal from the maternal abdominal wall ECG signal through an adaptive algorithm, and the adaptive algorithm includes an adaptive least mean square algorithm, Adaptive recursive least squares algorithm, least squares lattice adaptive algorithm.

此外,尽管在上文详细描述中提及了胎儿心电分离装置的若干单元,但是这种划分仅仅并非强制性的。实际上,根据本发明的实施方式,上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。同样,上文描述的一个单元的特征和功能也可以进一步划分为由多个单元来具体化。Furthermore, although several units of the fetal electrocardioseparation device are mentioned in the above detailed description, this division is merely not mandatory. Actually, according to the embodiment of the present invention, the features and functions of two or more units described above may be embodied in one unit. Likewise, the features and functions of one unit described above can also be further divided to be embodied by a plurality of units.

以下通过一个具体实施例来说明采用本发明的基于平滑窗的胎儿心电分离方法实现胎儿心电信号分离的效果。A specific example will be used below to illustrate the effect of using the method for separating fetal ECG signals based on the smoothing window of the present invention to separate fetal ECG signals.

参看图6所示,为该具体实施例采用的5路母体腹壁混合心电信号的波形图。如图6所示,每路信号的数据长度为2500个样本点,采样率250Hz。由图6可知,第1路腹壁信号为胎母心电信噪比最大的一路信号,即此路信号效果最好,第4路和第5路腹壁心电信号用肉眼已经无法识别出胎儿心电R波峰位置,此两路信号提取效果已不能为本发明的算法所用,故不在本发明的算法的考虑范围之内。Referring to FIG. 6 , it is a waveform diagram of 5 maternal abdominal wall mixed ECG signals used in this specific embodiment. As shown in Figure 6, the data length of each signal is 2500 sample points, and the sampling rate is 250Hz. It can be seen from Figure 6 that the first abdominal wall signal is the one with the largest fetal-maternal ECG signal-to-noise ratio, that is, this signal has the best effect, and the fourth and fifth abdominal wall ECG signals cannot be recognized by the naked eye. The position of the electric R wave peak, the extraction effect of these two signals can no longer be used by the algorithm of the present invention, so it is not within the scope of consideration of the algorithm of the present invention.

本实施例选择第2路信号作为算法验证数据,自适应滤波采用自适应递归最小二乘方算法,其实验结果分析参看图7。In this embodiment, the second channel signal is selected as the algorithm verification data, and the self-adaptive recursive least squares algorithm is used for the adaptive filtering. For the analysis of the experimental results, refer to FIG. 7 .

1、胎儿QRS波群远离母体QRS波群:1. The fetal QRS complex is far from the maternal QRS complex:

胎儿QRS波群远离母体QRS波群,如图7中标注a,由图7可知,这种情况下的胎儿QRS波分量被提取出,左右的母体QRS波分量被抑制掉。此种胎儿QRS波群和母体QRS波群的位置关系,是腹壁心电中最普遍的位置关系,也是最容易提取出的位置关系。The fetal QRS complex is far away from the maternal QRS complex, as marked a in Figure 7. It can be seen from Figure 7 that the fetal QRS complex components in this case are extracted, and the left and right maternal QRS complex components are suppressed. The positional relationship between the fetal QRS complex and the maternal QRS complex is the most common positional relationship in the abdominal wall ECG, and it is also the easiest positional relationship to extract.

2、胎儿QRS波群紧邻母体QRS波群:2. The fetal QRS complex is next to the maternal QRS complex:

胎儿QRS波群紧邻母体QRS波群,如图7中标注b所示,由图7可知,这种情况下的胎儿QRS波分量被提取出,紧邻的母体QRS波分量被抑制掉,并未对胎儿心电QRS波群的识别造成影响。此种胎儿QRS波群和母体QRS波群的位置关系,是腹壁心电中少数比例存在的位置关系,也是较难提取出的位置关系。Fetal QRS complexes are adjacent to maternal QRS complexes, as shown by b in Figure 7. It can be seen from Figure 7 that the fetal QRS complex components in this case are extracted, and the adjacent maternal QRS complex components are suppressed. The identification of fetal ECG QRS complex is affected. This positional relationship between the fetal QRS complex and the maternal QRS complex is a positional relationship that exists in a small proportion of abdominal wall ECG, and it is also a positional relationship that is difficult to extract.

3、胎儿QRS波群和母体QRS波群完全重合:3. Fetal QRS complexes completely overlap with maternal QRS complexes:

胎儿QRS波群完全和母体QRS波群重合,如图7中标注c所示,由图7可知,这种情况下的胎儿QRS波分量被提取出,重合的母体QRS波分量被抑制掉。此种胎儿QRS波群和母体QRS波群的位置关系,是腹壁心电中少数比例存在的位置关系,也是最难提取出的位置关系。The fetal QRS complex completely overlaps with the maternal QRS complex, as shown by c in Figure 7. It can be seen from Figure 7 that the fetal QRS complex components are extracted in this case, and the overlapping maternal QRS complex components are suppressed. This positional relationship between the fetal QRS complex and the maternal QRS complex is the positional relationship that exists in a small proportion of the abdominal wall ECG, and it is also the most difficult positional relationship to extract.

由以上结果分析可知,通过本发明的基于平滑窗的胎儿心电分离方法,胎儿QRS波群远离母体QRS波群、胎儿QRS波群紧邻母体QRS波群、胎儿QRS波群和母体QRS波群完全重合三种情况,均有非常理想的提取效果。From the analysis of the above results, it can be known that by the method for fetal electrocardiogram separation based on the smoothing window of the present invention, the fetal QRS complexes are far away from the maternal QRS complexes, the fetal QRS complexes are close to the maternal QRS complexes, and the fetal QRS complexes are completely separated from the maternal QRS complexes. The three situations overlap, all of which have very ideal extraction effects.

本发明提出的胎儿心电分离方法及装置,具有以下有益效果:The fetal electrocardiographic separation method and device proposed by the present invention have the following beneficial effects:

1、单导联:算法仅需要一路母体腹壁心电信号,即可完成基于自适应的胎儿信号分离。传统的采用基于自适应的胎儿心电信号分离算法,需要两路信号,即:一路腹壁心电信号(主输入信号),一路胸导联心电信号(参考输入信号),通过对主输入信号进行母体QRS波群间平滑处理,消除母体QRS波群间信号段的胎儿分量,对母体QRS波群进行幅值修正,从而估计出腹壁母体信号波形,信号被输入自适应滤波器,同时未处理的同步腹壁心电信号作为自适应滤波的期望信号,因此,本发明只需要一路腹壁心电信号,这将大大简化了系统的导联系统。1. Single lead: The algorithm only needs one maternal abdominal wall ECG signal to complete the separation of fetal signals based on self-adaptation. The traditional adaptive fetal ECG signal separation algorithm requires two signals, namely: one abdominal wall ECG signal (main input signal), one chest lead ECG signal (reference input signal), through the main input signal Perform smoothing between maternal QRS complexes, eliminate the fetal component of the signal segment between maternal QRS complexes, and correct the amplitude of the maternal QRS complexes, thereby estimating the maternal signal waveform of the abdominal wall. The signal is input into an adaptive filter and is not processed The synchronous abdominal wall ECG signal is used as the expected signal for adaptive filtering. Therefore, the present invention only needs one abdominal wall ECG signal, which greatly simplifies the lead system of the system.

2、分离效果好:通过平滑和母体QRS波波群幅值修正处理,估计母体腹壁心电信号,此方法避免了胸导联和腹壁导联中母体分量非相关所带来的母体心电抑制不彻底的影响,具有更好的分离效果。2. Good separation effect: through smoothing and maternal QRS wave group amplitude correction processing, the maternal abdominal wall ECG signal is estimated. This method avoids the maternal ECG inhibition caused by the non-correlation of the maternal component in the chest lead and abdominal wall lead Incomplete impact, with better separation effect.

3、算法设计思路简单直观,针对于异常的临床数据,易于查找原因和添加修正模块。3. The design idea of the algorithm is simple and intuitive. For abnormal clinical data, it is easy to find the cause and add a correction module.

本发明实施例还提供一种计算机可读程序,其中当执行所述程序时,所述程序使得计算机执行本发明实施例所述的胎儿心电分离方法。An embodiment of the present invention also provides a computer-readable program, wherein when the program is executed, the program causes the computer to execute the method for fetal electrocardiography separation described in the embodiment of the present invention.

本发明实施例还提供一种存储有计算机可读程序的存储介质,其中所述计算机可读程序使得计算机执行本发明实施例所述的胎儿心电分离方法。An embodiment of the present invention also provides a storage medium storing a computer-readable program, wherein the computer-readable program enables a computer to execute the method for fetal electrocardiography separation described in the embodiment of the present invention.

本发明以上的装置和方法可以由硬件实现,也可以由硬件结合软件实现。本发明涉及这样的计算机可读程序,当该程序被逻辑部件所执行时,能够使该逻辑部件实现上文所述的装置或构成部件,或使该逻辑部件实现上文所述的各种方法或步骤。本发明还涉及用于存储以上程序的存储介质,如硬盘、磁盘、光盘、DVD、flash存储器等。The above devices and methods of the present invention can be implemented by hardware, or by combining hardware and software. The present invention relates to such a computer-readable program that, when the program is executed by a logic component, enables the logic component to realize the above-mentioned device or constituent component, or enables the logic component to realize the above-mentioned various methods or steps. The present invention also relates to a storage medium for storing the above program, such as hard disk, magnetic disk, optical disk, DVD, flash memory and the like.

本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。Those skilled in the art should understand that the embodiments of the present invention may be provided as methods, systems, or computer program products. Accordingly, the present invention can take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) having computer-usable program code embodied therein.

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It should be understood that each procedure and/or block in the flowchart and/or block diagram, and a combination of procedures and/or blocks in the flowchart and/or block diagram can be realized by computer program instructions. These computer program instructions may be provided to a general purpose computer, special purpose computer, embedded processor, or processor of other programmable data processing equipment to produce a machine such that the instructions executed by the processor of the computer or other programmable data processing equipment produce a An apparatus for realizing the functions specified in one or more procedures of the flowchart and/or one or more blocks of the block diagram.

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing apparatus to operate in a specific manner, such that the instructions stored in the computer-readable memory produce an article of manufacture comprising instruction means, the instructions The device realizes the function specified in one or more procedures of the flowchart and/or one or more blocks of the block diagram.

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。These computer program instructions can also be loaded onto a computer or other programmable data processing device, causing a series of operational steps to be performed on the computer or other programmable device to produce a computer-implemented process, thereby The instructions provide steps for implementing the functions specified in the flow chart or blocks of the flowchart and/or the block or blocks of the block diagrams.

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。In the present invention, specific examples have been applied to explain the principles and implementation methods of the present invention, and the descriptions of the above examples are only used to help understand the method of the present invention and its core idea; meanwhile, for those of ordinary skill in the art, according to this The idea of the invention will have changes in the specific implementation and scope of application. To sum up, the contents of this specification should not be construed as limiting the present invention.

Claims (10)

  1. A kind of 1. Fetal ECG separation method, it is characterised in that including:
    QRS complex detection is carried out to the maternal abdominal electrocardiosignal received, determines the original position and QRS complex of QRS complex Between original position;
    The signal segment the QRS complex is smoothed using smoothing windows;
    Linear averaging is carried out to the QRS complex, generates average template;
    Parent QRS complex in the average template is related to parent QRS complex corresponding to the current cardiac electrical cycle Property, to the parent QRS complex Mobile state amplitude correction in the average template, generation maternal abdominal electrocardio estimation signal;
    Using maternal abdominal electrocardio estimation signal as reference signal, believe the maternal abdominal electrocardiosignal as expectation Number, by adaptive algorithm, Fetal ECG signal is isolated from the maternal abdominal electrocardiosignal.
  2. 2. Fetal ECG separation method according to claim 1, it is characterised in that using smoothing windows to the QRS complex Between signal segment be smoothed, including:
    Using regular length window, the signal segment the QRS complex is smoothed, the value root of the regular length window Determined according to signal sampling rate and Fetal ECG wave group length.
  3. 3. Fetal ECG separation method according to claim 1, it is characterised in that linearly put down to the QRS complex , average template is generated, including:
    From the maternal abdominal electrocardiosignal, parent electrocardio R ripple crest locations are extracted;
    With the centered about point of the R ripples crest location, linear averaging is carried out to the QRS complex, generates the average template.
  4. 4. Fetal ECG separation method according to claim 3, it is characterised in that according to the parent in the average template The correlation of QRS complex parent QRS complex corresponding with the currently cardiac electrical cycle, to the parent QRS in the average template Wave group Mobile state amplitude correction, generation maternal abdominal electrocardio estimation signal, including:
    It is related to the QRS complex fragment in the average template to calculate QRS complex fragment corresponding to the current cardiac electrical cycle Property, generate coefficient correlation;
    According to the coefficient correlation, QRS complex fragment is to the mother in the average template corresponding to the current cardiac electrical cycle Body QRS complex carries out dynamic amplitudes amendment, generation maternal abdominal electrocardio estimation signal.
  5. 5. Fetal ECG separation method according to claim 1, it is characterised in that by adaptive algorithm, from the mother Fetal ECG signal is isolated in body stomach wall electrocardiosignal, including:
    The adaptive algorithm includes adaptive least mean square algorithm, most adaptive RLS algorithm, a young waiter in a wineshop or an inn Multiplicative lattice shape adaptive algorithm.
  6. A kind of 6. Fetal ECG separator, it is characterised in that including:
    Wave group detection module, for carrying out QRS complex detection to the maternal abdominal electrocardiosignal received, determine QRS complex Original position between original position and QRS complex;
    Smoothing module, for being smoothed using smoothing windows the signal segment the QRS complex;
    Average treatment module, for carrying out linear averaging to the QRS complex, generate average template;
    Wave group correcting module, it is corresponding with the current cardiac electrical cycle for the parent QRS complex in the average template The correlation of parent QRS complex, to the parent QRS complex Mobile state amplitude correction in the average template, generate maternal abdominal Electrocardio estimates signal;
    Separation module, for the maternal abdominal electrocardio to be estimated into as reference signal, the maternal abdominal electrocardio is believed for signal Number desired signal is used as, by adaptive algorithm, Fetal ECG signal is isolated from the maternal abdominal electrocardiosignal.
  7. 7. Fetal ECG separator according to claim 6, it is characterised in that the smoothing module is using smoothly Window is smoothed the signal segment the QRS complex, is specifically included:
    Using regular length window, the signal segment the QRS complex is smoothed, the value root of the regular length window Determined according to signal sampling rate and Fetal ECG wave group length.
  8. 8. Fetal ECG separator according to claim 6, it is characterised in that the average treatment module is to described QRS complex carries out linear averaging, generates average template, including:
    From the maternal abdominal electrocardiosignal, parent electrocardio R ripple crest locations are extracted;
    With the centered about point of the R ripples crest location, linear averaging is carried out to the QRS complex, generates the average template.
  9. 9. Fetal ECG separator according to claim 8, it is characterised in that the wave group correcting module is according to The correlation of parent QRS complex parent QRS complex corresponding with the currently cardiac electrical cycle in average template, to described average Parent QRS complex Mobile state amplitude correction in template, generation maternal abdominal electrocardio estimation signal, including:
    Coefficient correlation computing unit, for calculating QRS complex fragment and the average template corresponding to the current cardiac electrical cycle In QRS complex fragment correlation, generate coefficient correlation;
    Dynamic corrections unit, for according to the coefficient correlation, QRS complex fragment to be to institute corresponding to the current cardiac electrical cycle State the parent QRS complex in average template and carry out dynamic amplitudes amendment, generation maternal abdominal electrocardio estimation signal.
  10. 10. Fetal ECG separator according to claim 6, it is characterised in that the separation module passes through adaptive Algorithm, Fetal ECG signal is isolated from the maternal abdominal electrocardiosignal, including:
    The adaptive algorithm includes adaptive least mean square algorithm, most adaptive RLS algorithm, a young waiter in a wineshop or an inn Multiplicative lattice shape adaptive algorithm.
CN201610639848.8A 2016-08-05 2016-08-05 A kind of Fetal ECG separation method and device Pending CN107684423A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610639848.8A CN107684423A (en) 2016-08-05 2016-08-05 A kind of Fetal ECG separation method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610639848.8A CN107684423A (en) 2016-08-05 2016-08-05 A kind of Fetal ECG separation method and device

Publications (1)

Publication Number Publication Date
CN107684423A true CN107684423A (en) 2018-02-13

Family

ID=61152069

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610639848.8A Pending CN107684423A (en) 2016-08-05 2016-08-05 A kind of Fetal ECG separation method and device

Country Status (1)

Country Link
CN (1) CN107684423A (en)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108836315A (en) * 2018-07-23 2018-11-20 智联时空(北京)科技有限公司 Intelligent waistband for Fetal ECG monitoring
CN112120688A (en) * 2019-06-25 2020-12-25 深圳市理邦精密仪器股份有限公司 Electrocardiosignal processing method, electrocardiosignal processing equipment and computer-readable storage medium
CN114119968A (en) * 2021-12-16 2022-03-01 安徽心之声医疗科技有限公司 Paper electrocardiogram signal region extraction method and system based on deep neural network
CN114767126A (en) * 2022-04-29 2022-07-22 广东粤港澳大湾区国家纳米科技创新研究院 Fetal electrocardiogram monitoring method and system based on marked interference signals

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1287729C (en) * 2001-05-29 2006-12-06 生殖健康技术公司 System for detection and analysis of material uterine, material and fetal cardiac and fetal brain activity
CN101513345A (en) * 2008-12-26 2009-08-26 华南理工大学 Quick blind source separating fetal electrocardioscanner and detection method
CN102204817A (en) * 2011-05-31 2011-10-05 重庆大学 Dual-lead device and method for extracting fetal electrocardiosignals
CN102525452A (en) * 2011-10-19 2012-07-04 广东工业大学 Single-channel fetal electrocardiogram blind separation device based on oblique projection and separation method
CN105305439A (en) * 2015-11-24 2016-02-03 华中科技大学 Probability dynamic power flow computing method and system in view of input variable correlation
CN105342594A (en) * 2015-12-04 2016-02-24 中国人民解放军重庆通信学院 Single-channel maternal abdominal wall fetal heart rate robust estimation method for home monitoring
CN105411577A (en) * 2015-12-30 2016-03-23 深圳先进技术研究院 Method and system for separating fetal ECG (electrocardiogram)
CN105640545A (en) * 2015-12-31 2016-06-08 深圳先进技术研究院 Fetal electrocardiosignal extraction method and device

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1287729C (en) * 2001-05-29 2006-12-06 生殖健康技术公司 System for detection and analysis of material uterine, material and fetal cardiac and fetal brain activity
CN101513345A (en) * 2008-12-26 2009-08-26 华南理工大学 Quick blind source separating fetal electrocardioscanner and detection method
CN102204817A (en) * 2011-05-31 2011-10-05 重庆大学 Dual-lead device and method for extracting fetal electrocardiosignals
CN102525452A (en) * 2011-10-19 2012-07-04 广东工业大学 Single-channel fetal electrocardiogram blind separation device based on oblique projection and separation method
CN105305439A (en) * 2015-11-24 2016-02-03 华中科技大学 Probability dynamic power flow computing method and system in view of input variable correlation
CN105342594A (en) * 2015-12-04 2016-02-24 中国人民解放军重庆通信学院 Single-channel maternal abdominal wall fetal heart rate robust estimation method for home monitoring
CN105411577A (en) * 2015-12-30 2016-03-23 深圳先进技术研究院 Method and system for separating fetal ECG (electrocardiogram)
CN105640545A (en) * 2015-12-31 2016-06-08 深圳先进技术研究院 Fetal electrocardiosignal extraction method and device

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
叶大田: "《基于匹配滤波的胎儿心电监护仪的研制》", 《中国医疗器械杂志》 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108836315A (en) * 2018-07-23 2018-11-20 智联时空(北京)科技有限公司 Intelligent waistband for Fetal ECG monitoring
CN112120688A (en) * 2019-06-25 2020-12-25 深圳市理邦精密仪器股份有限公司 Electrocardiosignal processing method, electrocardiosignal processing equipment and computer-readable storage medium
CN112120688B (en) * 2019-06-25 2024-08-09 深圳市理邦精密仪器股份有限公司 Electrocardiosignal processing method, equipment and computer readable storage medium
CN114119968A (en) * 2021-12-16 2022-03-01 安徽心之声医疗科技有限公司 Paper electrocardiogram signal region extraction method and system based on deep neural network
CN114767126A (en) * 2022-04-29 2022-07-22 广东粤港澳大湾区国家纳米科技创新研究院 Fetal electrocardiogram monitoring method and system based on marked interference signals
CN114767126B (en) * 2022-04-29 2025-04-11 广东粤港澳大湾区国家纳米科技创新研究院 A fetal electrocardiogram monitoring method and system based on marker interference signal

Similar Documents

Publication Publication Date Title
WO2018023697A1 (en) Fetal electrocardiosignal separation method and device
CN101972145B (en) Fetus electrocardio blind separation method based on relative sparsity of time domain of source signal
CN103549950B (en) Improved difference threshold detection algorithm for mobile ECG (electrocardiogram) monitoring
CN103156599B (en) Detection method of electrocardiosignal R characteristic waves
CN104173043B (en) Electrocardio data analysis method suitable for mobile platform
CN104887220B (en) A kind of method and system that Fetal ECG signal is extracted by stomach wall electrocardiosignal
CN103405227B (en) Double-layer morphological filter based electrocardiosignal preprocessing method
CN102626310A (en) Electrocardiogram signal feature detection algorithm based on wavelet transformation lifting and approximate envelope improving
CN106037694A (en) Continuous blood pressure measuring device based on pulse waves
Rahmati et al. A PCA/ICA based fetal ECG extraction from mother abdominal recordings by means of a novel data-driven approach to fetal ECG quality assessment
CN105640545A (en) Fetal electrocardiosignal extraction method and device
CN104305992B (en) A kind of interactive fast automatic extracting method of Fetal ECG
CN105411577A (en) Method and system for separating fetal ECG (electrocardiogram)
CN108272451A (en) A kind of QRS wave recognition methods based on improvement wavelet transformation
CN107684423A (en) A kind of Fetal ECG separation method and device
CN103054572B (en) Blind source separation method on basis of late potential signals of single-lead cardiac ventricle
Franchevska et al. The method and algorithm for detecting the fetal ECG signal in the presence of interference
CN102783945A (en) Fetal electrocardiogram signal extracting method based on wavelet threshold denoising
Di Maria et al. Extracting fetal heart beats from maternal abdominal recordings: selection of the optimal principal components
CN103750835A (en) Electrocardiosignal characteristic detection algorithm
Kahankova Influence of gestation age on the performance of adaptive systems for fetal ECG extraction
CN106108850A (en) The recognition methods of the interference data of ecg database and device
WO2018023698A1 (en) Fetal-electrocardiogram separation method and device
CN107684422A (en) A kind of Fetal ECG separation method and device
CN110179456A (en) Electrocardio Noise Identification model training and electrocardio noise detecting method, device

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180213