一种胎儿心电分离方法及装置与流程

本发明涉及生物医学工程技术领域，尤其涉及一种胎儿心电分离方法及装置，具体的讲是一种基于两级模板修正的胎儿心电分离方法及装置。

背景技术：

本部分旨在为权利要求书中陈述的本发明的实施方式提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。

胎儿心电信号(fetalelectrocardiogram，fecg)是胎儿宫内生理活动最主要的客观指标之一，通过分析胎儿心电信号，能够早期诊断胎儿宫内缺氧、宫内窘迫以及先天性心脏病等，有效地降低新生儿各类疾病的发病率和死亡率。因此，围产期监护仪的临床应用将是保证围产期各个阶段的胎儿安全，保障胎儿优生的有力措施之一。

在围产期，胎儿在母体腹部由几个不同的电导率的包裹层包围，如胎盘、羊水、母体子宫和母体腹部皮下组织等，其中羊水的导电能力最高，而胎脂的导电能力最低，此外孕妇的皮肤和皮下脂肪也有微弱的导电性，但是导电能力比肌肉组织约小十倍。在采集胎儿心电时，由于胎儿心电信号需要经过上述导电层，本来就很微弱的胎儿心电信号更容易受外界的干扰。这些不同的组织，形成了所谓的“人体容积导体”，这个容积导体不是一个稳定的导体，其导电性和几何形状随着整个妊娠期不断改变。因此，这将导致母体腹部心电混合信号具有如下特点：1)信噪比低，胎儿心电分量太弱，仅为10～50uv左右，而母亲心电信号则可达750uv，两者可差十倍以上。肌电噪声与其它干扰也很强，信号完全淹没于噪声之中。2)时域中胎儿qrs波与母体qrs波有10～30％的重叠。3)频域中胎儿频谱与母体的频谱大部分重叠。4)信号随机性很强，是非平稳过程，特别是由于胎儿的自由移动造成胎儿心电传输通道的时变性。

母体腹部心电混合信号的上述特点，为获得高精度的胎儿心电信号制造了比较大的障碍，因此从母体腹壁混合信号中提取出胎儿心电信号，以实现围产期胎儿生命信息的监护，是目前广泛研究和关注的问题。

目前，从母体腹壁混合信号中提取出胎儿心电信号的方法主要有以下几种：

1)匹配滤波：匹配滤波主要利用的是信号的自相关和互相关等统计特性，但是在实际环境中所要面临的情况远比假设情况要复杂的多，采用匹配滤波获得的胎儿心电信号不是很理想。

2)独立成分分析法(ica)：此方法需要多路腹壁信号，无法实现单导联胎儿心电信号的提取。

3)小波阈值去噪的胎儿心电信号提取方法滤波：利用小波变换结合先验信号知识能够获得效果不错的胎儿心电信号，但是提取过程中容易受到幅度和胎儿心电信号类似信号的干扰，往往导致胎儿心电信号中r波漏检和误检等情况产生，在计算胎儿心率时会引起误差，不利于于临床诊断。

4)神经网络法：此方法可以根据神经网络的预测误差检测心电信号的qrs复合波。pt段信号平缓，可以被较好的预测，而对于qrs复合波，相邻点差异较大，会产生预测误差，利用这个预测误差训练网络，对qrs复合波进行预测，预测时的人工神经网络具有多层结构，利用神经网络学习训练过程需要花费较长时间，而且要求训练样本具有一定的代表性，这在实际腹壁心电信号降噪应用中很难实现。

5)自适应滤波法：自适应滤波方法包括使用一个或几个参考通道训练自适应滤波消除母体心电图。现有的自适应滤波方法去除母体ecg成分，需要形态上类似于被淹没fecg信号的母体ecg参考通道，但是现有的自适应对消方法要求额外提供参考信号，也就是说，除了腹壁电极外，还需要配备额外的母体胸导连电极，自适应方法对消掉胸导联和腹壁导联中相关的母体分量，这意味着母体腹壁心电信号和母体胸导联心电信号需要相关，才能够被较好的对消掉，并且，母体心电信号的传导路径非常复杂，某些情况下，腹壁中的母体心电信号和胸导联中的母体心电信号相关性差。因此，利用现有的自适应滤波算法对消提取胎儿心电信号会有大量的母体心电泄露，对正确识别胎儿心电分量造成了很大干扰。并且，现有的自适应滤波对消无法实现单导联的胎儿心电信号提取。

技术实现要素：

本发明提出一种胎儿心电分离方法及装置，用以解决现有的胎儿心电分离算法计算复杂且无法实现单导联的问题。

为了达到上述目的，本发明实施例提供一种胎儿心电分离方法，包括：

以接收到的母体腹壁心电信号中的r波波峰位置为中心定位点，截取特定长度的心电信号片段；对当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段进行线性平均，生成平均模板；根据所述平均模板中的母体qrs波群与所述当前心电周期对应的母体qrs波群的相关性，对所述平均模板中的母体qrs波群行动态幅值修正，生成母体腹壁心电估计信号；根据所述母体腹壁心电估计信号，从所述母体腹壁心电信号中分离出胎儿心电信号qrs。

为了达到上述目的，本发明实施例还提供一种胎儿心电分离装置，包括：片段截取模块，用于以接收到的母体腹壁心电信号中的r波波峰位置为中心定位点，截取特定长度的心电信号片段；平均模板生成模块，用于对当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段进行线性平均，生成平均模板；波群修正模块，用于根据所述平均模板中的母体qrs波群与所述当前心电周期对应的母体qrs波群的相关性，对所述平均模板中的母体qrs波群进行动态幅值修正，生成母体腹壁心电估计信号；分离模块，用于根据所述母体腹壁心电估计信号，从所述母体腹壁心电信号中分离出胎儿心电信号qrs。

本发明提出的胎儿心电分离方法及装置，实行两级模板修正，对第一级估计后的母体腹壁心电信号进行了母体qrs波群修正，修正后的母体腹壁信号估计量，将腹壁混合信号和估计的参考信号进行匹配对消，因此本发明的实现只需要一路母体腹壁心电信号，故可以完成单导联胎儿心电信号分离；并且，本发明的胎儿分离算法设计思路简单直观，针对于异常的临床数据，易于查找原因和添加修正模块。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例的胎儿心电分离方法的系统结构示意图；

图2为本发明实施例的基于两级模板修正的腹壁母体心电估计机制示意图；

图3为本发明实施例的胎儿心电分离方法的处理流程图；

图4为本发明实施例的第一级腹壁母体心电估计机制的示意框图；

图5为本发明实施例的第二级腹壁母体心电估计机制示意图；

图6为本发明实施例的胎儿心电分离装置的结构示意图；

图7为图6所示实施例的波群修正模块103的结构示意图；

图8为本发明的具体实施例采用的5路母体腹壁混合心电信号的波形图；

图9为本发明的具体实施例的实验结果分析图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本领域技术技术人员知道，本发明的实施方式可以实现为一种系统、装置、设备、方法或计算机程序产品。因此，本公开可以具体实现为以下形式，即：完全的硬件、完全的软件(包括固件、驻留软件、微代码等)，或者硬件和软件结合的形式。

下面参考本发明的若干代表性实施方式，详细阐释本发明的原理和精神。

本发明针对从母体腹壁提取出的混合心电信号，利用平均法对腹壁母体心电信号进行估计，并且针对于母体qrs波群泄露会对正确识别胎儿心电分量造成很大影响的缺点，本发明对估计后的母体腹壁心电信号进行了母体qrs波群修正，生成修正后的母体腹壁信号估计量，将腹壁混合信号和母体腹壁信号估计量进行匹配对消，从而分离出胎儿心电信号。并且，本发明实现只需要一路母体腹壁心电信号，因此可以实现单导联胎儿心电信号分离。

图1为本发明实施例的胎儿心电分离方法的系统结构示意图。如图1所示，y(n)表示采集到的母体腹壁心电信号，f(n)为腹壁信号中的胎儿心电信号，x(n)为腹壁信号中的母体心电信号，为估计出的腹壁母体心电信号，即为分离出的胎儿心电信号。

本发明的核心是如何计算出腹壁母体心电估计信号，以实现母体分量的有效对消。在本发明中，采用两级模板修正的方法，如图2所示，腹壁心电信号经过第一级模板，生成腹壁母体心电信号的第一级估计，但是在腹壁胎儿心电分离中，由于母体r波波峰值是胎儿心电峰值的十几倍，qrs波群估计不准确，会造成大量的母体qrs分量泄露，泄露的分量对准确识别胎儿心电节拍位置造成很大的影响，针对于此问题，本发明中还设计了第二级模板模块，用于修正第一级母体心电估计信号中的qrs波群，以实现母体心电分量的有效对消。

图3为本发明实施例的胎儿心电分离方法的处理流程图。如图3所示，包括：

步骤s101，以接收到的母体腹壁心电信号中的r波波峰位置为中心定位点，截取特定长度的心电信号片段；

步骤s102，对当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段进行线性平均，生成平均模板；

步骤s103，根据所述平均模板中的母体qrs波群与所述当前心电周期对应的母体qrs波群的相关性，对所述平均模板中的母体qrs波群进行动态幅值修正，生成母体腹壁心电估计信号；

步骤s104，根据所述母体腹壁心电估计信号，从所述母体腹壁心电信号中分离出胎儿心电信号。

其中，步骤s101和s102对应于第一级模板的处理，生成腹壁母体心电信号的第一级估计，具体实施时，参看图4，为本发明实施例的第一句腹壁母体心电估计机制示意图。

在本发明实施例的步骤s101中，首先要以腹壁心电信号中r波波峰位置为中心定位点，截取特定长度的心电信号。

其中，确定母体腹壁心电信号中的r波波峰位置的方法为：对接收到的母体腹壁心电信号进行差分运算，以获取腹壁心电信号的斜率曲线，对斜率曲线做绝对滑动积分运算，记录每秒钟内的最大值，对前8秒内的最大值做平均，则r波判定阈值＝系数比值*前8秒内的最大值平均，对差分值大于r波判定阈值的段，进行r波波峰位置最大值查找，以确定最终的r波波峰位置点。

并且，在步骤s101中，截取特定长度的心电信号片段，n为截取的片段的长度，且为r波波峰间距离，γ的取值范围为(1.1，1.4)。截取片段的长度n的选择至关重要，γ的取值范围(1.1，1.4)是为了保证片段截取的起止点都落在心电信号的tp段，此段内信号平缓，并且片段重合部分将会被平均，避免了平均模板连接处信号较大的幅值落差。参看图4的(b)部分，在该实施例中，根据特定的片段长度分割了5个心电信号片段。

在本发明实施例的步骤s102中，对当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段进行线性平均，生成平均模板。参看图4的(c)部分，在该实施例中，五角星对应的为当前心电周期，取其前面的4个心电信号片段与五角星对应的心电周期的心电信号片段进行线性平均。

在步骤s102中，具体实现的方法包括：

将所述当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段的r波波峰位置对齐，进行线性平均，生成平均模板：

其中，y1(t)表示生成的平均模板，表示所述平均模板中的母体心电分量，表示所述平均模板中的胎儿心电分量。

由于分割片段的中心点是母体r波波峰，胎儿心电和母体心电来自于两个独立的信号源，即胎儿qrs波群分布在母体心电周期内的不同位置。

因此，平均模板中的胎儿心电分量可表示为：

由于母体心电幅值一般是胎儿心电幅值的十几倍，即：

因此：

也就是说，y1(t)是对腹壁母体心电信号的第一级的平均估计，估计出的信号中几乎不含有胎儿心电信号，如图4所示实施例的(d)部分。

在本发明实施例中，步骤s103对应于第二级模板的处理，用以对所述平均模板中的母体qrs波群进行动态幅值修正，生成最终的母体腹壁心电估计信号。

在本发明实施例中，能否有效的抑制掉母体qrs波群分量，对正确识别分离结果中的胎儿心电分量起着决定性作用。平均模板中的母体心电波形，是对较早的心电波形和最新波形的等权累加和线性平均，估计出的母体心电片段不具有动态性。因为，本发明的第二模板根据最新母体qrs波群和平均片段中母体qrs波群的相关性，对平均片段中的母体qrs波群进行动态修正。

具体实施时，参看图5，方法如下：

首先，计算所述当前心电周期对应的qrs波群片段与所述平均模板中的qrs波群片段的相关性，生成相关系数，具体步骤为：

(1)参看图5的(b)部分，通过母体心电r波检测提取当前心电周期中的母体qrs波群片段：y3＝[y3(1),y3(2),...,y3(p)]，其中p为母体qrs波群片段的长度；

(2)提取出平均模板中母体qrs波群片段：y2＝[y2(1),y2(2),...,y2(p)]；

(3)计算当前心电周期中的母体qrs波群片段与平均模板中母体qrs波群片段的相关系数r23：

其次，参看图5的(c)部分，根据所述相关系数，用所述当前心电周期对应的qrs波群片段对所述平均模板中的母体qrs波群进行动态幅值修正，所采用的表达式如下：

最后，生成经过第二级模板修正的母体腹壁心电估计信号，如图5的(d)部分所示。由于对第一级模板生成的母体腹壁心电估计信号进行了qrs波群修正，因此估计的更精确，可以有效对消掉母体心电分量，具有较好的胎儿心电分离效果。

在步骤s104中，根据所述母体腹壁心电估计信号，对消掉所述母体腹壁心电信号中的母体心电分量，以获得所述胎儿心电分量。

应当注意，尽管在附图中以特定顺序描述了本发明方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。

在介绍了本发明示例性实施方式的方法之后，接下来，参考图6对本发明示例性实施方式的胎儿心电分离装置进行介绍。该装置的实施可以参见上述方法的实施，重复之处不再赘述。以下所使用的术语“模块”和“单元”，可以是实现预定功能的软件和/或硬件。尽管以下实施例所描述的模块较佳地以软件来实现，但是硬件，或者软件和硬件的组合的实现也是可能并被构想的。

图6为本发明实施例的胎儿心电分离装置的结构示意图。如图6所示，包括：

片段截取模块101，用于以接收到的母体腹壁心电信号中的r波波峰位置为中心定位点，截取特定长度的心电信号片段；平均模板生成模块102，用于对当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段进行线性平均，生成平均模板；波群修正模块103，用于根据所述平均模板中的母体qrs波群与所述当前心电周期对应的母体qrs波群的相关性，对所述平均模板中的母体qrs波群进行动态幅值修正，生成母体腹壁心电估计信号；分离模块104，用于根据所述母体腹壁心电估计信号，从所述母体腹壁心电信号中分离出胎儿心电信号。

进一步地，在本发明实施例中，所述片段截取模块101以接收到的母体腹壁心电信号中的r波波峰位置为中心定位点，截取特定长度的心电信号片段，具体包括：

所述心电信号片段的特定长度为n，且其中γ的取值范围(1.1，1.4)。

进一步地，在本发明实施例中，所述平均模板生成模块102对当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段进行线性平均，生成平均模板，具体包括：

将所述当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段的r波波峰位置对齐，进行线性平均，生成平均模板：

其中，y1(t)表示生成的平均模板，表示所述平均模板中的母体心电分量，表示所述平均模板中的胎儿心电分量。

进一步地，在本发明实施例中，所述波群修正模块103根据所述平均模板中的母体qrs波群与所述当前心电周期对应的母体qrs波群的相关性，对所述平均模板中的母体qrs波群进行动态幅值修正，生成母体腹壁心电估计信号，如图7所示，具体包括：

相关系数计算单元1031，用于计算所述当前心电周期对应的qrs波群片段与所述平均模板中的qrs波群片段的相关性，生成相关系数；

动态修正单元1032，用于根据所述相关系数，用所述当前心电周期对应的qrs波群片段对所述平均模板中的母体qrs波群进行动态幅值修正，生成母体腹壁心电估计信号。

进一步地，在本发明实施例中，所述分离模块104根据所述母体腹壁心电估计信号，对消掉所述母体腹壁心电信号中的母体心电分量，以获得所述胎儿心电分量。

此外，尽管在上文详细描述中提及了胎儿心电分离装置的若干单元，但是这种划分仅仅并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。同样，上文描述的一个单元的特征和功能也可以进一步划分为由多个单元来具体化。

以下通过一个具体实施例来说明采用本发明的基于两级模板的胎儿心电分离方法实现胎儿心电信号分离的效果。

参看图8所示，为该具体实施例采用的5路母体腹壁混合心电信号的波形图。如图8所示，每路信号的数据长度为2500个样本点，采样率250hz。由图8可知，第1路腹壁信号为胎母心电信噪比最大的一路信号，即此路信号效果最好，第4路和第5路腹壁心电信号用肉眼已经无法识别出胎儿心电r波峰位置，此两路信号提取效果已不能为本发明的算法所用，故不在本发明的算法的考虑范围之内。

本实施例选择第2路信号作为算法验证数据，其实验结果分析参看图9。

1、胎儿qrs波群远离母体qrs波群：

胎儿qrs波群远离母体qrs波群，如图9中标注a，由图9可知，这种情况下的胎儿qrs波分量被提取出，左右的母体qrs波分量被抑制掉。此种胎儿qrs波群和母体qrs波群的位置关系，是腹壁心电中最普遍的位置关系，也是最容易提取出的位置关系。

2、胎儿qrs波群紧邻母体qrs波群：

胎儿qrs波群紧邻母体qrs波群，如图9中标注b所示，由图9可知，这种情况下的胎儿qrs波分量被提取出，紧邻的母体qrs波分量被抑制掉，并未对胎儿心电qrs波群的识别造成影响。此种胎儿qrs波群和母体qrs波群的位置关系，是腹壁心电中少数比例存在的位置关系，也是较难提取出的位置关系。

3、胎儿qrs波群和母体qrs波群完全重合：

胎儿qrs波群完全和母体qrs波群重合，如图9中标注c所示，由图9可知，这种情况下的胎儿qrs波分量被提取出，重合的母体qrs波分量被抑制掉。此种胎儿qrs波群和母体qrs波群的位置关系，是腹壁心电中少数比例存在的位置关系，也是最难提取出的位置关系。

由以上结果分析可知，通过本发明的基于两级模板的胎儿心电分离方法，胎儿qrs波群远离母体qrs波群、胎儿qrs波群紧邻母体qrs波群、胎儿qrs波群和母体qrs波群完全重合三种情况，均有非常理想的提取效果。

本发明提出的胎儿心电分离方法及装置，实行两级模板修正，对第一级估计后的母体腹壁心电信号进行了母体qrs波群修正，修正后的母体腹壁信号估计量，将腹壁混合信号和估计的参考信号进行匹配对消，因此本发明的实现只需要一路母体腹壁心电信号，故可以完成单导联胎儿心电信号分离；并且，本发明的胎儿分离算法设计思路简单直观，针对于异常的临床数据，易于查找原因和添加修正模块。

本发明实施例还提供一种计算机可读程序，其中当执行所述程序时，所述程序使得计算机执行本发明实施例所述的胎儿心电分离方法。

本发明实施例还提供一种存储有计算机可读程序的存储介质，其中所述计算机可读程序使得计算机执行本发明实施例所述的胎儿心电分离方法。

本发明以上的装置和方法可以由硬件实现，也可以由硬件结合软件实现。本发明涉及这样的计算机可读程序，当该程序被逻辑部件所执行时，能够使该逻辑部件实现上文所述的装置或构成部件，或使该逻辑部件实现上文所述的各种方法或步骤。本发明还涉及用于存储以上程序的存储介质，如硬盘、磁盘、光盘、dvd、flash存储器等。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。