一种时域自适应加窗的胎心音降噪方法技术

本发明专利技术提供了一种时域自适应加窗的胎心音降噪方法。该方法通过对胎心音数据的包络进行自适应加窗，保留了包络上有用胎心音信号，剔除了大部分干扰信号。该方法可抑制与心音信号同频带的噪声，提高信噪比，并且可补偿强噪声背景下被淹没的微弱心音信号。当噪声干扰造成心率曲线断续时，该降噪方法可使心率曲线更加连续。

A time domain adaptive windowing algorithm for reducing fetal heart sounds

The present invention provides a time domain adaptive windowing method for reducing sound of fetal heart sound. In this method, an adaptive window is added to the envelope of the fetal heart sound data, which preserves the useful fetal heart sound signals and eliminates most of the interfering signals. The method can suppress the noise of the same frequency band as the heart sound signal, improve the signal to noise ratio, and compensate the faint heart sound signal under the strong noise background. The noise reduction method can make the heart rate curve more continuous when the noise disturbance causes the heart rate curve to be interrupted.

【技术实现步骤摘要】

本专利技术涉及音频信号处理领域，尤其涉及一种针对微弱胎心音的时域降噪方法。
技术介绍
胎心率是胎儿健康监护的重要指标，胎心率曲线直观地显示了胎儿的心脏功能，并且与胎儿的供氧状况和脐带缠绕等信息紧密相关。目前，最常规的胎心率检测方法是采用超声波胎心仪，超声波胎心仪主动发射超声波作用于胎儿心脏，利用超声波的多普勒效应，经过信号处理计算出胎心率。超声波存在安全剂量问题，医学上不建议高频率、长时间地使用超声波胎心仪，这导致孕妇无法随时监控胎儿的健康状况。针对这种问题，更加理想的胎儿健康监护手段是采用无源式收集胎心音的检测装备，这种检测装备不发射超声波，被动式接收胎儿的心跳声信号，对胎儿和孕妇都十分安全，可长时间连续使用。无源式胎心音检测装备采用高灵敏度的声信号采集模块，在孕妇腹部采集胎心音，然后对胎心音做信号处理，计算胎心率并绘制曲线图。胎心音是一种微弱的低频声信号，易受环境噪声、母体内噪声等干扰，采集的胎心音数据中不可避免地含有各种低频噪声干扰，严重的噪声干扰会影响心率值的计算。在实现本专利技术的过程中，申请人发现常规的频域滤波方法可滤除非胎心音频段的噪声干扰，在一定程度上提高信噪比，却无法滤除与胎心音信号同频段的噪声，同频段噪声干扰很严重时，无法计算出心率值，导致心率曲线断断续续。此时需要一种可滤除与心音信号同频段噪声的方法，进一步提高信噪比，改善胎心率曲线断断续续现象。
技术实现思路
(一)要解决的技术问题鉴于上述技术问题，本专利技术提供了一种时域自适应加窗的胎心音降噪方法，以实现提高信噪比并改善胎心率曲线断断续续现象的效果。(二)技术方案本专利技术一种时域自适应加窗的胎心音降噪方法包括：步骤A，对第i秒原始胎心音数据di(n)进行滤波，滤除胎心音频段外的干扰信号，得到胎心音滤波数据fi(n)，其中，n＝1，2，…，Ⅳ，Ⅳ为1秒总采样点数；步骤B，对滤波数据fi(n)取包络，得到包络数据ei(n)；步骤C，根据前M秒胎心率判定是否开启时域自适应加窗，若开启加窗，则进入步骤D，若不开启加窗，则跳至步骤E；步骤D，时域自适应加窗具体包括：子步骤D1：根据上1秒胎心率值，生成当前1秒两个矩形窗之间的参考间距L及间距的变化长度dH和dL；子步骤D2：设置宽度为w的矩形窗，在当前1秒的包络上移动，移动步长为v，计算每个起始位置为lk处矩形窗内包络的总能量pk，其中，k＝1，2，…，K，K为总位置点数；子步骤D3：将当前1秒包络数据划分成Q段，从pk中找出每段中能量最大值pq及其相应位置lq，q＝1，2，…，Q；子步骤D4：从Q段最大值pq中找出满足条件的两个最大值pmax1和pmax2，依次记录其位置lmax1和lmax2；子步骤D5：以lmax1和lmax2为加窗的起始点，分别在包络上加宽度为w的矩形窗，将在[lmax1+w，lmax2]范围内的包络置零，以及lmax1位置前和在lmax2+w位置后的lmax2-lmax1-w距离内的包络置零。步骤E，根据包络峰值间隔计算当前1秒心率值。(三)有益效果从上述技术方案可以看出，本专利技术一种时域自适应加窗的心音降噪方法具有以下有益效果：(1)抑制与心音信号同频带的噪声，提高信噪比；(2)补偿强噪声背景下被淹没的微弱心音信号，使心率曲线更加连续。附图说明图1为根据本专利技术实施例时域自适应加窗的胎心音降噪方法的流程图。图2为实测胎心音数据1的不加窗处理的心率曲线与自适应加窗处理后的心率曲线对比；图3为实测胎心音数据2的不加窗处理的心率曲线与自适应加窗处理后的心率曲线对比；具体实施方式为使本专利技术的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本专利技术进一步详细说明。需要说明的是，在附图或说明书描述中，相似或相同的部分都使用相同的图号。附图中未绘示或描述的实现方式，为所属
中普通技术人员所知的形式。另外，虽然本文可提供包含特定值的参数的示范，但应了解，参数无需确切等于相应的值，而是可在可接受的误差容限或设计约束内近似于相应的值。实施例中提到的方向用语，例如“上”、“下”、“前”、“后”、“左”、“右”等，仅是参考附图的方向。因此，使用的方向用语是用来说明并非用来限制本专利技术的保护范围。本专利技术一种时域自适应加窗的胎心音降噪方法，该方法假定胎心音包络数据上至少存在两次心跳产生的两个心音峰，在胎心率连续的情况下开启自适应加窗，通过一些列筛选和判决，找到两个心音包络峰，在心音包络峰上加矩形窗，两个矩形窗中间包络及矩形窗外侧一定范围内的包络置零，从而保留包络上有用胎心音信号，剔除大部分干扰信号。在本专利技术的一个示例性实施例中，提供了对实测胎心音进行时域自适应加窗的胎心音降噪方法的演示。图1是根据本专利技术实施例时域自适应加窗的胎心音降噪方法的流程图。请参照图1，本实施例包括：步骤A，对第i秒原始胎心音数据di(n)进行滤波，滤除胎心音频段外的干扰信号，得到胎心音滤波数据fi(n)，其中，n＝1，2，…，N，Ⅳ为1秒总采样点数；步骤B，对滤波数据fi(n)取包络，得到包络数据ei(n)；步骤C，判断前M秒是否连续计算出在正常范围内的胎心率，若是，则开启时域自适应加窗，进入步骤D，若不是，则不开启加窗，跳至步骤E：步骤D，在包络数据ei(n)进行自适应加窗，输出新包络数据，具体包括：子步骤D1：根据上1秒胎心率值，生成当前1秒两个矩形窗之间的参考间距L及间距的变化长度dH和dL；子步骤D1中，参考间距L、间距的变化长度dH和dL可通过以下公式计算：其中，Ⅳ为1秒总采样点数，f0为上一秒胎心率，胎心率为f0表示每分钟f0次心跳，只为胎心率变化值。这里需要注意的是，心跳的频率一般不会突变，胎心率变化值不宜设置过大。子步骤D2：设置宽度为w的矩形窗，在当前1秒的包络上移动，移动步长为v，计算每个起始位置为lk处矩形窗内包络的总能量pk，其中，k＝1，2，…，K，K为总位置点数；子步骤D2中需要注意的是，建议矩形窗的宽度w小于1/2L，若矩形窗宽度太宽，降噪效果可能不佳，同时，移动步长为v不宜过大。子步骤D3：将当前1秒包络数据划分成Q段，从pk中找出每段中能量最大值pq及其相应位置lq，q＝1，2，…，Q；子步骤D4：从Q段最大值pq中找出满足条件的两个最大值pmax1和pmax2，依次记录其位置lmax1和lmax2；子步骤D4选择两个最大值的必需满足以下三个条件：条件a：lmax2-lmax1的值在[L-dL，L+dH]范围内；条件b：lmax1在设定的起始范围内；条件c：pmax1与pmax2的比值不超过2。子步骤D5：以lmax1和lmax2为加窗的起始点，分别在包络上加宽度为w的矩形窗，将在[lmax1+w，lmax2]范围内的包络置零，以及lmax1位置前和在lmax2+w位置后的lmax2-lmax1-w距离内的包络置零。子步骤D5中需要注意的是，如当前1秒包络上存在三个心跳包络峰，在对其中两个包络峰加矩形窗后，并对矩形窗之间及两侧一定范围内包络置零后，第三个心跳包络峰仍然保留了下来，且其附近的噪声也得到一定程度抑制。步骤E，根据包络峰值间隔计算当前1秒心率值。图2至图3给出了两组实测胎心音数据的不加窗处理的心率曲线与自适应加窗处理后的心率曲线对比图，可以看出，由于噪本文档来自技高网...

【技术保护点】
一种时域自适应加窗的胎心音降噪方法，其特征在于，包括：步骤A，对第i秒原始胎心音数据di(n)进行滤波，滤除胎心音频段外的干扰信号，得到胎心音滤波数据fi(n)，其中，n＝1，2，…，Ⅳ，Ⅳ为1秒总采样点数；步骤B，对滤波数据fi(n)取包络，得到包络数据ei(n)；步骤C，根据前M秒胎心率值判定是否开启时域自适应加窗，若开启加窗，则进入步骤D，若不开启加窗，则跳至步骤E；步骤D，在包络数据ei(n)上进行自适应加窗，输出新包络数据：步骤E，根据包络峰值间隔计算当前1秒心率值。

【技术特征摘要】
1.一种时域自适应加窗的胎心音降噪方法，其特征在于，包括：步骤A，对第i秒原始胎心音数据di(n)进行滤波，滤除胎心音频段外的干扰信号，得到胎心音滤波数据fi(n)，其中，n＝1，2，…，Ⅳ，Ⅳ为1秒总采样点数；步骤B，对滤波数据fi(n)取包络，得到包络数据ei(n)；步骤C，根据前M秒胎心率值判定是否开启时域自适应加窗，若开启加窗，则进入步骤D，若不开启加窗，则跳至步骤E；步骤D，在包络数据ei(n)上进行自适应加窗，输出新包络数据：步骤E，根据包络峰值间隔计算当前1秒心率值。2.根据权利要求1所述的方法，其特征在于，所述步骤C开启自适应加窗的条件是：前M秒连续计算出在正常范围内的胎心率。3.根据权利要求1所述的方法，其特征在于，所述步骤D具体包括：子步骤D1：根据上1秒胎心率值，生成当前1秒两个矩形窗之间的参考间距L及间距的变化长度dH和dL；子步骤D2：令矩形窗宽度为w，在当前1秒的包络上移动，移动步长为v，计算每个起始位置为lk处的矩形窗内包络的总能量pk，其中，k＝1，2，…，K，K为总位置点数；子步骤D3：将当前1秒包络数据划分成Q段，从pk中找出每段中能量最大值pq及其相应位置...

【专利技术属性】
技术研发人员：周营，饶丽婷，
申请(专利权)人：苏州萌动医疗科技有限公司，
类型：发明
国别省市：江苏;32