[发明专利]基于时频分析的多普勒胎心音的瞬时心率检测方法及装置

说明书

技术领域

本发明涉及医用监护技术领域，更具体地，涉及一种基于时频分析的多普勒胎心音的瞬时心率检测方法及装置。

背景技术

胎心监护是通过监测胎心率的变化来评价宫内安危的一种常用的胎儿监护方法。在围产阶段对母体内的胎儿进行监护，可以在了解胎儿发育状况的同时，大大地减少因缺氧、缺血而出现的窘迫，降低胎儿出生缺陷，甚至胎儿死亡的情况的发生。出生缺陷目前已成为影响我国人口质量的一个严重问题。因此，密切关注胎儿心率的变化情况对提高生育质量有着重要的意义。早至19世纪初就有产科医师通过胎心听诊估计胎儿宫内状况，随着超声多普勒技术的发展，产时电子胎心监护（electronicfetal monitoring，EFM）已成为目前最常用的胎儿监护方法。超声多普勒测量方法又是目前最常用的胎儿心率测量方法。

但是，由于从超声波换能器上检测的超声多普勒声音混杂有幅度大，分布广的噪声干扰，以及胎儿身体在母体内运动时，声音信号强弱变化明显。在时域和频域上，这些干扰信号混叠在一起对胎儿心音信号的瞬时心率计算造成很大的影响。因此，研究如何准确、有效地测量母体内胎儿心音的瞬时心率具有重要的理论价值和临床意义。

对胎心的监护和胎儿瞬时心率的研究在很早就已经开始了，并且采用的处理方法有很多种，总结起来，可分为以下几种：

（1）基于匹配滤波法计算胎儿心率：该方法的基本思想是用以往多次得到的母体心电作为模板，抵消腹部采样信号中的母体心电成分，将胎儿心电分离出来。由于腹部信号与模板相减要求有较高的精度，则需在模板计算、相位、幅值修正上采取多种措施来保证减去母体心电的准确性，这是基于心电模式的方法。

（2）采用自相关的方法计算胎儿心率；总所周知，相关方法是从加性噪声中提取已知波形，特别是对于确定性的周期信号是一种好方法。用自相关方法提取胎儿心率信号，效果不够好，主要是因为胎儿心率是重复性信号，但不是确定性的周期信号，并且胎儿心音信号的波形也不是固定不变的，周期与波形均带有随机性变化，因而造成自相关函数的波形不易检测；这是基于心音模式的方法。

正常心脏有四个心音：第一(S1)、第二(S2)、第三(S3)、第四(S4)心音，但多数情况下只能听到第一和第二心音。第一心音的出现标志着心室收缩期的开始；第二心音的出现标志心室缩张期的开始。从出现第一心音到出现第二心音的期间定义为心室收缩期；从出现第二心音到下一心动周期中出现第一心音的期间定义为心室舒张期。在一个心动周期中，心音的主要成分为：第一心音、收缩期、第二心音、舒张期，通过它们可以完整的表述出心音的时间特性。对于一般人来说，正常情况下心脏的收缩期短于舒张期。胎儿的平均心率为120-160次/分，心动周期约为0.5(秒)，其中收缩期大约占0.2(秒)，舒张期大约占0.3(秒)。也就是说，一般人的心音信号中，S1音到S2音在时域上的间隔大约为0.2(秒)。

由于胎儿心音信号属于非平稳信号，常规傅里叶变换方法不能刻画任一时刻的频率成分，致使无法对其进行全面的分析。时频分析是非平稳信号分析的一个有力工具，该方法能够将一维时域信号变换到二维的时频平面，提供了时间域和频率域的联合分布信息，清楚的描述了信号频率随时间变换的关系。

发明内容

本发明的主要目的在于克服现有胎儿心率检测方法的不足，提供一种基于时频分析的多普勒胎心音的瞬时心率检测方法，该检测方法将胎儿心音的时间域和频率域的分布信息联合起来，再利用了心音信号的先验信息——心音信号中S1音和S2音在观测信号时域上的间隔，从而检测胎心音的瞬时心率。

为解决上述技术问题，本发明的技术方案如下：

一种基于时频分析的多普勒胎心音的瞬时心率检测方法，包括如下步骤：

S1.信号预处理：对采集到的多普勒胎心音进行带通滤波，带通滤波器的通带频率为f_L到f_H；

S2.时频分析：对于经过步骤S1预处理的多普勒胎心音进行时频分析，得到时频图；

S3.选取特征频段和模板：从时频图中选取信号中的特征频段，特征频段为f_cL到f_cH；从时频图选取包含S1音和S2音特征的时频块，所述时频块的时间间隔为0.2(秒)＜t₀＜0.5(秒)；

S4.求互相关函数：对特征频段和模板求互相关函数，根据互相关函数的结果绘出相关曲线；

S5.用阈值检测法求互相关曲线的峰值；

S6.对峰值做差分，求瞬时心率值；根据瞬时心率值绘制出瞬时心率图。