自动定位并修复电子听诊器削波失真波形的方法.pdf

《自动定位并修复电子听诊器削波失真波形的方法.pdf》由会员分享，可在线阅读，更多相关《自动定位并修复电子听诊器削波失真波形的方法.pdf（13页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 201910063222.0 (22)申请日 2019.01.23 (71)申请人 苏州美糯爱医疗科技有限公司 地址 215000 江苏省苏州市高新区科技城 科灵路78号05号楼-5-501 (72)发明人 蔡盛盛胡南徐兴国周宁 (74)专利代理机构 上海精晟知识产权代理有限 公司 31253 代理人 冯子玲 (51)Int.Cl. A61B 7/04(2006.01) (54)发明名称 一种自动定位并修复电子听诊器削波失真 波形的方法 (57)摘要 本发明公开一种自动定位并修复。

2、电子听诊 器削波失真波形的方法， 具有以下步骤： 读取听 诊信号； 对信号x求差分， 得到差分信号xd； 经门 限比较以及端点配对， 确定发生削波失真的K个 不交迭的时间区间， 按时间先后顺序排序， 并按 先进入后退出的方式形成堆栈； 若K0， 则无削 波失真区域， 不再需要修复， 若K1， 确定对堆栈 插值数据时间点集合K； 利用时间点集合K及 其每一个时间点所对应的x上的值x(K)， 经插 值拟合出第K段削波失真区域中的值以替代x中 原本的削波失真值； 移除堆栈顶部的时间区间； 输出无削波失真的听诊信号x。 本发明能够对听 诊器信号的削波失真区域进行自动检测与定位， 可自动修复削波失真区域。

3、的信号以及可自动修 复多个距离较近的削波失真区域。 权利要求书2页 说明书6页 附图4页 CN 109893163 A 2019.06.18 CN 109893163 A 1.一种自动定位并修复电子听诊器削波失真波形的方法， 其特征在于： 具有以下步骤： 步骤1、 读取缓存中时长为N的听诊信号采样序列x(n),n1,2,.,N， 并将其表示为向 量形式x； 步骤2、 对原始信号x求差分， 得到差分信号xd， 所述计算原始信号的差分信号的表达式 为xd(n)x(n)-x(n-1)， 其中， xd(1)0； 步骤3、 根据差分信号xd， 经门限比较以及端点配对， 确定发生削波失真的K个不交迭的 时。

4、间区间(K0)， 按时间先后顺序依次为： nK ,begin,nK ,end,nK-1 ,begin,nK-1 ,end,., n1,begin,n1,end， 按后出现的削波失真区间先进入后退出的方式形成堆栈； 步骤4、 若K0， 该段数据已无削波失真区域， 不再需要修复， 直接进入步骤7， 若K1， 确定对堆栈中最早的削波失真区域nK,begin,nK,end进行修复所需要的插值数据时间点集合 K； 步骤5、 利用时间点集合K及其每一个时间点所对应的x上的值x(K)， 经Hermite插 值， 拟合出第K段削波失真区域中的值以替代x中原本的削波失真值， 从而更新x； 步骤6、 移除堆栈顶部。

5、的时间区间nK,begin,nK,end， 令KK-1， 并返回步骤4。 步骤7、 输出无削波失真的听诊信号x。 2.根据权利要求1所述的一种自动定位并修复电子听诊器削波失真波形的方法， 其特 征在于： 所述步骤3中形成堆栈， 具体包括以下步骤： S1、 找出差分信号xd中绝对值大于阈值 max(|x|)的数值所对应的时间点， 形成集合 n1,n2,.,nL， 其中， 0 1为一个预设常数； S2、 从集合n1,n2,.,nL中找到相邻的xd上对应值符号相反的时间点两两进行配对， 确定其分别为某一削波失真区间的起始点nk,begin和结束点nk,end； 若在整个数据起始或者结 束附近找到一个。

6、未能配对的单个时间点， 则结合前一段或者后一段数据对其进行额外配 对； S3、 根据两两配对结果， 若共得到K个不交迭的时间区间以时间先后分别为nK,begin, nK,end,nK-1,begin,nK-1,end,.,n1,begin,n1,end的顺序， 按后出现的削波失真区间先进入堆 栈的方式形成堆栈。 3.根据权利要求1所述的一种自动定位并修复电子听诊器削波失真波形的方法， 其特 征在于： 所述步骤4中所需要的插值数据时间点集合K获取方法， 具体包括以下步骤： H1、 从时间点nK,begin出发， 逆着时间轴方向找到从nK,begin开始的x上第3个过零点所对应 的时间点nK,le。

7、ft3， 若x上对应该点不等于0， 则取最接近0的值所对应的时间点； H2、 从时间点nK,end出发， 顺着时间轴方向找到从nK,end开始的x上第1个过零点所对应的 时间点nK,right1， 若x上对应该点不恰好等于0， 则取最接近0的值所对应的时间点； H3、 确定修复削波失真区域nK,begin,nK,end所需要的数据时间点集合K为nK,left3, nK,begin-PnK,end+P,nK,right1， 其中， P为考虑过度按压时听诊器输出转换到饱和状态的转 换时间所设置的预留点数， 一般取值为110的非负整数。 4.根据权利要求1所述的一种自动定位并修复电子听诊器削波失真波。

8、形的方法， 其特 征在于： 所述步骤5中更新x的方法， 具体包括以下步骤： F1、 若时间集合K中共有Tk+1个时间点：估算每个时间点上的导数 权利要求书 1/2 页 2 CN 109893163 A 2 x (ni),i0,1,.,Tk； F1、 若时间集合K中共有Tk+1个时间点：估算每个时间点上数据点 x(ni)的导数x (ni),i0,1,.,Tk； F2、 对需要拟合修复的任意一点x(m)， 其中mnK,begin-P+1,nK,end+P-1是待拟合区域 上的时间点， 利用K中Tk+1个时间点上的数据点x(ni)及其导数x (ni),i0,1,.,Tk， 求 得其Hermite插值。

9、： F3、 更新x中第K个待修复区间上的值： x(m)xH(m),mnK,begin-P+1,nK,end+P-1。 权利要求书 2/2 页 3 CN 109893163 A 3 一种自动定位并修复电子听诊器削波失真波形的方法 技术领域 0001 本发明属于听诊器技术领域， 涉及到一种自动定位并修复电子听诊器削波失真波 形的方法。 背景技术 0002 听诊是内科医生第一时间了解患者病情的重要技术之一， 然而长期以来该项技术 都受到就诊地点、 医技水平等因素的制约， 随着物联网(IoT)技术的蓬勃发展， 各种类型的 电子听诊器不断涌现， 为病人的心肺音数据的实时监控、 自动转录、 云端诊疗与智能。

10、诊断带 来了可能。 0003 电子听诊器通过换能器(麦克风)将微弱的生理声信号(如心音、 肺音等)转化为电 信号， 目前电子听诊器所使用的一般都是灵敏度高、 频率特性优越的压电式麦克风， 谐振腔 上的薄膜收集声音振动导致晶体形变最终转化为电流。 然而在实际应用中， 由于听诊时的 过度按压导致电子听诊器输出电流快速达到饱和， 从而使得所记录的听诊信号发生削波失 真， 如图1所示的现象时有发生： 如儿科诊室由于幼儿患者的配合程度低， 在开展听诊过程 中可能会发生过度按压； 又如在患者自主使用电子听诊器将听诊数据上传进行云端诊疗的 过程中， 由于操作不当也很容易导致出现削波失真， 削波失真的出现， 。

11、一方面使得电子听诊 器的使用体验变差， 另一方面影响了信号质量进而影响后续心音定位、 心肺音自动诊断的 效果。 0004 对于按压导致的信号削波失真问题， 现有技术的缺点具体体现在以下方面： a.目 前尚无信号削波失真区域的自动定位方法； b.目前信号削波失真区域的自动修复方法尚不 完善： 主要是解决该问题的重要前提削波失真区域的自动定位的缺失， 同时如何能够 根据上下文信号内容更好地恢复出因削波导致的信号缺失部分仍是需要考虑的问题； c.目 前尚无对多个距离较近的削波失真区域进行修复的技术方案。 0005 特别的， 对于旨在实现高精度的胎儿实时监测、 心/肺功能智能评估、 心/肺疾病自 动诊。

12、断等电子听诊器的人工智能远景目标来说， 能够自动定位并修复听诊信号中的削波失 真， 是这些人工智能功能得以实现的重要前提之一， 因此， 电子听诊器信号发生削波失真是 电子听诊器亟待解决的共性问题。 0006 目前已有的与电子听诊器相关的专利的关注点包括设备的硬件采集系统、 传输系 统、 外观、 信号预处理(包括降噪、 心音定位、 心肺音分离等)、 信号智能分析(胎心监测、 基于 心音的心脏疾病智能诊断、 基于肺音的呼吸道疾病智能诊断)等等， 但尚无专利针对信号削 波失真问题提出自动修复解决方案。 0007 且论文 “Emmanouilidou D,Mccollum E D,Park D E,e。

13、t al.Computerized Lung Sound Screening for Pediatric Auscultation in Noisy Field EnvironmentsJ .IEEE Transactions on Biomedical Engineering,2018,65(7):1564-1574.” 中提出了一 种利用三次样条插值进行削波失真修复的方法。 然而该文献中的方法存在以下问题： (1)未 提出任何削波失真区域自动检测定位方法， 因此该方法得以实施的前提是要先人为选定削 说明书 1/6 页 4 CN 109893163 A 4 波失真区域， 不符合实际应用的要求。

14、； (2)该方法使用的三次样条插值拟合出来的曲线只能 保证在连接点的连续性， 不能保证其光滑性； 整条拟合曲线受所有插值点控制， 如果插值点 中有任何一个改动， 则整条曲线都受影响， 这对于听诊数据的修复而言是不符合其实际特 性的， 特别是对于呼吸道疾病的听诊而言其特异性的呼吸附加音的持续时间往往很短(如 单个湿啰音的持续时长在20毫秒左右)， 在这种情况下三次样条插值拟合会带来较大的误 差； (3)未提出当有距离较近的两段或者多段削波失真区域时的修复策略， 而如果直接对它 们同时进行修复可能因为待拟合点数过多导致欠拟合。 发明内容 0008 本发明的目的在于提供的自动定位并修复电子听诊器削波。

15、失真波形的方法， 解决 了现有技术无法对听诊器削波失真区域进行自动检测与定位， 且无法对削波失真区域进行 自动修复以及无法修复削波失真区域的信号的问题。 0009 本发明的目的可以通过以下技术方案实现： 0010 一种自动定位并修复电子听诊器削波失真波形的方法， 具有以下步骤： 0011 步骤1、 读取缓存中时长为N的听诊信号采样序列x(n),n1,2,.,N， 并将其表示 为向量形式x； 0012 步骤2、 对原始信号x求差分， 得到差分信号xd， 所述计算原始信号的差分信号的表 达式为xd(n)x(n)-x(n-1)， 其中， xd(1)0； 0013 步骤3、 根据差分信号xd， 经门限。

16、比较以及端点配对， 确定发生削波失真的K个不交 迭的时间区间(K0)， 按时间先后顺序依次为： nK,begin,nK,end,nK-1,begin,nK-1,end,., n1,begin,n1,end， 按后出现的削波失真区间先进入后退出的方式形成堆栈； 0014 步骤4、 若K0， 该段数据已无削波失真区域， 不再需要修复， 直接进入步骤7， 若K 1， 确定对堆栈中最早的削波失真区域nK,begin,nK,end进行修复所需要的插值数据时间点 集合K； 0015 步骤5、 利用时间点集合K及其每一个时间点所对应的x上的值x(K)， 经Hermite 插值， 拟合出第K段削波失真区域中的。

17、值以替代x中原本的削波失真值， 从而更新x； 0016 步骤6、 移除堆栈顶部的时间区间nK,begin,nK,end， 令KK-1， 并返回步骤4。 0017 步骤7、 输出无削波失真的听诊信号x。 0018 进一步地， 所述步骤3中形成堆栈， 具体包括以下步骤： 0019 S1、 找出差分信号xd中绝对值大于阈值 max(|x|)的数值所对应的时间点， 形成 集合n1,n2,.,nL， 其中， 0 1为一个预设常数； 0020 S2、 从集合n1,n2,.,nL中找到相邻的xd上对应值符号相反的时间点两两进行 配对， 确定其分别为某一削波失真区间的起始点nk,begin和结束点nk,end。

18、； 若在整个数据起始 或者结束附近找到一个未能配对的单个时间点， 则结合前一段或者后一段数据对其进行额 外配对； 0021 S3、 根据两两配对结果， 若共得到K个不交迭的时间区间以时间先后分别为 nK,begin,nK,end,nK-1,begin,nK-1,end,.,n1,begin,n1,end的顺序， 按后出现的削波失真区间 先进入堆栈的方式形成堆栈。 0022 进一步地， 所述步骤4中所需要的插值数据时间点集合K获取方法， 具体包括以 说明书 2/6 页 5 CN 109893163 A 5 下步骤： 0023 H1、 从时间点nK,begin出发， 逆着时间轴方向找到从nK,be。

19、gin开始的x上第3个过零点所 对应的时间点nK,left3， 若x上对应该点不等于0， 则取最接近0的值所对应的时间点； 0024 H2、 从时间点nK,end出发， 顺着时间轴方向找到从nK,end开始的x上第1个过零点所对 应的时间点nK,right1， 若x上对应该点不恰好等于0， 则取最接近0的值所对应的时间点； 0025 H3、 确定修复削波失真区域nK ,begin,nK ,end所需要的数据时间点集合K为 nK,left3,nK,begin-PnK,end+P,nK,right1， 其中， P为考虑过度按压时听诊器输出转换到饱和 状态的转换时间所设置的预留点数， 一般取值为11。

20、0的非负整数。 0026 进一步地， 所述步骤5中更新x的方法， 具体包括以下步骤： 0027F1、 若时间集合K中共有Tk+1个时间点： n0,n1,估算每个时间点上数据点 x(ni)的导数x (ni),i0,1,.,Tk； 0028 F2、 对需要拟合修复的任意一点x(m)， 其中mnK,begin-P+1,nK,end+P-1是待拟合 区域上的时间点， 利用K中Tk+1个时间点上的数据点x(ni)及其导数x (ni),i0,1,., Tk， 求得其Hermite插值： 0029 0030 F3、 更新x中第K个待修复区间上的值： x(m)xH(m),mnK,begin-P+1,nK,en。

21、d+P-1。 0031 本发明的有益效果： 0032 本发明提供的自动定位并修复电子听诊器削波失真波形的方法， 通过对信号求差 分、 与阈值比较、 配对， 以实现对听诊器信号的削波失真区域进行自动检测与定位； 通过根 据削波失真区域的自动定位结合过零点选择修复所需的插值数据时间点集合， 并通过 Hermite插值可自动修复削波失真区域的信号； 通过采用堆栈的 “先进后出” 特点， 可依次实 现多个削波失真区域的自动修复， 并能够解决削波失真区域相距很近时的修复问题， 具有 定位与修复的速度快， 准确性高的特点。 附图说明 0033 为了更清楚地说明本发明实施例的技术方案， 下面将对实施例描述所。

22、需要使用的 附图作简单地介绍， 显而易见地， 下面描述中的附图仅仅是本发明的一些实施例， 对于本领 域普通技术人员来讲， 在不付出创造性劳动的前提下， 还可以根据这些附图获得其他的附 图。 0034 图1为现有技术中听诊信号发生削波失真的示意图； 0035 图2为本发明中一种自动定位并修复电子听诊器削波失真波形的方法的示意图； 0036 图3为本发明实施例1中修复的示意图； 0037 图4为本发明实施例2中修复的示意图； 说明书 3/6 页 6 CN 109893163 A 6 具体实施方式 0038 下面将结合本发明实施例中的附图， 对本发明实施例中的技术方案进行清楚、 完 整地描述， 显然。

23、， 所描述的实施例仅仅是本发明一部分实施例， 而不是全部的实施例。 基于 本发明中的实施例， 本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它 实施例， 都属于本发明保护的范围。 0039 如图1所示， 一种自动定位并修复电子听诊器削波失真波形的方法， 具有以下步 骤： 0040 步骤1、 读取缓存中时长为N的听诊信号采样序列x(n),n1,2,.,N， 并将其表示 为向量形式x； 0041 步骤2、 对原始信号x求差分， 得到差分信号xd， 所述计算原始信号的差分信号的表 达式为xd(n)x(n)-x(n-1)， 其中， xd(1)0； 0042 步骤3、 根据差分信号xd， 经。

24、门限比较以及端点配对， 确定发生削波失真的K个不交 迭的时间区间(K0)， 按时间先后顺序依次为： nK,begin,nK,end,nK-1,begin,nK-1,end,., n1,begin,n1,end， 按后出现的削波失真区间先进入后退出的方式形成堆栈； 0043 步骤3中形成堆栈， 具体包括以下步骤： 0044 S1、 找出差分信号xd中绝对值大于阈值 max(|x|)的数值所对应的时间点， 形成 集合n1,n2,.,nL， 其中， 0 1为一个预设常数； 0045 S2、 从集合n1,n2,.,nL中找到相邻的xd上对应值符号相反的时间点两两进行 配对， 确定其分别为某一削波失真区。

25、间的起始点nk,begin和结束点nk,end； 若在整个数据起始 或者结束附近找到一个未能配对的单个时间点， 则结合前一段或者后一段数据对其进行额 外配对； 0046 S3、 根据两两配对结果， 若共得到K个不交迭的时间区间以时间先后分别为 nK,begin,nK,end,nK-1,begin,nK-1,end,.,n1,begin,n1,end的顺序， 按后出现的削波失真区间 先进入堆栈的方式形成堆栈。 0047 步骤4、 若K0， 该段数据已无削波失真区域， 不再需要修复， 直接进入步骤7， 若K 1， 确定对堆栈中最早的削波失真区域nK,begin,nK,end进行修复所需要的插值数据。

26、时间点 集合K； 0048 步骤4中若K1确定对堆栈中最早的削波失真区域nK,begin,nK,end进行修复所需 要的插值数据时间点集合K获取方法， 具体包括以下步骤： 0049 H1、 从时间点nK,begin出发， 逆着时间轴方向找到从nK,begin开始的x上第3个过零点所 对应的时间点nK,left3， 若x上对应该点不等于0， 则取最接近0的值所对应的时间点； 0050 H2、 从时间点nK,end出发， 顺着时间轴方向找到从nK,end开始的x上第1个过零点所对 应的时间点nK,right1， 若x上对应该点不恰好等于0， 则取最接近0的值所对应的时间点； 0051 H3、 确定。

27、修复削波失真区域nK ,begin,nK ,end所需要的数据时间点集合K为 nK,left3,nK,begin-PnK,end+P,nK,right1， 其中， P为考虑过度按压时听诊器输出转换到饱和 状态的转换时间所设置的预留点数， 一般取值为110的非负整数。 0052 步骤5、 利用时间点集合K及其每一个时间点所对应的x上的值x(K)， 经Hermite 插值， 拟合出第K段削波失真区域中的值以替代x中原本的削波失真值， 从而更新x； 0053 更新x的方法， 包括以下步骤： 说明书 4/6 页 7 CN 109893163 A 7 0054F1、 若时间集合K中共有Tk+1个时间点：。

28、 n0,n1,估算每个时间点上数据点 x(ni)的导数x (ni),i0,1,.,Tk； 0055 F2、 对需要拟合修复的任意一点x(m)， 其中mnK,begin-P+1,nK,end+P-1是待拟合 区域上的时间点， 利用K中Tk+1个数据点x(ni)及其导数x (ni),i0,1,.,Tk， 求得其 Hermite插值： 0056 0057 F3、 更新x中第K个待修复区间上的值： x(m)xH(m),mnK,begin-P+1,nK,end+P-1。 0058 步骤6、 移除堆栈顶部的时间区间nK,begin,nK,end， 令KK-1， 并返回步骤4。 0059 步骤7、 输出无削。

29、波失真的听诊信号x。 0060 实施例1： 0061 读取一段小儿肺炎患者的肺音听诊数据x， 其中包含两处削波失真区域， 对其进行 削波失真区域的自动定位与修复， 数据采样率为4KHz， 时长为2.5秒。 0062 首先， 对该段数据求差分， 差分信号xd如图3(a)所示， 可见在两处削波失真区域的 开始与结束位置都有一对极性相反的幅度很大的峰值。 0063 然后， 找出差分信号xd中绝对值大于阈值0.8max(|x|)的数值所对应的时间点， 从中找到相邻的xd上对应值符号相反的时间点两两进行配对， 可检测出两处削波失真区 域， 并确定其各自的起始点和终点， 见图3(b)所示； 0064 最后。

30、， 将这两个定位出的削波失真区间按逆时顺序放入堆栈， 利用堆栈操作依次 确定其插值点集合， 并通过Hermite插值依次进行修复， 并完成信号更新， 最终修复后的信 号如图3(c)所示， 其中考虑过度按压时听诊器输出转换到饱和状态的转换时间所设置的预 留点数P5。 0065 实施例2： 0066 读取一段正常人心音听诊数据x， 数据中共有5处削波失真区域， 其中， 有两处削波 失真区域离得很近， 相距0.0136秒， 对其进行削波失真区域的自动定位与修复， 数据采样率 为4KHz， 时长为1.5秒。 0067 首先， 对该段数据求差分， 差分信号xd如图4(a)所示， 可见在5处削波失真区域的。

31、 开始与结束位置都有一对极性相反的幅度很大的峰值。 0068 然后， 找出差分信号xd中绝对值大于阈值0.8max(|x|)的数值所对应的时间点， 从中找到相邻的xd上对应值符号相反的时间点两两进行配对， 可检测出全部5处削波失真 区域， 并确定其各自的起始点和终点， 见图4(b)所示； 0069 最后， 将这5个定位出的削波失真区间按逆时顺序放入堆栈， 利用堆栈操作依次确 定其插值点集合并通过Hermite插值依次进行修复， 并完成信号更新， 最终修复后的信号如 图4(c)所示， 其中， 考虑过度按压时听诊器输出转换到饱和状态的转换时间所设置的预留 点数P5， 可见即便是两处相距很近的削波失。

32、真区域， 利用本发明所提的方法仍然能够实 现较为精确的修复。 说明书 5/6 页 8 CN 109893163 A 8 0070 以上内容仅仅是对本发明的构思所作的举例和说明， 所属本技术领域的技术人员 对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代， 只要不偏离发明 的构思或者超越本权利要求书所定义的范围， 均应属于本发明的保护范围。 说明书 6/6 页 9 CN 109893163 A 9 图1 说明书附图 1/4 页 10 CN 109893163 A 10 图2 说明书附图 2/4 页 11 CN 109893163 A 11 图3 说明书附图 3/4 页 12 CN 109893163 A 12 图4 说明书附图 4/4 页 13 CN 109893163 A 13 。