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

本发明涉及生物医学工程技术领域，尤其涉及一种胎儿心电分离方法及装置，具体的讲是一种自适应单导联胎儿心电分离方法及装置。

背景技术：

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

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

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

上述的这些特点，为获得高精度的胎儿心电信号制造了比较大的障碍，因此从母体腹壁信号中提取出胎儿心电信号，以实现围产期胎儿生命信息的监护，是目前广泛研究和关注的问题。目前，对此问题的解决主要采用以下几种方法：

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

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

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

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

针对上述问题，大量研究中将自适应方法用于胎儿心电分离，自适应滤波主要通过自主学习达到适应条件变化的效果。但是传统的自适应对消方法要求额外提供参考信号，即，除了腹壁电极外，还需要配备额外的母体胸导连电极，自适应方法对消掉胸导联和腹壁导联中相关的母体分量，这意味着母体腹壁心电信号和母体胸导联心电信号需要相关，才能够被较好的对消掉，并且，母体心电信号的传导路径非常复杂，某些情况下，腹壁中的母体心电信号和胸导联中的母体心电信号相关性差。因此，利用自适应对消提取胎儿心电信号会有大量的母体心电泄露，对正确识别胎儿心电分量造成了很大干扰。此外，自适应对消仍然无法实现单导联的胎儿心电信号提取。

技术实现要素：

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

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

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

本发明提出的胎儿心电分离方法及装置，只需要一路腹壁心电信号，实现了单导联；并且，经过多次叠加、加权后的递归模板基本不再含有胎儿心电分量，自适应滤波算法消除了参考信号与主输入信号中相关的部分，保留不相关的部分，因此腹壁心电信号中的胎儿心电信号分离可以被最大化的提取出来；本发明的胎儿分离算法设计思路简单直观，针对于异常的临床数据，易于查找原因和添加修正模块。

附图说明

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

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

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

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

图4为图3所示实施例的步骤s102的处理流程图；

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

图6为图5所示实施例的递归模板生成模块102的结构示意图；

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

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

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

具体实施方式

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

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

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

本发明针对从母体腹壁提取出的混合心电信号，利用平均法对腹壁母体心电信号进行估计，针对于母体qrs波群泄露，对正确识别胎儿心电分量造成很大影响，本发明对估计后的母体腹壁心电信号进行了母体qrs波群修正，修正后的母体腹壁信号估计量，作为自适应算法中的参考信号，由于该方法的实现只需要一路母体腹壁心电信号，因此，本发明可以使基于自适应的胎儿心电信号分离方案在单导联系统下实现。

图1为本发明实施例的胎儿心电分离方法的系统结构示意图。如图1所示，y(n)表示采集到的母体腹壁心电信号，f(n)为腹壁信号中的胎儿心电信号，x(n)为腹壁信号中的母体心电信号，y(n)为自适应滤波器的期望信号，x'(n)为通过腹壁母体心电估计机制估计的腹壁母体心电，自适应滤波器对输入信号不断的进行“学习”和“跟踪”，使输出信号最大限度的逼近腹壁心电信号，从而使误差信号最小，即为分离出的胎儿心电信号。

也就是说，基于自适应的胎儿心电分离算法，需要两路输入信号，即：一路母体腹壁心电信号(主输入信号)，一路估计出的腹壁母体心电信号(参考输入信号)，这相对于现有技术中的需要接入一路胸导联心电信号作为参考输入信号，实现了单导联。

本发明的核心是如何计算出腹壁母体心电估计信号，以作为自适应算法的参考输入信号。在本发明中，如图2所示，为本发明的腹壁母体心电估计机制的示意框图。母体腹壁信号经过r波波峰点检测，确定母体腹壁心电信号中的r波波峰位置，并根据r波波峰位置采用特定长度对母体腹壁信号进行分割，分割后的信号片段进行线性加权平均后，生成递归模板，递归模板基本不含有胎儿心电分量；然后提取当前心电周期的qrs波群对递归模板中的qrs波群做修正，生成母体腹壁心电估计信号。

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

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

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

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

步骤s104，将所述母体腹壁心电估计信号作为参考信号，将所述母体腹壁心电信号作为期望信号，通过自适应算法，从所述母体腹壁心电信号中分离出胎儿心电信号。

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

并且，在步骤s101中，截取特定长度的心电信号片段，n为截取的片段的长度，且的取值范围为(1.1，1.4)。以母体r波波峰位置为定位点，对r波波峰左右各取固定长度，形成r-r～2r-r长度的混合心电模板，对当r-r间隔前的固定数目的模板求平均模板。其中，截取片段长度n的选择至关重要，的取值范围(1.1，1.4)是为了保证片段截取的起止点都落在心电信号的tp段，此段内信号平缓，并且片段重合部分将会被平均，避免了递归模板连接处信号较大的幅值落差。

在本发明实施例的步骤s102中，对当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段进行线性加权平均，生成递归模板，如图4所示，具体包括：

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

对当前心电周期前面的m-1个连续的心电信号片段进行线性平均，生成前m-1个心电信号片段的平均模板，具体方法为：

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

其中，y1(t)表示生成的前m-1个心电信号片段的平均模板，表示所述平均模板中的母体心电分量，表示所述平均模板中的胎儿心电分量。

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

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

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

因此：

也就是说，y1(t)是前m-1个心电信号片段的平均模板，估计出的信号中几乎不含有胎儿心电信号

步骤s1022，根据所述当前心电周期对应的心电信号片段与所述平均模板获得递归模板＝a*平均模板+b*当前心电周期对应的心电信号片段，a，b为设定的加权系数。加权系数a，b可以由人为经验设定。

经过多次叠加、加权后的递归模板基本不再含有胎儿心电分量。

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

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

具体实施时，方法如下：

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

(1)提取当前心电周期中的母体qrs波群片段：y3＝[y3(1),y3(2),...,y3(p)]，其中p为母体qrs波群片段的长度；

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

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

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

所采用的表达式如下：

由于对递归模板中的母体腹壁心电估计信号进行了qrs波群修正，因此估计的更精确，具有较好的胎儿心电分离效果。

在步骤s104中，腹部电极取得主输入信号，将估计出的母体腹壁心电信号作为输入信号，此信号中基本不含有胎儿心电信号，只包含欲消除的噪声和母体心电信号；腹壁心电信号为期望信号，通过自适应算法，消除了输入信号与期望信号中相关的部分，保留不相关的部分，因此腹壁心电信号中的胎儿心电信号分离可以被最大化的提取出来。

本发明实施例中，自适应算法可以为自适应的最小均方(lms)算法、自适应的递归最小二乘方(rls)算法、最小二乘格形(lsl)自适应算法等。

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

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

图5为本发明实施例的胎儿心电分离装置的结构示意图。如图5所示，包括：片段截取模块101，用于以接收到的母体腹壁心电信号中的r波波峰位置为中心定位点，截取特定长度的心电信号片段；递归模板生成模块102，用于对当前心电周期对应的心电信号片段以及前面m-1个连续的心电信号片段进行线性加权平均，生成递归模板；波群修正模块103，用于根据所述递归模板中的母体rqs波群与所述当前心电周期对应的母体qrs波群的相关性，对所述递归模板中的母体qrs波群进行动态幅值修正，生成母体腹壁心电估计信号；分离模块104，用于将所述母体腹壁心电估计信号作为参考信号，将所述母体腹壁心电信号作为期望信号，通过自适应算法，从所述母体腹壁心电信号中分离出胎儿心电信号。

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

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

平均模板生成单元1021，对当前心电周期前面的m-1个连续的心电信号片段进行线性平均，生成前m-1个心电信号片段的平均模板；

加权单元1022，用于根据所述当前心电周期对应的心电信号片段与所述平均模板获得递归模板＝a*平均模板+b*当前心电周期对应的心电信号片段，a，b为设定的加权系数。

进一步地，在本发明实施例中，所述波群修正模块103根据所述递归模板中的母体rqs波群与所述当前心电周期对应的母体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波群完全重合三种情况，均有非常理想的提取效果。

本发明提出的胎儿心电分离方法及装置，具有以下有益效果：

1、单导联：本发明仅需要一路母体腹壁心电信号，即可完成基于自适应的胎儿信号分离。传统的采用基于自适应的胎儿心电信号分离算法，需要两路信号，即：一路腹壁心电信号(主输入信号)，一路胸导联心电信号(参考输入信号)，通过对主输入信号进行模板递归生成参考信号，因此，本发明设计中只需要一路腹壁心电信号，这将大大简化了系统的导联系统。

2、分离效果好：由于母体心电分量幅值是胎儿心电分量幅值的十几倍，因此经过多次叠加、加权后的递归模板基本不再含有胎儿心电分量；自适应滤波算法消除了参考信号与主输入信号中相关的部分，保留不相关的部分，因此腹壁心电信号中的胎儿心电信号分离可以被最大化的提取出来。

3、算法设计思路简单直观，针对于异常的临床数据，易于查找原因和添加修正模块。

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

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

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

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

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

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

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

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