CN111714151B - 一种心音的分析与标注方法和系统 - Google Patents

一种心音的分析与标注方法和系统 Download PDF

Info

Publication number
CN111714151B
CN111714151B CN202010598567.9A CN202010598567A CN111714151B CN 111714151 B CN111714151 B CN 111714151B CN 202010598567 A CN202010598567 A CN 202010598567A CN 111714151 B CN111714151 B CN 111714151B
Authority
CN
China
Prior art keywords
heart sound
heart
map
marking
playing
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010598567.9A
Other languages
English (en)
Other versions
CN111714151A (zh
Inventor
孙锟
王景涛
周欣
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Airdoc Technology Co Ltd
Original Assignee
Beijing Airdoc Technology Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Airdoc Technology Co Ltd filed Critical Beijing Airdoc Technology Co Ltd
Priority to CN202010598567.9A priority Critical patent/CN111714151B/zh
Publication of CN111714151A publication Critical patent/CN111714151A/zh
Application granted granted Critical
Publication of CN111714151B publication Critical patent/CN111714151B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B7/00Instruments for auscultation
    • A61B7/02Stethoscopes
    • A61B7/04Electric stethoscopes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/72Signal processing specially adapted for physiological signals or for diagnostic purposes
    • A61B5/7235Details of waveform analysis
    • A61B5/7264Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems
    • A61B5/7267Classification of physiological signals or data, e.g. using neural networks, statistical classifiers, expert systems or fuzzy systems involving training the classification device
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/03Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 characterised by the type of extracted parameters
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10LSPEECH ANALYSIS OR SYNTHESIS; SPEECH RECOGNITION; SPEECH OR VOICE PROCESSING; SPEECH OR AUDIO CODING OR DECODING
    • G10L25/00Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00
    • G10L25/48Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use
    • G10L25/51Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination
    • G10L25/66Speech or voice analysis techniques not restricted to a single one of groups G10L15/00 - G10L21/00 specially adapted for particular use for comparison or discrimination for extracting parameters related to health condition

Abstract

本发明提供了一种心音的分析与标注方法和系统,其方法包括:根据导入的心音信号样本转换得到心音图谱,心音图谱包括心音图和心音频响图;播放心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内。本发明通过导入的心音样本,将心音样本转换成心音图和心音频响图,通过观察心音频响图和心音图的时域变化趋势和情况,判断心音的响度,观察心音响度颜色的划分,区分正常心音和异常心音,配合听取导入的音频文件,通过倍速调节、放大时间轴,更精准的判断出异常心音出现的区间,以及响度的变化,并且通过对异常心音的区间进行标注,导出形成样本数据库,有益于后期的心音计算以及患者病症数据库的研究分析。

Description

一种心音的分析与标注方法和系统
技术领域
本发明涉及医疗测试技术领域领域,尤指一种心音的分析与标注方法和系统。
背景技术
心音是由于心肌收缩、心脏瓣膜关闭和血液撞击心室壁、大动脉壁等引起的振动所产生的,是人体最重要的生理信号之一,当人体出现一些疾病时,心音就会产生一些重要的病理信息,而这些病理信息在许多疾病中都有特征体现,因此通过观察心音信号的变化可以对症状的判断提供一定的依据。
在现有的技术中,心音的分类识别功能大多数只是区分正常和异常心音的信号,而且在心音的分段过程中,通常采取计算单一阈值来进行的分段,这种计算方法容易实现,但得出的计算结果准确性相对来说不是很准确,对外界杂音的抗干扰能力不强,从而会导致心音的分段出现严重的误差,进而影响患者的症状判断结果。
因此,如何以直观、简洁的方式辅助医师判断心音的正异常情况是亟需解决的问题。
发明内容
本发明的目的是提供一种心音的分析与标注方法和系统,实现以直观、简洁的方式辅助医师判断心音的正异常情况。
本发明提供的技术方案如下:
本发明提供一种心音的分析与标注方法,包括步骤:
根据导入的心音信号样本转换得到心音图谱,所述心音图谱包括心音图和心音频响图;
播放所述心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内。
进一步的,所述根据导入的心音信号样本转换得到心音图谱,所述心音图谱包括心音图和心音频响图之前包括步骤:
采集获取多个听诊区的心音音频文件;
将所述心音音频文件与对应患者身份进行绑定,且将所述心音音频文件转换为统一音频格式得到待导入的心音信号样本。
进一步的,所述根据导入的心音信号样本转换得到心音图谱,所述心音图谱包括心音图和心音频响图包括步骤:
提取所导入心音信号样本的包络特征以及特征参数;所述特征参数包括:频率、第一心音、第二心音、收缩期以及舒张期的持续时间;
根据包络特征以及特征参数自动分割得到五个瓣膜区,后根据所述瓣膜区、包络特征以及特征参数分割得到各瓣膜区对应的多个心动周期;
根据心动周期进行区间划分并生成心音图谱;所述心音图谱包括心动周期、第一心音、第二心音的区间位置。
进一步的,所述播放所述心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内包括步骤:
同进度播放各区间位置对应的心音图谱和心音信号样本;
在获取到触发标记信息后弹窗显示心音评价项及其对应的标记选项按钮;
获取医生根据心音的响度变化针对各心音评价项所输入的标记信息,根据所述标记信息将对应区间范围进行标注;所述标记信息包括医生所选中的标记选项按钮。
进一步的,还包括步骤:
在播放期间根据获取到的播放调节指令同步调整图谱滚动播放状态和音频语音播放状态。
进一步的,所述播放所述心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内之后包括步骤:
获取到标记结束指令以停止标注功能后,将标注后的心音图谱导出。
本发明还提供一种心音的分析与标注系统,包括:
导入转换模块,用于根据导入的心音信号样本转换得到心音图谱,所述心音图谱包括心音图和心音频响图;
处理模块,用于播放所述心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内。
进一步的,还包括:
采集模块,用于采集获取多个听诊区的心音音频文件;
所述处理模块,还用于将所述心音音频文件与对应患者身份进行绑定,且将所述心音音频文件转换为统一音频格式得到待导入的心音信号样本。
进一步的,所述导入转换模块包括:
提取单元,用于提取所导入心音信号样本的包络特征以及特征参数;所述特征参数包括:频率、第一心音、第二心音、收缩期以及舒张期的持续时间;
分割单元,用于根据包络特征以及特征参数自动分割得到五个瓣膜区,后根据所述瓣膜区、包络特征以及特征参数分割得到各瓣膜区对应的多个心动周期;
生成单元,用于根据心动周期进行区间划分并生成心音图谱;所述心音图谱包括心动周期、第一心音、第二心音的区间位置。
进一步的,所述处理模块包括:
播放单元,同进度播放各区间位置对应的心音图谱和心音信号样本;
显示单元,在获取到触发标记信息后弹窗显示心音评价项及其对应的标记选项按钮;
获取单元,用于获取医生根据心音的响度变化针对各心音评价项所输入的标记信息;还用于获取播放调节指令;
标注单元,根据所述标记信息将对应区间范围进行标注;所述标记信息包括医生所选中的标记选项按钮;
调节单元,用于在播放期间根据获取到的播放调节指令同步调整图谱滚动播放状态和音频语音播放状态。
通过本发明提供的一种心音的分析与标注方法和系统,能够以直观、简洁的方式辅助医师判断心音的正异常情况。
附图说明
下面将以明确易懂的方式,结合附图说明优选实施方式,对一种心音的分析与标注方法和系统的上述特性、技术特征、优点及其实现方式予以进一步说明。
图1是本发明一种心音的分析与标注方法的一个实施例的流程图;
图2是本发明一种心音的分析与标注方法的心音图的示意图;
图3是本发明本发明一种心音的分析与标注方法的心音频响图的示意图;
图4是本发明一种心音的分析与标注方法的一个实施例的流程图;
图5是本发明一种心音的分析与标注方法的心音频响图的X轴放大示意图;
图6是本发明一种心音的分析与标注方法的标注样本界面图。
具体实施方式
以下描述中,为了说明而不是为了限定,提出了诸如特定系统结构、技术之类的具体细节,以便透彻理解本申请实施例。然而,本领域的技术人员应当清楚,在没有这些具体细节的其他实施例中也可以实现本申请。在其他情况中,省略对众所周知的系统、装置、电路以及方法的详细说明,以免不必要的细节妨碍本申请的描述。
应当理解,当在本说明书和所附权利要求书中使用时,术语“包括”指示所述描述特征、整体、步骤、操作、元素和/或组件的存在,但并不排除一个或多个其他特征、整体、步骤、操作、元素、组件和/或集合的存在或添加。
为使图面简洁,各图中只示意性地表示出了与本发明相关的部分,它们并不代表其作为产品的实际结构。另外,以使图面简洁便于理解,在有些图中具有相同结构或功能的部件,仅示意性地绘示了其中的一个,或仅标出了其中的一个。在本文中,“一个”不仅表示“仅此一个”,也可以表示“多于一个”的情形。
还应当进一步理解,在本申请说明书和所附权利要求书中使用的术语“和/或”是指相关联列出的项中的一个或多个的任何组合以及所有可能组合,并且包括这些组合。
另外,在本申请的描述中,术语“第一”、“第二”等仅用于区分描述,而不能理解为指示或暗示相对重要性。
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对照附图说明本发明的具体实施方式。显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图,并获得其他的实施方式。
本发明的一个实施例,如图1所示,一种心音的分析与标注方法,包括:
S100根据导入的心音信号样本转换得到心音图谱,心音图谱包括心音图和心音频响图;
具体的,由于人体心脏的收缩与舒张,血液在心压作用下在主动脉内流动,均能产生机械振动,这种振动传递到人体表面产生一种声音,通常称为心音。心音信号经过换能器转换成电讯号,再经滤波、放大,形成一组随时间变化的波形图,这就是心音图。如图2所示,心音图是一种二维图,其中横坐标表示时间,纵坐标表示电压幅度值,通过对心音信号建立时域-电压幅度值关系得到心音图以表征心音的时域特性。
另外,如图3所示,心音频响图是一种三维图,其中横坐标表示时间,纵坐标表示频域,颜色柱表示声音的响度。通过对心音信号建立时域-频谱能量关系,进一步建立频谱能量与声音响度的关系,从而得到对应的心音频响图以表征心音的时频与响度特性。
S200播放心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内。
具体的,通过系统显示界面随着时间的流逝,滚动播放心音图谱,在滚动播放的过程中,医生可以观察心音图谱进行分析判断得到心音分析结果,其中,心音分析结果包括但是不限于判断各区间的心音质量等级(或者分数),各区间对应的心音是正常心音还是异常心音,以及异常心音的类型,进而根据判断结果输入标记信息,以便将心音分析结果在对应区间位置处进行标注说明。
本实施例中,在医生观察心音图和/或心音频响图的过程中,可以对心音图和/或心音频响图进行标注,从而进入标注样本界面,这样以便医生后期对患者心音异常区间的审核,此外,后期也可界面显示心音图谱上不同的标注内容,以直观、简洁的方式辅助医师判断心音的正异常情况。
本发明的一个实施例,如图4所示,一种心音的分析与标注方法,包括:
S010采集获取多个听诊区的心音音频文件;
具体的,听诊区选择为五个,分别为:二尖瓣听诊区、肺动脉瓣听诊区、主动脉瓣听诊区、主动脉瓣第二听诊区、三尖瓣听诊区,并且采集各听诊区的心音信号时,采用电子听诊器分别依次采集五个听诊区的心音,避免最后呈现的心音图谱是连续的,分别依次采集五个听诊区的心音能够使得分割出的心动周期间隔明显,有利于利于医生观察,标注。
S020将心音音频文件与对应患者身份进行绑定,且将心音音频文件转换为统一音频格式得到待导入的心音信号样本;
具体的,每一份心音音频文件都对应一位患者,将心音音频文件与患者信息(包括但是不限于患者的身份证号码、电话号码或者就诊卡卡号)绑定在一起,避免出现心音分析结果与患者身份不一致的情况,从而大大降低误诊率。将采集到的心音音频文件整理成系统所支持且一致的音频格式,这样,将心音音频文件转换为统一音频格式的心音信号样本后,导入系统进行播放分析时便于系统识别,使得系统能够快速识别并播放心音信号样本,减少识别转码的时间,进而间接降低分析心音的等待时间。
S110提取所导入心音信号样本的包络特征以及特征参数;特征参数包括:频率、第一心音、第二心音、收缩期以及舒张期的持续时间;
具体的,通过心音分割流程提取包络特征以及特征参数:对心音信号样本利用离散小波(sym6)去噪,然后对去噪信号加窗计算归一化平均香农能量来表征信号的包络特征,后利用双门限法结合信号的能量分布特征寻找心音段的边界。另外,特征参数中的频率是对原始时域信号即心音信号样本进行分段处理,然后用小波、短时傅里叶变换或者Wigner-Ville分布等方法将心音的从时域搬移到时频空间上,纵坐标的数值即表示该空间的频域维度。其中包络特征利用加窗求平均香农能量方法来提取包络,此为现有技术在此不再一一赘述。
心音信号的包络特征能够比心音信号更好反应心音的特征,包络特征包括心音信号样本的有效时域信息,这样,通过包络特征能够进行心音各个成分的识别,大大提升检测精确度,减少工作量和计算量,简化心音分析识别的流程,进而大大提高分析识别效率。
S120根据包络特征以及特征参数自动分割得到五个瓣膜区,后根据瓣膜区、包络特征以及特征参数分割得到各瓣膜区对应的多个心动周期;
具体的,五个瓣膜区包括二尖瓣区、肺动脉瓣区、主动脉瓣区、主动脉瓣第二区、三尖瓣区。根据心音信号的包络特征以及特征参数先得到五个瓣膜区,在五个瓣膜区的基础上,根据心音信号的包络特征以及特征参数分割出心动周期。
一个心动周期中有第一心音、第二心音、第三心音和第四心音。第一心音发生于心缩初期,标志心室收缩开始,是由于心室肌收缩,房室瓣突然关闭及随后射血入动脉等引起的振动而形成。第二心音发生于心舒期开始,由于动脉瓣迅速关闭,血流冲击大动脉根部及心室内壁振动而形成。通过系统自动识别心音,分别将五个瓣膜区的多个心动周期分割出来,心动周期包括收缩期和舒张期,即心脏每收缩和舒张一次构成一个心动周期。通常以心房开始收缩作为一个心动周期的起点,以心房结束舒张作为一个心动周期的终点。
S130根据心动周期进行区间划分并生成心音图谱;心音图谱包括心动周期、第一心音、第二心音的区间位置;
具体的,建立心音频率图谱,心音频率图谱包括:心音频响图和心音图。根据包络特征、特征参数以及五个瓣膜区自动分割出对应心动周期。心动周期时程的长短与心率有关,以便医生通过观察心动周期中的第一心音、第二心音、收缩期和舒张期的情况,来判断患者的心血管系统是否正常,从而使得医生在观察心音图谱的过程中对异常心音进行标注。心动周期的划分是划分成第一心音(S1)和第二心音(S2),也就是根据上文提到的心音分割流程提取包络特征,计算第一心音(The First Heart Sound,以下简称S1)或第二心音(TheSecond Heart Sound,以下简称S2)期的持续时间,根据收缩期大都小于舒张期,将收缩期前的心音成分确定为S1,将舒张期前的心音成分确定为S2,从而可以确定心音图谱的第一心音和第二心音的区间位置。
当然,由于异常心音导致心动周期划分不准确,不能及时的进行更改调整,从而容易导致在后续对异常心音进行标注的过程中,由于心动周期划分不准确从而增加了工作难度。当实际心音信号出现的区间位置与分割出的心动周期的区间位置出现偏差时,通过该系统选中要重设的心动周期,系统会自动选中原始的心动周期范围,然后拖动滚动条,自由调节心动周期范围,通过选中重设范围进行保存,以此实现心动周期的范围调整,以此来辅助医生准确的观察、判断正常心音与异常心音出现的区间。
S210同进度播放各区间位置对应的心音图谱和心音信号样本;
具体的,系统能够根据心音音频播放的进度,与心音图谱中相对应的心音强度变化同步同进度显示。此外,系统能够实现心音频响图与心音图的切换,进行切换的过程中,心音图与心音频响图会同步同进度显示,并且心音图谱的滚动播放进度与心音信号样本的语音播放进度系统,以便于医生根据声图结合的方式对比病理特征,分析患者心音是否异常。示例性的,系统滚动播放第一个心动区间的第一心音区间所对应的心音图的同时,语音播放第一个心动区间的第一心音区间所对应的心音音频,此时,如果医生需求切换心音频响图进一步查看确定分析时,系统切换滚动播放的心音频响图过程中,心音图与心音频响图会同步同进度显示第一个心动区间的第一心音区间对应的波形变化,同时,语音播放第一个心动区间的第一心音区间所对应的心音音频。
通过该系统的频响/时能将心音频响图切换至心音图,上述功能在心音图中同样适用,通过观察心音图的变化趋势和情况可以看出心音增强或减弱,结合心音频响图,取长补短,可以更准确地记录正常心音和异常心音。即观察心音图的变化趋势和情况得到心音信号的频率信息,观察频率的大小,确定病理杂音分布的区间,通过来回切换频响/时能的方式,对心音图中对应的心音异常区间和与之对应的心音频响图中的心音异常区间相对比,进一步地得到很多与心音频响图相对应的病理特征,辅助医生准确的做出诊断。
S215在播放期间根据获取到的播放调节指令同步调整图谱滚动播放状态和音频语音播放状态;
具体的,播放调节指令可以包括三种类型的指令,一是即针对图谱滚动播放状态,也针对音频语音播放状态的指令,二是只针对图谱滚动播放状态的指令,三是只针对音频语音播放状态的指令。
第一种类型的播放调节指令包括播放时间调节指令,这样系统可通过拖动时间轴来调节包揽区域,系统对心音音频的播放区间和心音图谱的播放区间进行调节,其中,包揽与未包揽的区域由颜色的深浅来划分。通过该系统放大时间轴进行心音音频文件的播放时,播放的进度会与心音图中相对应的心音强度变化同步显示,在判断心动周期时,更好的通过观察所选区域内的心音变化,来判断心音是否异常和是否存在杂音。
第一种类型的播放调节指令还包括播放倍速调节指令,这样系统可对心音音频的播放倍速和心音图谱中相对应的心音强度变化的播放倍速进行调节。
第一种类型的播放调节指令还包括循环播放指令,这样系统可对心音音频和心音图谱的播放次数进行调节。
第二类型的播放调节指令包括播放调节指令包括缩放调节指令,这样系统可对通过对心音图谱的X轴、Y轴和颜色柱进行缩放,更加清楚的观察心音的频率大小、响度强弱和变化情况。如图5所示,医生在进行标注的过程中,直接通过缩放坐标轴读取心音图谱,通过频响/时能进行心音频响图和心音图的切换,初步的判断正常心音和异常心音,并且根据该系统已经设定好的心音标注的类型进行标注。对于不太确定或者需要更加进一步地了解心音的正异常情况时,通过播放心音音频,自由的拖动进度条,从而确定需要听取的心音音频的范围,再通过暂停播放或者循环播放的方式重复听取心音音频,进而确定心音的正异常情况,完成标注,在减少了医生诊断时间的过程中,又提高了心音标注的准确性。
第三种类型的播放调节指令包括播放调节指令包括音量倍率调节指令,这样系统可对心音音频的播放音量进行调节。通过重播或者循环重复播放音频文件,配合调节播放倍速和音量大小,进而通过听取心音的音频文件和观察心音频响图,准确地判断出异常心音,并对其异常心音区间进行标注。
通过系统自由的调节音频的播放倍速,用于辅助医生更精准的来判断患者的心音异常部分,同时通过对心音图的X轴和Y轴进行缩放,更加清楚的观察心音的频率大小、响度强弱和变化情况,用于提高诊断过程的效率以及准确性,并且还能通过该系统对心音的音频播放的音量大小进行倍率调节,结合除该系统以外的调节音量大小的播放器来使用,用于辅助判断心音异常,提高工作效率。
S220在获取到触发标记信息后弹窗显示心音评价项及其对应的标记选项按钮;
具体的,医生通过上述方式观察心音图谱上各个心动周期的波形变化,以及听取各个心动周期对应心音音频后,如果确定某一区间(包括心动周期、第一心音、第二心音)时,那么就点击选中“暂停开始标注”,或者确定整段具有标记需求时,就点击选中“整段标注”,以便界面切换显示心音评价项及其对应的标记选项按钮,使得医生进行标记。
心音评价项包括心音信号质量、第一心音异常类型,第二信号异常类型、是否闻及第三心音和是否有杂音等。不同的心音评价项具有对应的标记选项按钮。例如系统按照现有的心音信号质量评估标准进行划分:心音信号质量有六个等级,还例如第一心音和第二心音的异常类型有三类。上述只是示例,当然,系统能够根据用户不同的标注需求,及时的更新心音评价项及其对应的标记选项按钮。
示例性的,如图6所示,在医生观察心音频率图谱的过程中,对心动周期、第一心音、第二心音、以及整段进行标注,从而进入标注样本界面。
系统的界面显示心音评价项即其对应的标记选项按钮:
按心音质量分为六级:很差(一级)、较差(二级)、一般(三级)、较好(四级)、很好(五级)、极好(六级)。
按第一心音异常分为三类:第一心音增强、第一心音减弱、第一心音分裂。
按第二心音异常分为三类:第二心音增强、第二心音减弱、第二心音分裂。
还对判断心音是否闻及第三心音以及是否有杂音的结果进行备注。
其中,正异常心音信号区分:通过上述心音分割流程提取包络特征以及特征参数得到心音各成分后,划分出各个心音周期(一般一个完整的心音周期包括第一心音、第二心音、第三心音和第四心音四个成分,其中S1与S2成分相较而言比较明显,因此方法中采用S1起始段开始到下一次S1起始段算一个完整的心音周期),并结合专家标注系统进行心音的标注。为了便于心音快速分类识别,这里构建心音信号的特征矢量,其中每个心音周期的特征矢量包括S1的间期T1(即从S1结束到S2开始的时间为T1)和S2的间期T2(即从S2结束到S1开始的时间为T2),S1持续时间T3,S2持续时间T4,S1波峰强度值P1,S2波峰强度峰值P2,S1能量值E1,S2能量值E2,收缩期与舒张期的持续时间T2+T3,T1+T4,频域特征:2~10HZ,10Hz~20Hz,20Hz~30Hz,30~40Hz,40~60Hz,60~150Hz,150~200Hz,200Hz以上能量占比以及梅尔倒谱特征等。根据特征矢量与标注结果训练SVM分类器。优化参数并得到训练模型后即可进行正,异常心音信号的区分。系统可自行、智能识别各个心音周期的心音是否异常的判断,并弹窗提示医生,医生根据提示结果再结合自身观察心音图谱和听取心音音频的结果进行综合判断。
S230获取医生根据心音的响度变化针对各心音评价项所输入的标记信息,根据标记信息将对应区间范围进行标注;标记信息包括医生所选中的标记选项按钮;
具体的,系统根据导入的心音文件,自动将心音的响度,通过颜色的划分表现出来,并且医生通过观察心音频响图来区分正常心音和异常心音。具体为系统根据五个听诊区分别分割出心动周期、第一心音和第二心音的区间以及心音信号的频率大小,根据系统划分的五个瓣膜区作为听诊区,通过医生对心音频响图观察到频谱随时域变化的趋势和情况,同时根据心音频响图的区域像素的饱和度,明度,色度来辨别心音响度的大小,用于区分正常心音和异常心音,然后医生在心音频响图的各心动周期的对应不同区间进行标注,当然也可同时在第一心音以及第二心音异常的区间内进行标注,这样记录存档后,便于后续辅助医生准确的观察、判断正常心音与异常心音出现的生理区间、频率特征及强度。
目前对于心音信号的分段都是通过计算单一阈值来进行的,虽然很容易去实现,但是结果准确性不高,抗干扰能力不强,并且当心音信号中出现一段噪音大于所选定的阈值时,就会出现分段不准确的情况,从而影响最终的结果。本发明通过系统显示心音频率图谱,便于更精准、细致的分辨出正常心音与异常心音。
然后,医生分别通过观察心音图谱的五个瓣膜区的心动周期分布情况,心音响度的强弱,从而判断出异常心音分布的区间,并对其进行标注,即医生通过心音频响图(或心音图)可以观察到时域变化趋势和情况,同时可以通过颜色的划分来判断心音的响度,用于区分正常心音和异常心音,通过对正常心音和异常心音的标注,从而用于医生根据心音的正异常情况进行诊断,后期能根据医生在诊断过程中不同的标注需求,通过该系统对标注内容进行更新,以此来辅助医生减少诊断的时间,提高工作效率。
S300获取到标记结束指令以停止标注功能后,将标注后的心音图谱导出。
具体的,通过系统导出标注结果即标注后的心音图谱,即医生在标注完以后,通过点击选中系统的标注结束功能,从而将标注的结果导出,形成样本数据库,用于后期的心音计算以及患者病症数据库的研究分析。
本发明系统通过导入的心音样本,将音频文件转换成心音图和心音频响图,用于记录正常心音和异常心音,通过观察心音频响图的时域变化趋势和情况,判断心音的响度,观察心音响度颜色的划分,区分正常心音和异常心音,配合听取导入的音频文件,通过倍速调节,以及放大时间轴,更精准的判断出异常心音出现的区间,以及响度的变化,并且通过该系统对其异常心音的区间进行标注,导出,形成样本数据库,用于后期的心音计算以及患者病症数据库的研究分析。
本发明的一个实施例,一种心音的分析与标注系统包括:
导入转换模块,用于根据导入的心音信号样本转换得到心音图谱,心音图谱包括心音图和心音频响图;
处理模块,用于播放心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内。
具体的,本实施例是上述方法实施例对应的系统实施例,具体效果参见上述系统实施例,在此不再一一赘述。
基于前述实施例,还包括:
采集模块,用于采集获取多个听诊区的心音音频文件;
处理模块,还用于将心音音频文件与对应患者身份进行绑定,且将心音音频文件转换为统一音频格式得到待导入的心音信号样本。
具体的,本实施例是上述方法实施例对应的系统实施例,具体效果参见上述系统实施例,在此不再一一赘述。
基于前述实施例,导入转换模块包括:
提取单元,用于提取所导入心音信号样本的包络特征以及特征参数;特征参数包括:频率、第一心音、第二心音、收缩期以及舒张期的持续时间;
分割单元,用于根据包络特征以及特征参数自动分割得到五个瓣膜区,后根据瓣膜区、包络特征以及特征参数分割得到各瓣膜区对应的多个心动周期;
生成单元,用于根据心动周期进行区间划分并生成心音图谱;心音图谱包括心动周期、第一心音、第二心音的区间位置。
具体的,本实施例是上述方法实施例对应的系统实施例,具体效果参见上述系统实施例,在此不再一一赘述。
基于前述实施例,所述处理模块包括:
播放单元,同进度播放各区间位置对应的心音图谱和心音信号样本;
显示单元,在获取到触发标记信息后弹窗显示心音评价项及其对应的标记选项按钮;
获取单元,用于获取医生根据心音的响度变化针对各心音评价项所输入的标记信息;还用于获取播放调节指令;
标注单元,根据所述标记信息将对应区间范围进行标注;所述标记信息包括医生所选中的标记选项按钮;
调节单元,用于在播放期间根据获取到的播放调节指令同步调整图谱滚动播放状态和音频语音播放状态。
具体的,本实施例是上述方法实施例对应的系统实施例,具体效果参见上述系统实施例,在此不再一一赘述。
基于前述实施例,处理模块,还用于获取到标记结束指令以停止标注功能后,将标注后的心音图谱导出。
具体的,本实施例是上述方法实施例对应的系统实施例,具体效果参见上述系统实施例,在此不再一一赘述。
所属领域的技术人员可以清楚地了解到,为了描述的方便和简洁,仅以上述各程序模块的划分进行举例说明,实际应用中,可以根据需要而将上述功能分配由不同的程序模块完成,即将所述装置的内部结构划分成不同的程序单元或模块,以完成以上描述的全部或者部分功能。实施例中的各程序模块可以集成在一个处理单元中,也可是各个单元单独物理存在,也可以两个或两个以上单元集成在一个处理单元中,上述集成的单元既可以采用硬件的形式实现,也可以采用软件程序单元的形式实现。另外,各程序模块的具体名称也只是为了便于相互区分,并不用于限制本申请的保护范围。
在上述实施例中,对各个实施例的描述都各有侧重,某个实施例中没有详细描述或记载的部分,可以参见其他实施例的相关描述。
本领域普通技术人员可以意识到,结合本文中所公开的实施例描述的各示例的单元及算法步骤,能够以电子硬件、或者计算机软件和电子硬件的结合来实现。这些功能究竟以硬件还是软件来执行,取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本申请的范围。
应当说明的是,上述实施例均可根据需要自由组合。以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。

Claims (6)

1.一种心音的分析与标注方法,其特征在于,包括步骤:
根据导入的心音信号样本转换得到心音图谱,所述心音图谱包括心音图和心音频响图;
播放所述心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内;
所述根据导入的心音信号样本转换得到心音图谱,所述心音图谱包括心音图和心音频响图包括步骤:
提取所导入心音信号样本的包络特征以及特征参数;所述特征参数包括:频率、第一心音、第二心音、收缩期以及舒张期的持续时间;
根据包络特征以及特征参数自动分割得到五个瓣膜区,后根据所述瓣膜区、包络特征以及特征参数分割得到各瓣膜区对应的多个心动周期;
根据心动周期进行区间划分并生成心音图谱;所述心音图谱包括心动周期、第一心音、第二心音的区间位置;
所述播放所述心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内包括步骤:
同进度播放各区间位置对应的心音图谱和心音信号样本;
在获取到触发标记信息后弹窗显示心音评价项及其对应的标记选项按钮;
获取医生根据心音的响度变化针对各心音评价项所输入的标记信息,根据所述标记信息将对应区间范围进行标注;所述标记信息包括医生所选中的标记选项按钮。
2.根据权利要求1所述的心音的分析与标注方法,其特征在于,所述根据导入的心音信号样本转换得到心音图谱,所述心音图谱包括心音图和心音频响图之前包括步骤:
采集获取多个听诊区的心音音频文件;
将所述心音音频文件与对应患者身份进行绑定,且将所述心音音频文件转换为统一音频格式得到待导入的心音信号样本。
3.根据权利要求2所述的心音的分析与标注方法,其特征在于,还包括步骤:
在播放期间根据获取到的播放调节指令同步调整图谱滚动播放状态和音频语音播放状态。
4.根据权利要求1-3任一项所述的心音的分析与标注方法,其特征在于,所述播放所述心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内之后包括步骤:
获取到标记结束指令以停止标注功能后,将标注后的心音图谱导出。
5.一种心音的分析与标注系统,其特征在于,包括:
导入转换模块,用于根据导入的心音信号样本转换得到心音图谱,所述心音图谱包括心音图和心音频响图;
处理模块,用于播放所述心音图谱并获取输入的标记信息,将心音分析结果标注在对应区间范围内;
所述导入转换模块包括:提取单元,用于提取所导入心音信号样本的包络特征以及特征参数;所述特征参数包括:频率、第一心音、第二心音、收缩期以及舒张期的持续时间;
分割单元,用于根据包络特征以及特征参数自动分割得到五个瓣膜区,后根据所述瓣膜区、包络特征以及特征参数分割得到各瓣膜区对应的多个心动周期;
生成单元,用于根据心动周期进行区间划分并生成心音图谱;所述心音图谱包括心动周期、第一心音、第二心音的区间位置;
所述处理模块包括:
播放单元,同进度播放各区间位置对应的心音图谱和心音信号样本;
显示单元,在获取到触发标记信息后弹窗显示心音评价项及其对应的标记选项按钮;
获取单元,用于获取医生根据心音的响度变化针对各心音评价项所输入的标记信息;还用于获取播放调节指令;
标注单元,根据所述标记信息将对应区间范围进行标注;所述标记信息包括医生所选中的标记选项按钮;
调节单元,用于在播放期间根据获取到的播放调节指令同步调整图谱滚动播放状态和音频语音播放状态。
6.根据权利要求5所述的心音的分析与标注系统,其特征在于,还包括:
采集模块,用于采集获取多个听诊区的心音音频文件;
所述处理模块,还用于将所述心音音频文件与对应患者身份进行绑定,且将所述心音音频文件转换为统一音频格式得到待导入的心音信号样本。
CN202010598567.9A 2020-06-28 2020-06-28 一种心音的分析与标注方法和系统 Active CN111714151B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010598567.9A CN111714151B (zh) 2020-06-28 2020-06-28 一种心音的分析与标注方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010598567.9A CN111714151B (zh) 2020-06-28 2020-06-28 一种心音的分析与标注方法和系统

Publications (2)

Publication Number Publication Date
CN111714151A CN111714151A (zh) 2020-09-29
CN111714151B true CN111714151B (zh) 2023-04-14

Family

ID=72569375

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010598567.9A Active CN111714151B (zh) 2020-06-28 2020-06-28 一种心音的分析与标注方法和系统

Country Status (1)

Country Link
CN (1) CN111714151B (zh)

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5025809A (en) * 1989-11-28 1991-06-25 Cardionics, Inc. Recording, digital stethoscope for identifying PCG signatures
EP2014235A1 (en) * 2007-07-09 2009-01-14 Medicalgorithmics Sp. Z.O.O. Method, device and system for cardio-accoustic signal analysis
CN102271589A (zh) * 2008-12-30 2011-12-07 皇家飞利浦电子股份有限公司 用于处理心音信号的方法和系统
CN102283670B (zh) * 2011-07-08 2013-07-03 西华大学 一种心音信号量化分析方法及装置
CN108523920B (zh) * 2018-02-26 2023-08-22 河南善仁医疗科技有限公司 基于多路心音的心脏杂音类型的识别方法
US20190362022A1 (en) * 2018-05-25 2019-11-28 Risto Haukioja Audio file labeling process for building datasets at scale

Also Published As

Publication number Publication date
CN111714151A (zh) 2020-09-29

Similar Documents

Publication Publication Date Title
Varghees et al. Effective heart sound segmentation and murmur classification using empirical wavelet transform and instantaneous phase for electronic stethoscope
CN108647565B (zh) 一种基于深度学习模型对心电信号分类的数据预处理方法
Thiyagaraja et al. A novel heart-mobile interface for detection and classification of heart sounds
JP4256683B2 (ja) 心臓診断システム
US7458939B2 (en) Procedure for extracting information from a heart sound signal
US20040260188A1 (en) Automated auscultation system
US20110257548A1 (en) Method and system for processing heart sound signals
Lee et al. Comparison between short time Fourier and wavelet transform for feature extraction of heart sound
Haghighi-Mood et al. A sub-band energy tracking algorithm for heart sound segmentation
Javed et al. A signal processing module for the analysis of heart sounds and heart murmurs
CN113076846B (zh) 心音分类识别方法及系统
Brusco et al. Development of an intelligent PDA-based wearable digital phonocardiograph
Debbal et al. Analysis of the second heart sound using continuous wavelet transform
Mondal et al. Boundary estimation of cardiac events S1 and S2 based on Hilbert transform and adaptive thresholding approach
CN106037792B (zh) 一种心音波形数据处理和图形显示方法及其听诊系统
CN112494000B (zh) 一种基于多路心音的冠状动脉狭窄可视化量化方法及设备
CN111714151B (zh) 一种心音的分析与标注方法和系统
Ghaffari et al. Phonocardiography signal processing for automatic diagnosis of ventricular septal defect in newborns and children
Kamarulafizam et al. Classification of heart sound based on multipoint auscultation system
Tan et al. Real time analysis of fetal phonocardiography
Behbahani A hybrid algorithm for heart sounds segmentation based on phonocardiogram
WO2016207951A1 (ja) シャント音解析装置、シャント音解析方法、コンピュータプログラム及び記録媒体
Ari et al. On a robust algorithm for heart sound segmentation
CN110811673A (zh) 基于概率神经网络模型的心音分析系统
US20220280132A1 (en) Medical device and program

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
TA01 Transfer of patent application right

Effective date of registration: 20210617

Address after: 200050 room E-11, 17th floor, 8 Xingyi Road, Changning District, Shanghai

Applicant after: Sound smart technology (Shanghai) Co.,Ltd.

Applicant after: XIN HUA HOSPITAL AFFILIATED TO SHANGHAI JIAO TONG University SCHOOL OF MEDICINE

Address before: 200050 room 5v-5396, building 2, 753 Yuyuan Road, Changning District, Shanghai

Applicant before: Sound smart technology (Shanghai) Co.,Ltd.

TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20220801

Address after: 100080 room 21, floor 4, building 2, yard a 2, West Third Ring North Road, Haidian District, Beijing

Applicant after: Beijing Yingtong Technology Development Co.,Ltd.

Address before: 200050 room E-11, 17th floor, 8 Xingyi Road, Changning District, Shanghai

Applicant before: Sound smart technology (Shanghai) Co.,Ltd.

Applicant before: XIN HUA HOSPITAL AFFILIATED TO SHANGHAI JIAO TONG University SCHOOL OF MEDICINE

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant