一种基于非负盲分离胎儿心电瞬时心率识别方法及系统的制作方法
【专利摘要】本发明公开一种基于非负盲分离胎儿心电瞬时心率识别方法及系统，该方法包括：对胎儿心电图机采集得到的心电信号进行预处理，小波去噪后利用相干平均法去掉母亲心电成分获取胎儿心电信号，将胎儿心电信号变换到时频域，得到Choi-Williams谱和Born-Jordan谱，再变为非负频谱；利用非负矩阵分解进行盲分离，得到胎儿心电信号的非负时频分量；求瞬时心率，对所求胎儿心电信号的非负时频分量进行极大值检测，设定阈值；采用区间遍历法求信号区间的最大值，求出心电曲线的所有峰值点；将求得心电曲线的所有峰值坐标进行换算，依次进行前向差分，由瞬时心率公式求得胎儿的瞬时心率。本发明的识别具有快速简单，准确率高的特点。
【专利说明】一种基于非负盲分离胎儿心电瞬时心率识别方法及系统
【技术领域】
[0001] 本发明涉及医疗器械领域，更具体地，涉及一种基于非负盲分离胎儿心电瞬时心 率识别方法及系统。
【背景技术】
[0002] 胎儿心电图是一种重要的胎儿监护手段，是一种检测胎儿在子宫内健康状况的重 要方法，能够反映胎儿心脏的全貌。健康的胎儿随子宫内环境的不同时刻发生着微妙的变 化，而且胎儿心率的变化是循环系统和中枢神经系统机能调节的表现，因此在围产期对胎 儿进行胎心率检测可以了解胎儿在子宫内的重要监护指标。然而由于采集的胎儿心电信号 幅度小，信号的噪声比低且大都混杂有幅度大，分布广的噪声干扰，例如常见的工频干扰、 呼吸、肌电等引起的噪声，其中最主要的干扰是母体心电信号（MaternalECG，MECG)，它给 医学诊断带来极大的困难。在时域和频域上，这些干扰信号混叠在一起对胎儿心音信号的 瞬时心率计算造成很大的影响。因此，研宄如何准确、便捷，有效地测量围产期胎心电的瞬 时心率具有及其重要的实用价值和临床意义。
[0003] 目前已经有许多采用其它技术来提取胎儿心电信号（FECG)的报道，如相干平均、 自相关和互相关、自适应滤波、基于奇异值分解和主分量分析等。但这些技术的提取方法受 很多限制，或计算过于复杂，或需要更多的人工干预。
[0004] 非负矩阵分解（NonnegativeMatrixFactorization，NMF)是近十年来随着盲信 号分离技术的发展而出现的一种新的信号处理技术。由于它的算法简单，计算量小，适用范 围广，已逐渐成为信号处理领域的热点之一。
【发明内容】
[0005] 为了克服上述现有技术的不足，本发明首先提出一种基于非负盲分离胎儿心电瞬 时心率识别方法，针对由胎儿心电图机采集得到的胎儿心电信号，本发明能快速简便的计 算出胎儿的瞬时心率，更好的服务于临床胎儿监护。
[0006] 本发明的又一目的是提出一种基于非负盲分离胎儿心电瞬时心率识别系统。
[0007] 为了实现上述目的，本发明的技术方案为：
[0008] 一种基于非负盲分离胎儿心电瞬时心率识别方法，包括以下步骤：
[0009] SI.对胎儿心电图机采集得到的心电信号进行预处理，其预处理包括小波去噪、 利用相干平均法去掉心电信号中母亲心电成分获取胎儿心电信号，并其将胎儿心电信 号变换到时频域，得到Choi-Williams谱和Born-Jordan谱，再将Choi-Williams谱和 Born-Jordan谱变为非负频谱；
[0010] S2.利用非负矩阵分解进行盲分离，对胎儿心电信号的Choi_Wi11iams谱和 Born-Jordan谱进行处理，得到胎儿心电信号的非负时频分量；
[0011] S3.求取瞬时心率，包括峰值点检测和差分法求瞬时心率，其中峰值点检测是对所 求胎儿心电信号的非负时频分量进行极大值检测，设定阈值；采用区间遍历法求心电信号
[0012] 进一步的，所述步骤Sl中的小波去噪的具体实现方式为：采用Daubechies小波， 小波基为db4,进行2级分解，取小波变换后的细节系数第二级信号xdl进行阈值去噪操作， 去噪的阈值取信号幅度绝对值最大值的〇. 06?0. 08倍；
[0013] 利用相干平均法去掉心电信号中母亲心电成分的具体实现方式为：根据母亲心电 信号成分的模板将心电信号中的母亲心电信号依次减去，所得心电信号为预处理后的胎儿 心电信号；
[0014] 所述信号变换到时频域的具体实现方式为：利用Choi-Williams时频变换和 Born-Jordan时频变换对胎儿心电信号xdl进行处理得到Choi-Williams谱、Born-Jordan 谱，然后分别与它们的共轭复数相乘即可得到非负Choi-Williams谱和Born-Jordan谱， 记非负Choi-Williams谱为X1，非负Born-Jordan谱为X2，并将xJPX2合并为矩阵X=
[X1,X2] 〇
[0015] 进一步的，所述步骤S2中利用非负矩阵分解进行盲分离的方式为：
【权利要求】
1. 一种基于非负盲分离胎儿心电瞬时心率识别方法，其特征在于，包括以下步骤：
51. 对胎儿心电图机采集得到的心电信号进行预处理，其预处理过程包括小波去噪、 利用相干平均法去掉心电信号中母亲心电成分获取胎儿心电信号，并其将胎儿心电信 号变换到时频域，得到Choi-Williams谱和Born-Jordan谱，再将Choi-Williams谱和 Born-Jordan谱变为非负频谱；
52. 利用非负矩阵分解进行盲分离，对胎儿心电信号的Choi-Williams谱和 Born-Jordan谱进行处理，得到胎儿心电信号的非负时频分量；
53. 求取瞬时心率，包括峰值点检测和差分法求瞬时心率，其中峰值点检测是对所求胎 儿心电信号的非负时频分量进行极大值检测，设定阈值；采用区间遍历法求心电信号区间 的最大值，再求出心电曲线的所有峰值点；差分法求瞬时心率是将所求得心电曲线的所有 峰值坐标进行换算，然后依次进行前向差分，再根据瞬时心率公式求得胎儿的瞬时心率。
2. 根据权利要求1所述的基于非负盲分离胎儿心电瞬时心率识别方法，其特征在于， 所述步骤S1中的小波去噪的具体实现方式为：采用Daubechies小波，小波基为db4,进行 2级分解，取小波变换后的细节系数第二级信号xdl进行阈值去噪操作，去噪的阈值取信号 幅度绝对值最大值的〇. 06?0. 08倍； 利用相干平均法去掉心电信号中母亲心电成分的具体实现方式为：根据母亲心电信号 成分的模板将心电信号中的母亲心电信号依次减去，所得心电信号为预处理后的胎儿心电 信号； 所述信号变换到时频域的具体实现方式为：利用Choi-Williams时频变换和Born-Jordan时频变换对胎儿心电信号xdl进行处理得到Choi-Williams谱、Born-Jordan 谱，然后分别与它们的共轭复数相乘即可得到非负Choi-Williams谱和Born-Jordan谱， 记非负Choi-Williams谱为Xi，非负Born-Jordan谱为x2，并将xJPx2合并为矩阵X= [x1；x2] 〇
3. 根据权利要求2所述的基于非负盲分离胎儿心电瞬时心率识别方法，其特征在于， 所述步骤S3中峰值点检测的具体方式为：对胎儿心电图非负时频分量进行分解和取值，得 到独立非负时频分量1和独立非负时频分量2,然后对独立非负时频分量1进行极大值检 测，对结果进行观测，设定对应的阈值；采用区间遍历法求信号区间的最大值，再求出该曲 线的所有峰值点；该步骤中采用阈值检测，阈值的定义如下：
其中，P为阈值，对取得的结果进行峰值观测，y取〇. 3，为所有极值点ei，i= 1，…，n的幅值的平均值。
4. 根据权利要求3所述的基于非负盲分离胎儿心电瞬时心率识别方法，其特征在于， 所述步骤S3中差分法求瞬时心率的具体方式为：将所求得曲线的所有峰值坐标进行换算， 然后依次进行前向差分，再根据瞬时心率公式求得胎儿的瞬时心率； 其中瞬时心率为：
At每两相邻峰值的时间间隔，即所求的u为瞬时心率。
5. -种基于非负盲分离胎儿心电瞬时心率识别系统，其特征在于，包括： 预处理模块，用于对胎儿心电图机采集得到的心电信号进行预处理，其预处理过程包 括小波去噪、利用相干平均法去掉心电信号中母亲心电成分获取胎儿心电信号，并其将胎 儿心电信号变换到时频域，得到Choi-Williams谱和Born-Jordan谱，再将Choi-Williams 谱和Born-Jordan谱变为非负频谱； 盲分离模块，采用非负矩阵分解进行盲分离，是对胎儿心电信号的Choi-Williams谱 和Born-Jordan谱进行处理，得到胎儿心电信号的非负时频分量； 瞬时心率获取模块，用于求取瞬时心率，具体为：包括峰值点检测和差分法求瞬时心 率，其中峰值点检测是对所求胎儿心电信号的非负时频分量进行极大值检测，设定阈值；采 用区间遍历法求心电信号区间的最大值，再求出心电曲线的所有峰值点；差分法求瞬时心 率是将所求得心电曲线的所有峰值坐标进行换算，然后依次进行前向差分，再根据瞬时心 率公式求得胎儿的瞬时心率。
6. 根据权利要求5所述的基于非负盲分离胎儿心电瞬时心率识别系统，其特征在于， 所述预处理模块中的小波去噪的具体实现方式为：采用Daubechies小波，小波基为db4,进 行2级分解，取小波变换后的细节系数第二级信号xdl进行阈值去噪操作，去噪的阈值取信 号幅度绝对值最大值的〇. 06?0. 08倍； 利用相干平均法去掉心电信号中母亲心电成分的具体实现方式为：根据母亲心电信号 成分的模板将心电信号中的母亲心电信号依次减去，所得心电信号为预处理后的胎儿心电 信号； 所述信号变换到时频域的具体实现方式为：利用Choi-Williams时频变换和Born-Jordan时频变换对胎儿心电信号xdl进行处理得到Choi-Williams谱、Born-Jordan 谱，然后分别与它们的共轭复数相乘即可得到非负Choi-Williams谱和Born-Jordan谱， 记非负Choi-Williams谱为Xi，非负Born-Jordan谱为x2，并将xJPx2合并为矩阵X= [x1；x2] 〇
7. 根据权利要求6所述的基于非负盲分离胎儿心电瞬时心率识别系统，其特征在于， 所述瞬时心率获取模块中峰值点检测的具体方式为：对胎儿心电图非负时频分量进行分解 和取值，得到独立非负时频分量1和独立非负时频分量2,然后对独立非负时频分量1进行 极大值检测，对结果进行观测，设定对应的阈值；采用区间遍历法求信号区间的最大值，再 求出该曲线的所有峰值点；该步骤中采用阈值检测，阈值的定义如下：
其中，P为阈值，对取得的结果进行峰值观测，y取〇. 3，
为所有极值点ei，i= 1，…，n的幅值的平均值。
8. 根据权利要求7所述的基于非负盲分离胎儿心电瞬时心率识别系统，其特征在于， 所述瞬时心率获取模块中差分法求瞬时心率的具体方式为：将所求得曲线的所有峰值坐标 进行换算，然后依次进行前向差分，再根据瞬时心率公式求得胎儿的瞬时心率； 其中瞬时心率为：
At每两相邻峰值的时间间隔，即所求的u为瞬时心率。
【文档编号】A61B5/0402GK104473631SQ201410777741
【公开日】2015年4月1日 申请日期:2014年12月12日 优先权日:2014年12月12日
【发明者】康嘉文, 谢侃, 谢胜利, 蔡坤, 王旭 申请人:广东工业大学