CN104688213A - 一种心音信号分段方法 - Google Patents
一种心音信号分段方法 Download PDFInfo
- Publication number
- CN104688213A CN104688213A CN201510081867.9A CN201510081867A CN104688213A CN 104688213 A CN104688213 A CN 104688213A CN 201510081867 A CN201510081867 A CN 201510081867A CN 104688213 A CN104688213 A CN 104688213A
- Authority
- CN
- China
- Prior art keywords
- spacing
- envelope
- cycle
- heart sound
- less
- 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.)
- Granted
Links
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/02—Detecting, measuring or recording pulse, heart rate, blood pressure or blood flow; Combined pulse/heart-rate/blood pressure determination; Evaluating a cardiovascular condition not otherwise provided for, e.g. using combinations of techniques provided for in this group with electrocardiography or electroauscultation; Heart catheters for measuring blood pressure
- A61B5/024—Detecting, measuring or recording pulse rate or heart rate
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7203—Signal processing specially adapted for physiological signals or for diagnostic purposes for noise prevention, reduction or removal
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/72—Signal processing specially adapted for physiological signals or for diagnostic purposes
- A61B5/7221—Determining signal validity, reliability or quality
Abstract
本发明公开了一种心音信号分段方法,其特征在于,包括以下步骤:第一步:预处理;第二步:校正;第三步:筛选。本发明的心音信号分段方法不仅步骤简单,而且便于使用,本发明通过包含多层次的校正和筛选的分段方法,来确保能够获得正确的心音分段,从而为心率检测提供更为便捷的方法,因此适合推广使用。
Description
技术领域
本发明涉及一种心率检测方法,尤其涉及一种心音信号分段方法。
背景技术
心音是指由心肌收缩、心脏瓣膜关闭和血液撞击心室壁、大动脉壁等引起的振动所产生的声音。它可在胸壁一定部位用听诊器听取;现在也可以用电子听诊器听取,从而转化为数字信号,处理后进行识别。
心音信号作为一种不稳定的周期循环信号,对其进行分段从而获得单独的周期,主要有两个目的:1.能够获得心率;2.以周期为基本单位进行识别、诊断。因此,对心音信号进行分段一直是心音信号处理中的必需环节。
作为心音信号处理的整个流程的第一个环节,分段的质量直接关系到能否输出正确的心率和能否进行正确的诊断。但是由于心音信号的来源以及采集方式,实际所获得的心音信号通常是高噪声、不稳定的,这给分段带来了极大难度。而目前已有的相关方法存在以下问题:
1.在处理高噪声的信号时,由于噪声幅值较大抬高阈值引起漏检或者心音信号包络之间夹杂了大量噪声包络,则会出现大量错误分段;
2.处理一些疾病信号或者特殊的信号时,由于其心音包络形状大小非典型或者心音信号的幅值处于不稳定中,会出现错误;
3.一些采用自相关的分段,计算量较大,移植到移动平台时运算时间过长。因此在实际应用中,一般只能应用于低噪声的理想信号,或者对采集的环境和设备要求苛刻。
仅仅通过简单的校正,在高噪声下无法获得正确的分段。因此构造一种包含多层次的校正和筛选的分段方法,来确保能够获得正确的心音分段显得尤为重要。
发明内容
本发明的目的在于提供一种心音信号分段方法,以期望构造一种包含多层次的校正和筛选的分段方法,来确保能够获得正确的心音分段。
本发明的目的通过下述技术方案实现:
1.一种心音信号分段方法,包括以下步骤:
第一步:预处理
(1)获得心音信号,对获得的心音信号进行滤波;
(2)计算提取信号的THres阀值,其中peak为经过第一步的步骤(1)所取的波峰幅值,N为波峰数量,T为系数;T取决于信号长度与信号波峰的离散度;T=0.2*(1-(40/t)*σ),其中σ为peak的标准差,t为信号时长(5<t<40秒);
(3)使用THres阈值提取心音信号,大于THres阈值的保留,小于THres阈值的去除;
第二步:校正
(1)合并间距小于硬阈值的包络
通过interval(i)=envelop(i+1)–envelop(i)计算相邻的两个包络之间的间距,其中envelop为包络起点位置,i为包络序号;
根据间距和频率Fs计算时长,interval(i)/Fs,若时长小于0.2秒,则将构成该间距的两个包络合并;
(2)去除导致间距小于软阈值的包络
计算合并后两个包络之间的间距,然后根据间距均值计算阈值:
Threshold=K*Σinterval/n,其中,interval为合并后的间距,n为间距数量,K为系数(0<K<0.5);
若某一间距小于该阈值,则判断它为错误间距;
判断该间距相邻的两个间距的大小关系:若前一个间距大于后一个间距,则去掉组成该错误间距的前一个包络;反之,则去掉组成该错误间距的后一个包络;
(3)检查漏检
计算经过合并、去除之后的包络间距;
判断相邻两个包络间距的大小是否大致相等:若大致相等,说明可能出现第二心音漏检,则保留形成这两个间距的三个包络;若大小有明显差异,则根据大小关系,区分包络属于第一心音还是第二心音,只保留第一心音的包络;
第三步筛选
(1)根据所获得的第一心音包络,计算周期:
cycle(i)=(S1envelop(i+1)–S1envelop(i))/Fs,其中S1envelop为第一心音包络的起点,i为包络序列;
(2)将大于2秒或者小于0.3秒的周期判断为错误的周期,并去除;
(3)把剩下的周期按长短排序,获得周期中位数,只保留长短在周期中位数附近的周期;
(4)所剩下周期的起点和终点即为每个分段的起点和终点。
进一步的技术方案是,第一步的步骤2中计算提取信号的THres阀值时当相邻的两个波峰之间的间距小于0.05秒时,只取较大的波峰幅值。
本发明较现有技术相比,具有以下优点及有益效果:
本发明的心音信号分段方法不仅步骤简单,而且便于使用,本发明通过包含多层次的校正和筛选的分段方法,来确保能够获得正确的心音分段,从而为心率检测提供更为便捷的方法。
具体实施方式
下面结合实施例对本发明作进一步地详细说明:
实施例
本发明的一种心音信号分段方法,通过对一定时长的心音信号(通常大于10秒)进行分段;每个分段为一心音周期。本发明分段方法通过多次校正,最后筛选出正确的心音周期,用于计算心率与自动诊断。具体,本发明的分段方法分为三个步骤,包括预处理、校正和筛选。
第一步:预处理
(1)首先获得心音信号,并对获得的心音信号进行滤波。
(2)计算提取信号的THres阀值,当相邻的两个波峰之间的间距小于0.05秒时,只取较大的波峰幅值。所述的THres阀值通过以下计算公式进行计算:
其中peak为经过第一步的步骤(1)所取的波峰幅值,N为波峰数量,T为系数;T取决于信号长度与信号波峰的离散度;T=0.2*(1-(40/t)*σ),其中σ为peak的标准差,t为信号时长(5<t<40秒)。
(3)使用THres阈值提取心音信号,大于THres阈值的保留,小于THres阈值的去除。
第二步:校正
(1)合并间距小于硬阈值的包络
通过interval(i)=envelop(i+1)–envelop(i)计算相邻的两个包络之间的间距,其中envelop为包络起点位置,i为包络序号。然后根据间距和频率Fs计算时长,interval(i)/Fs,若时长小于0.2秒,则将构成该间距的两个包络合并。
(2)去除导致间距小于软阈值的包络
计算合并后两个包络之间的间距,然后根据间距均值计算阈值:Threshold=K*Σinterval/n,其中,interval为合并后的间距,n为间距数量,K为系数(0<K<0.5)。若某一间距小于该阈值,则判断它为错误间距。判断该间距相邻的两个间距的大小关系:若前一个间距大于后一个间距,则去掉组成该错误间距的前一个包络;反之,则去掉组成该错误间距的后一个包络。
(3)检查漏检
计算经过合并、去除之后的包络间距;判断相邻两个包络间距的大小是否大致相等:若大致相等,说明可能出现第二心音漏检,则保留形成这两个间距的三个包络;若大小有明显差异,则根据大小关系,区分包络属于第一心音还是第二心音,只保留第一心音的包络。
第三步筛选
(1)根据所获得的第一心音包络,计算周期:
cycle(i)=(S1envelop(i+1)–S1envelop(i))/Fs,其中S1envelop为第一心音包络的起点,i为包络序列。
(2)将大于2秒或者小于0.3秒的周期判断为错误的周期,并去除。
(3)把剩下的周期按长短排序,获得周期中位数,只保留长短在周期中位数附近的周期。
(4)所剩下周期的起点和终点即为每个分段的起点和终点。
本发明给出了一种完全基于时域上的判决条件,对心音进行正确分段的方法。一个较小的阈值过滤包络,解决了高噪声背景和不稳定的信号所带来的问题,尽可能保留了心音信号包络,但同时带入了大量噪声包络。根据心音信号本身的特点,设立了诸多判断条件,来区分噪声或者舍去错误的结果,最终达到获得正确分段结果的目的。本方法的计算量小,非常适用于移植到各种平台的实际应用。
以上所述仅为本发明的较佳实施例而已,并不用以限制本发明,凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明的保护范围之内。
Claims (2)
1.一种心音信号分段方法,其特征在于,包括以下步骤:
第一步:预处理
(1)获得心音信号,对获得的心音信号进行滤波;
(2)计算提取信号的THres阀值,其中peak为经过第一步的步骤(1)所取的波峰幅值,N为波峰数量,T为系数;T取决于信号长度与信号波峰的离散度;T=0.2*(1-(40/t)*σ),其中σ为peak的标准差,t为信号时长(5<t<40秒);
(3)使用THres阈值提取心音信号,大于THres阈值的保留,小于THres阈值的去除;
第二步:校正
(1)合并间距小于硬阈值的包络
通过interval(i)=envelop(i+1)–envelop(i)计算相邻的两个包络之间的间距,其中envelop为包络起点位置,i为包络序号;
根据间距和频率Fs计算时长,interval(i)/Fs,若时长小于0.2秒,则将构成该间距的两个包络合并;
(2)去除导致间距小于软阈值的包络
计算合并后两个包络之间的间距,然后根据间距均值计算阈值:
其中,interval为合并后的间距,n为间距数量,K为系数(0<K<0.5);
若某一间距小于该阈值,则判断它为错误间距;
判断该间距相邻的两个间距的大小关系:若前一个间距大于后一个间距,则去掉组成该错误间距的前一个包络;反之,则去掉组成该错误间距的后一个包络;
(3)检查漏检
计算经过合并、去除之后的包络间距;
判断相邻两个包络间距的大小是否大致相等:若大致相等,说明可能出现第二心音漏检,则保留形成这两个间距的三个包络;若大小有明显差异,则根据大小关系,区分包络属于第一心音还是第二心音,只保留第一心音的包络;
第三步筛选
(1)根据所获得的第一心音包络,计算周期:
cycle(i)=(S1envelop(i+1)–S1envelop(i))/Fs,其中S1envelop为第一心音包络的起点,i为包络序列;
(2)将大于2秒或者小于0.3秒的周期判断为错误的周期,并去除;
(3)把剩下的周期按长短排序,获得周期中位数,只保留长短在周期中位数附近的周期;
(4)所剩下周期的起点和终点即为每个分段的起点和终点。
2.按照权利要求1所述的一种心音信号分段方法,其特征在于:第一步的步骤2中计算提取信号的THres阀值时当相邻的两个波峰之间的间距小于0.05秒时,只取较大的波峰幅值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510081867.9A CN104688213B (zh) | 2015-02-15 | 2015-02-15 | 一种心音信号分段方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510081867.9A CN104688213B (zh) | 2015-02-15 | 2015-02-15 | 一种心音信号分段方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104688213A true CN104688213A (zh) | 2015-06-10 |
CN104688213B CN104688213B (zh) | 2017-02-15 |
Family
ID=53336085
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510081867.9A Active CN104688213B (zh) | 2015-02-15 | 2015-02-15 | 一种心音信号分段方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104688213B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105212960A (zh) * | 2015-08-19 | 2016-01-06 | 四川长虹电器股份有限公司 | 心音信号质量评估方法 |
CN105232023A (zh) * | 2015-10-10 | 2016-01-13 | 四川长虹电器股份有限公司 | 一种胎心音去噪的方法 |
CN106037704A (zh) * | 2016-05-19 | 2016-10-26 | 四川长虹电器股份有限公司 | 一种心音心率计算方法 |
CN107822597A (zh) * | 2017-10-17 | 2018-03-23 | 四川长虹电器股份有限公司 | 基于智能平台心脏活力指数计算及图像化显示系统与方法 |
CN110123371A (zh) * | 2019-05-30 | 2019-08-16 | 四川长虹电器股份有限公司 | 第二心音部分缺失心音信号的心率计算方法 |
CN110123370A (zh) * | 2019-05-30 | 2019-08-16 | 四川长虹电器股份有限公司 | 第一心音部分遗漏心音信号的心率计算方法 |
CN110151156A (zh) * | 2019-04-07 | 2019-08-23 | 西安电子科技大学 | 一种基于窗内自动加权平均算法的胎心率提取方法及系统 |
CN110368020A (zh) * | 2019-07-03 | 2019-10-25 | 北京必安必恒科技发展有限公司 | 一种心音信号预处理方法及装置 |
CN111528900A (zh) * | 2020-05-21 | 2020-08-14 | 广东工业大学 | 基于巴特沃斯滤波器与香农熵法的心音分段方法和装置 |
CN112185545A (zh) * | 2020-09-17 | 2021-01-05 | 苏州美糯爱医疗科技有限公司 | 一种心音段完整性判断模型及心音自适应分段方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4446873A (en) * | 1981-03-06 | 1984-05-08 | Siemens Gammasonics, Inc. | Method and apparatus for detecting heart sounds |
CN102271589A (zh) * | 2008-12-30 | 2011-12-07 | 皇家飞利浦电子股份有限公司 | 用于处理心音信号的方法和系统 |
CN102283670A (zh) * | 2011-07-08 | 2011-12-21 | 西华大学 | 一种心音信号量化分析方法及装置 |
US20130211273A1 (en) * | 2012-02-12 | 2013-08-15 | National Taiwan University | Method and apparatus for heart rate measurement |
CN103948398A (zh) * | 2014-04-04 | 2014-07-30 | 杭州电子科技大学 | 适用于Android系统的心音定位分段方法 |
CN103393435B (zh) * | 2013-07-31 | 2015-01-07 | 深圳市理邦精密仪器股份有限公司 | 一种胎心音信号包络的取得方法及装置 |
-
2015
- 2015-02-15 CN CN201510081867.9A patent/CN104688213B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4446873A (en) * | 1981-03-06 | 1984-05-08 | Siemens Gammasonics, Inc. | Method and apparatus for detecting heart sounds |
CN102271589A (zh) * | 2008-12-30 | 2011-12-07 | 皇家飞利浦电子股份有限公司 | 用于处理心音信号的方法和系统 |
CN102283670A (zh) * | 2011-07-08 | 2011-12-21 | 西华大学 | 一种心音信号量化分析方法及装置 |
US20130211273A1 (en) * | 2012-02-12 | 2013-08-15 | National Taiwan University | Method and apparatus for heart rate measurement |
CN103393435B (zh) * | 2013-07-31 | 2015-01-07 | 深圳市理邦精密仪器股份有限公司 | 一种胎心音信号包络的取得方法及装置 |
CN103948398A (zh) * | 2014-04-04 | 2014-07-30 | 杭州电子科技大学 | 适用于Android系统的心音定位分段方法 |
Non-Patent Citations (4)
Title |
---|
于倩 等: "《基于简单度的包络提取算法在心音分段中的应用》", 《生物医学工程研究》 * |
张磊邦 等: "《心音信号的预处理与包络提取算法研究》", 《生物医学工程学杂志》 * |
赵治栋 等: "《心音自动分段算法研究》", 《航天医学与医学工程》 * |
陈萌辉 等: "《基于信号包络及短时过零率的心音分段算法》", 《北京生物医学工程》 * |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105212960A (zh) * | 2015-08-19 | 2016-01-06 | 四川长虹电器股份有限公司 | 心音信号质量评估方法 |
CN105212960B (zh) * | 2015-08-19 | 2018-03-30 | 四川长虹电器股份有限公司 | 心音信号质量评估方法 |
CN105232023A (zh) * | 2015-10-10 | 2016-01-13 | 四川长虹电器股份有限公司 | 一种胎心音去噪的方法 |
CN106037704A (zh) * | 2016-05-19 | 2016-10-26 | 四川长虹电器股份有限公司 | 一种心音心率计算方法 |
CN107822597A (zh) * | 2017-10-17 | 2018-03-23 | 四川长虹电器股份有限公司 | 基于智能平台心脏活力指数计算及图像化显示系统与方法 |
CN110151156A (zh) * | 2019-04-07 | 2019-08-23 | 西安电子科技大学 | 一种基于窗内自动加权平均算法的胎心率提取方法及系统 |
CN110123370A (zh) * | 2019-05-30 | 2019-08-16 | 四川长虹电器股份有限公司 | 第一心音部分遗漏心音信号的心率计算方法 |
CN110123371A (zh) * | 2019-05-30 | 2019-08-16 | 四川长虹电器股份有限公司 | 第二心音部分缺失心音信号的心率计算方法 |
CN110123371B (zh) * | 2019-05-30 | 2021-06-29 | 四川长虹电器股份有限公司 | 第二心音部分缺失心音信号的心率计算方法 |
CN110123370B (zh) * | 2019-05-30 | 2021-06-29 | 四川长虹电器股份有限公司 | 第一心音部分遗漏心音信号的心率计算方法 |
CN110368020A (zh) * | 2019-07-03 | 2019-10-25 | 北京必安必恒科技发展有限公司 | 一种心音信号预处理方法及装置 |
CN110368020B (zh) * | 2019-07-03 | 2022-03-08 | 北京必安必恒科技发展有限公司 | 一种心音信号预处理方法及装置 |
CN111528900A (zh) * | 2020-05-21 | 2020-08-14 | 广东工业大学 | 基于巴特沃斯滤波器与香农熵法的心音分段方法和装置 |
CN112185545A (zh) * | 2020-09-17 | 2021-01-05 | 苏州美糯爱医疗科技有限公司 | 一种心音段完整性判断模型及心音自适应分段方法 |
CN112185545B (zh) * | 2020-09-17 | 2024-03-12 | 苏州美糯爱医疗科技有限公司 | 一种心音段完整性判断模型及心音自适应分段方法 |
Also Published As
Publication number | Publication date |
---|---|
CN104688213B (zh) | 2017-02-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104688213A (zh) | 一种心音信号分段方法 | |
WO2019119518A1 (zh) | 房颤信号的识别方法、装置和设备 | |
WO2017202292A1 (zh) | 一种回声时延跟踪方法及装置 | |
CN102160296B (zh) | 双端通话检测方法及装置 | |
CN104161510A (zh) | 一种多级导联心电信号qrs波形识别方法 | |
CN103584854A (zh) | 心电信号r波的提取方法 | |
CN108814590B (zh) | 一种心电qrs波群的检测方法及其心电分析方法 | |
CN105448303A (zh) | 语音信号的处理方法和装置 | |
CN105030233B (zh) | 一种心电信号st段的识别方法 | |
CN108814591B (zh) | 一种心电qrs波群宽度的检测方法及其心电分析方法 | |
CN104905785A (zh) | 室性/室上性早搏检测方法和系统 | |
CN107361764B (zh) | 一种心电信号特征波形r波的快速提取方法 | |
CN107274913A (zh) | 一种声音识别方法及装置 | |
EP3413310B1 (en) | Acoustic meaningful signal detection in wind noise | |
CN110853677B (zh) | 歌曲的鼓声节拍识别方法、装置、终端和非临时性计算机可读存储介质 | |
CN115935144A (zh) | 一种运检维护数据去噪重构方法 | |
CN104382589B (zh) | 基于部分按段重采样的胎儿心电图分离提取方法 | |
WO2020228107A1 (zh) | 一种音频修复方法、设备及可读存储介质 | |
CN106571138B (zh) | 一种信号端点的检测方法、检测装置及检测设备 | |
CN115902391A (zh) | 基于形态学滤波和特征值分解的动态门限信号检测方法 | |
CN110123371B (zh) | 第二心音部分缺失心音信号的心率计算方法 | |
CN107548007B (zh) | 一种音频信号采集设备的检测方法及装置 | |
CN109100441B (zh) | 一种液相色谱曲线去除脉冲干扰的方法 | |
CN103278702B (zh) | 提取pace波的智能检测单元、方法及系统 | |
CN105997036A (zh) | 基于心音信号时域特征的心率计算方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |