《一种胎儿心音瞬时心率检测识别方法及装置.pdf》由会员分享，可在线阅读，更多相关《一种胎儿心音瞬时心率检测识别方法及装置.pdf（7页完整版）》请在专利查询网上搜索。
1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201610928161.6 (22)申请日 2016.10.31 (71)申请人 广东工业大学 地址 510062 广东省广州市越秀区东风东 路729号 (72)发明人 资道利 谢侃 谢胜利 蔡坤 (74)专利代理机构 广东广信君达律师事务所 44329 代理人 杨晓松 (51)Int.Cl. A61B 5/00(2006.01) A61B 5/024(2006.01) (54)发明名称 一种胎儿心音瞬时心率检测识别方法及装 置 (57)摘要 本发明公开了一种胎儿心音瞬时心率检。
2、测 识别方法及装置， 本发明技术方案首先通过仪器 采集胎儿心音信号作为原始混合信号， 然后进行 滤波和归一化处理而得到胎儿心音信号作为预 处理信号， 再通过提取胎儿心音信号的特征以建 立并训练混合高斯模型， 通过训练完毕的混合高 斯模型对待测定的胎儿心音信号进行识别， 最终 得到胎儿心音瞬时心率。 本发明的技术方案相对 于现有技术检测识别方法， 能够快速便捷并准确 地检测识别出胎儿心音瞬时心率， 更好地服务于 临床胎儿监护。 权利要求书1页 说明书4页 附图1页 CN 106491085 A 2017.03.15 CN 106491085 A 1.一种胎儿心音瞬时心率检测识别方法， 其特征在于。
3、， 包括以下步骤： S1.通过仪器对胎儿心音信号进行采集并作为原始混合信号； S2.对采集到的所述原始混合信号进行滤波处理， 并作归一化处理， 而得到胎儿心音信 号作为预处理信号； S3.根据所述预处理信号， 提取所述胎儿心音信号的特征值； S4.将提取出来的所述特征值建立并训练混合高斯模型； S5.通过训练完毕的所述混合高斯模型对待测定的所述胎儿心音信号进行识别， 最终 得出所述胎儿心音信号瞬时心率。 2.如权利要求1所述的胎儿心音瞬时心率检测识别方法， 其特征在于， 所述步骤S4包括 以下步骤： S41.对所述混合高斯模型进行定义； S42.利用最大似然法对所述混合高斯模型的各个变量的参数。
4、进行估计； S43.直至所述混合高斯模型的参数收敛为止。 3.一种胎儿心音瞬时心率检测识别装置， 其特特征在于， 所述胎儿心音瞬时心率检测 识别装置使用如权利1所述的胎儿心音瞬时心率检测识别方法。 权 利 要 求 书 1/1 页 2 CN 106491085 A 2 一种胎儿心音瞬时心率检测识别方法及装置 技术领域 0001 本发明涉及胎心音信号领域， 特别涉及一种胎儿心音瞬时心率检测识别方法及装 置。 背景技术 0002 众所周知， 通过对胎儿的监护， 可以判断胎儿在子宫内的健康状况并作出合理的 预防与治疗。 早在19世纪， 通过胎儿听诊以估计胎儿在子宫内的健康状况的方法已经被产 科应用， 。
5、随着超声波多普勒技术的发展， 特别是电子胎儿监护已经成为目前最常见的胎儿 监护方法。 0003 现代医学证明， 由于胎盘脐带因素、 孕妇自身因素等都会造成胎儿宫内缺氧， 统称 为胎儿窘迫。 而通过对胎儿监护， 可及时发现胎儿的健康状况， 并及时进行治疗， 能够最大 限度降低先天缺陷儿的出生率， 从而提高我国的人口质量。 因此， 胎儿心率监测对于我国大 力提倡的优生优育有着极其重要的意义。 0004 现有技术， 由于硬件采集到的胎儿信号中， 通常包含胎儿信号、 母亲信号以及其他 大量的干扰等信号， 并且胎儿在母亲体内运动的声音信号强弱变化幅度较大， 而这些干扰 信号在时域和频域上混叠在一起， 对。
6、胎儿心音信号的瞬时心率的计算造成很大影响。 因此， 如何准确、 便捷、 有效监测胎儿心音的瞬时心率具有极其重要的应用价值和临床意义。 0005 现有技术中关于胎儿心音的瞬时心率的检测识别方法， 大致可以分为以下两类： 0006 1)基于心电模式的匹配滤波法计算胎儿心音心率。 其基本思路是， 首先将胎儿心 音电信号分离出来， 然后对胎儿心音电信号进行R波定位， 最后求出胎儿心音的瞬时心率。 0007 由于对胎儿心音电信号进行R波定位前的匹配滤波中， 对模板的要求比较高， 混合 信号与模板相减要求有较高的精度， 因此在选取模板、 阈值、 相位、 幅度值修正上需要采取 多种措施以保证分离后的信号达到。
7、一定准确性， 从而使得整个过程比较复杂， 算法存在较 多冗余。 0008 2)基于胎儿心音模式的自相关方法计算胎儿心音心率。 自相关方法是从加性噪声 中提取已知滤波， 特别是对于确定性的周期信号时一种好方法， 但是自相关提取胎儿心率 信号， 效果不够好， 主要是因为胎儿心音信号具有周期性， 这是不确定性的周期信号， 并且 胎儿心音信号波形也是受到胎儿本身以及周边环境的影响， 使得周期不是固定不变的， 因 此周期与波形均带有随机性变化造成的自相关函数波形难以检测。 0009 因此， 对胎心心音瞬时心率快速并准确的检测识别， 成为急需解决的技术问题。 发明内容 0010 本发明的主要目的是提出一种。
8、胎儿心音瞬时心率检测识别方法及装置， 旨在能快 速便捷并准确地检测识别出胎儿心音瞬时心率， 更好地服务于临床胎儿监护。 0011 为实现上述目的， 本发明提出的一种胎儿心音瞬时心率检测识别方法， 包括以下 步骤： 说 明 书 1/4 页 3 CN 106491085 A 3 0012 S1.通过仪器对胎儿心音信号进行采集并作为原始混合信号； 0013 S2.对采集到的所述原始混合信号进行滤波处理， 并作归一化处理， 而得到胎儿心 音信号作为预处理信号； 0014 S3.根据所述预处理信号， 提取所述胎儿心音信号的特征值； 0015 S4.将提取出来的所述特征值建立并训练混合高斯模型； 0016。
9、 S5.通过训练完毕的所述混合高斯模型对待测定的所述胎儿心音信号进行识别， 最终得出所述胎儿心音信号的瞬时心率。 0017 优选地， 所述步骤S4包括以下步骤： 0018 S41.对所述混合高斯模型进行定义； 0019 S42.利用最大似然法对所述混合高斯模型的各个变量的参数进行估计； 0020 S43.直至所述混合高斯模型的参数收敛为止。 0021 本发明还提出一种胎儿心音瞬时心率检测识别装置， 所述胎儿心音瞬时心率检测 识别装置使用胎儿心音瞬时心率检测识别方法。 0022 本发明技术方案首先通过仪器采集胎儿心音信号作为原始混合信号， 然后进行滤 波和归一化处理而得到胎儿心音信号作为预处理信。
10、号， 再通过提取胎儿心音信号的特征以 建立并训练混合高斯模型， 通过训练完毕的混合高斯模型对待测定的胎儿心音信号进行识 别， 最终得到胎儿心音瞬时心率。 本发明的技术方案相对于现有技术检测识别方法， 能够快 速便捷并准确地检测识别出胎儿心音瞬时心率， 更好地服务于临床胎儿监护。 附图说明 0023 为了更清楚地说明本发明实施例或现有技术中的技术方案， 下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍， 显而易见地， 下面描述中的附图仅仅是本 发明的一些实施例， 对于本领域普通技术人员来讲， 在不付出创造性劳动的前提下， 还可以 根据这些附图示出的结构获得其他的附图。 0024 图1。
11、为本发明胎儿心音瞬时心率检测识别方法原理图。 0025 本发明目的的实现、 功能特点及优点将结合实施例， 参照附图做进一步说明。 具体实施方式 0026 下面将结合本发明实施例中的附图， 对本发明实施例中的技术方案进行清楚、 完 整地描述， 显然， 所描述的实施例仅仅是本发明的一部分实施例， 而不是全部的实施例。 基 于本发明中的实施例， 本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其 他实施例， 都属于本发明保护的范围。 0027 需要说明， 若本发明实施例中有涉及方向性指示(诸如上、 下、 左、 右、 前、 后)， 则该方向性指示仅用于解释在某一特定姿态(如附图所示)下各部件之。
12、间的相对位置关系、 运动情况等， 如果该特定姿态发生改变时， 则该方向性指示也相应地随之改变。 0028 另外， 若本发明实施例中有涉及 “第一” 、“第二” 等的描述， 则该 “第一” 、“第二” 等 的描述仅用于描述目的， 而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技 术特征的数量。 由此， 限定有 “第一” 、“第二” 的特征可以明示或者隐含地包括至少一个该特 征。 另外， 各个实施例之间的技术方案可以相互结合， 但是必须是以本领域普通技术人员能 说 明 书 2/4 页 4 CN 106491085 A 4 够实现为基础， 当技术方案的结合出现相互矛盾或无法实现时应当认为这种。
13、技术方案的结 合不存在， 也不在本发明要求的保护范围之内。 0029 参见图1， 本发明提出一种胎儿心音瞬时心率检测识别方法， 包括以下步骤： 0030 第一步， 首先通过仪器对胎儿心音信号进行采集。 通过常规技术手段和硬件， 如胎 心仪等设备装置对胎儿进行检测， 并采集到相应信号作为原始混合信号。 因为对胎儿心音 信号检测过程中， 通常包含母亲心音信号、 其他设备等工频干扰等信号， 因此可设定S(n)表 示原始混合信号， y(n)表示胎儿心音信号， x(n)表示干扰信号， 包括母亲心音信号以及工频 干扰等信号， 最终公式表示为： 0031 S(n)y(n)+x(n) 0032 第二步， 然后。
14、对采集的原始混合信号进行滤波， 归一化处理得到胎儿心音信号作 为预处理信号。 通过自适应滤波技术， 对原始混合信号S(n)进行数字滤波处理， 去除干扰信 号x(n)， 从而可以得到比较单纯的胎儿心音信号y(n)， 然后采用巴特沃斯二阶低通滤波器 得到胎儿心音的预处理信号。 优选地， 本发明中巴特沃斯二阶低通滤波器的截止频率设置 为fc881hz。 0033 第三步， 根据预处理信号， 提取胎儿特征值序列d(n)。 采用Mel差分倒谱系数法提 取胎儿心音信号的特征值， 采用上述方法不仅可以反映胎儿心音信号的静态特性， 而且在 动态变化过程中得到很好的体现， 其中Mel差分倒谱次数的计算公式如下：。
15、 0034 0035 其中， 式中的参数c， d表示一帧胎儿心音信号序列， k为常数。 通过上述公式求出的 Mel差分倒谱系数d(n)即为胎儿心音信号的特征值， 用于后续的模式训练和识别。 0036 第四步， 然后将上述步骤公式求出的Mel差分倒谱系数d(n)即为胎儿心音信号的 特征值， 建立混合高斯模型并进行相应混合高斯模型的训练， 包括以下步骤： 0037 第一， 对混合高斯模型进行定义， 公式为： 0038 0039 式中的参数k为模型的个数， k为第k个混合高斯模型的权重， p(x|k)为第k个混合 高斯模型概率密度， 其均值为 k， 方差为 k。 对此概率密度相当于要求出 k、 k、。
16、 k各个变量。 当求出P(x)的表达式之后， 求和时的各项的结果就分别代表样本x属于各个变量的概率。 0040 第二， 利用最大似然法估计混合高斯模型参数 k、 k、 k各个变量的参数。 其中， 最 大似然法就是使样本点在估计的概率密度函数上的概率值最大， 由于概率值一般都很少， 当N很大的时候， 连乘的结果非常小， 容易造成浮点数下溢， 所以通常采取对数取值方式， 将 目标改写成： 0041 0042 上式的求解时， 通常采用EM求解算法， 将求解过程分解为两个步骤。 第一步， 假设 已经知道混合高斯模型的各个参数， 然后去估计各个参数在混合高斯模型的权值； 第二步， 基于估算的权值， 然后。
17、去确定混合高斯模型的各个参数的数值。 通过重复上述的两个步骤， 使得混合高斯模型的各个参数的波动逐渐变小， 以近似达到极值， 具体步骤如下： 说 明 书 3/4 页 5 CN 106491085 A 5 0043 1)E步骤： 0044 对于第i个样本xi， 由第k个mod el生成的概率为： 0045 0046 在这一步骤中， 假设混合高斯模型的参数为已知， 如由上述迭代或由初始值决定。 0047 2)M步骤： 0048 由上述的E步骤得到每个点的i(k)后， 对于样本xi的i(k)xi的值是由第k个混 合高斯模型的参数， 用i(k)xi(i1N)的数据去作相应的参数估计， 公式如下： 00。
18、49 0050 0051 0052 3)通过重复E步骤和M步骤， 直至所述混合高斯模型的参数收敛为止。 0053 第五步， 通过训练好的混合高斯模型来识别待测的胎儿心音信号属于哪一类的心 音， 从而根据判断得到心音信号的瞬时心率而作为待测胎儿心音的瞬时心率。 0054 本发明还提出一种胎儿心音瞬时心率检测识别装置， 该胎儿心音瞬时心率检测识 别装置使用了上述一种胎儿心音瞬时心率检测识别方法。 0055 以上所述仅为本发明的优选实施例， 并非因此限制本发明的专利范围， 凡是在本 发明的构思下， 利用本发明说明书及附图内容所作的等效结构变换， 或直接/间接运用在其 他相关的技术领域均包括在本发明的专利保护范围内。 说 明 书 4/4 页 6 CN 106491085 A 6 图1 说 明 书 附 图 1/1 页 7 CN 106491085 A 7 。