CN102178522A - 一种母亲及胎儿心电信号qrs波中r波的检测定位方法 - Google Patents

Abstract

本发明公开了母亲及胎儿心电信号QRS波中R波的检测定位方法。前端滤波：这里的预处理过程依次采用了低通、高通以及微分器滤波器；小波处理过程：根据母亲心电信号和胎儿心电信号的各自医学特征，分布的频段、幅度大小，采用小波分解过程，使母亲心电信号的R波和胎儿心电信号的R波相对更加明显；自适应阈值R波判决，Rdetect算法中采用浮动的自适应阈值，首先判决母亲心电信号的R波出现位置，接着经过前端处理，使用该算法对胎儿心电信号中QRS中R波特征点出现的位置进行定位，经后处理及信号调整，完成在混合信号上检测到母亲心电及胎儿心电的R波位置。本发明对于采集到的混合心电信号，能很好的达到分别检测的目的。

Description

一种母亲及胎儿心电信号QRS波中R波的检测定位方法

技术领域

本发明属于生物医学工程技术领域，尤其涉及一种母亲及胎儿心电信号QRS波中R波的检测定位方法。

背景技术

心电图检测是20世纪建立起来的并广泛应用于临床诊断和检测的重大技术才成果之一。随着科学的发展和技术的不断进步，尤其计算机技术的日益完善，心电图检测技术的发展步入了更高的层次。QRS复合波(QRS Complex)(参考附图1)检测是心电图(Electrocardiogram，ECG)信号分析中的一个关键问题，只有在QRS波特征(峰值点和起始点)确定以后才能计算基线，检测P波、T波及ST段等参数，是心电自动分析的基础。但由于心电信号在生理意义上的多样性及人与人之间的差异性，再加上存在实时心电信号中的漂移、干扰、噪声，所以QRS波检测也一直是心电信号实时自动分析的难点。

一般来说，与心电信号中其他波形相比，R波具有较高的幅值；另外，由对心电信号的频谱分析可知，QRS波的中心频带在17Hz左右(该频率也被称为QRS波的特征频率)，带宽约为10Hz，而T波、P波、基线漂移等的频带均在此频带的低端外。以上两点是QRS波区别与其他波形的显著特点。除此之外，QRS波的面积、持续期、R波的斜率也是QRS波的重要特点。基于QRS的这些特点，目前出现了很多新的算法。例如基于人工神经网络(Artificial Neural Network，ANN)的算法、遗传算法、滤波器组、希尔伯特变换、匹配滤波、自适应滤波器、小波变换以及奇异值检测等算法。这些算法对于QRS波的提取，都能有较好的效果，也被引申出很多其他的优良算法。

然而，这些算法所应用的前提是，检测单一心电信号的R波位置，或者说检测混合心电信号中幅度大的心电信号分量的R波。而对于其他微弱的、幅度很小的心电分量的R波，则无能为力。当然，在很多情况下，我们都只需要检测出主要心电信号的R波位置，至于其他微弱的分量，都不是很重要。

但是，由于现在社会对母婴健康和安全的要求日益提高，胎儿监护成为了一项很具重要意义的研究课题。胎儿心电图提供了有关胎儿健康的信息，是反映胎儿在孕期的生长和健康状况的一项重要指标。通过对围产期的胎儿心电信号提取和分析，可以确定胎儿心率，分析胎儿心脏功能参数，判断胎儿是否患有窘迫、心律紊乱、酸中毒等新生儿疾病，以便及早诊断宫内缺氧及其它宫内异常，及时采取适当治疗措施，保证孕妇和胎儿的安全。这对提高围产儿监护质量具有非常重要的意义，从深层次，也说明了产前对胎儿进行监护是十分必要的。

目前，获取胎儿心电的方式，有无侵入式的，还有侵入式这两种。无侵入式方式，是通过在孕妇体表放置电极，一般是孕妇的胸部、腹部等位置，通过电极采集得到混合的心电信号；混合的心电信号里头，主要有母亲的心电信号，胎儿的心电信号，以及其它各类干扰噪声，这里母亲心电信号的幅度一般是胎儿心电信号幅度的上百倍。整个过程，我们的初衷是得到胎儿心电信号的各类信息，以作为判断胎儿健康状况的参考。而获取胎儿的QRS波也显得很是必需。如之前所言，母亲心电信号频谱与胎儿心电信号频谱存在重叠，而且母亲心电信号的幅度远大于胎儿心电信号幅度，如果采用之前提出的这些算法，所能检测的是母亲的QRS波，而对于胎儿的QRS波则检测不出来。之前提到的那些算法在这种情况下，是有很大局限性的。

发明内容

本发明的目的在于克服上述现有技术的缺点和不足，提供一种母亲及胎儿心电信号QRS波中R波的检测定位方法，本发明针对无侵入采集得到的混合心电信号，能同时定位检测出母亲和胎儿心电信号中的R波位置，为后续胎儿心电信号各类信息的获取提供支持，以为更好的进行胎儿监护。

本发明通过下述技术方案实现：

一种母亲及胎儿心电信号QRS波中R波的检测定位方法，包括如下步骤：

(1)前端滤波处理，包括低通滤波、高通滤波、微分器以及平方操作；

(2)小波分解，对于经过步骤(1)中前端滤波处理的混合心电信号，根据信号的特征，需要经过小波分解过程，以更加突出胎儿心电信号R波的份量；

(3)对步骤(2)处理后的信号采用Rdetect算法检测出混合心电信号中母亲心电信号QRS波中R波特征点出现的位置；

(4)对经过上述处理的混合心电信号，进行前端处理，并保留信号的辅助信息；

(5)采用Rdetect算法检测出经过前端处理后的信号中，胎儿心电信号QRS波中R波特征点出现的位置；

(6)后处理，利用前端处理中的辅助信息，恢复信号；

(7)处理过程中，会引入延迟；根据这些延迟调整信号。

上述步骤(2)中，采用小波一层分解，取细节部分，便于后续胎儿心电信号QRS波中R波位置的检测。

上述步骤(3)中，经过滤波以及小波分解后，采用Rdetect算法，其中Rdetect算法采用自适应阈值的方式，定位检测到母亲心电的R波位置，Rdetect算法参数的选择如下：

初始阈值：threshold＝param×max{y(n)}，这里param取0.7，也可以取0.7相近的一些数值；

自适应阈值：threshold＝param1×threshold+param1×maxy，其中，maxy是在检测R波中，出现的幅值最大值，param1取值为0.3～0.5。

上述步骤(4)中，所述前端处理中，采用的方法是去除母亲R波位置出前后一定范围的数据点，并把经过这一操作的信号进行重新组合得到新的信号，在进行前端处理的过程中，记录信号的辅助信息，即母亲心电信号R波出现的位置点，该R波的数目以及每个R波出现位置处去除的点数；根据心电信号准周期特性，每个R波位置处去除的点数近似相等，记为POINTS，POINTS的经验值是：

POINTS＝round(Fs/10)

其中，Fs是采样频率，round(x)表示对数值x取四舍五入的整数值。

上述步骤(5)采用Rdetect算法，检测出胎儿心电信号的R波位置，这里的输入信号是步骤(4)中重新组合的信号，所以经过Rdetect算法处理后，步骤(6)后处理中需要根据步骤(4)中母亲R波点的位置，加上减去的信号采样点数。

上述步骤(7)中，需要考虑之前处理过程中引入的延迟，计算滤波处理的延迟周期，即对应延迟的采样点数，后续的处理步骤中，是把信号从时域变换到小波域。

与现有技术相比本发明的优点和效果在于：

(1)本方法利用了心电信号的准周期特性，对于从孕妇腹部采集到的混合心电信号，相比于其他方法，除了能够准确地定位检测出母体心电信号R波位置，还能检测出其中幅度相对很微弱的胎儿心电信号的R波位置；

(2)本方法采用了浮动的阈值，阈值选取方式上实现简单；在检测过程中，对于输入的混合心电信号，能立即定位到R波点的位置，不存在学习调整过程；

(3)本方法对于孕妇腹部采集到的混合心电信号，定位检测母亲与胎儿心电中的R波位置准确率很高；

(4)本方法中的Rdetect算法实现简单，运行速度很快，便于软件和硬件实现；

附图说明

图1是现有心电信号中的QRS波形；

图2是现有通用的QRS波检测的算法结构；

图3是本发明母亲及胎儿心电信号QRS波中R波的检测定位方法流程示意框图；

图4是本发明使用到的低通高通滤波器的幅频与相频相应；

图5是本发明使用到的微分滤波器的幅频与相频相应；

图6是原始测试数据与滤波处理后的心电信号波形；

图7是本发明经过前端滤波处理后的数据再经过非线性滤波处理后的波形；

图8是图7的放大波形；

图9是图7中数据经过一级小波分解后的近似部分与细节部分波形图；

图10是Rdetect算法流程框图，用于定位检测出母亲心电信号中R波的位置；

图11是本发明胎儿心电信号R波定位的算法结构；

图12是针对Xsig(n)进行Rdetect算法处理后的波形，包括胎儿以及母亲心电的R波定位；

图13是图12中前端处理后得到的数据波形；

图14是图12中第五路信号放大的局部波形图；

图15是经过信号调整后在原始测试数据上标记的结果图；

图16是原始测试数据放大后的波形图。

具体实施方式

下面结合具体实施例对本发明作进一步具体详细描述，但本发明的实施方式不限于此，对于未特别注明的工艺参数，可参照常规技术进行。

实施例

图1是现有心电信号中的QRS波形；图2是现有通用的QRS波检测的算法结构；

如图3所示，本发明母亲及胎儿心电信号QRS波中R波的检测定位方法，包括如下步骤：

(1)前端滤波处理，包括低通滤波、高通滤波、微分器以及平方操作；

(2)小波分解，对于经过步骤(1)中前端滤波处理的混合心电信号，根据信号的特征，需要经过小波分解过程，以更加突出胎儿心电信号R波的份量；

(3)对步骤(2)处理后的信号采用Rdetect算法检测出混合心电信号中母亲心电信号QRS波中R波特征点出现的位置；

(4)对经过上述处理的混合心电信号，进行前端处理，并保留信号的辅助信息；

(5)采用Rdetect算法检测出经过前端处理后的信号中，胎儿心电信号QRS波中R波特征点出现的位置；

(6)后处理，利用前端处理中的辅助信息，恢复信号；

(7)处理过程中，会引入延迟；根据这些延迟调整信号。

上述步骤(2)中，采用小波一层分解，取细节部分，便于后续胎儿心电信号QRS波中R波位置的检测。

上述步骤(3)中，经过滤波以及小波分解后，采用Rdetect算法，其中Rdetect算法采用自适应阈值的方式，定位检测到母亲心电的R波位置，Rdetect算法参数的选择如下：

初始阈值：threshold＝param×max{y(n)}，这里param取0.7，也可以取0.7相近的一些数值；

自适应阈值：threshold＝param1×threshold+param1×maxy，其中，maxy是在检测R波中，出现的幅值最大值，param1取值为0.4，也可以取与0.4很相近的其他数值。

上述步骤(4)中，所述前端处理中，采用的方法是去除母亲R波位置出前后一定范围的数据点，并把经过这一操作的信号进行重新组合得到新的信号；在进行前端处理的过程中，需记录信号的辅助信息，即母亲心电信号R波出现的位置点，该R波的数目以及每个R波出现位置处去除的点数；根据心电信号准周期特性，每个R波位置处去除的点数近似相等，记为POINTS，POINTS的经验值是：

POINTS＝round(Fs/10)

其中，Fs是采样频率，round(x)表示对数值x取四舍五入的整数值。

上述步骤(5)采用Rdetect算法，检测出胎儿心电信号的R波位置，这里的输入信号是步骤(4)中重新组合的信号，所以经过Rdetect算法处理后，步骤(6)后处理中需要根据步骤(4)中母亲R波点的位置，加上减去的信号采样点数。

上述步骤(7)中，需要考虑之前处理过程中引入的延迟，计算滤波处理的延迟周期，即对应延迟的采样点数，后续的处理步骤中，是把信号从时域变换到小波域。

下面结合原理具体说明本发明的母亲及胎儿心电信号QRS波中R波的检测定位方法，

(1)对于典型的QRS波，其频谱分量一般分布在10Hz到25Hz，中心频率大致在17Hz。因此，大多QRS检测算法都会采用滤波器预处理过程，滤除其他信号分量，这些信号分量主要包括P波、T波、基线漂移、50Hz的电极线干扰，以及其他非耦合噪声等。P波、T波以及基线漂移其频带分布在较低频段，滤除它们，一般采用数字高通滤波器；而对于非耦合的噪声，对其进行抑制，一般选用低通滤波器。高通与低通滤波器组合一起，也就形成了我们的所说的带通滤波器，其截至频率大约为10Hz以及25Hz。

在本发明中，采用Hamilton和Tompkins在A real-time QRS Detection Algorithm(IEEE)中提出的滤波器模型。主要有三类线性滤波器：低通、高通、微分器。

低通滤波器的差分方程式：

y(nT)＝2y(nT-T)-y(nT-2T)+x(nT)-2x(nT-6T)+x(nT-12T)   (1)

这里T是采样周期，n是整数。其系统传输函数是

$H\left(z\right)=\frac{{(1-z^{-6})}^2}{{(1-z^{-1})}^2}---\left(2\right)$

对应的幅频响应以及相频响应如附图4(a)所示。

高通滤波器的差分方程：

y(nT)＝32x(nT-16T)-[y(nT-T)+x(nT)-x(nT-32T)]         (3)

其系统传输函数为

$H\left(z\right)=\frac{(-1+32z^{--16}+z^{-32})}{(1+z^{-1})}---\left(4\right)$

对应的幅频相频响应描绘在图4(b)中。

微分器滤波器的差分方程：

y(nT)＝(2x(nT)+x(nT-T)-x(nT-3T)-2x(nT-4T))/8         (5)

对应的系统传输函数为：

H(z)＝(1/8)(2+z^-1-z^-3-z^-4)                         (6)

对应的幅频相频响应显示在图5中。这里，本发明中采用的心电数据是通过在母亲胸部放置3根电极，在腹部放置5根电极采集得到，采样频率为250Hz，采样时间10s(当然这些参数是可以改变的)。这里选用其中一路混合心电数据来测试本发明的方法。图6显示了采用上述前端滤波处理后混合心电信号的波形。Ecg是原始的混合心电数据，a是心电信号经过低通后的波形，b是经过高通处理得到的结果，c是采用微分滤波器滤除后获得的波形。

这里前端滤波处理采用的是平方操作，对上述处理后的数据，其每一点进行平方操作，把数据点值都变为正值，以为后续的处理。图7为平方操作后的波形图。

(2)对图7进行放大，显示在图8中。从该图8中可以观测得到的是，主要的分量是母亲心电，胎儿的心电很微弱。通过算法检测出母亲的R波位置是不困难的，因为从图8可以看到，母亲R波位置上的幅度非常的之大。但目前，我们关键的是确定胎儿的R波位置。把胎儿的R波检测看做是奇异点检测，想到被称之为“数字显微镜”的小波分解。经过实验发现，小波分解是很恰当的。而且实验也发现，采用db6小波，进行1级分解就可以取得很好的效果。图9为小波分解处理后的波形。其中，a1是近似部分，而d1是细节部分。从心电信号特征分析可知，母亲心电信号的频率比胎儿心电信号频率稍低，而且部分噪声也分布在低频部分。图9显示，a1近似部分保留了图8中信号的大部分能量，波形很相似，胎儿心电依然不是很明显；而相对来说，d1细节部分，胎儿心电R波位置较明显，那些幅度较小的部位就是胎儿R波出现的位置。根据这种特性，后续处理中，我们选用的就是细节部分d1。

(3)自适应阈值R波判决。这里需要判决母亲心电信号的R波位置，以及胎儿心电R波的位置。目前，有较多的算法也是采用自适应阈值的方法来检测出母亲心电信号R波的位置，但由于之前预处理阶段所采用的方法不同，这里的阈值选择方案也相应的有些变化，

这里，针对检测母亲R波位置。如图10的算式结构。首先根据预处理后的数据集{x(n)}，选取约十分之一(这个数值可调，可以根据具体情况而变，只要选取的数据点数不至于过小就可，这里选择：的数据点，进行类似微分操作，

y(n)＝2x(n)+x(n-1)-x(n-3)-2x(n-4)     (7)

并根据

threshold＝0.7×max{y(n)}             (8)

获得初始的阈值。之后遍历数据x(n)，与初始阈值进行比较；在比较的过程之中，确定第一个R波的位置，并同时根据

threshold＝0.4×threshold+0.4×maxy   (9)

(其中，maxy是在检测R波中，出现的幅值最大值。)更新阈值，以为后续R波的检测。整个过程在图10中有很详细的操作流程。

以上是就母亲心电R波检测所采用的算法(这里称之为Rdetect算法)，而对于胎儿心电R波位置的定位，所采用的方法和以上的很相似，但是在检测之前，需要进行一定的预处理。由于在图9中的d1图形中，母亲的R波位置波形幅度明显相对于胎儿R波位置幅度大得多，而且周围的点数值也较大，这里根据采用频率，在以上检测到的母亲R波位置的前后减去POINTS个点，其中，POINTS的经验值是

POINTS＝round(Fs/10)                (10)

其中，Fs是采样频率，round(x)表示的意思是对数值x取四舍五入的整数值。再对新的到的数据集{xn(n)}，进行Rdetect算法处理，搜寻得到胎儿R波的位置点。这里，为确保搜寻结果的正确，在减去数据点时，需要记录母亲R波位置点处的信息，以便恢复之用。图11重复描述了整个过程。

还有非常重要的一点需要说明的是，在预处理阶段，我们使用了三类滤波器，这三类滤波器在处理时，都会引入延迟，式(2)的低通滤波器延迟6个采样周期，式(4)高通滤波器引入16个采样周期，而微分滤波器会引入2个采样周期。加起来，总共24个采样周期。而在后续的R波检测过程中，使用的数据集是滤波器处理后的数据，所以，在原始的混合心电数据上，需要把延迟考虑进去，以准确的确定母亲以及胎儿R波的位置。

(4)实验测试结果，在图12中，第一路信号是小波分解后的细节部分的数据d1，仍然是2500个数据点，这里把数据d1记为Xsig(n)，其中n的1，2，L。第二路是对Xsig(n)进行Rdetect算法处理得到的结果，也就是母亲心电R波位置的标记，红×位置处标记的就是母亲心电信号R波出现的位置，在第三路中把Xsig(n)和红×绘制到一起，从这路信号波形可以看出，标记的非常的准确。

图13是对信号Xsig(n)进行图12中所示的前端预处理后的数据波形，这里记图13的数据为Fsig(n)，n＝1，2，L。图12中第4路信号是对信号Xsig(n)进行图11所示算法处理后得到的R波位置标记图。也可以说是，对Fsig(n)进行Rdetect算法处理，并进行图11中的后处理得到的胎儿心电信号中R波位置的标记。

图14是对图12中第五路信号放大观察的局部波形，从实验现象来看，对胎儿心电的R波的定位，也是非常的准确的。

图15是考虑延迟后，在原始混合心电数据标记的结果。另外，试验中，我们对其他七路采集到的混合心电信号进行了同样的操作，其实验结果列在下表1中。从表1中可以看到，检测的效果接近100％，几乎没有遗漏的情况。

表1：8路混合心电信号的检测结果

在表1中，第1至5路信号是电极从腹部采集得到，而第6、7、8路信号是电极贴在孕妇胸部采集得到，其中，MR表示母亲心电的R波数目，FR表示胎儿R波的数目，TR是实际的R波数目，而DR是检测到的R波数目，errR表示相对误差率，

$\mathrm{errR}=\frac{|\mathrm{DR}-\mathrm{TR}|}{\mathrm{TR}}\×100\%---\left(12\right)$

图16是各路采集信号放大后的波形。从波形上也可以看出，从腹部采集的心电数据中，胎儿的成分比较明显，而在胸部采集到的数据，几乎很难看到胎儿的成分，这也是造成表1中6、7、8路数据误差率较大的主要原因。而第四路信号虽说是从腹部采集的到，但是，从图16看到，噪声干扰过大，造成了误检。

如上所述便可较好地实现本发明。

上述实施例为本发明较佳的实施方式，但本发明的实施方式并不受上述实施例的限制，其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化，均应为等效的置换方式，都包含在本发明的保护范围之内。

Claims (6)

1.一种母亲及胎儿心电信号QRS波中R波的检测定位方法，其特征在于包括如下步骤：

(1)前端滤波处理，包括低通滤波、高通滤波、微分器以及平方操作；

(2)小波分解，对于经过步骤(1)中前端滤波处理的混合心电信号，根据信号的特征，需要经过小波分解过程，以更加突出胎儿心电信号R波的份量；

(3)对步骤(2)处理后的信号采用Rdetect算法检测出混合心电信号中母亲心电信号QRS波中R波特征点出现的位置；

(4)对经过上述处理的混合心电信号，进行前端处理，并保留信号的辅助信息；

(5)采用Rdetect算法检测出经过前端处理后的信号中，胎儿心电信号QRS波中R波特征点出现的位置；

(6)后处理，利用前端处理中的辅助信息，恢复信号；

(7)处理过程中，会引入延迟；根据这些延迟调整信号。

2.根据权利要求1所述的检测定位方法，其特征在于上述步骤(2)中，采用小波一层分解，取细节部分，便于后续胎儿心电信号QRS波中R波位置的检测。

3.根据权利要求1所述的检测定位方法，其特征在于上述步骤(3)中，经过滤波以及小波分解后，采用Rdetect算法，其中Rdetect算法采用自适应阈值的方式，定位检测到母亲心电的R波位置，Rdetect算法参数的选择如下：

初始阈值：threshold＝param×max{y(n)}，这里param取0.7；

自适应阈值：threshold＝param1×threshold+param1×maxy，其中，maxy是在检测R波中，出现的幅值最大值，param1取值为0.4。

4.根据权利要求1所述的检测定位方法，其特征在于上述步骤(4)中，所述前端处理中，采用的方法是去除母亲R波位置出前后一定范围的数据点，并把经过这一操作的信号进行重新组合得到新的信号，在进行前端处理的过程中，记录信号的辅助信息，即母亲心电信号R波出现的位置点，该R波的数目以及每个R波出现位置处去除的点数；根据心电信号准周期特性，每个R波位置处去除的点数近似相等，记为POINTS，POINTS的经验值是：

POINTS＝round(Fs/10)

其中，Fs是采样频率，round(x)表示对数值x取四舍五入的整数值。

5.根据权利要求1所述的检测定位方法，其特征在于上述步骤(5)采用Rdetect算法，检测出胎儿心电信号的R波位置，这里的输入信号是步骤(4)中重新组合的信号，所以经过Rdetect算法处理后，步骤(6)后处理中需要根据步骤(4)中母亲R波点的位置，加上减去的信号采样点数。

6.根据权利要求1所述的检测定位方法，其特征在于上述步骤(7)中，需要考虑之前处理过程中引入的延迟，计算滤波处理的延迟周期，即对应延迟的采样点数，后续的处理步骤中，是把信号从时域变换到小波域。

Images