CN110491504A - 一种心音信号医学指标数据的获取方法 - Google Patents

Abstract

本发明公开了一种心音信号医学指标数据的获取方法，使用双阈值限定来进行心音信号分段，解决了单一阈值分段存在的抗干扰能力差的问题，提高了心音信号分段的准确性，为后续可以准确的提取出心音信号的医学指标数据提供了保障。同时，在分段前对采集的心音信号进行尖峰信号消除，采用取对数，用线性的低通滤波器进行低通滤波得获得一个非常平滑的心音信号包络，这在对心音信号进行分段，确定心音中的各项成分是非常有帮助的。最后，根据心音信号的分段，得到心音信号的第一心音、第二心音、收缩期以及舒张期的持续时间，然后根据这些持续时间，就可以提取出心音信号的各项医学指标数据用作后续诊断使用。

Description

一种心音信号医学指标数据的获取方法

技术领域

本发明属于医疗测试技术领域，特别是，更为具体地讲，涉及一种心音信号医学指标数据的获取方法。

背景技术

心音信号是由于心脏有规律的机械振动产生的，经由心胸传导系统传导到人体皮肤表面的声音信号。对于心音信号的听诊(简称心音听诊)可以诊断出心脏是否存在某些异常情况，它较心电更早的发现心脏问题，从而把握住早期的治疗时机。同时心音听诊具有便捷、无创、经济性好的特点，所以心音听诊在临床诊断中扮演着相当重要的角色。

但传统的心音听诊同样存在着一定的局限性，诊断的准确性往往取决于听诊大夫的医疗知识水平、经验丰富程度，存在很多的主观因素。因此，实现心音信号的自动化诊断，可以排除不稳定的主观因素，更加客观地获取心音信号医学指标数据，对及时发现心脏是否存在问题具有深刻的研究意义。

识别心音信号中第一心音、第二心音以及收缩期、舒张期的确切位置，对心音进行分段是实现心音信号自动化诊断的必要前提，同时也使得利用计算机对正常/异常心音进行分类成为了可能。

对于心音信号的分段，传统、经典的方法是对心音信号进行预处理，获取心音的某一种包络，再基于单一阈值来对心音信号进行分段。虽然易于实现，但是准确性不高，抗干扰能力不强。当心音信号中某一个噪声成分大于所选定的阈值时，往往会出现分段错误的情况。

发明内容

本发明的目的在于克服现有技术的不足，提出一种心音信号医学指标数据的获取方法，以提高分段抗干扰能力。

为实现上述发明目的，本发明心音信号医学指标数据的获取方法，其特征在于，包括以下步骤：

(1)、心音信号采集与预处理

首先使用电子听诊器以采样频率F_s采集人体心音信号，得到离散的心音信号x₀(n)，然后对所采集到的离散心音信号x₀(n)进行预处理：

1.1)、分别使用截止频率为25Hz四阶巴特沃斯高通滤波器和截止频率为400Hz四阶巴特沃斯低通滤波器，对离散心音信号x₀(n)进行滤波处理，得到滤波后的心音信号x(n)；

1.2)、对滤波后的心音信号x(n)进行尖峰消除处理，消除心音信号x(n)中的尖峰；

1.3)、对消除尖峰的心音信号x(n)进行归一化处理，得到归一化的心音信号x(n)；

(2)、心音信号包络特征提取

对步骤(1)获得的归一化心音信号x(n)的绝对值取对数，得到心音信号ln|x(n)|，并通过一个线性的低通滤波器对心音信号ln|x(n)|进行低通滤波，得到除去心音信号中高频成分的信号L[ln|x(n)|]，将信号L[ln|x(n)|]作为e的指数，得到心音信号的同态包络(简称心音信号包络)Env_homomorphic(n)；

(3)、心音信号分段

3.1)、设定两个阈值thr₁与thr₂，其中：

thr₁＝para₁*max(Env_homomorphic(n))，para₁∈[0.1，0.4]

thr₂＝para₂*mean(Env_homomorphic(n))，para₂∈[0，0.1]

max()函数为求最大值函数，mean()函数为求平均值函数；

3.2)、对步骤(2)提取到的心音信号包络Env_homomorphic(n)，按照时间顺序找出所有幅值大于阈值thr₁包络段的第一个点，并依次存入数组LargerThanThr₁中；

3.3)、从数组LargerThanThr₁中依次逐个取点，从该点前后两个方向，心音信号包络Env_homomorphic(n)找出小于thr₂的第一个点，前向找到的点依次存入数组thr2end，后向找到的点依次存入数组thr2start；

3.4)、计算时间间隔T1、T2：

T1＝thr2start[2]-thr2end[1]

T2＝thr2start[3]-thr2end[2]

其中，thr2start[2]为数组thr2start的第二个点对应的时刻点，thr2start[3]为数组thr2start的第二个点对应的时刻点，thr2end[1]为数组thr2end的第一个点对应的时刻点，thr2end[2]为数组thr2end的第二个点对应的时刻点；

如果T1<T2，则数组thr2start的奇数项为第一心音的起始点，偶数项为第二心音的起始点，数组thr2end数组奇数项为第一心音的终止点，偶数项为第二心音的终止点；

如果T1>T2，则数组thr2start的奇数项为第二心音的起始点，偶数项为第一心音的起始点，数组thr2end奇数项为第二心音的终止点，偶数项为第一心音的终止点；

3.5)、记录下第一、第二心音的起始与终止点，并根据第一、第二心音的起始与终止点得到心音信号中第一心音、第二心音以及收缩期、舒张期的确切位置，从而完成心音信号的分段工作；

(4)、心音信号医学指标数据提取

根据心音信号的分段，得到心音信号的第一心音、第二心音、收缩期以及舒张期的持续时间。

本发明的目的是这样实现的。

本发明心音信号医学指标数据的获取方法，使用双阈值限定来进行心音信号分段，解决了单一阈值分段存在的抗干扰能力差的问题，提高了心音信号分段的准确性，为后续可以准确的提取出心音信号的医学指标数据提供了保障。同时，在分段前对采集的心音信号进行尖峰信号消除，采用取对数，用线性的低通滤波器进行低通滤波得获得一个非常平滑的心音信号包络，这在对心音信号进行分段，确定心音中的各项成分是非常有帮助的。最后，根据心音信号的分段，得到心音信号的第一心音、第二心音、收缩期以及舒张期的持续时间，然后根据这些持续时间，就可以提取出心音信号的各项医学指标数据用作后续诊断使用。如为了实现心音信号的自动辅助诊断，我们可以从测试心音样本信号提取出的医学指标与健康、正常的心音医学指标范围进行比对，找出该心音信号是否正常以及可能存在异常的地方，给出参考的辅助诊断结果。本发明能够辅助医护人员或被测用户对于心音的健康情况作出初步的评估，方便、及时、快速的找出心音可能存在的问题，具有较高的准确率。

附图说明

图1是本发明心音信号医学指标数据的获取方法一种具体实施方式流程图；

图2是一个心音信号提取同态包络前后对照图

图3是心音信号的分段结果图；

图4某心音信号得出的诊断结果图。

具体实施方式

下面结合附图对本发明的具体实施方式进行描述，以便本领域的技术人员更好地理解本发明。需要特别提醒注意的是，在以下的描述中，当已知功能和设计的详细描述也许会淡化本发明的主要内容时，这些描述在这里将被忽略。

图1是本发明心音信号医学指标数据的获取方法一种具体实施方式流程图。

在本实施例中，如图1所述，本发明心音信号医学指标数据的获取方法包括以下步骤：

一、心音信号采集与预处理S1

首先使用电子听诊器以采样频率F_s采集人体心音信号，得到离散的心音信号x₀(n)，然后对所采集到的离散心音信号x₀(n)进行预处理。

1、巴特沃斯带通滤波

分别使用截止频率为25Hz四阶巴特沃斯高通滤波器和截止频率为400Hz四阶巴特沃斯低通滤波器，对采集的离散心音信号x₀(n)进行滤波处理，得到滤波后的心音信号x(n)。

其中，截止频率为25Hz的四阶巴特沃斯高通滤波器传递函数为：

截止频率为400Hz的四阶巴特沃斯低通滤波器传递函数为：

2、尖峰消除算法

对滤波后的心音信号x(n)进行尖峰消除处理，在本实施例中，具体尖峰消除方法如下：

2.1以500ms的时长将心音信号x(n)分割成若干个心音片段，则每一片段的心音信号可表示为x_m(n)，m为心音片段序列号；

2.2找出各心音片段中的最大绝对幅度(即幅值的绝对值最大的点)，记作f_m(t)：

f_m(t)＝max[|x_m(n)|]；

其中，m表示的是第m个心音片段，t表示的是该心音片段的最大绝对幅度对应的横坐标即时刻；

2.3如果存在某一即第i个心音片段的最大绝对幅度f_i(t)，大于其他各片段心音信号最大绝对幅度的中位数的三倍，则执行以下子步骤：

2.3.1选定包含f_i(t)的第i段心音片段x_i(n)；

2.3.2在心音片段x_i(n)中，最大绝对幅度f_i(t)的位置坐标(t，f_i(t))被认为是噪声尖峰的顶点；

2.3.3将噪声尖峰的起始点N_start定义为在最大绝对幅度点(t，f_i(t))之前的最后一个过零点(N_start，0)；

2.3.4将噪声尖峰的终止点N_end定义为在最大绝对幅度点(t，f_i(t))之后的第一个过零点(N_end，0)；

2.3.5将噪声尖峰的起止范围内视作一个完整的噪声尖峰，并将范围内的幅值全部置0，即x(N_start:N_end)＝0；

2.4返回继续处理，直到大于其他各片段心音信号最大绝对幅度的中位数的三倍心音片段全部处理完毕。

3、归一化处理

由于采样得到的心音信号强度大小不一，不利于后期对心音信号的研究分析，故采用归一化的方法，来解决心音信号强度不一致的问题，将心音信号的强度转化到[-1，1]的范围内。归一化后的心音信号可以表示为下式：

其中，＝表示赋值，即归一化后再赋值给x(n)；

二、心音信号包络特征提取S2

心音信号的包络可以更好地反映出心音信号的特征，因此在对心音信号进行上述过程的预处理之后，就可以通过提取心音信号包络特征，以便对心音信号进行更加充分的研究和分析。

提取心音信号的同态包络是一种有效的获取心音图幅度包络的方法。同态包络的优点在于它可以获得一个非常平滑的心音信号的包络，这在对心音信号进行分段，确定心音中的各项成分是非常有帮助的。

在本实施例中，同态包络的提取是通过同态滤波器完成的。将心音信号x(n)看作是一个缓慢变化的控制信号幅值的成分a(n)和一个快速变化的代表信号振荡部分的成分o(n)这两部分乘积的产物，故心音信号可以这样来表示：

x(n)＝a(n)·o(n) a(n)＞0；

为了将表达式从乘号转换成加号，对心音信号x(n)的绝对值取对数，即：

ln|x(n)|＝ln|a(n)|+ln|o(n)|；

因此心音信号x(n)由之前的乘法变成了加法，并且可以通过一个线性的低通滤波器对心音信号ln|x(n)|进行低通滤波，得到除去心音信号中高频成分的信号L[ln|x(n)|]，即：

L(ln|x(n)|)＝L(ln|a(n)|)+L(ln|o(n)|)≈L(ln|a(n)|)

其中，L代表的是一个低通滤波器，这里使用的一阶巴特沃斯低通滤波器，截止频率是8Hz。经过低通滤波之后，将信号L[ln|x(n)|]作为e的指数，得到心音信号的同态包络即心音信号包络)Env_homomorphic(n)：

Env_homomorphic(n)＝exp(L(ln|x(n)|))≈a(n)

最终a(n)即为一个心音信号的同态包络，如图2所示。上面一张为原始即预处理后的心音信号，其幅度为相对值，下面一张为心音信号包络，其幅度也为相对值。从图2可以看出，本发明很好地提取了心音信号的同态包络。

三、心音信号分段S3

对于提取到的心音信号包络，一种比较传统经典的方法是基于单阈值的心音信号分段方法。但采用单阈值进行分段的方法存在一个缺点就是抗干扰噪声能力差，当心音信号中某一个噪声成分大于所选定的阈值时，单阈值分段的方法会将干扰噪声按正常心音信号进行处理，从而造成心音信号分段出现错误。

本发明采用一种基于双阈值的心音信号分段方法，即设定两个阈值，其中一个阈值用于消除干扰噪声，另一个阈值用于对心音进行分段。具体分段方法的步骤如下：

1、设定两个阈值thr₁与thr₂，其中：

thr₁＝para₁*max(Env_homomorphic(n))，para₁∈[0.1，0.4]

thr₂＝para₂*mean(Env_homomorphic(n))，para₂∈[0，0.1]

max()函数为求最大值函数，mean()函数为求平均值函数。

2、对提取到的心音信号包络Env_homomorphic(n)，按照时间顺序找出所有幅值大于阈值thr₁包络段的第一个点，并依次存入数组LargerThanThr₁中；

3、从数组LargerThanThr₁中依次逐个取点，从该点前后两个方向，心音信号包络Env_homomorphic(n)找出小于thr₂的第一个点，前向找到的点依次存入数组thr2end，后向找到的点依次存入数组thr2start；

4、由于心音的特点是舒张期的时间间隔大于收缩期时间间隔，基于这一特性，计算时间间隔T1、T2：

T1＝thr2start[2]-thr2end[1]

T2＝thr2start[3]-thr2end[2]

其中，thr2start[2]为数组thr2start的第二个点对应的时刻点，thr2start[3]为数组thr2start的第二个点对应的时刻点，thr2end[1]为数组thr2end的第一个点对应的时刻点，thr2end[2]为数组thr2end的第二个点对应的时刻点

如果T1<T2，则数组thr2start的奇数项为第一心音的起始点，偶数项为第二心音的起始点，数组thr2end数组奇数项为第一心音的终止点，偶数项为第二心音的终止点；

如果T1>T2，则数组thr2start的奇数项为第二心音的起始点，偶数项为第一心音的起始点，数组thr2end奇数项为第二心音的终止点，偶数项为第一心音的终止点；

5、记录下第一、第二心音的起始与终止点，并根据第一、第二心音的起始与终止点得到心音信号中第一心音、第二心音以及收缩期、舒张期的确切位置，从而完成心音信号的分段工作，心音信号的分段结果如图3所示，图3中标注的虚线是第一心音的起始和终止点，标注的实线是第二心音的起始和终止点。由图3可以看出，应用本发明提出的心音分段方法可以准确的划分第一心音和第二心音的起始终止点(位置)。

四、心音信号医学指标数据提取

在完成对一个心音信号进行分段的工作之后，就可以清楚的知道心音信号的第一心音、第二心音、收缩期以及舒张期的持续时间等各项数据，即可以提取出心音信号的各项医学指标数据用作后续诊断使用。

在本实施例中，提取出心音信号的5项医学指标数据，具体医学指标类别如下：

1、第一心音持续时间

经过分段后的心音信号可以得知每个心动周期的第一心音持续时间，求所有第一心音持续时间的平均值作为该心音信号的第一心音持续时间meanS1：

其中，k表示该心音信号中共有k个第一心音出现，tS1_j表示第j个第一心音所持续的时间。

2、第二心音持续时间

经过分段后的心音信号可以得知每个心动周期的第二心音持续时间，求所有第二心音的平均值作为该心音信号样本的第二心音持续时间meanS2：

其中，k表示该心音信号中共有k个第二心音出现，tS2_j表示第j个第一心音所持续的时间。

3、心动周期间隔时间：

心动周期间隔是指从第一心音开始，经理收缩期，第二心音，直到舒张期结束这样一个周期所持续的时间，因此心动周期间隔平均值meanS11可表示为：

其中，k表示该心音信号中共有k个心动周期，tS1_j表示第j个第一心音所持续的时间，tSys_j表示第j个收缩期所持续的时间，tS2_j表示第j个第二心音所持续的时间，tDia_j表示第j个舒张期所持续的时间。

4、第一二心音间隔

第一二心音间隔是指从第一心音发生开始到第二心音结束所持续的时间，因此第一二心音间隔平均值meanS12可表示为：

其中，k表示该心音信号中共有k个心动周期，tS1_j表示第j个第一心音所持续的时间，tSys_j表示第j个收缩期所持续的时间，tS2_j表示第j个第二心音所持续的时间。

5、心率

心率是指心脏每分钟跳动的次数。正常人在静息状态下心率的范围为60～100次每分钟，低于60次/分钟叫做心率过缓，高于100次/每分钟叫做心率过速。心率值hr的计算可用下式来完成：

五、基于心音医学指标的心音辅助诊断

相较于传统的临床听诊，基于本发明得到心音信号医学指标数据作为诊断心音是否正常的依据，蒋提高心音诊断的准确率。

对于一个获取到的心音信号样本，在具体诊断中，进行以下处理：

(1)心音信号预处理、包络提取，心音分段处理

通过对一个未知心音信号进行预处理、包络提取，心音分段处理，我们可以得知该心音的构成，为下面进一步诊断打好基础。

(2)心音医学指标数据提取

有了上一步对心音信号的正确分段作为前提，即可提取到该心音信号的第一心音、第二心音持续时间等5项医学指标，此步骤为现有技术。

(3)医学指标范围比对

针对提取到的某一心音信号的5项医学指标，与健康、正常心音的医学指标范围进行比对，判断心音的各项医学指标是否正常。

其中，第一心音持续时间正常值的范围为0.10秒～0.16秒，小于0.10秒认为第一心音狭窄，大于0.16秒认为第一心音肥厚；第二心音持续时间正常值的范围为0.08秒～0.12秒，小于0.08秒认为第二心音狭窄，大于0.12秒认为第二心音肥厚；心动周期间隔的正常值为0.50秒～1.20秒，不在此范围内则被认为心动周期间隔可能存在异常；一二心音间隔的正常值为0.30秒～0.50秒，不在此范围内则被认为一二心音间隔可能存在异常。

(4)得出诊断结论

根据比对的结果，给出参考的诊断结论，指出存在问题的地方，提醒被诊断者及时检查、尽快就医，诊断结果图如图4所示。

六、结论

基于本发明方法，可以较准确的完成对心音信号进行自动分段、正异常诊断工作，能够辅助医护人员或普通用户对于心音的健康情况作出初步的评估，一定程度上减轻了医疗工作者的工作量，同时也有利于心脏病患者及时发现自身的问题并及时进行相应的治疗工作，具有显著的社会价值。

尽管上面对本发明说明性的具体实施方式进行了描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。

Claims (3)

1.一种心音信号医学指标数据的获取方法，其特征在于，包括以下步骤：

(1)、心音信号采集与预处理

首先使用电子听诊器以采样频率F_s采集人体心音信号，得到离散的心音信号x₀(n)，然后对所采集到的离散心音信号x₀(n)进行预处理：

1.1)、分别使用截止频率为25Hz四阶巴特沃斯高通滤波器和截止频率为400Hz四阶巴特沃斯低通滤波器，对离散心音信号x₀(n)进行滤波处理，得到滤波后的心音x(n)；

1.2)、对滤波后的心音信号x(n)进行尖峰消除处理，消除心音信号x(n)中的尖峰；

1.3)、对消除尖峰的心音信号x(n)进行归一化处理，得到归一化的心音信号x(n)；

(2)、心音信号包络特征提取

对步骤(1)获得的归一化心音信号x(n)的绝对值取对数，得到心音信号ln|x(n)|，并通过一个线性的低通滤波器对心音信号ln|x(n)|进行低通滤波，得到除去心音信号中高频成分的信号L[ln|x(n)|]，将信号L[ln|x(n)|]作为e的指数，得到心音信号的同态包络(简称心音信号包络)Env_homomorphic(n)；

(3)、心音信号分段

3.1)、设定两个阈值thr₁与thr₂，其中：

thr₁＝para₁*max(Env_homomorphic(n)),para₁∈[0.1,0.4]

thr₂＝para₂*mean(Env_homomorphic(n)),para₂∈[0,0.1]

max()函数为求最大值函数，mean()函数为求平均值函数；

3.2)、对步骤(2)提取到的心音信号包络Env_homomorphic(n)，按照时间顺序找出所有幅值大于阈值thr₁包络段的第一个点，并依次存入数组LargerThanThr₁中；

3.3)、从数组LargerThanThr₁中依次逐个取点，从该点前后两个方向，心音信号包络Env_homomorphic(n)找出小于thr₂的第一个点，前向找到的点依次存入数组thr2end，后向找到的点依次存入数组thr2start；

3.4)、计算时间间隔T₁、T₂：

T1＝thr2start[2]-thr2end[1]

T2＝thr2start[3]-thr2end[2]

其中，thr2start[2]为数组thr2start的第二个点对应的时刻点，thr2start[3]为数组thr2start的第二个点对应的时刻点，thr2end[1]为数组thr2end的第一个点对应的时刻点，thr2end[2]为数组thr2end的第二个点对应的时刻点；

如果T1<T2，则数组thr2start的奇数项为第一心音的起始点，偶数项为第二心音的起始点，数组thr2end数组奇数项为第一心音的终止点，偶数项为第二心音的终止点；

如果T1>T2，则数组thr2start的奇数项为第二心音的起始点，偶数项为第一心音的起始点，数组thr2end奇数项为第二心音的终止点，偶数项为第一心音的终止点；

3.5)、记录下第一、第二心音的起始与终止点，并根据第一、第二心音的起始与终止点得到心音信号中第一心音、第二心音以及收缩期、舒张期的确切位置，从而完成心音信号的分段工作；

(4)、心音信号医学指标数据提取

根据心音信号的分段，得到心音信号的第一心音、第二心音、收缩期以及舒张期的持续时间。

2.根据权利要求1所述的获取方法，其特征在于，在步骤(1)中，所述的截止频率为25Hz的四阶巴特沃斯高通滤波器传递函数为：

截止频率为400Hz的四阶巴特沃斯低通滤波器传递函数为：

3.根据权利要求1所述的获取方法，其特征在于，在步骤(1)中，所述的尖峰消除为：

1.1)以500ms的时长将心音信号x(n)分割成若干个心音片段，则每一片段的心音信号可表示为x_m(n)，m为心音片段序列号；

1.2)找出各心音片段中的最大绝对幅度(即幅值的绝对值最大的点)，记作f_m(t)：

f_m(t)＝max[|x_m(n)|]；

其中，m表示的是第m个心音片段，t表示的是该心音片段的最大绝对幅度对应的横坐标即时刻；

1.3)如果存在某一即第i个心音片段的最大绝对幅度f_i(t)，大于其他各片段心音信号最大绝对幅度的中位数的三倍，则执行以下子步骤：

1.3.1)选定包含f_i(t)的第i段心音片段x_i(n)；

1.3.2)在心音片段x_i(n)中，最大绝对幅度f_i(t)的位置坐标(t,f_i(t))被认为是噪声尖峰的顶点；

1.3.3)将噪声尖峰的起始点N_start定义为在最大绝对幅度点(t,f_i(t))之前的最后一个过零点(N_start,0)；

1.3.4)将噪声尖峰的终止点N_end定义为在最大绝对幅度点(t，f_i(t))之后的第一个过零点(N_end，0)；

1.3.5)将噪声尖峰的起止范围内视作一个完整的噪声尖峰，并将范围内的幅值全部置0，即x(N_start:N_end)＝0；

1.4)返回继续处理，直到大于其他各片段心音信号最大绝对幅度的中位数的三倍心音片段全部处理完毕。