专利名称：一种结合心音的超声频谱包络峰识别的方法
技术领域：
本发明涉及一种多普勒超声频谱的分析和计算方法，包括了多普勒频谱的自动包
络、在心音指导下的准确有效包络峰的识别和提取、有效峰的心功能参数的计算，属于医学 信号与图像处理技术领域。
背景技术：
超声多普勒心排量检测技术是监测心输出量的无创方法，这种方法与已有技术 (Swan-Ganz导管、PiCC0连续心排量监测)相比，具有直观、准确、无创、灵敏及重复性好等 特点。采用相关技术的现有设备包括心功能彩超、黑白超、小型便携超声检测设备等。
超声多普勒心排量检测技术实时地以动态图像显示心排量的变化，具有直观的效 果。并且超声多普勒是一项无创检测技术，对人体副作用很小。超声多普勒检测技术已比较 成熟，准确度得到国际广泛认可。其灵敏度也相当高，可观察到药物作用对心排量的变化。
心音信号是人体最重要的生理信号之一，它含有关于心脏各个部分如心房、心室 及各个瓣膜功能状态的大量的生理病理信息，在评价心血管系统功能和疾病方面具有重要 的实用价值。心音还可以用来做信号同步，例如国内外部分动脉硬化检测仪加入了心音信 号，给人体四肢脉搏信号作参考。心音信号无论是自身的多项参数还是为超声心功能信号 作同步都具有重要意义。 目前还没有一种方法将心音信号与超声多普勒心排量检测技术结合进行分析和 计算的方法。
发明内容
本发明的目的是提供一种将心音信号与超声多普勒心排量检测技术结合来识别 超声频谱包络峰的方法。 为了达到上述目的，本发明的技术方案是提供了一种结合心音的超声频谱包络峰
识别的方法，其特征在于，步骤为 步骤1、结合心音的超声多普勒频谱分析 步骤1. 1、在超声探头采集信号的同时采集心音信号，对于实时信号，每次从超声 硬件模块得到一段新的IQ信号超声数据，在计算机中计算多普勒频谱，过程是先对IQ信号 做傅立叶变换，然后计算幅值，得到一维的频谱数组； 步骤1. 2、对频谱数组做包络计算，得到当前频谱数组的包络点，为了时间域上包 络的连贯性，需将当前频谱数组的包络点结合时间域上前5个包络点做平滑处理，得到新 的平滑包络点； 步骤1.3、在处理超声数据的同时处理心音数据，在频域上使用起始频率约为
180Hz的布莱克曼窗高通滤波器将采集到的心音数据做滤波，得到心音波形； 步骤1. 4、对通过步骤1. 2得到的至少50个平滑包络点作为一个点集做波峰识别，
通过值比较先找到点集中的最大值点，若最大值点前的所有点使用最小二乘法直线拟合后所得直线斜率大于l，最大值后的点使用最小二乘直线拟合所得直线斜率大于-l，则认为 识别到波峰，如未识别到波峰则返回步骤1. 2计算下一个时间域上的频谱数组的平滑包络 点的点集，如识别到波峰，则进入下一步； 步骤1. 5、对步骤1. 4识别到的超声频谱的波峰和步骤1. 3得到的心音波形的波峰
做时间差计算，采用峰起点做时间差计算，或采用峰顶点做时间差计算，若最终得到的时间
差小于0. 2ms范围，则认为该识别到的超声频谱的波峰是有效峰，否则认为是伪峰； 步骤2、计算参数 步骤2. 1 、计算血液流速公式为V = ^fe; 其中，c是超声在人体中的声速，fd是超声发射频率和接收频率的频移，f是发射 频率，e为探头与血管的角度； 步骤2. 2、计算心搏距离心搏距离由步骤1. 5得到的超声频谱的有效峰的积分面
积(速度时间积分)得到，公式为 其中，s为心搏距离，^为波峰起始时刻，t2为波峰结束时刻，v(t)为速度随时间 变化函数，在图像上反映即为包络；
步骤2. 3、计算心排量，公式为
co = s hr xsa ; 其中，hr为心率，由60秒和单个心音波形的波峰的时间比值得到，xsa是人体动脉 血管的截面积，co为输出的心排量。 同时采用多普勒超声探头，及心音探头，在人体相关部位同步采集心排量信号和 心音信号，两者信号输入到计算节点做处理。计算节点将心排量信号和心音信号相结合自 动分析出有效峰并包络，然后计算相关心排量参数和心音参数。
与现有技术相比，本发明具有以下优点 1、超声多普勒信号与心音信号同步采集，通过在时间轴上峰位置的比较，可以有 效提取心排量频谱峰，滤除因探头摩擦或人体运动引起的伪峰，提高包络识别和心排量参 数计算的准确度； 2、通过在时间轴上比较超声多普勒频谱峰和心音波峰的时差分析，可以有效地判 别超声探头放置位置的正确与否，并可以根据医学仪器屏幕上反映出的时差变化来轻微调 整探头的方向，该方法直观清晰，结合超声仪器常用的声音判别方向法一起使用，更好地指 导使用者对心排量检测的操作； 3、心音在心血管、心房、心室、大血管及各个瓣膜功能状态方面具有重要的参考价 值，结合心音与心排量的参数来综合分析人体心脏功能具有相当重要的医学价值；
4、多普勒超声和心音同样都是无创检测方法，正常的检测对人体基本无副作用， 并且操作简单方便，亦可做监护使用。
图1为识别到有效波峰时所显示的多普勒频谱及包络图像；
图2为未识别到有效波峰时的图像； 图3a为频谱数组的值变化图； 图3b为频谱数组值做累积和之后得到的图。
具体实施例方式
以下结合实施例来具体说明本发明。
实施例 本发明提供了的一种结合心音的超声频谱包络峰识别的方法，其步骤为
步骤1、结合心音的超声多普勒频谱分析 步骤1. 1、在超声探头采集信号的同时采集心音信号。对于实时信号，每次从超声 硬件模块得到一段新的IQ信号超声数据，在计算机中计算多普勒频谱，过程是先对IQ信 号做傅立叶变换，然后计算幅值，得到一维的频谱数组。在该步骤中调试设备的过程为将 心音探头放置人体左胸心脏上方位置，观察仪器屏幕心音信号的波形(如图2上方波形所 示)，待心音波形趋于稳定状态时，固定心音探头位置，具体可采用沙包压住心音探头的方 法来固定。将多普勒超声探头置于人体脖子下方锁骨间的下凹处，中医称之为天突穴。观 察超声多普勒频谱，根据包络线形状调整超声探头的位置，同时关注扬声器发出的由超声 血流信号转换的声音，直至出现清晰且稳定的声音。此处观察频谱图像的波峰及包络线，直 至出现稳定的频谱波峰和包络线，如图1所示。 步骤1. 2、对频谱数组做包络计算，假设频谱数组的值变化如图3a所示，X轴为数 组点索引，Y轴为值，将数组值做累积和得图3b中的曲线，则该频谱数组的包络点为图3b中 在直线上方的曲线的拐点A，得到当前频谱数组的包络点后，为了时间域上包络的连贯性， 需将当前频谱数组的包络点结合时间域上前5个包络点做高斯平滑处理，得到新的平滑包 络点，高斯平滑窗大小为6，与包络点数一致； 步骤1.3、在处理超声数据的同时处理心音数据，在频域上使用起始频率约为 180Hz的布莱克曼窗高通滤波器将采集到的心音数据做滤波，得到如图l及图2上方所示心 音波形； 步骤1. 4、对通过步骤1. 2得到的至少50个平滑包络点作为一个点集做波峰识别， 通过值比较先找到点集中的最大值点，若最大值点前的所有点使用最小二乘法直线拟合后 所得直线斜率大于l，最大值后的点使用最小二乘直线拟合所得直线斜率大于-l，则认为 识别到波峰，如未识别到波峰则返回步骤1. 2计算下一个时间域上的频谱数组的平滑包络 点的点集，如识别到波峰，则进入下一步； 步骤1. 5、对步骤1. 4识别到的超声频谱的波峰和步骤1. 3得到的心音波形的波峰
做时间差计算，采用峰起点或峰顶点做时间差计算，若采用峰顶点做时间差计算，则对心音
波形的波峰需要采用与对超声频谱的波峰相同的方法进行包络计算，若最终得到的时间差
小于0. 2ms范围，则认为该识别到的超声频谱的波峰是有效峰，否则认为是伪峰； 步骤2、计算参数 步骤2. 1 、计算血液流速公式为 V = ^ ;
5
其中，c是超声在人体中的声速，fd是超声发射频率和接收频率的频移，f是发射 频率，e为探头与血管的角度； 步骤2. 2、计算心搏距离心搏距离由步骤1. 5得到的超声频谱的有效峰的积分面 积(速度时间积分)得到，公式为 = f2 ; 其中，s为心搏距离，^为波峰起始时刻，t2为波峰结束时刻，v(t)为速度随时间
变化函数，在图像上反映即为包络；
步骤2. 3、计算心排量，公式为
co 二 s hr xsa ; 其中，hr为心率，由60秒和单个心音波形的波峰的时间比值得到，xsa是人体动脉 血管的截面积，co为输出的心排量。
权利要求
一种结合心音的超声频谱包络峰识别的方法，其特征在于，步骤为步骤1、结合心音的超声多普勒频谱分析步骤1.1、在超声探头采集信号的同时采集心音信号，对于实时信号，每次从超声硬件模块得到一段新的IQ信号超声数据，在计算机中计算多普勒频谱，过程是先对IQ信号做傅立叶变换，然后计算幅值，得到一维的频谱数组；步骤1.2、对频谱数组做包络计算，得到当前频谱数组的包络点，为了时间域上包络的连贯性，需将当前频谱数组的包络点结合时间域上前至少5个包络点做平滑处理，得到新的平滑包络点；步骤1.3、在处理超声数据的同时处理心音数据，在频域上使用起始频率约为180Hz的布莱克曼窗高通滤波器将采集到的心音数据做滤波，得到心音波形；步骤1.4、对通过步骤1.2得到的至少50个平滑包络点作为一个点集做波峰识别，通过值比较先找到点集中的最大值点，若最大值点前的所有点使用最小二乘法直线拟合后所得直线斜率大于1，最大值后的点使用最小二乘直线拟合所得直线斜率大于-1，则认为识别到波峰，如未识别到波峰则返回步骤1.2计算下一个时间域上的频谱数组的平滑包络点的点集，如识别到波峰，则进入下一步；步骤1.5、对步骤1.4识别到的超声频谱的波峰和步骤1.3得到的心音波形的波峰做时间差计算，采用峰起点或峰顶点做时间差计算，若最终得到的时间差小于0.2ms范围，则认为该识别到的超声频谱的波峰是有效峰，否则认为是伪峰；步骤2、计算参数步骤2.1、计算血液流速公式为 <mrow><mi>v</mi><mo>=</mo><mfrac> <mrow><msub> <mi>f</mi> <mi>d</mi></msub><mo>·</mo><mi>c</mi> </mrow> <mrow><mn>2</mn><mi>f</mi><mo>·</mo><mi>cos</mi><mi>θ</mi> </mrow></mfrac><mo>;</mo> </mrow>其中，c是超声在人体中的声速，fd是超声发射频率和接收频率的频移，f是发射频率，θ为探头与血管的角度；步骤2.2、计算心搏距离心搏距离由步骤1.5得到的超声频谱的有效峰的积分面积(速度时间积分)得到，公式为 <mrow><mi>s</mi><mo>=</mo><msubsup> <mo>∫</mo> <msub><mi>t</mi><mn>1</mn> </msub> <msub><mi>t</mi><mn>2</mn> </msub></msubsup><mi>v</mi><mrow> <mo>(</mo> <mi>t</mi> <mo>)</mo></mrow><mi>dt</mi><mo>;</mo> </mrow>其中，s为心搏距离，t1为波峰起始时刻，t2为波峰结束时刻，v(t)为速度随时间变化函数，在图像上反映即为包络；步骤2.3、计算心排量，公式为co＝s·hr·xsa；其中，hr为心率，由60秒和单个心音波形的波峰的时间比值得到，xsa是人体动脉血管的截面积，co为输出的心排量。
2.如权利要求1所述的一种结合心音的超声频谱包络峰识别的方法，其特征在于，步 骤1. 5中若采用峰顶点做时间差计算，则对心音波形的波峰需要采用与对超声频谱的波峰 相同的方法进行包络计算。
全文摘要
本发明提供了一种结合心音的超声频谱包络峰识别的方法，其步骤为同时采用多普勒超声探头及心音探头，在人体相关部位同步采集心排量信号和心音信号，两者信号输入到计算机做处理。计算机将心排量信号和心音信号相结合自动分析出有效峰并包络，然后计算心排量参数和心音参数。本发明能更加有效的提取真实峰，同时通过对有效包络峰的计算得出能够反映心功能的相关参数。本发明亦能够有效的得出频谱包络峰的多项参数，能够真实的反映人体心功能的基本生理状况和变化。
文档编号A61B8/00GK101779966SQ201010113860
公开日2010年7月21日 申请日期2010年2月25日 优先权日2010年2月25日
发明者任春剑, 方志强, 马扬喜 申请人:上海三埃弗电子有限公司