CN116421196A - 一种基于母体心电信号的胎儿心率提取方法 - Google Patents

Abstract

本发明的一种基于母体心电信号的胎儿心率提取方法，属于医疗器械技术领域，包括以下步骤：S1：给孕妇佩戴听诊式胎儿监护仪，将信号分解成由高频到低频，对收集到的孕妇腹部所有频段的心电信号作初步处理；S2：将获得的胎儿心电信号再做降噪处理。本方法通过对孕妇腹部电信号进行预处理，使用改进的假阳性和假阴性校正算法对初步定位结果进行校正，消除母体心电信号干扰，减弱相关性，提取孕妇腹部电信号中的胎儿心电信号信息，然后重构胎儿心电信号,从而建立一个新的多通道胎儿心电信号提取方法，整个方法简便且易于实现，能够有效地在母体腹壁混合心电信号中提取出胎儿心电信号，为医生对胎儿疾病的诊疗提供保证。

Description

一种基于母体心电信号的胎儿心率提取方法

技术领域

本发明涉及医疗器械领域，具体讲是一种基于母体心电信号的胎儿心率提取方法。

背景技术

胎儿心电信号是各种心电信号中重要且常见的一种，目前临床上主要通过监测心音、心动和心电图的变化来诊断胎儿在子宫内的发育情况。其中，胎儿心电图能最好地体现胎儿心脏活动，对胎儿发育过程中出现的问题更敏感，能为临床诊断提供更可靠的依据。但是胎儿心电信号非常微弱，从母体体表提取胎儿心电信号会受到各种噪声的严重干扰，最主要的为母体心电信号、50Hz工频干扰以及基线漂移，这对后期的分析处理带来了很大的不便。

虽然测量孕妇腹部电信号的方式更加简便,但由于胎儿心电信号比较微弱,与母体心电在频谱上有重叠,并且孕妇腹部的胎儿心电信号还受到各种噪声的污染,噪声包括胎儿大脑活动、母亲心电干扰、子宫收缩干扰、肌电干扰、50Hz/60Hz电力线干扰和运动伪迹干扰等，想要替代常用的多普勒超声检测技术仍然还具有很大的挑战性。

发明内容

针对上述技术的不足，本发明公开一种基于母体心电信号的胎儿心率提取方法，本发明通过对提取母体腹部信号的预处理，对胎儿QRS波群实现定位、校正，通过去噪和对母体心电信号信息消除，再基于周期性心电信号方法的信息提取和胎儿心电信号重构，最后相加得到提取的最终胎儿心电信号。

为了实现上述技术效果，本发明采用以下技术方案：

一种基于母体心电信号的胎儿心率提取方法，包括以下步骤：

S1：给孕妇佩戴听诊式胎儿监护仪，通过其胎心率探头和宫缩探头提取孕妇腹部的心电信号，通过迭代筛选过程，将信号分解成由高频到低频，对收集到的孕妇腹部所有频段的心电信号作初步处理；

S2：将初步获取的孕妇腹部电信号减去经过微调处理后的最佳匹配的母体心电信号模板信号，获得胎儿心电信号，将获得的胎儿心电信号再做降噪处理；

S3：依据胎儿心电和母体心电的基频和幅值差异，消除频段信号中母体心电的QRS波群干扰,利用胎儿心电的QRS波群中显著的频段信号，检测胎儿心电的特征周期；

S4：采用准周期分量提取胎儿心电，然后根据周期性修缮重构胎儿心电，将多通道腹壁混合心电信号分别进行分解，得到各通道各自心电频段，将各通道信号中相同阶数的心电频段进行提取，得到胎儿心电的各阶心电频段，最后相加得到提取的最终胎儿心电；

S5：在胎儿心电信号中根据初步定位的胎儿QRS波群位置对应的幅值判断胎儿心电信号的极性；

当极性判断为正时，则在初步提取的R峰位置周围0.08s的局部区间内使用局部最大值函数重新确定胎儿R峰的位置；

根据胎儿RR间期序列的中间部分计算RR间期均值RRm作为整个RR间期序列的平均值；

利用每个RR间期序列的间期长度特征对初步定位的胎儿QRS波群进行校正。

S6：基于聚类分析提取的胎儿QRS波群的可行性分析过程，从预处理的胎儿心电信号中选取特征作为未分类的数据集，采用欧几里得距离作为相似性度量来进行聚类分析；

S7：根据预处理后的胎儿心电信号的极性,另选取相应的极值对(极大-极小对)作为未聚类的数据集，进行聚类分析，最后对两种聚类结果使用交叉阑值的方法进行处理。

进一步的，所述步骤S1中初步处理方法包括：

S1.1使用小波软阈值降噪算法去除单通道孕妇腹部电信号的高频噪声；

S1.2使用小波分解去除8Hz以下的基线和运动伪迹干扰,去除的低频成分包括母体心电信号中的P、T波；

S1.3使用50Hz陷波器去除电力线干扰。

进一步的，所述步骤S2中降噪处理方法包括：

S2.1对胎儿心电信号经过小波软阈值降噪算法去除高频噪声；

S2.2使用零相位低通滤波器提取低频成分。

进一步的，所述步骤S3中消除信号干扰方法包括：

S3.1对预处理后的胎儿心电信号进行导数滤波，然后对信号导数的幅度平方处理；

S3.2使用三次曲线函数增强平方处理后的信号中的对应胎儿QRS波群的部分；

S3.3对增强后的信号选取合适的闯值,提取出胎儿QRS波群对应的子区间；

S3.4针对信号导数的幅度平方处理后的信号分别在提取的子区间中寻找最大值对应的位置,即是初步提取的胎儿QRS波群位置；

S3.5对胎儿QRS波群的假阳性和假阴性进行校正。

进一步的，所述步骤S3中胎儿心电的特征周期获取方法包括：

假设由m个混合信号构成的矩阵为X＝(x₁.x₂,...x_m,)^T,其中,x为列向量表示的信号,是由n个相互独立的源信号s₁,s₂,...,s_n,线性叠加而成的，且s为信号值构成的列向量。设S＝(s₁,s₂,...,s_n)^T,则有:

X＝AS

式中，A为m*n阶系数矩阵,称为混合矩阵。需要找到一个正交矩阵W＝(w1,w2,...,w_n)(其中w为列向量),使得:

Y＝W^TS

为源信号S的估计，则考虑Y＝(y₁,y₂,...y_n,)^T中的第i个信号为目标信号胎儿心电),具有准周期特性,其特征周期为T。

进一步的，所述步骤S5中校正方法包括：

定义R_i为胎儿的第i个R波的位置，i的取值为1,2，...,N,N为初步提取的胎儿R波的个数，则胎儿心电信号中第i个RR间期序列为：

RR_i＝R_i+1-R_i

所有R峰幅值的绝对值的均值定义为Ra，对胎儿QRS波群的校正过程为：

S5.1若RR_i≤0.55RRm,则移除R_i+1；

S5.2若R_i位置在原信号对应的幅值小于0.4Ra,则该位置的R峰是误检,需要移除R_i；

S5.3根据假阳性校正结果更新RRm和Ra，进行假阴性校正，从起始位置到第一个胎儿R峰的间期定义为

RR_head＝R₁-1；

最后一个胎儿R峰到结束位置的间期定义为

RR_tail＝length(fECG)-R_N。

进一步的，所述校正方法对起始间期RR_head或RR_tail的处理方法如下：

如果RR_head或RR_tail＞1.5RRm,起始位置缺失一个胎儿R峰或R谷,则在指定区间[R₁-1.3RRm,R_i-1.3RRm+0.6RRm]内根据信号极性寻找最值作为漏检的胎儿R峰或R谷；

如果0.75RRm＜RR_head或RR_tail＜1.55RRm，该间期可能错过一个R峰或R谷,指定搜索区间为[1，0.5RRm]，搜索区间的最值为Ramp，在满足判断条件abs(Ramp)＞0.6Ra时,则此区间有漏检胎儿R峰或R谷，反之则RR_head或RR_tail间期正常。

进一步的，所述校正方法对RR间期序列剩余部分RR_i的处理方法如下：

若RRi＞1.55RRm，则判断为在该RR_i间期缺失至少一个R峰,反之则表示该区间表现正常；

对于缺失R峰的RR_i间期，根据:

m＝round(RR_i/RRm)-1

判断该区间胎儿R峰缺失的个数,该区间指定对应的m个片段,在每个片段中搜索缺失的对应的R峰,每个片段的指定搜索区间为

其中,n决定每个片段的搜索区间的大小,最大值为6,即在每个片段中逐渐扩展搜索区间寻找一个大于指定阀值的局部最值,其对应的位置为漏检的胎儿R峰位置。

进一步的，所述步骤S6聚类分析方法包括：

S6.1根据预处理的胎儿心电信号的极大值和极小值均值的绝对值判断信号的极性；

S6.2如果胎儿心电信号的极性为正，则选择极大值作为未聚类的数据集；

S6.3如果胎儿心电信号的极性为负,则选择极小值作为未聚类的数据集。

进一步的，所述步骤S7中交叉阈值分析与处理过程包括：

S7.1分析过程：基于极值特征的聚类分析结果被定义为E_i(i＝1,2,...,m),其中m是被选为胎儿QRS波群类别极值的个数,并计算对应位置极值的均值为EM，以极值对为聚类特征的聚类结果被定义为P_i(i＝1,2,...,n),其中n是被选为胎儿QRS波群类别的极值对的个数,并计算对应位置极值对的特征值的均值为PM；

S7.2处理过程：若P_i(i＝1,2,...,n)对应位置的胎儿R峰极值大于0.5M，则该极值对对应位置为初步提取到的胎儿QRS波群位置；若E_i(i＝1,2,...,m)对应位置的极值对的特征值大于0.02PM，则该极值对应的位置为提取到的胎儿QRS波群位置。

本发明通过改进在此提供一种基于母体心电信号的胎儿心率提取方法，与现有技术相比，具有如下改进及优点：

本发明一种基于母体心电信号的胎儿心率提取方法，通过对孕妇腹部电信号进行预处理，使用基于聚类分析算法初步定位母体QRS波群,使用改进的假阳性和假阴性校正算法对初步定位结果进行校正，然后从预处理后的信号中移除母体心电信号模板信号后得到的剩余量就是胎儿心电信号，消除母体心电信号干扰，减弱相关性，提取孕妇腹部电信号中的胎儿心电信号信息，然后重构胎儿心电信号,从而建立一个新的多通道胎儿心电信号提取方法，整个方法简便且易于实现，能够有效地在母体腹壁混合心电信号中提取出胎儿心电信号，为医生对胎儿疾病的诊疗提供保证。

附图说明

下面结合附图和实施例对本发明作进一步解释:

图1为本发明的步骤示意图；

图2为本发明的初步处理方法步骤示意图；

图3为本发明的降噪处理方法步骤示意图；

图4为本发明的消除信号干扰方法步骤示意图；

图5为本发明的聚类分析方法步骤示意图。

具体实施方式

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

本发明通过改进在此提供一种基于母体心电信号的胎儿心率提取方法，如图1-2所示，本发明一种基于母体心电信号的胎儿心率提取方法第一步S1：给孕妇佩戴听诊式胎儿监护仪，通过其胎心率探头和宫缩探头提取孕妇腹部的心电信号，通过迭代筛选过程，将信号分解成由高频到低频，对收集到的孕妇腹部所有频段的心电信号作初步处理，初步处理方法包括：S1.1使用小波软阈值降噪算法去除单通道孕妇腹部电信号的高频噪声；S1.2使用小波分解去除8Hz以下的基线和运动伪迹干扰,去除的低频成分包括母体心电信号中的P、T波；S1.3使用50Hz陷波器去除电力线干扰；

如图3所示，本发明一种基于母体心电信号的胎儿心率提取方法第二步S2：将初步获取的孕妇腹部电信号减去经过微调处理后的最佳匹配的母体心电信号模板信号，获得胎儿心电信号，将获得的胎儿心电信号再做降噪处理，降噪处理方法包括：S2.1对胎儿心电信号经过小波软阈值降噪算法去除高频噪声；S2.2使用零相位低通滤波器提取低频成分；

如图4所示，本发明一种基于母体心电信号的胎儿心率提取方法第三步S3：依据胎儿心电和母体心电的基频和幅值差异，消除频段信号中母体心电的QRS波群干扰,利用胎儿心电的QRS波群中显著的频段信号，检测胎儿心电的特征周期，消除信号干扰方法包括：S3.1对预处理后的胎儿心电信号进行导数滤波，然后对信号导数的幅度平方处理；S3.2使用三次曲线函数增强平方处理后的信号中的对应胎儿QRS波群的部分；S3.3对增强后的信号选取合适的闯值,提取出胎儿QRS波群对应的子区间；S3.4针对信号导数的幅度平方处理后的信号分别在提取的子区间中寻找最大值对应的位置,即是初步提取的胎儿QRS波群位置；S3.5对胎儿QRS波群的假阳性和假阴性进行校正；

特征周期获取方法包括：

假设由m个混合信号构成的矩阵为X＝(x₁.x₂,...x_m,)^T,其中,x为列向量表示的信号,是由n个相互独立的源信号s₁,s₂,...,s_n,线性叠加而成的，且s为信号值构成的列向量。设S＝(s₁,s₂,...,s_n)^T,则有:

X＝AS

式中，A为m*n阶系数矩阵,称为混合矩阵。需要找到一个正交矩阵W＝(w1,w2,...,w_n)(其中w为列向量),使得:

Y＝W^TS

为源信号S的估计，则考虑Y＝(y₁,y₂,...y_n,)^T中的第i个信号为目标信号胎儿心电),具有准周期特性,其特征周期为T。

本发明一种基于母体心电信号的胎儿心率提取方法第四步S4：采用准周期分量提取胎儿心电，然后根据周期性修缮重构胎儿心电，将多通道腹壁混合心电信号分别进行分解，得到各通道各自心电频段，将各通道信号中相同阶数的心电频段进行提取，得到胎儿心电的各阶心电频段，最后相加得到提取的最终胎儿心电；

本发明一种基于母体心电信号的胎儿心率提取方法第五步S5：在胎儿心电信号中根据初步定位的胎儿QRS波群位置对应的幅值判断胎儿心电信号的极性；当极性判断为正时，则在初步提取的R峰位置周围0.08s的局部区间内使用局部最大值函数重新确定胎儿R峰的位置；根据胎儿RR间期序列的中间部分计算RR间期均值RRm作为整个RR间期序列的平均值；利用每个RR间期序列的间期长度特征对初步定位的胎儿QRS波群进行校正；校正方法包括：定义R_i为胎儿的第i个R波的位置，i的取值为1,2，...,N,N为初步提取的胎儿R波的个数，则胎儿心电信号中第i个RR间期序列为：

RR_i＝R_i+1-R_i

所有R峰幅值的绝对值的均值定义为Ra，对胎儿QRS波群的校正过程为：S5.1若RR_i≤0.55RRm,则移除R_i+1；S5.2若R_i位置在原信号对应的幅值小于0.4Ra,则该位置的R峰是误检,需要移除R_i；S5.3根据假阳性校正结果更新RRm和Ra，进行假阴性校正，从起始位置到第一个胎儿R峰的间期定义为

RR_head＝R₁-1；

最后一个胎儿R峰到结束位置的间期定义为

RR_tail＝length(fECG)-R_N；

校正方法对起始间期RR_head或RR_tail的处理方法：如果RR_head或RR_tail＞1.5RRm,起始位置缺失一个胎儿R峰或R谷,则在指定区间[R₁-1.3RRm,R_i-1.3RRm+0.6RRm]内根据信号极性寻找最值作为漏检的胎儿R峰或R谷；如果0.75RRm＜RR_head或RR_tail＜1.55RRm，该间期可能错过一个R峰或R谷,指定搜索区间为[1，0.5RRm]，搜索区间的最值为Ramp，在满足判断条件abs(Ramp)＞0.6Ra时,则此区间有漏检胎儿R峰或R谷，反之则RR_head或RR_tail间期正常；校正方法对RR间期序列剩余部分RR_i的处理方法如下：若RRi＞1.55RRm，则判断为在该RR_i间期缺失至少一个R峰,反之则表示该区间表现正常；对于缺失R峰的RR_i间期，根据:

m＝round(RR_i/RRm)-1

判断该区间胎儿R峰缺失的个数,该区间指定对应的m个片段,在每个片段中搜索缺失的对应的R峰,每个片段的指定搜索区间为

其中,n决定每个片段的搜索区间的大小,最大值为6,即在每个片段中逐渐扩展搜索区间寻找一个大于指定阀值的局部最值,其对应的位置为漏检的胎儿R峰位置。

如图5所示，本发明一种基于母体心电信号的胎儿心率提取方法第六步S6：基于聚类分析提取的胎儿QRS波群的可行性分析过程，从预处理的胎儿心电信号中选取特征作为未分类的数据集，采用欧几里得距离作为相似性度量来进行聚类分析，聚类分析方法包括：S6.1根据预处理的胎儿心电信号的极大值和极小值均值的绝对值判断信号的极性；S6.2如果胎儿心电信号的极性为正，则选择极大值作为未聚类的数据集；S6.3如果胎儿心电信号的极性为负,则选择极小值作为未聚类的数据集。

本发明一种基于母体心电信号的胎儿心率提取方法第七步S7：根据预处理后的胎儿心电信号的极性,另选取相应的极值对(极大-极小对)作为未聚类的数据集，进行聚类分析，最后对两种聚类结果使用交叉阑值的方法进行处理，S7.1分析过程：基于极值特征的聚类分析结果被定义为E_i(i＝1,2,...,m),其中m是被选为胎儿QRS波群类别极值的个数,并计算对应位置极值的均值为EM，以极值对为聚类特征的聚类结果被定义为P_i(i＝1,2,...,n),其中n是被选为胎儿QRS波群类别的极值对的个数,并计算对应位置极值对的特征值的均值为PM；S7.2处理过程：若P_i(i＝1,2,...,n)对应位置的胎儿R峰极值大于0.5M，则该极值对对应位置为初步提取到的胎儿QRS波群位置；若E_i(i＝1,2,...,m)对应位置的极值对的特征值大于0.02PM，则该极值对应的位置为提取到的胎儿QRS波群位置。

工作原理：首先对孕妇腹部电信号进行预处理，使用基于聚类分析算法初步定位母体QRS波群,使用改进的假阳性和假阴性校正算法对初步定位结果进行校正，然后从预处理后的信号中移除母体心电信号模板信号后得到的剩余量就是胎儿心电信号，消除母体心电信号干扰，减弱相关性，提取孕妇腹部电信号中的胎儿心电信号信息，然后重构生成胎儿心电信号。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。

Claims (10)

1.一种基于母体心电信号的胎儿心率提取方法，其特征在于，包括以下步骤：

S1：给孕妇佩戴听诊式胎儿监护仪，通过其胎心率探头和宫缩探头提取孕妇腹部的心电信号，通过迭代筛选过程，将信号分解成由高频到低频，对收集到的孕妇腹部所有频段的心电信号作初步处理；

S2：将初步获取的孕妇腹部电信号减去经过微调处理后的最佳匹配的母体心电信号模板信号，获得胎儿心电信号，将获得的胎儿心电信号再做降噪处理；

S3：依据胎儿心电和母体心电的基频和幅值差异，消除频段信号中母体心电的QRS波群干扰,利用胎儿心电的QRS波群中显著的频段信号，检测胎儿心电的特征周期；

S4：采用准周期分量提取胎儿心电，然后根据周期性修缮重构胎儿心电，将多通道腹壁混合心电信号分别进行分解，得到各通道各自心电频段，将各通道信号中相同阶数的心电频段进行提取，得到胎儿心电的各阶心电频段，最后相加得到提取的最终胎儿心电；

S5：在胎儿心电信号中根据初步定位的胎儿QRS波群位置对应的幅值判断胎儿心电信号的极性；

当极性判断为正时，则在初步提取的R峰位置周围0.08s的局部区间内使用局部最大值函数重新确定胎儿R峰的位置；

根据胎儿RR间期序列的中间部分计算RR间期均值RRm作为整个RR间期序列的平均值；

利用每个RR间期序列的间期长度特征对初步定位的胎儿QRS波群进行校正。

S6：基于聚类分析提取的胎儿QRS波群的可行性分析过程，从预处理的胎儿心电信号中选取特征作为未分类的数据集，采用欧几里得距离作为相似性度量来进行聚类分析；

S7：根据预处理后的胎儿心电信号的极性,另选取相应的极值对(极大-极小对)作为未聚类的数据集，进行聚类分析，最后对两种聚类结果使用交叉阑值的方法进行处理。

2.根据权利要求1所述的一种基于母体心电信号的胎儿心率提取方法，其特征在于：所述步骤S1中初步处理方法包括：

S1.1使用小波软阈值降噪算法去除单通道孕妇腹部电信号的高频噪声；

S1.2使用小波分解去除8Hz以下的基线和运动伪迹干扰,去除的低频成分包括母体心电信号中的P、T波；

S1.3使用50Hz陷波器去除电力线干扰。

3.根据权利要求1所述的一种基于母体心电信号的胎儿心率提取方法，其特征在于：所述步骤S2中降噪处理方法包括：

S2.1对胎儿心电信号经过小波软阈值降噪算法去除高频噪声；

S2.2使用零相位低通滤波器提取低频成分。

4.根据权利要求1所述的一种基于母体心电信号的胎儿心率提取方法，其特征在于：所述步骤S3中消除信号干扰方法包括：

S3.1对预处理后的胎儿心电信号进行导数滤波，然后对信号导数的幅度平方处理；

S3.2使用三次曲线函数增强平方处理后的信号中的对应胎儿QRS波群的部分；

S3.3对增强后的信号选取合适的闯值,提取出胎儿QRS波群对应的子区间；

S3.4针对信号导数的幅度平方处理后的信号分别在提取的子区间中寻找最大值对应的位置,即是初步提取的胎儿QRS波群位置；

S3.5对胎儿QRS波群的假阳性和假阴性进行校正。

5.根据权利要求4所述的一种基于母体心电信号的胎儿心率提取方法，其特征在于：所述步骤S3中胎儿心电的特征周期获取方法包括：

假设由m个混合信号构成的矩阵为X＝(x₁.x₂,...x_m,)^T,其中,x为列向量表示的信号,是由n个相互独立的源信号s₁,s₂,...,s_n,线性叠加而成的，且s为信号值构成的列向量。设S＝(s₁,s₂,...,s_n)^T,则有:

X＝AS

式中，A为m*n阶系数矩阵,称为混合矩阵。需要找到一个正交矩阵W＝(w1,w2,...,w_n)(其中w为列向量),使得:

Y＝W^TS

为源信号S的估计，则考虑Y＝(y₁,y₂,...y_n,)^T中的第i个信号为目标信号胎儿心电),具有准周期特性,其特征周期为T。

6.根据权利要求1所述的一种基于母体心电信号的胎儿心率提取方法，其特征在于：所述步骤S5中校正方法包括：

定义R_i为胎儿的第i个R波的位置，i的取值为1,2，...,N,N为初步提取的胎儿R波的个数，则胎儿心电信号中第i个RR间期序列为：

RR_i＝R_i+1-R_i

所有R峰幅值的绝对值的均值定义为Ra，对胎儿QRS波群的校正过程为：

S5.1若RR_i≤0.55RRm,则移除R_i+1；

S5.2若R_i位置在原信号对应的幅值小于0.4Ra,则该位置的R峰是误检,需要移除R_i；

S5.3根据假阳性校正结果更新RRm和Ra，进行假阴性校正，从起始位置到第一个胎儿R峰的间期定义为

RR_head＝R₁-1；

最后一个胎儿R峰到结束位置的间期定义为

RR_tail＝length(fECG)-R_N。

7.根据权利要求6所述的一种基于母体心电信号的胎儿心率提取方法，其特征在于：所述校正方法对起始间期RR_head或RR_tail的处理方法如下：

如果RR_head或RR_tail＞1.5RRm,起始位置缺失一个胎儿R峰或R谷,则在指定区间[R₁-1.3RRm,R_i-1.3RRm+0.6RRm]内根据信号极性寻找最值作为漏检的胎儿R峰或R谷；

如果0.75RRm＜RR_head或RR_tail＜1.55RRm，该间期可能错过一个R峰或R谷,指定搜索区间为[1，0.5RRm]，搜索区间的最值为Ramp，在满足判断条件abs(Ramp)＞0.6Ra时,则此区间有漏检胎儿R峰或R谷，反之则RR_head或RR_tail间期正常。

8.根据权利要求6所述的一种基于母体心电信号的胎儿心率提取方法，其特征在于：所述校正方法对RR间期序列剩余部分RR_i的处理方法如下：

若RRi＞1.55RRm，则判断为在该RR_i间期缺失至少一个R峰,反之则表示该区间表现正常；

对于缺失R峰的RR_i间期，根据:

m＝round(RR_i/RRm)-1

判断该区间胎儿R峰缺失的个数,该区间指定对应的m个片段,在每个片段中搜索缺失的对应的R峰,每个片段的指定搜索区间为

其中,n决定每个片段的搜索区间的大小,最大值为6,即在每个片段中逐渐扩展搜索区间寻找一个大于指定阀值的局部最值,其对应的位置为漏检的胎儿R峰位置。

9.根据权利要求1所述的一种基于母体心电信号的胎儿心率提取方法，其特征在于：所述步骤S6聚类分析方法包括：

S6.1根据预处理的胎儿心电信号的极大值和极小值均值的绝对值判断信号的极性；

S6.2如果胎儿心电信号的极性为正，则选择极大值作为未聚类的数据集；

S6.3如果胎儿心电信号的极性为负,则选择极小值作为未聚类的数据集。

10.根据权利要求1所述的一种基于母体心电信号的胎儿心率提取方法，其特征在于：所述步骤S7中交叉阈值分析与处理过程包括：

S7.1分析过程：基于极值特征的聚类分析结果被定义为E_i(i＝1,2,...,m),其中m是被选为胎儿QRS波群类别极值的个数,并计算对应位置极值的均值为EM，以极值对为聚类特征的聚类结果被定义为P_i(i＝1,2,...,n),其中n是被选为胎儿QRS波群类别的极值对的个数,并计算对应位置极值对的特征值的均值为PM；

S7.2处理过程：若P_i(i＝1,2,...,n)对应位置的胎儿R峰极值大于0.5M，则该极值对对应位置为初步提取到的胎儿QRS波群位置；若E_i(i＝1,2,...,m)对应位置的极值对的特征值大于0.02PM，则该极值对应的位置为提取到的胎儿QRS波群位置。

Images