CN1907230A - 频谱多普勒血流速度的自动检测方法 - Google Patents
频谱多普勒血流速度的自动检测方法 Download PDFInfo
- Publication number
- CN1907230A CN1907230A CNA2005100363720A CN200510036372A CN1907230A CN 1907230 A CN1907230 A CN 1907230A CN A2005100363720 A CNA2005100363720 A CN A2005100363720A CN 200510036372 A CN200510036372 A CN 200510036372A CN 1907230 A CN1907230 A CN 1907230A
- Authority
- CN
- China
- Prior art keywords
- noise
- frequency
- blood flow
- spectral line
- threshold value
- 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
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/06—Measuring blood flow
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F1/00—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow
- G01F1/66—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters
- G01F1/663—Measuring the volume flow or mass flow of fluid or fluent solid material wherein the fluid passes through a meter in a continuous flow by measuring frequency, phase shift or propagation time of electromagnetic or other waves, e.g. using ultrasonic flowmeters by measuring Doppler frequency shift
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/52017—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging
- G01S7/52023—Details of receivers
- G01S7/52025—Details of receivers for pulse systems
- G01S7/52026—Extracting wanted echo signals
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/08—Detecting organic movements or changes, e.g. tumours, cysts, swellings
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B8/00—Diagnosis using ultrasonic, sonic or infrasonic waves
- A61B8/48—Diagnostic techniques
- A61B8/488—Diagnostic techniques involving Doppler signals
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
- G01S15/8906—Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
- G01S15/8979—Combined Doppler and pulse-echo imaging systems
Abstract
一种频谱多普勒血流速度的自动检测方法,用于超声系统对血流速度的测量过程,包括步骤:射频超声回波信号经过解调、滤波和模数转换后获得流体的多普勒信号;对所述多普勒信号进行频谱分析,获得该多普勒信号随时间变化的各功率谱线;确定阈值;根据所述阈值和当前功率谱线来确定该功率谱线对应的频偏参数或血流速度;其中,所述阈值与前一基本时间段内的平均噪声水平相关,并自动根据所述平均噪声水平被更新调整。所述方法确定频偏参数前,还可以包括一判断该谱线是否有效的过程及一平滑所述功率谱的过程。采用本发明方法,可以使阈值自动收敛,从而便于系统准确并稳健地进行包络检测,自动实时显示参数计算结果。
Description
技术领域 本发明涉及超声技术,特别涉及利用频谱多普勒效应测量流体速度的超声技术,尤其是涉及超声诊断仪利用包络检测来获得血流速度的方法。
背景技术 在超声成像系统中,多普勒效应被广泛用于血流检测。例如在应用超声多普勒技术检查心脏与动静脉的超声系统中,需要从多普勒频谱图中提取相关参数来评价心脏与血管的血流动力学状态。该参数包括:峰值速度,舒张最低速度,阻力指数,收缩/舒张指数等,这些参数的确定许多都依赖于血流的峰值速度和平均速度。而对血流峰值速度和平均速度的检测可以转换成对超声回波最大频偏和平均频偏的检测。其原理是,超声换能器发射超声波到受试者体内,并接收从受试者体内返回的超声回波;由于血管中运动红细胞对超声波的散射作用,使得超声回波的频率发生偏移,而该频率偏移大小与发射超声波的中心频率及红细胞的运动速度相关,故只要测量超声回波频率偏移,就可以计算获知红细胞的运动情况。
图1所示为传统超声多普勒频谱分析系统处理信号的基本流程。超声回波信号经波束合成后形成射频回波信号,又经解调模块分解成两路分量,一路是同相位分量I(In-phasecomponent),另一路是正交分量Q(quadrate component)。在脉冲多普勒信号处理系统中,系统分别对所述两路I、Q分量作距离选通,即在特定的时间段内分别对两路分量信号进行累加(累加的时间段和脉冲多普勒发射脉冲的长度由操作者根据实际情况选择);而连续波多普勒信号处理系统则不需该距离选通过程。所述I、Q两路分量(对脉冲多普勒信号处理系统则指经过距离累加后的两路分量)再分别经壁滤波(高通滤波)处理,滤除由静止或者慢速运动组织引起的杂波,得到以红细胞运动引起的回波为主的两路分量信号,并被送往谱估计模块。所述谱估计模块一般用快速傅立叶变换(FFT)来估算功率谱,该快速傅立叶变换的点数可以是128点,也可以是256点。由于估算出来的功率谱动态范围太大,系统需要将每次估算出来的功率谱进行压缩处理,使之落入灰度显示范围,并在屏幕上的多普勒频谱图上显示出对应时刻和速度(频偏)的功率谱强度。
因血流峰值速度的检测对应于多普勒最大频移的检测,所以早期利用超声多普勒频谱图来手动检测血流峰值速度的方法包括步骤:先保存几个心动周期的功率谱图;接着操作者根据显示屏上的功率谱图来手动标记谱的峰值;最后计算机根据所述标记的峰值来计算出各种参数。这种手动检测的缺点是:操作者标记峰值速度比较单调费时,重复性差,估计精度低;并且在检测时,为了标记峰值速度,操作者需要中断采集多普勒信号,从而无法实时估计。
因而,图1的系统还须包括自动包络检测模块,可以自动跟踪血流峰值速度和平均速度与时间相关的变化,并在多普勒频谱图上实时显示。该自动检测可以在谱估计之后对估算的功率谱进行,也可以如图1所示对压缩后的功率谱进行。这样医生可以实时自动获取血流的峰值速度,平均速度及其它相关参数。
从原理上来说,检测多普勒血流峰值速度似乎比较简单,只要检测最大多普勒频移就可以了。实际上该检测受两个主要因素影响:一是声谱的内在变宽,因为估算功率谱的数据点有限(如128或256点),所以估算出的频谱带宽比理想情况要宽,在整个截止频率带内均有功率谱分布;该因素对多普勒最大频移的影响不好定量判断。二是多普勒信号本身含有的噪声(noise),在多普勒频谱图中总包含有信号和噪声,检测最大频移即找出二者之间的转折点,但频谱由噪声到信号的转折点不是很明显,多普勒频移在一个心动周期的心脏射血期间变化很快,但在其他期间则变化缓慢;而且,对特定频率,频谱图的信噪比也随时间而变化。为此,先后有很多自动检测血流峰值速度和平均速度的方法被提出来力图克服上述影响。
为了能准确并稳健地估算出多普勒频谱图的包络,需要区分频谱图上多普勒谱信号和噪声的边界,所以方法之一是设定阈值,将频谱图上大于阈值的视为信号,反之视为噪声。该阈值可以被设定为所有信号和噪声之和的某个固定百分比,在信噪比(signal to noise rate;SNR)比较高的条件下,这种方法效果比较好。临床应用中,这种方法受SNR和带宽的影响很大。在SNR较低的情况下,固定百分比设定的噪声阈值会低于实际噪声水平;而在SNR较高的情况下,固定百分比设定的噪声阈值会高于实际噪声水平。因为该方法对噪声的估计是通过平均靠近截止频率附近谱线来实现的,所以在SNR较低时,估计峰值有正向偏移;在SNR较高时,估计的峰值有负向偏移。由于噪声具有一定的随机性,这种方法的平均阈值会随着每一根谱线发生变化,导致估计的阈值水平随之而变,而实际上平均噪声水平相对来说会比较平稳,所以这种会自动包络跟踪的方法可能导致结果不可预料。
在文章“Comparison of four digital maximum frequency estimators for Dopplerultrasound.”,Ultrasound in Med.& Biol.,Vol.14,No.5,pp.355-363,(1988)中,LarryY.L.等人对四种估算多普勒频率最大频偏的方法作出比较,并提出一种百分比法的改进方法(称之为混合方法)。该方法首先计算每一根谱线的频谱积分曲线,再对该积分谱线进行分析。通常情况下,因频谱图能量主要集中在频率较低部分,故频谱积分曲线在频率较低部分变化较快,频率较高部分变化较缓。该混合方法寻找预定直线和所述频谱积分曲线的交叉点,并将该交叉点对应的频率视为最大频偏。其中所述直线的斜率和噪声水平有关,而估计噪声采用了同百分比方法类似的方法,即平均截止频率附近的谱线。
在文章“Comparison of the performance of three maximum Doppler frequencyestimators coupled with different spectral estimation methods.”,Ultrasound in Med.& Biol.,Vol.20,No.7,pp.629-638,(1994)中,K.Marasek等人对混合方法进行改进,提出了几何方法,也是先计算出每一根谱线的频谱积分曲线,再对每条积分谱线进行分析。该方法和所述混合方法不同的是对最大频偏点的确认法:它设计一条直线,计算积分谱线上每点到该直线的距离,认为距离最短的那个点对应的频率即最大频偏。
在文章“The performance of three maximum frequency envelope detection algorithmsfor Doppler signals.”,J.vasc.Invest.1:126-134,(1995)中,R.Moraes等人对几何方法进行改进,提出改进的几何方法:设计直线使经过坐标原点,积分曲线和直线之间的垂直距离最大的那个点对应的频率即最大频偏。为了防止信号较弱时出现不可估计的误差,该方法还增加使用一经验阈值。若信号弱于所述阈值,则不进行最大频偏检测,而直接将最大频偏值设为壁滤波截止频率。
在美国专利US 5,287,753公开的技术方案中,Routh等人介绍了一种自适应阈值包络检测方法。它的基本思想是通过对比多普勒功率谱强度和一设定阈值,来确定出最大频偏。在该方法中,阈值能够根据每个心动周期的信噪比来自适应调整,与前几种方法最大的不同点是,该方法假设在一个心动周期内,噪声水平和平均信噪比比较稳定,因此根据心动周期来确定平均噪声水平和平均信噪比。而前几种方法中阈值会随着每一根谱线剧烈变化,由此估计的阈值水平也会随着每一根谱线剧烈变化。
上述现有技术的主要不足在于:混合方法以频谱图中最大频偏应该小于截止频率为假设前提。在该前提下,对声谱图上边缘部分平均估计噪声才会有意义;这样当截止频率小于最大频偏时,该方法将失效。同时因最大频偏估计结果受实际估计噪声和积分曲线的影响很大,当噪声阈值设定较小,最大频偏估计会偏小,反之会偏大;故该方法受实际SNR和准确估计噪声的影响较大。
几何方法中因最大频偏估计是一种有偏估计,当截止频率小于最大频偏时,该方法有效性会降低。该方法没有直接估计噪声水平,但受积分曲线影响较大:在SNR较小时,积分曲线不能够保证在频率低时变化较快,频率高时变化较慢,此时估计出的最大频率的误差较大。
改进的几何方法中,因最大频率估计是永远小于真实值的有偏估计,和几何方法类似,当最大频偏和截止频率较接近时,估计误差会变大。虽然该方法增加一阈值来降低弱信号时的误判,但该阈值的选择具有经验性,在不同信噪比条件下,选择差异性较大。同几何方法一样,在信号能量较强,但SNR较小时,积分曲线不能够保证在频率低时变化较快,频率高时变化较慢,此时估计出的最大频率的误差较大。
自适应阈值方法中,因噪声阈值由平均信噪比得到,而信噪比受信号强度影响较大,当噪声水平估计不正确时将会影响到最大频偏估计。另外,因阈值和前一心动周期包络检测结果有关,当多普勒频谱图比较稳定时,自适应阈值调整并不能够保证阈值能被收敛到一个稳定值。在确定包络时,采用简单的阈值判断法,由于FFT计算出来的功率谱起伏较大,也影响了包络检测的稳健性。
发明内容 本发明要解决的技术问题是针对上述现有技术的不足,而提出一种频谱多普勒血流速度的自动检测方法,能减小信号信噪比和带宽的影响,准确并且稳健地进行包络检测,从而有利于血流参数被连续自动地准确计算。
为解决上述技术问题,本发明的基本构思为:在自适应阈值方法的基础上,改用基于噪声水平来确定阈值水平,这样当多普勒频谱图稳定时阈值也能够稳定收敛,从而使包络检测更加稳健。此外,还可以对FFT计算出来的功率谱进行平滑处理,以减少功率谱起伏太大引起的包络检测错误。
作为实现本发明构思的技术方案是,提供一种频谱多普勒血流速度的自动检测方法,用于超声系统对血流速度的测量过程,包括步骤:
A.射频超声回波信号经过解调、滤波和模数转换后获得流体的多普勒信号;
B.对所述多普勒信号进行频谱分析,获得该多普勒信号随时间变化的各功率谱线;
C.确定阈值;
D.根据所述阈值和当前功率谱线来确定该功率谱线对应的频偏参数或血流速度,直至测量过程结束或各功率谱线均被处理完毕;
尤其是,所述阈值与前一基本时间段内的平均噪声水平相关,除了在系统启动测试时由系统经验确定外,均自动根据所述平均噪声水平来更新调整。
上述方案中,所述阈值的更新调整是以一预定时间间隔进行的;所述预定时间间隔包括但不限于一个心动周期。
上述方案中,所述阈值更新调整的计算方法是:假设第i-1个时间段的平均噪声水平是NOISEi-1,第i个时间段的阈值是Ti,则当前心动周期的平均噪声水平NOISEi可自动更新为
其中l(t,f)是功率谱密度,t是时间,f是多普勒频率;tis是第i个时间段的起始,
tie是第i个时间段的结束;fo是声谱图上的截止频率,K是遗忘因子;
及第i+1个时间段的阈值Ti+1可自动更新为Ti+1=M*NOISEi 其中M是一个经验常数。
上述方案中,所述步骤D中每确定一功率谱线对应的频偏参数或血流速度前,包括一判断该谱线是否有效的过程;若判断有效,则进行频偏计算,否则直接确定最大频偏。
上述方案中,所述步骤D中每确定一功率谱线对应的频偏参数或血流速度前,还包括一采用滤波方法来平滑所述功率谱的过程。
上述方案中,所述步骤D中的频偏参数包括最大频率,是基于所述阈值,利用谱线上的多个点来共同确定的。
采用上述各技术方案,可以用简单的阈值更新方法使阈值自动收敛,从而准确并稳健地进行包络检测,还能减小信号信噪比和带宽的影响,有利于血流参数的准确计算。
附图说明 图1是传统超声多普勒频谱分析系统处理信号的基本流程图
图2是本发明包络检测软件流程图
图3是本发明最大频率和平均频率检测流程图
图4是平滑前的功率谱示意图
图5是平滑后的功率谱示意图
图6是最大频率判断示意图
图7是平均噪声估计示意图
具体实施方式 下面,结合附图所示之最佳实施例进一步阐述本发明。
本发明频谱多普勒血流速度的自动检测方法,基于图1,用于超声系统(因属于现有技术,不在此赘述)对血流速度的测量过程,包括步骤:
A.射频超声回波信号经过解调、滤波和模数转换后获得流体的多普勒信号;
B.对所述多普勒信号进行频谱分析,获得该多普勒信号随时间变化的各功率谱线;
C.确定阈值;
D.根据所述阈值和当前功率谱线来确定该功率谱线对应的频偏参数或血流速度,直至测量过程结束或各功率谱线均被处理完毕;
其中,所述阈值与前一基本时间段内的平均噪声水平相关,除了在系统启动测试时由系统经验确定外,均自动根据所述平均噪声水平被更新调整。
本发明对所述阈值的更新调整是以一预定时间间隔进行的;所述预定时间间隔包括但不限于一个心动周期。这是根据一个基本时间段内噪声水平和阈值水平可以保持不变的前提而作的,主要因为:一,背景噪声主要来源于超声仪器系统,并且受血流速度的影响很小,不同病人的SNR可能相差非常大,但是当系统参数确定后,噪声水平应该比较一致。二,噪声水平在不同的心动周期有差异,当波形检查比较稳定时,噪声水平才会相对比较稳定。三,噪声分布符合一定的模型,为了简化计算,假设最大噪声能量仅仅和噪声平均水平相关。所以,对于当前心动周期内的当前谱线,所述平均噪声和所述阈值水平Ti均不变;但当谱线进入下一个心动周期,则所述的平均噪声和阈值水平将发生变化。
系统计算其它血流参数的基础是包括在步骤D中确定的最大频率,图2给出了该最大频率偏移的检测过程:系统在开始时,可以根据增益水平和系统参数来初始化平均噪声和阈值;在后续处理中,对FFT估算出来的每一条功率谱谱线逐条处理,分别就各谱线分别计算正向和负向血流引起的最大频偏和平均频偏。
其中,平均噪声估计如图7所示(所述预定时间间隔以一个心动周期为例)。横轴表示时间,纵轴表示频率偏移,图中实线表示计算出来的最大频偏,中线以上的实线表示正向血流引起的最大频偏,以下的实线表示负向血流引起的最大频偏。虚线处的时间表示一个心动周期的结束,此刻需要重新计算平均噪声水平和阈值水平,则计算噪声的噪声区域是最大频偏以外的区域。假设第i-1个心动周期的平均噪声水平是NOISEi-1,假设Ti是第i次心动周期的阈值,则当前心动周期的平均噪声水平是NOISEi可自动更新为
其中l(t,f)是谱密度,t是时间,f是多普勒频率。tis是第i个心动周期的起始,tie是第i个心动周期的结束。fo是声谱图上的截止频率,K是遗忘因子,可根据经验设为0.6~0.9。而第i+1个时间段的阈值Ti+1因假设仅与平均噪声水平NOISEi相关,可自动更新为
Ti+1=M*NOISEi
其中M是一个经验常数。这种阈值更新方法计算简单,而且计算收敛,在信噪比较差的情况下还能稳定检测出包络。
为了排除超声探头从病人身上挪开或者超声探头未对准身体检测区域的情况,所述步骤D中每确定一功率谱线对应的频偏参数或血流速度之前,包括一判断该谱线是否有效的过程;若判断有效,则进行频偏计算,否则直接确定最大频偏。
具体是,假设前一基本时间段的平均噪声是NOISEi,首先寻找当前谱线中最强峰值P(f)max,然后比较P(f)max和平均噪声水平NOISEi,来判断当前谱线是否有效:
若P(f)max小于NOISEi,则该谱线被认为不含任何信号,不必进行最大频率和平均频率检测,可以直接将最大频率和平均频率设定为壁滤波的截止频率;
若有连续的若干条(比如5条)谱线都出现这种情况,则将平均噪声水平NOISEi减半,直到该NOISEi被减少到一个设定的阈值;这种情况被认为当前接收到的只是噪声,而没有信号;
若P(f)max在某个确定的范围内,比如KL*NOISEi和KH*NOISEi之间(其中KL和KH是根据经验选择的系数),则认为当前采集到的多普勒信号信噪比太小,最大频率的检测不是很稳健,当前检测到的最大频率和平均频率与上次检测到的相同;
若P(f)max大于KH*NOISEi,则判断该谱线为有效谱线,可以进行频偏检测。
在这里,所述基本时间段表示一段时间,在本发明实施例中可以被设定但不限为一个心动周期。
最大频偏和平均频偏的检测流程可以如图3所示。一般采用FFT变换估算出来的功率谱会有非常大的起伏,如图4所示,采用谱平均方法可以部分减弱该起伏,但所述起伏还是较大,并且不利于后续包络检测的稳健进行。噪声的起伏会使得最大频偏估计偏高,而信号的起伏会使得最大频偏估计偏低。因此,本发明实施例在检测最大频率之前,还包括一对FFT计算出来的功率谱进行平滑预处理的过程。该平滑处理可以采用简单的滤波。可以在对谱线执行是否有效的所述判断之前进行,也可以在该判断之后进行。不妨设原始计算出来的功率谱是P(fk),则经过平滑处理后的功率谱Paverage(fk)可以表示为
Paverage(fi)=a2P(fi-2)+a1P(fi-1)+a0P(fi)+a1P(fi+1)+a2P(fi+2)
其中a2,a1,a0是平滑滤波系数。图5所示是经过平滑后的功率谱图,与图4平滑前的谱图对比可见,功率谱起伏被较大消除,从而平滑了许多。
如图3所示,检测最大频率时,可以以零频率为基础将谱线分成上下两个部分分别进行处理。其中上半部分对应于正向频率偏移,下半部分对应于负向频率偏移,这样能分别检测出正向频偏和负向频偏。在正向和反向血流同时存在的情况下,这种双边检测的结果更加有效。
本发明实施例中,最大频率是系统基于所述阈值,利用谱线上的多个点来共同稳健地确定的。假设当前心电周期对应的阈值是Ti,图6示意了最大频偏的判断过程。横坐标表示频率变化,纵坐标表示功率谱密度,曲线表示平滑后的功率谱。设需要检测的是正向频率的功率谱。在频率较高部分,功率谱的主要成分是功率谱估计噪声和系统本身噪声;在频率较低部分,功率谱的主要成分是血流运动引起的频率偏移。最大频偏判断需要判断出从噪声到信号的转折点。具体如:从声谱图的边缘开始寻找第一个满足预定条件的点,则该点所对应的频率为检测到的最大频偏;若找不到该点,则认为最大频偏小于壁滤波截止频率。
如图6所示,以所述多个点包括依次连续的四个点为例,第一点对应最大频率应当满足两个条件:1)第一点及第二点的功率谱密度均大于当前心电周期的对应阈值Ti;2)第三点和第四点均大于一预定值,该预定值可以是第一点和第二点功率谱密度的平均值与心电周期对应阈值Ti之和。通过这两条件能够部分消除多普勒频谱图上斑点噪声对包络检测的影响,使包络检测的结果更加稳健。本实施例并不限定本发明方法多点判定的条件,采用其它条件的类似多点判定方法也在本发明的保护范围之内。
检测出最大频偏fmax后,系统可以进一步计算平均频偏fmean,如
由上述计算出的各谱线最大频偏或平均频偏,可以进一步得到最大频偏或平均频偏的包络曲线,本发明图2实施例还包括所述包络曲线的平滑处理过程,目的是消除包络毛刺,使检测出来的包络更加平滑。可以采用中值滤波方法或多点平滑方法。中值滤波是从连续几根谱线计算出来的包络中,取中间值作为当前谱线的包络值。多点平滑即对连续几根谱线计算出来的包络做平滑处理,将这个平滑结果作为最后包络的结果。
本发明方法经过试验验证,在各种信噪比情况下,都能够实时稳健地计算出血流最大速度和平均速度。
Claims (9)
1.一种频谱多普勒血流速度的自动检测方法,用于超声系统对血流速度的测量过程,包括步骤:
A.射频超声回波信号经过解调、滤波和模数转换后获得流体的多普勒信号;
B.对所述多普勒信号进行频谱分析,获得该多普勒信号随时间变化的各功率谱线;
C.确定阈值;
D.根据所述阈值和当前功率谱线来确定该功率谱线对应的频偏参数或血流速度,直至测量过程结束或各功率谱线均被处理完毕;
其特征在于,
所述阈值与前一基本时间段内的平均噪声水平相关,除了在系统启动测试时由系统经验确定外,均自动根据所述平均噪声水平被更新调整。
2.根据权利要求1所述频谱多普勒血流速度的自动检测方法,其特征在于:
所述阈值的更新调整是以一预定时间间隔进行的;所述预定时间间隔包括但不限于一个心动周期。
3.根据权利要求1或2所述频谱多普勒血流速度的自动检测方法,其特征在于,
所述阈值更新调整的计算方法是:假设第i-1个时间段的平均噪声水平是NOISEi-1,第i个时间段的阈值是Ti,则当前心动周期的平均噪声水平NOISEi可自动更新为
其中l(t,f)是功率谱密度,t是时间,f是多普勒频率;tis是第i个时间段的起始,tie是第i个时间段的结束;fo是声谱图上的截止频率,K是遗忘因子;及第i+1个时间段的阈值Ti+1可自动更新为
Ti+1=M*NOISEi
其中M是一个经验常数。
4.根据权利要求3所述频谱多普勒血流速度的自动检测方法,其特征在于:
所述遗忘因子K根据经验设为0.6~0.9。
5.根据权利要求1所述频谱多普勒血流速度的自动检测方法,其特征在于:
所述步骤D中每确定一功率谱线对应的频偏参数或血流速度之前,包括一判断该谱线是否有效的过程;若判断有效,则进行频偏计算,否则直接确定最大频偏。
6.根据权利要求5所述频谱多普勒血流速度的自动检测方法,其特征在于:
所述谱线是否有效的判断方法是,将当前谱线中最强峰值P(f)max与前一基本时间段的平均噪声水平NOISEi进行比较,
若P(f)max小于NOISEi,则认为该谱线不含信号,直接将最大频率和平均频率设定为壁滤波的截止频率;
若连续若干次P(f)max小于NOISEi,则将所述平均噪声水平NOISEi减半,直至一个设定值,并认为只接收到噪声;
若P(f)max在KL*NOISEi和KH*NOISEi之间,其中KL和KH根据经验选择的系数,则认为当前信号信噪比太小,将最大频率和平均频率设定为与上一谱线相同的检测结果;
若P(f)max大于KH*NOISEi,则判断该谱线为有效谱线。
7.根据权利要求1或5所述频谱多普勒血流速度的自动检测方法,其特征在于:
所述步骤D中每确定一功率谱线对应的频偏参数或血流速度之前,还包括一采用滤波方法来平滑所述功率谱的过程。
8.根据权利要求1所述频谱多普勒血流速度的自动检测方法,其特征在于:
所述步骤D中的频偏参数包括最大频率,是基于所述阈值,利用谱线上的多个点来共同确定的。
9.根据权利要求8所述频谱多普勒血流速度的自动检测方法,其特征在于:
所述多个点包括依次连续的四个点,第一点对应最大频率当且仅当满足条件:
1)第一点及第二点的功率谱密度均大于当前所述阈值;
2)第三点和第四点均大于一预定值,该预定值是第一点和第二点功率谱密度的平均值与当前所述阈值之和。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100363720A CN100496409C (zh) | 2005-08-02 | 2005-08-02 | 频谱多普勒血流速度的自动检测方法 |
US11/317,848 US7798968B2 (en) | 2005-08-02 | 2005-12-23 | Automatic detection system and method of spectral Doppler blood flow velocity |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2005100363720A CN100496409C (zh) | 2005-08-02 | 2005-08-02 | 频谱多普勒血流速度的自动检测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1907230A true CN1907230A (zh) | 2007-02-07 |
CN100496409C CN100496409C (zh) | 2009-06-10 |
Family
ID=37698578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2005100363720A Active CN100496409C (zh) | 2005-08-02 | 2005-08-02 | 频谱多普勒血流速度的自动检测方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US7798968B2 (zh) |
CN (1) | CN100496409C (zh) |
Cited By (24)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101301212B (zh) * | 2007-05-11 | 2011-08-17 | 深圳迈瑞生物医疗电子股份有限公司 | 实时估计多普勒参数的方法及装置 |
CN102238915A (zh) * | 2008-12-02 | 2011-11-09 | 株式会社东芝 | 超声波诊断装置、多普勒测量装置和多普勒测量方法 |
CN101461720B (zh) * | 2007-12-18 | 2012-01-25 | 深圳迈瑞生物医疗电子股份有限公司 | 基于频谱多普勒的运动速度测量范围的调整方法及装置 |
CN101461719B (zh) * | 2007-12-18 | 2012-02-01 | 深圳迈瑞生物医疗电子股份有限公司 | 一种多普勒扫描变换实现方法和装置 |
CN102370499A (zh) * | 2010-08-26 | 2012-03-14 | 深圳迈瑞生物医疗电子股份有限公司 | 多普勒图像、b型图像和彩色血流图像同时显示的方法和系统 |
CN101766496B (zh) * | 2008-12-31 | 2012-08-29 | 深圳迈瑞生物医疗电子股份有限公司 | 一种噪声估计方法和一键优化方法及其系统 |
CN103505246A (zh) * | 2012-06-18 | 2014-01-15 | 深圳市蓝韵实业有限公司 | 多普勒参数实时自动标记的方法 |
CN103841898A (zh) * | 2011-09-30 | 2014-06-04 | 皇家飞利浦有限公司 | 具有随样本体积移动的自动动态多普勒流设置的超声系统 |
CN105662472A (zh) * | 2016-01-13 | 2016-06-15 | 北京悦琦创通科技有限公司 | 生成多普勒频谱图的方法和设备 |
CN105997147A (zh) * | 2016-07-28 | 2016-10-12 | 飞依诺科技(苏州)有限公司 | 一种超声波脉冲多普勒成像方法及装置 |
CN107389787A (zh) * | 2017-07-28 | 2017-11-24 | 电子科技大学 | 一种基于频域的气侵检测信号处理方法 |
CN108186050A (zh) * | 2018-01-03 | 2018-06-22 | 声泰特(成都)科技有限公司 | 一种基于超声通道数据的多普勒血流速度成像方法和系统 |
CN108680786A (zh) * | 2018-04-11 | 2018-10-19 | 东南大学 | 一种脉冲信号频域自适应滤波包络提取方法 |
CN109171810A (zh) * | 2018-09-21 | 2019-01-11 | 沈阳华龄智能科技产业有限公司 | 一种护理监护系统及方法 |
CN109414245A (zh) * | 2016-09-30 | 2019-03-01 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流运动谱的显示方法及其超声成像系统 |
CN109470890A (zh) * | 2018-09-26 | 2019-03-15 | 沈畅 | 一种超声仪表的动态监测修正算法 |
CN110337273A (zh) * | 2017-11-14 | 2019-10-15 | 深圳市理邦精密仪器股份有限公司 | 用于超声成像的多分辨率判别分析的系统与方法 |
CN110786883A (zh) * | 2019-11-06 | 2020-02-14 | 深圳开立生物医疗科技股份有限公司 | 频谱包络线的生成方法、装置、超声设备及存储介质 |
CN113066083A (zh) * | 2021-04-25 | 2021-07-02 | 青岛海信医疗设备股份有限公司 | 流体的多普勒参数确定方法和电子设备 |
CN113221828A (zh) * | 2021-05-31 | 2021-08-06 | 煤炭科学研究总院 | 一种爆破震动响应信号的去噪方法、装置及电子设备 |
WO2021189671A1 (zh) * | 2020-03-26 | 2021-09-30 | 深圳开立生物医疗科技股份有限公司 | 超声多普勒血流成像方法、装置、设备及可读存储介质 |
CN114041826A (zh) * | 2021-11-11 | 2022-02-15 | 苏州国科昂卓医疗科技有限公司 | 椎骨钉道血管成像方法、系统及椎骨钉道创建引导装置 |
CN114414836A (zh) * | 2021-11-29 | 2022-04-29 | 北京信息科技大学 | 一种地下管道流速测量方法 |
CN116299318A (zh) * | 2023-05-18 | 2023-06-23 | 成都凯天电子股份有限公司 | 基于点云密度的直升机大气数据的去噪方法 |
Families Citing this family (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7695439B2 (en) * | 2005-08-22 | 2010-04-13 | Siemens Medical Solutions Usa, Inc. | Automated identification of cardiac events with medical ultrasound |
EP1982211A1 (en) * | 2006-01-27 | 2008-10-22 | Koninklijke Philips Electronics N.V. | Automatic ultrasonic doppler measurements |
US20140336510A1 (en) * | 2013-05-08 | 2014-11-13 | Siemens Medical Solutions Usa, Inc. | Enhancement in Diagnostic Ultrasound Spectral Doppler Imaging |
WO2015033469A1 (ja) * | 2013-09-09 | 2015-03-12 | パイオニア株式会社 | 流速検出装置及び流速検出方法 |
EP2989986B1 (en) * | 2014-09-01 | 2019-12-18 | Samsung Medison Co., Ltd. | Ultrasound diagnosis apparatus and method of operating the same |
US20160377717A1 (en) * | 2015-06-29 | 2016-12-29 | Edan Instruments, Inc. | Systems and methods for adaptive sampling of doppler spectrum |
US20170119356A1 (en) * | 2015-10-30 | 2017-05-04 | General Electric Company | Methods and systems for a velocity threshold ultrasound image |
JP6291557B2 (ja) * | 2016-10-26 | 2018-03-14 | パイオニア株式会社 | 流体評価装置及び方法 |
EP3323350A3 (en) * | 2016-11-16 | 2018-06-20 | Huntleigh Technology Limited | Ultrasound doppler monitoring system |
JP2018105879A (ja) * | 2018-02-09 | 2018-07-05 | パイオニア株式会社 | 流体評価装置及び方法 |
US11779311B2 (en) * | 2018-09-14 | 2023-10-10 | Fujifilm Sonosite, Inc. | Method and apparatus for performing spectral doppler imaging |
JP2019168465A (ja) * | 2019-05-28 | 2019-10-03 | パイオニア株式会社 | 流体評価装置及び方法 |
JP7343342B2 (ja) * | 2019-09-25 | 2023-09-12 | キヤノンメディカルシステムズ株式会社 | 超音波診断装置、及び画像処理装置 |
JP2021067696A (ja) * | 2020-12-28 | 2021-04-30 | パイオニア株式会社 | 流体評価装置及び方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS61170442A (ja) * | 1985-01-23 | 1986-08-01 | 松下電器産業株式会社 | 超音波ドツプラ血流装置 |
US5287753A (en) * | 1992-05-02 | 1994-02-22 | Advanced Technology Laboratories, Inc. | Continuous display of peak and mean blood flow velocities |
NO942222D0 (no) * | 1994-06-14 | 1994-06-14 | Vingmed Sound As | Fremgangsmåte ved bestemmelse av hastighet/tid-spektrum ved blodströmning |
WO1996021196A1 (en) * | 1994-12-30 | 1996-07-11 | Acuson Corporation | Adaptive temporal filtering to enhance fluid flow or tissue motion imaging |
US5634465A (en) * | 1995-06-09 | 1997-06-03 | Advanced Technology Laboratories, Inc. | Continuous display of cardiac blood flow information |
US5935074A (en) * | 1997-10-06 | 1999-08-10 | General Electric Company | Method and apparatus for automatic tracing of Doppler time-velocity waveform envelope |
JP2000005177A (ja) | 1998-06-17 | 2000-01-11 | Matsushita Electric Ind Co Ltd | 超音波ドプラ血流計 |
US6464641B1 (en) * | 1998-12-01 | 2002-10-15 | Ge Medical Systems Global Technology Company Llc | Method and apparatus for automatic vessel tracking in ultrasound imaging |
US6418343B1 (en) * | 1999-10-01 | 2002-07-09 | Cardiac Pacemakers, Inc. | Method and apparatus for adjusting the sensing threshold of a cardiac rhythm management device |
KR100381875B1 (ko) * | 2000-07-08 | 2003-04-26 | 주식회사 메디슨 | 도플러 효과를 이용하여 혈류 속도를 측정하기 위한초음파 진단 장치 및 방법 |
JP4245428B2 (ja) * | 2003-07-29 | 2009-03-25 | 東芝医用システムエンジニアリング株式会社 | 超音波ドプラ診断装置及び画像データ生成方法 |
-
2005
- 2005-08-02 CN CNB2005100363720A patent/CN100496409C/zh active Active
- 2005-12-23 US US11/317,848 patent/US7798968B2/en active Active
Cited By (36)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101301212B (zh) * | 2007-05-11 | 2011-08-17 | 深圳迈瑞生物医疗电子股份有限公司 | 实时估计多普勒参数的方法及装置 |
CN101461720B (zh) * | 2007-12-18 | 2012-01-25 | 深圳迈瑞生物医疗电子股份有限公司 | 基于频谱多普勒的运动速度测量范围的调整方法及装置 |
CN101461719B (zh) * | 2007-12-18 | 2012-02-01 | 深圳迈瑞生物医疗电子股份有限公司 | 一种多普勒扫描变换实现方法和装置 |
CN102238915B (zh) * | 2008-12-02 | 2014-12-03 | 株式会社东芝 | 超声波诊断装置、多普勒测量装置和多普勒测量方法 |
CN102238915A (zh) * | 2008-12-02 | 2011-11-09 | 株式会社东芝 | 超声波诊断装置、多普勒测量装置和多普勒测量方法 |
CN101766496B (zh) * | 2008-12-31 | 2012-08-29 | 深圳迈瑞生物医疗电子股份有限公司 | 一种噪声估计方法和一键优化方法及其系统 |
CN102370499A (zh) * | 2010-08-26 | 2012-03-14 | 深圳迈瑞生物医疗电子股份有限公司 | 多普勒图像、b型图像和彩色血流图像同时显示的方法和系统 |
US10166006B2 (en) | 2011-09-30 | 2019-01-01 | Koninklijke Philips N.V. | Ultrasound system with dynamically automated doppler flow settings as a sample volume is moved |
CN103841898A (zh) * | 2011-09-30 | 2014-06-04 | 皇家飞利浦有限公司 | 具有随样本体积移动的自动动态多普勒流设置的超声系统 |
CN103841898B (zh) * | 2011-09-30 | 2016-12-21 | 皇家飞利浦有限公司 | 具有随样本体积移动的自动动态多普勒流设置的超声系统 |
US11344282B2 (en) | 2011-09-30 | 2022-05-31 | Koninklijke Philips N.V. | Ultrasound system with dynamically automated doppler flow settings as a sample volume is moved |
CN103505246A (zh) * | 2012-06-18 | 2014-01-15 | 深圳市蓝韵实业有限公司 | 多普勒参数实时自动标记的方法 |
CN105662472A (zh) * | 2016-01-13 | 2016-06-15 | 北京悦琦创通科技有限公司 | 生成多普勒频谱图的方法和设备 |
CN105997147A (zh) * | 2016-07-28 | 2016-10-12 | 飞依诺科技(苏州)有限公司 | 一种超声波脉冲多普勒成像方法及装置 |
CN109414245A (zh) * | 2016-09-30 | 2019-03-01 | 深圳迈瑞生物医疗电子股份有限公司 | 超声血流运动谱的显示方法及其超声成像系统 |
CN107389787B (zh) * | 2017-07-28 | 2019-09-24 | 电子科技大学 | 一种基于频域的气侵检测信号处理方法 |
CN107389787A (zh) * | 2017-07-28 | 2017-11-24 | 电子科技大学 | 一种基于频域的气侵检测信号处理方法 |
CN110337273A (zh) * | 2017-11-14 | 2019-10-15 | 深圳市理邦精密仪器股份有限公司 | 用于超声成像的多分辨率判别分析的系统与方法 |
CN108186050B (zh) * | 2018-01-03 | 2021-02-23 | 声泰特(成都)科技有限公司 | 一种基于超声通道数据的多普勒血流速度成像方法和系统 |
CN108186050A (zh) * | 2018-01-03 | 2018-06-22 | 声泰特(成都)科技有限公司 | 一种基于超声通道数据的多普勒血流速度成像方法和系统 |
CN108680786A (zh) * | 2018-04-11 | 2018-10-19 | 东南大学 | 一种脉冲信号频域自适应滤波包络提取方法 |
CN109171810A (zh) * | 2018-09-21 | 2019-01-11 | 沈阳华龄智能科技产业有限公司 | 一种护理监护系统及方法 |
CN109470890A (zh) * | 2018-09-26 | 2019-03-15 | 沈畅 | 一种超声仪表的动态监测修正算法 |
CN109470890B (zh) * | 2018-09-26 | 2020-07-31 | 沈畅 | 一种超声仪表的动态监测修正算法 |
CN110786883B (zh) * | 2019-11-06 | 2022-05-27 | 深圳开立生物医疗科技股份有限公司 | 频谱包络线的生成方法、装置、超声设备及存储介质 |
CN110786883A (zh) * | 2019-11-06 | 2020-02-14 | 深圳开立生物医疗科技股份有限公司 | 频谱包络线的生成方法、装置、超声设备及存储介质 |
WO2021189671A1 (zh) * | 2020-03-26 | 2021-09-30 | 深圳开立生物医疗科技股份有限公司 | 超声多普勒血流成像方法、装置、设备及可读存储介质 |
CN113066083A (zh) * | 2021-04-25 | 2021-07-02 | 青岛海信医疗设备股份有限公司 | 流体的多普勒参数确定方法和电子设备 |
CN113221828A (zh) * | 2021-05-31 | 2021-08-06 | 煤炭科学研究总院 | 一种爆破震动响应信号的去噪方法、装置及电子设备 |
CN113221828B (zh) * | 2021-05-31 | 2022-03-08 | 煤炭科学研究总院 | 一种爆破震动响应信号的去噪方法、装置及电子设备 |
CN114041826B (zh) * | 2021-11-11 | 2024-01-30 | 苏州国科昂卓医疗科技有限公司 | 椎骨钉道血管成像方法、系统及椎骨钉道创建引导装置 |
CN114041826A (zh) * | 2021-11-11 | 2022-02-15 | 苏州国科昂卓医疗科技有限公司 | 椎骨钉道血管成像方法、系统及椎骨钉道创建引导装置 |
CN114414836A (zh) * | 2021-11-29 | 2022-04-29 | 北京信息科技大学 | 一种地下管道流速测量方法 |
CN114414836B (zh) * | 2021-11-29 | 2023-09-19 | 北京信息科技大学 | 一种地下管道流速测量方法 |
CN116299318B (zh) * | 2023-05-18 | 2023-08-11 | 成都凯天电子股份有限公司 | 基于点云密度的直升机大气数据的去噪方法 |
CN116299318A (zh) * | 2023-05-18 | 2023-06-23 | 成都凯天电子股份有限公司 | 基于点云密度的直升机大气数据的去噪方法 |
Also Published As
Publication number | Publication date |
---|---|
CN100496409C (zh) | 2009-06-10 |
US7798968B2 (en) | 2010-09-21 |
US20070043294A1 (en) | 2007-02-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1907230A (zh) | 频谱多普勒血流速度的自动检测方法 | |
CN1875887A (zh) | 声谱图包络曲线提取方法 | |
US7535797B2 (en) | High-resolution ultrasound displacement measurement apparatus and method | |
US5935074A (en) | Method and apparatus for automatic tracing of Doppler time-velocity waveform envelope | |
US8050476B2 (en) | Heart rate demodulation of periodic movement in ultrasound data for object identification | |
CN101138507B (zh) | 多普勒血流声音信号的处理方法及装置 | |
US20100137717A1 (en) | Automatic Flow Tracking System and Method | |
KR102025328B1 (ko) | 평면파 합성을 이용한 초음파 벡터 도플러 영상의 생성 장치 및 방법 | |
CN101766496B (zh) | 一种噪声估计方法和一键优化方法及其系统 | |
JP3581676B2 (ja) | 血流速度測定向け超音波診断装置 | |
US20210290204A1 (en) | Medical imaging system and method | |
JP2013532570A (ja) | 胎児モニタリング中に超音波信号の不明確性を削減する方法 | |
CN112515704B (zh) | 一种基于超声的血管硬度测量方法 | |
CN110879389A (zh) | 基于多基地ir-uwb生物雷达信号的多人体目标识别定位方法 | |
CN105708494B (zh) | 一种基于超声的血压测量方法 | |
CN110879388A (zh) | 基于ir-uwb生物雷达信号的人与动物非接触探测区分方法 | |
CN104665877A (zh) | 颈动脉血管局部脉搏波传播速度测量方法 | |
CN103505246A (zh) | 多普勒参数实时自动标记的方法 | |
CN109471113B (zh) | 基于相位法的多波束声呐海底地形测量质量实时评估方法 | |
EP2034901A1 (en) | Global and local detection of blood vessel elasticity | |
US20230061869A1 (en) | System and methods for beamforming sound speed selection | |
CN101889877A (zh) | 基于希尔伯特变换-基频估计的平均骨小梁间距估计方法 | |
US9737283B2 (en) | Ultrasound measurement apparatus and ultrasound measurement method | |
CN103494615A (zh) | 最大频率包络曲线的计算方法和装置 | |
Hasegawa et al. | Reduction of influence of decrease in signal-to-noise ratio in measurement of change in thickness of arterial wall due to heartbeat |
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 | ||
EE01 | Entry into force of recordation of patent licensing contract | ||
EE01 | Entry into force of recordation of patent licensing contract |
Application publication date: 20070207 Assignee: Shenzhen Mindray Animal Medical Technology Co.,Ltd. Assignor: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd. Contract record no.: X2022440020009 Denomination of invention: An Automatic Detection Method of Spectral Doppler Blood Velocity Granted publication date: 20090610 License type: Common License Record date: 20220804 |