CN101646391B - 用于超声波成像的方法和装置 - Google Patents

用于超声波成像的方法和装置 Download PDF

Info

Publication number
CN101646391B
CN101646391B CN2008800102351A CN200880010235A CN101646391B CN 101646391 B CN101646391 B CN 101646391B CN 2008800102351 A CN2008800102351 A CN 2008800102351A CN 200880010235 A CN200880010235 A CN 200880010235A CN 101646391 B CN101646391 B CN 101646391B
Authority
CN
China
Prior art keywords
signal
power
image
gain control
type
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.)
Expired - Fee Related
Application number
CN2008800102351A
Other languages
English (en)
Other versions
CN101646391A (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.)
Hitachi Ltd
Original Assignee
Arockar K K
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 Arockar K K filed Critical Arockar K K
Publication of CN101646391A publication Critical patent/CN101646391A/zh
Application granted granted Critical
Publication of CN101646391B publication Critical patent/CN101646391B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52023Details of receivers
    • G01S7/52033Gain control of receivers
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/06Measuring blood flow
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
    • G01S15/88Sonar systems specially adapted for specific applications
    • G01S15/89Sonar systems specially adapted for specific applications for mapping or imaging
    • G01S15/8906Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques
    • G01S15/8979Combined Doppler and pulse-echo imaging systems
    • G01S15/8981Discriminating between fixed and moving objects or between objects moving at different speeds, e.g. wall clutter filter
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO 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/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/52Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
    • G01S7/52017Details 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/52077Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 particularly adapted to short-range imaging with means for elimination of unwanted signals, e.g. noise or interference

Landscapes

  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Biophysics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Acoustics & Sound (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • General Physics & Mathematics (AREA)
  • Medical Informatics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Hematology (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

血流信息用于抑制在血管超声波B型图像中的噪声表现。血流信号通过流体探测器获得。经过高通滤波后,只有血管内腔中的流体信号功率得到保留,而来自固定组织区域的信号功率受到抑制。流体信号成分用于计算流体成分参数,该流体成分参数用于生成减小在B型图像中的噪声的增益控制信号α。

Description

用于超声波成像的方法和装置
技术领域
本发明主要涉及超声波成像领域。本发明的实施例尤其涉及在对血管进行成像时用于减少噪声的方法和系统。
背景技术
超声波被用于对各种内部身体结构例如像心脏和肝脏、血管等器官以及孕妇体内的胎儿进行成像。B型成像是用于对血管进行成像的技术。在血管中,血液通常展示出比血管壁或者周边组织低的回波功率,从而产生高的反差。血管壁的表层或者内膜在诊断诸如动脉硬化、狭窄化或者冠状堵塞等心血管疾病方面特别有利。在诊断中通常需要测量和利用内膜厚度。然而由于不同的原因,通常在血管内腔中会出现杂波噪声,这使得内膜厚度的测量变得困难。
因而有必要减少在血管图像中存在的噪声。
发明内容
本发明的发明人发现可能需要一种利用血流信息减小血管B型图像中的杂波噪声的系统和方法。
本发明的一个方案提供了一种当通过在每一位置对血管传送数次超声波信号并接收数次自血管返回的超声波信号而对血管成像时通过返回的超声波信号抑制杂波噪声表现的方法。所述方法包括:对返回的超声波图像信号进行处理以产生B型图像输出;对返回的超声波信号进行解调以产生多普勒信号;对多普勒信号进行壁滤波以仅让信号的流体成分通过;计算流体成分参数;基于流体成分参数生成增益控制信号α;以及以此增益控制信号α控制B型图像处理器的输出,其中,在B型图像处理器的输出中的噪声表现得以抑制。
本发明的另一方案提供一种当通过在每一位置对血管传送数次超声波信号并接收数次自血管返回的超声波信号而对血管成像时通过返回的超声波信号抑制杂波噪声表现的方法。根据本发明这一方案的方法包括对返回的超声波图像信号进行处理以产生B型图像输出;对返回的超声波图像信号进行壁滤波以仅使信号中的流体成分通过;计算流体成分参数;基于流体成分参数生成增益控制信号α;以及以此增益控制信号α控制B型图像处理器的输出,其中,在B型图像处理器的输出中的噪声表现得以抑制。
本发明的另一方案提供一种当通过在每一位置对血管传送数次超声波信号并接收数次自血管返回的超声波信号而血管成像时通过返回的超声波信号抑制杂波噪声表现的系统。根据本发明这一方案的系统包括接收器,其配置成接收返回的超声波图像信号并输出接收到的信号;连结到所述接收器的B型图像处理器,其配置成从接收到的信号输出经B型处理的图像;连结到接收器的多普勒流体探测器,其配置成对返回的超声波图像信号进行解调以产生多普勒信号,该多普勒探测器具有壁滤波器(wall filter),该壁滤波器配置成对多普勒信号进行滤波并只输出其中的流体成分;连结到多普勒流体探测器的增益控制发生器,其配置成计算流体成分参数并基于流体成分生成增益控制信号α;以及连结到所述增益控制发生器和所述B型图像处理器的信号组合器,该信号组合器配置成以所述增益控制信号α对所述经B型处理的图像输出进行修改,其中,在经B型处理的图像中的噪声表现得以抑制。
本发明的另一方案提供一种当通过在每一位置对血管传送数次超声波信号并接收数次自血管返回的超声波信号而对血管成像时通过返回的超声波信号抑制杂波噪声表现的系统。根据本发明这一方案的系统包括接收器,其配置成接收返回的超声波图像信号并输出接收到的信号;连结到所述接收器的B型图像处理器,其配置成从接收到的信号输出经B型处理的图像;连结到所述接收器的射频流体探测器,具有壁滤波器的该流体探测器配置成对返回的超声波射频信号进行滤波,并只输出接收到的信号的流体成分;连结到射频流体探测器的增益控制发生器,其配置成计算流体成分参数并基于流体参数生成增益控制信号α;以及连结到所述增益控制发生器和所述B型图像处理器的信号组合器,该信号组合器被设置成以所述增益控制信号α对所述经B型处理的图像输出进行修改,其中,在经B型处理的图像中的噪声表现得以抑制。
在附图和以下说明中将对本发明的一个或多个实施例的细节加以阐明。通过说明和图解并通过权利要求,本发明的其他特征、目的及优点将一目了然。
附图说明
图1是示例性的呈现有杂波噪声的血管的B型图像。
图2是示例性的血流图像。
图3是示例性的在抑制了杂波噪声之后的B型图像。
图4是用于利用多普勒流体处理的B型图像的杂波噪声抑制系统的示例性系统图。
图5是用于利用射频流体处理的B型图像的杂波噪声抑制系统的示例性系统图。
图6是示例性的多普勒流体探测器。
图7是示例性的射频流体探测器。
图8是配置成查找表(look-up table)的示例性增益控制信号发生器。
图9是示例性的增益控制信号曲线。
图10是配置成查找表的示例性增益控制信号发生器。
图11是示例性的增益控制信号曲线。
图12是配置成数字信号处理器的示例性增益控制信号发生器。
图13是配置成数字信号处理器的示例性增益控制信号发生器。
图14是利用多普勒流体处理对B型图像杂波进行抑制的示例性方法。
图15是利用射频流体处理对B型图像杂波进行抑制的示例性方法。
具体实施方式
下面参照附图对本发明的实施例进行说明,其中,相同元件全部以相同数字表示。在对本发明的实施例进行详细解释之前,应理解,本发明并不限于在以下说明中提到的或在附图中描绘的具体例子中的应用。本发明能够以各种应用以及不同的方式运用或实施于其它的实施例。同时应理解,此间所用的措辞及术语只是为说明的目的而不应被认为是限制性的。此间所用的“包括(including)”,“由…组成(comprising)1”,或者“具有(having)”以及这些词的变形用于表示包括其后所列的细项及其等价物以及附加的细项。术语“装配的(mounted)”,“连接的(connected)”,和“连结的(coupled)”被广泛使用,其包括直接和间接的装配、连接和连结。另外,“连接的”和“连结的”并不限于物理或机械的连接或连结。
应注意,本发明不限于任何此间提到的或者在附图中间接表示的特别的软件语言。本领域的普通技术人员将会理解有多种可选软件语言可以用于实施本发明。也应当理解,其中有些部件和细项是当作硬件元件描绘和说明的,这是本领域内的普遍做法。然而,本领域的普通技术人员在阅读本详细说明的基础上将会理解,在至少一个实施例中,方法和系统中的部件可以用软件或硬件实现。
超声波通过超声波换能器被传送进入人体内,以对不同的器官、血管或孕妇体内的胎儿进行成像。超声波被组织内的散射体所散射,而被散射的超声波返回到换能器。接收波束形成器产生超声波束,后信息处理器从返回的超声波信号的振幅生成组织的图像作为B型图像。
由于血管显示了病人心血管的状况,因而常被成像。为了诊断经常需要测量并利用内膜厚度。然而,内膜的图像常常由于不同原因导致的噪声而模糊不清。血流信息通常利用彩色多普勒和频谱多普勒技术获得。
彩色多普勒是二维成像技术,通常用于通过把超声波发送到血管中并对被移动的血红细胞所散射的超声波进行探测来对血液成像。其由与B型图像相似的多束波组成。为探测流体速度,彩色多普勒在每一位置传送数次超声波信号以探测运动。为生成二维流体图像,传送位置以亚毫米或者约一个超声波波长的数量级变化。传送位置的变化需重复约100次才能覆盖几厘米以生成二维流体图像。对于相位阵列换能器或扇形图像格式,传送方向改变一个小角度,例如大约0.5至1.0度。这需要重复约100次才能覆盖扇形图像的大约90度范围。对于每一传送位置或方向,超声波都需传送数次。接收到的成束的射频超声波信号要经过正交解调,产生复合多普勒I-Q信号。
多普勒I-Q信号可能包括血流信号成分以及固定组织信号成分。典型地,固定组织信号成分要比血流信号成分大30至40分贝。因此,为了精确探测血流,必须将固定组织信号成分去除掉。
高通滤波器被应用于接收到的来自数次传送的多普勒信号,以仅获得流体信号成分,由于该种滤波器去除了血管壁噪声,因而通常也被称作壁滤波器。高通滤波器去除固定信号成分,且其只允许流体信号成分通过。高通滤波器的其中一种形式可以是信号减法或者2-tap FIR,其中,随后传送的多普勒信号被从在先传送的多普勒信号中减去。在接收到的连续传送的多普勒I-Q信号之间相位差显示了血流。另外,被高通滤波器滤去的多普勒I-Q信号的功率显示了血流的存在。
本发明的系统和方法通过利用上述血流信号降低了B型信号振幅而减少了杂波噪声。即使在血管内腔中存在杂波噪声,在内腔的相同位置上通常也会存在血流。血管中的杂波噪声使内膜模糊不清,并使得内膜厚度的测量变得困难。血流成分参数诸如振幅a,功率a2或者提高到ab(其中b为一实数)的功率,可通过计算得到,并用于生成增益控制信号α,以在同一图像位置处抑制B型图像处理器的输出信号。流体成分参数可以是总量或者平均量。抑制的杂波噪声量与所用的流体成分参数之间可以是成比例的或者预定的关系。
图1显示了从典型的B型图像处理器输出的血管101图像,其中,血管101具有近侧壁103和远侧壁105。示出的杂波噪声107靠近近侧壁103。图2显示了同一血管的血流201图像,其被用于减小图1所示的B型图像的增益。图3显示的是具有清晰血管内腔301的合成图像。
图4显示了根据本发明的超声波系统。图14显示了阐明本方法的流程图。超声波信号从由发送器406驱动的超声波探头402通过发送/接收开关404传送。接收器408通过开关404接收从探头402接收到的超声波信号,并对所述信号进行处理(步骤1405)。经处理的信号被连结至流体探测器410、B型图像处理器412(步骤1410)以及多普勒频谱处理器424。多普勒频谱处理器对从接收器408输出的信号的多普勒频谱进行计算,并以多普勒频谱模式或者以与B型图像和/或彩色流体图像相结合的模式将此多普勒频谱输出到扫描转换器420。
流体探测器410探测血流,计算并输出流体成分参数至增益控制信号发生器418,所述流体成分参数可以是振幅a,功率a2,提高到ab的功率或者这些值的组合。增益控制信号发生器418根据流体成分参数生成增益控制信号α,并将此增益控制信号α输出至信号组合器414,所述信号组合器414可以是例如乘法器或者可变增益放大器。
图6显示了多普勒流体探测器410的示意图。该多普勒流体探测器包括壁滤波器(即高通滤波器)602、速度计算器604、功率计算器606和方差计算器606。壁滤波器接收来自接收器408的经解调的I-Q信号(步骤1415)。高通滤波器的截止频率或者拐点频率可由操作者调节,其阻止低频率的固定组织信号成分通过,而只允许较高频率的流体信号成分通过。高通滤波器的结构可以是有限脉冲响应(FIR,finite impulse response)滤波器、无限脉冲响应(IIR,infinite impulse response)滤波器、多项式滤波器、回归线滤波器或者其他结构类型(步骤1420)。流体信号成分被连结至速度计算器604、功率计算器606以及方差计算器608。
速度计算器604计算血流速度,血流速度被输出至扫描转换器420,该扫描转换器420将速度信号转换成电视光栅扫描速度图像。然后,此速度图像以常规的彩色流体模式显示在显示监视器422上。方差计算器608计算方差或者湍流指标,该方差或者湍流指标被输出至扫描转换器420,所述扫描转换器420将方差信号转换成电视光栅扫描方差图像。然后,此方差图像以常规的彩色流体模式显示在显示监视器422上。功率计算器606计算流体成分参数:信号振幅a,功率a2,提高到ab的功率以及这些值的组合。流体成分参数通常利用下述公式由取样的经高通滤波所得的多普勒I-Q复合信号算得:
a 2 = Σ i = 1 N z i z i * - - - ( 1 )
其中zi是经高通滤波的复合多普勒信号,i表示一个超声波传送序列中的第i个成分,N是离散时间内经过高通滤波的信号样本的个数。*号代表复共轭性。经过滤输出的数目通常小于发送/接收信号的数目。
功率计算器606可以是DSP、FPGA、ASIC或者离散部件,例如乘法器、加法器、除法器以及绝对值计算器。
其他的流体成分参数可以通过下列公式由功率a2获得:
a = a 2 - - - ( 2 ) 或者
a b = ( a 2 ) b - - - ( 3 )
信号振幅a以及提高到ab的功率也可以通过下面的公式得到:
a = Σ i = 1 N | z i | - - - ( 4 ) 或者
a b = Σ i = 1 N | z i | b - - - ( 5 )
通过式(4)和式(5)得到的流体成分参数不同于通过式(2)和式(3)得到的值,因为两者处理步骤的次序不相同(步骤1425)。上述对信号振幅a,功率a2,提高到ab的功率这样的流体成分参数的计算反映的是总值。全部流体成分参数可以通过额外的计算,即通过将每个参数值除以样本数N而得到平均值的方法标准化。
振幅a、功率a2、提高到ab的功率或者这些值的组合这样的流体成分参数被连结至增益控制发生器418,并被用于生成增益控制信号α。为了增加信噪比(SNR),可以将流体成分参数组合成例如 d 1 a + d 2 a 2 + Σ i d i a b ( i ) 这样的组合流体成分参数,其中d1、d2和di是代表加权系数的实数,i第i个成分b(i)同样也是实数。组合多于一个的流体成分参数通常会增加信噪比,并减少无关联的噪声。 d 1 a + d 2 a 2 + Σ i d i a b ( i ) 可以在功率计算器606中计算得到,然后也输出给增益控制发生器418中。功率计算器606也可以以彩色流体模式输出给扫描转换器420。
增益控制发生器418可以是例如如图12中所示的数字信号处理器(DSP)1210,或者是如图8中所示的查找表(LUT)810,或者FPGA、ASIC或者诸如乘法器和加法器等离散部件。增益控制发生器418使用流体成分参数作为输入。图9显示了LUT 810或者DSP的信号处理能力的响应曲线。如果输入小,则增益高。如果输入大,则增益将归零或者接近于零。图9所示的曲线是示例性的。可以使用其他预定曲线来抑制B型图像处理器412的信号输出(步骤1430)。增益控制发生器418的输出是增益控制信号α,并且该增益控制信号α被连结到信号组合器414。
B型图像处理器412生成B型图像,并将图像输出给信号组合器414。图像被输出至扫描转换器420,该扫描转换器将图像信号转换成电视光栅扫描图像。然后,图像被显示在显示监视器422上。B型信号增益由增益控制信号α所控制。例如,如果流体成分参数大,信号组合器414的增益就小,那么如果血管内腔中存在杂波噪声,则此杂波噪声将被抑制。由于血流只存在于发生杂波噪声的血管内腔中,所以整体的增益控制就去除了杂波噪声。如果流体成分参数低(在组织区可能是这种情况),则信号组合器414的增益大,产生普通亮度的B型组织图像。这样,在血流区域中的杂波噪声就被流体信号减小了(步骤1435)。
在一个供替换的实施例中,不用彩色多普勒过程,而是利用互相关或者时移技术来在射频信号水平上探测血流。为探测流速,在每个位置传送数次超声波以探测运动。为生成二维流体图像,传送位置以亚毫米或者约一个超声波波长的数量级变化。传送位置的变化需重复约100次才能覆盖几厘米以生成二维流体图像。对于相位阵列换能器或扇形图像格式,传送方向改变一个小角度,例如大约0.5至1.0度。这需要重复约100次才能覆盖扇形图像的大约90度范围。对于每一传送位置或方向,超声波都需传送数次。高通滤波被应用于成束形成的射频信号以去除组织信号,并使血流信号通过。血流成分参数诸如振幅c,功率c2或者提高到cb的功率(其中b为一实数),或者这些值的组合,可通过计算得到,并用于产生增益控制信号α,以在同一图像位置处抑制B型图像处理器的输出信号。流体成分参数可以是总量或者是平均量。
图5、图7、图10、图11和图15显示了该种系统和方法。发送器506将超声波信号通过开关504发送至超声波探头502。超声波被传送到可能包括血管的人体对象。然后超声波返回到超声波探头502上,并被该探头转换成电子信号(步骤1505)。返回的超声波信号通过开关504被连结至接收器508。接收器对该信号进行处理,然后将射频信号输出至射频流体探测器510、B型图像处理器512(步骤1510)和多普勒频谱处理器524。多普勒频谱处理器计算来自从接收器508输出的信号的多普勒频谱,并以多普勒模式或者以与B型图像和/或彩色流体图像相结合的模式将此多普勒频谱输出到扫描转换器520。B型图像处理器512对射频信号进行处理,并输出B型图像信号。射频流体探测器510对射频信号进行处理,并将振幅c、功率c2、功率的b次幂cb或者这些值的组合这样的流体成分参数输出到增益控制信号发生器518。
图7显示了射频流体探测器510的示意图。该射频流体探测器包括高通滤波器702、速度计算器704、功率计算器706和方差计算器708。高通滤波器的其中一种形式可以是信号减法或者2-tap FIR,其中,随后传送的接收到的射频信号被从在先传送的射频信号中减去。这在两个或者更多个传送序列中进行高通滤波。在数次传送中进行高通滤波的滤波器去除了血管壁噪声,因而被称作壁滤波器。高通滤波器可使用在更高序列的FIR、IIR或者其他滤波器类型中从两个或者更多个传送中接收到的射频信号。其截止频率可由操作者调节的高通滤波器阻止低频的固定组织信号成分,只允许从接收器508中接收到的射频信号中的较高频率的流体信号成分通过(步骤1520),所述高通滤波器与速度计算器704、功率计算器706和方差计算器708相连结。
速度计算器704由经高通滤波的射频信号计算血流速度。速度被输出至扫描转换器520,由该扫描转换器将速度信号转换成电视光栅扫描速度图像。然后,此速度图像以常规的彩色流体模式显示在显示监视器522上。方差计算器708由经高通滤波的射频信号计算方差作为湍流指标。该方差被输出至扫描转换器520,由该扫描转换器将方差信号转换成电视光栅扫描方差图像。然后,将此方差图像以常规的彩色流体模式显示在显示监视器522上。功率计算器706利用下述公式由取样射频信号计算流体成分参数:振幅c、功率c2或者提高到cb的功率:
c = Σ i = 1 N | x i | - - - ( 6 )
c 2 = Σ i = 1 N x i 2 - - - ( 7 ) 以及
c b = Σ i = 1 N | x i | b - - - ( 8 )
其中,xi是经高通滤波的射频信号样本,而N是经高通滤波的信号样本的数目。
功率计算器706可以是DSP、FPGA、ASIC或者是例如乘法器、加法器、除法器以及绝对值计算器的离散部件。
可以选择计算这些值中的其中一个,而其他值可按以下方法从该首先计算得到的值获得。例如,可首先计算功率c2,而振幅c和提高到cb的功率可用下列公式由功率c2获得:
c 2 = Σ i = 1 N x i 2 - - - ( 9 )
c = c 2 - - - ( 10 ) 以及
c b = ( c 2 ) b - - - ( 11 )
通过式(10)和式(11)得到的流体成分参数不同于通过式(6)和式(8)得到的值,因为两者处理步骤的次序不相同(步骤1525)。上述信号振幅c、功率c2、提高到cb的功率这样的流体成分参数的计算反映的是总值。全部流体成分参数可以通过额外的计算,即通过将每个参数值除以样本数N而得到平均值的方法标准化。
增益控制发生器518接收流体成分参数:振幅c、功率c2、提高到cb的功率或者这些值的组合,并生成增益控制信号α。为增加信噪比(SNR),流体成分参数可以被组合成例如 d 1 c + d 2 c 2 + Σ i d i c b ( i ) 的组合流体成分参数,其中d1、d2和di是代表加权系数的实数,第i个成分b(i)同样也是实数。组合多于一个的流体成分参数通常会增加信噪比,并减少无关联的噪声。 d 1 c + d 2 c 2 + Σ i d i c b ( i ) 可以在功率计算器706中计算得到,然后也输出给增益控制发生器518中。功率计算器706也可以以彩色流体模式输出给扫描转换器。
增益控制发生器518可以是例如如图13中所示的数字信号处理器(DSP)1310,或者是如图10中所示的查找表(LUT),或者FPGA、ASIC或者诸如乘法器和加法器的离散部件。图11显示了DSP 1310或LUT 1010的示例响应。如果流体成分参数的值小,则增益控制信号α大,而高值可能产生接近于0的增益。图11所示的曲线是示例性的。可以使用其他预定曲线来抑制B型图像处理器512输出的信号(步骤1530)。增益控制发生器518输出增益控制信号α,该增益控制信号α被连结至信号组合器514。
信号组合器514可以是乘法器或者可变增益放大器。信号组合器514将B型图像与增益控制信号α相乘,并将B型图像输出至扫描转换器520,而扫描转换器则将图像信号转换成电视光栅扫描图像。然后,图像被显示在显示监视器522上(步骤1535)。
射频流体探测系统和方法可使用宽带超声波信号,因而能比多普勒流体探测器提供更高的空间解析度。
至此已描述了本发明的一个或多个实施例。然而应理解,在不脱离本发明精神和范围的前提下可以做出各种修改。相应地,其他实施例也被限定在所附权利要求的范围之中。

Claims (12)

1.一种当通过在每一位置对血管传送数次超声波信号并接收数次自血管返回的超声波信号而对血管成像时通过返回的超声波信号抑制杂波噪声表现的方法,该方法包括:
对返回的超声波图像信号进行处理以产生B型图像输出;
对返回的超声波信号进行解调以产生多普勒信号;
对多普勒信号进行壁滤波以仅让信号的流体成分通过;
计算流体成分参数;
基于流体成分参数生成增益控制信号α;以及
以所述增益控制信号α控制B型图像处理器的输出,其中,在B型图像处理器输出中的噪声表现得以抑制,并且其中
计算流体成分参数进一步包括:
由经壁滤波的多普勒信号计算振幅a、功率a2和提高到ab(i)的功率,这里幂b(i)是实数;
将振幅a、功率a2和提高到ab(i)的功率乘以相应的权重d1、d2和di,这里d1、d2和第i个分量di是实数;以及
对加权的振幅d1a、功率d2a2和提高到diab(i)的功率求和以作为组合流体成分参数
2.根据权利要求1所述的方法,其中,壁滤波进一步包括接收多普勒信号并利用高通滤波器响应对接收到的多普勒信号进行处理。
3.根据权利要求1所述的方法,其中,所述增益控制信号α是基于流体成分参数的预定曲线响应。
4.一种当通过在每一位置对血管传送数次超声波信号并接收数次自血管返回的超声波信号而对血管成像时通过返回的超声波信号抑制杂波噪声表现的方法,该方法包括:
对返回的超声波图像信号进行处理以产生B型图像输出;
对返回的超声波图像信号进行壁滤波以仅让信号的流体成分通过;
计算流体成分参数;
基于流体成分参数生成增益控制信号α;以及
以所述增益控制信号α控制B型图像处理器的输出,其中,在B型图像处理器输出中的噪声表现得以抑制,并且其中
计算流体成分参数进一步包括:
由经壁滤波的信号计算振幅c、功率c2和提高到cb(i)的功率,其中,幂b(i)是实数;
将振幅c、功率c2和提高到cb(i)的功率乘以相应的权重d1、d2和di,这里d1、d2和第i个分量di是实数;以及
对加权的振幅d1c、功率d2c2和提高到dicb(i)的功率求和以作为组合流体成分参数
Figure FSB00000522062800021
5.根据权利要求4所述的方法,其中,壁滤波进一步包括接收射频信号并利用高通滤波器响应对接收到的射频信号进行处理。
6.根据权利要求4所述的方法,其中,所述增益控制信号α是基于流体成分参数的预定曲线响应。
7.一种当通过在每一位置对血管传送数次超声波信号并接收数次自血管返回的超声波信号而对血管成像时通过返回的超声波信号抑制杂波噪声表现的系统,该系统包括:
接收器,其配置成接收返回的超声波图像信号并输出接收到的信号;
连结到接收器的B型图像处理器,其配置成从接收到的信号输出经B型处理的图像;
连结到接收器的多普勒流体探测器,其配置成对返回的超声波图像信号进行解调以产生多普勒信号,所述多普勒流体探测器具有壁滤波器,所述壁滤波器配置成对多普勒信号进行滤波并只输出所述多普勒信号的流体成分;
连结到多普勒流体探测器的增益控制发生器,其配置成计算流体成分参数并基于流体成分生成增益控制信号α;以及
连结到所述增益控制发生器和所述B型图像处理器的信号组合器,所述信号组合器配置成以所述增益控制信号α对经B型处理的图像输出进行修改,其中,在经B型处理的图像中的噪声表现得以抑制,并且其中
所述流体成分参数是由经壁滤波的多普勒信号得到的包括有振幅a、功率a2和提高到ab(i)的功率的组合流体成分参数这里幂b(i)是实数,振幅a、功率a2和提高到ab(i)的功率乘以相应的权重d1、d2和di并一起求和,这里d1、d2和第i个分量di是实数。
8.根据权利要求7所述的系统,其中所述壁滤波器响应是高通滤波器响应。
9.根据权利要求7所述的系统,其中,所述增益控制信号α是基于流体成分的预定曲线响应。
10.一种当通过在每一位置对血管传送数次超声波信号并接收数次自血管返回的超声波信号而对血管成像时通过返回的超声波信号抑制杂波噪声表现的系统,该系统包括:
接收器,其配置成接收返回的超声波图像信号并输出接收到的信号;
连结到接收器的B型图像处理器,其配置成从接收到的信号输出经B型处理的图像;
连结到接收器的射频流体探测器,具有壁滤波器的所述流体探测器配置成对返回的超声波射频信号进行滤波,并只输出接收到的信号的流体成分;
连结到射频流体探测器的增益控制发生器,其配置成计算流体成分参数并基于流体成分生成增益控制信号α;以及
连结到所述增益控制发生器和所述B型图像处理器的信号组合器,所述信号组合器配置成以所述增益控制信号α对所述经B型处理的图像输出进行修改,其中,在经B型处理的图像中的噪声表现得以抑制,并且其中
所述流体成分参数是由经壁滤波的射频信号得到的包括有振幅c、功率c2和提高到cb(i)的功率的组合流体成分参数
Figure FSB00000522062800031
这里幂b(i)是实数,振幅c、功率c2和提高到cb(i)的功率乘以相应的权重d1、d2和di并一起求和,这里d1、d2和第i个分量di是实数。
11.根据权利要求10所述的系统,其中所述壁滤波器响应是高通滤波器响应。
12.根据权利要求10所述的系统,其中,所述增益控制信号α是基于流体成分的预定曲线响应。
CN2008800102351A 2007-03-29 2008-03-28 用于超声波成像的方法和装置 Expired - Fee Related CN101646391B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US92063907P 2007-03-29 2007-03-29
US60/920,639 2007-03-29
US11/926,206 2007-10-29
US11/926,206 US8771190B2 (en) 2007-03-29 2007-10-29 Methods and apparatus for ultrasound imaging
PCT/JP2008/056727 WO2008123595A1 (en) 2007-03-29 2008-03-28 Methods and apparatus for ultrasound imaging

Publications (2)

Publication Number Publication Date
CN101646391A CN101646391A (zh) 2010-02-10
CN101646391B true CN101646391B (zh) 2011-08-24

Family

ID=39795580

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008800102351A Expired - Fee Related CN101646391B (zh) 2007-03-29 2008-03-28 用于超声波成像的方法和装置

Country Status (5)

Country Link
US (1) US8771190B2 (zh)
EP (1) EP2136713A4 (zh)
JP (1) JP5256210B2 (zh)
CN (1) CN101646391B (zh)
WO (1) WO2008123595A1 (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CA2785655A1 (en) * 2009-12-29 2011-07-07 Boston Scientific Scimed, Inc. Systems and methods for multi-frequency imaging of patient tissue using intravascular ultrasound imaging systems
JP5570877B2 (ja) * 2010-06-04 2014-08-13 株式会社東芝 超音波診断装置
JP5909498B2 (ja) * 2010-11-12 2016-04-26 ボストン サイエンティフィック サイムド,インコーポレイテッドBoston Scientific Scimed,Inc. カテーテルベースの撮像器及び同撮像器にて使用されるコンピュータ可読媒体
KR101313220B1 (ko) * 2010-11-23 2013-09-30 삼성메디슨 주식회사 특성 곡선 정보에 기초하여 컬러 도플러 모드 영상을 제공하는 초음파 시스템 및 방법
DE102011000090B4 (de) * 2011-01-11 2016-03-24 Leica Microsystems Cms Gmbh Verfahren und Einrichtung zur rastermikroskopischen Abbildung eines Objektes
JP6188594B2 (ja) * 2013-01-23 2017-08-30 東芝メディカルシステムズ株式会社 超音波診断装置、画像処理装置及び画像処理方法
CN103454640A (zh) * 2013-09-13 2013-12-18 北京东方惠尔图像技术有限公司 一种超声波成像装置和成像方法
WO2017035838A1 (zh) * 2015-09-06 2017-03-09 深圳迈瑞生物医疗电子股份有限公司 超声灰阶成像系统及方法
US11602328B2 (en) * 2015-12-30 2023-03-14 B-K Medical Aps Ultrasound flow imaging
EP4129194A1 (en) * 2021-08-03 2023-02-08 Koninklijke Philips N.V. Method and system for ultrasound parameter identification

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1135712A (zh) * 1994-09-22 1996-11-13 勃勒柯研究有限公司 超声波谱反衬成像
EP1529491A1 (en) * 2003-04-03 2005-05-11 Matsushita Electric Industrial Co., Ltd. Ultrasonograph and method for controlling ultrasonograph

Family Cites Families (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2604081A1 (fr) * 1986-09-19 1988-03-25 Labo Electronique Physique Dispositif d'exploration par echographie ultrasonore d'organes en mouvement et d'ecoulements sanguins
NO944736D0 (no) * 1994-12-07 1994-12-07 Vingmed Sound As Fremgangsmåte for bestemmelse av blodhastighet
US5871447A (en) * 1996-11-07 1999-02-16 Acuson Corporation Doppler energy-related parameters in an ultrasound imaging system
US5961460A (en) * 1997-04-11 1999-10-05 Acuson Corporation Ultrasound imaging enhancement methods and systems
JP4334032B2 (ja) 1998-02-23 2009-09-16 株式会社東芝 超音波診断装置
US6036643A (en) * 1998-05-14 2000-03-14 Advanced Technology Laboratories, Inc. Ultrasonic harmonic doppler imaging
JP4346147B2 (ja) * 1999-04-05 2009-10-21 株式会社東芝 超音波診断装置および超音波診断装置の作動方法
US6322505B1 (en) * 1999-06-08 2001-11-27 Acuson Corporation Medical diagnostic ultrasound system and method for post processing
JP2001079004A (ja) 1999-09-14 2001-03-27 Aloka Co Ltd 超音波診断装置
JP2002034987A (ja) 2000-07-19 2002-02-05 Ge Medical Systems Global Technology Co Llc Bモード画像生成方法および超音波診断装置
US6424275B1 (en) * 2001-01-18 2002-07-23 Scott R. Velazquez Linearity error compensator
US6599248B1 (en) * 2001-03-20 2003-07-29 Aloka Method and apparatus for ultrasound diagnostic imaging
KR100459616B1 (ko) * 2002-04-17 2004-12-04 주식회사 메디슨 도플러 효과를 이용하여 인체 조직 속도를 측정하기 위한초음파 진단 장치 및 방법
US6679843B2 (en) * 2002-06-25 2004-01-20 Siemens Medical Solutions Usa , Inc. Adaptive ultrasound image fusion

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1135712A (zh) * 1994-09-22 1996-11-13 勃勒柯研究有限公司 超声波谱反衬成像
EP1529491A1 (en) * 2003-04-03 2005-05-11 Matsushita Electric Industrial Co., Ltd. Ultrasonograph and method for controlling ultrasonograph

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JP特开2002-34987A 2002.02.05

Also Published As

Publication number Publication date
WO2008123595A1 (en) 2008-10-16
JP2010522582A (ja) 2010-07-08
EP2136713A1 (en) 2009-12-30
US20080242982A1 (en) 2008-10-02
CN101646391A (zh) 2010-02-10
US8771190B2 (en) 2014-07-08
EP2136713A4 (en) 2012-02-22
JP5256210B2 (ja) 2013-08-07

Similar Documents

Publication Publication Date Title
CN101646391B (zh) 用于超声波成像的方法和装置
US5971927A (en) Ultrasonic diagnostic apparatus for obtaining blood data
US5910118A (en) Removal of clutter components in ultrasound color-doppler imaging
KR20090042183A (ko) 합성 영상을 이용한 조직 도플러 영상 생성 장치 및 방법
KR20130076031A (ko) 결정 데이터에 기초하여 벡터 도플러 영상을 제공하는 초음파 시스템 및 방법
US6599248B1 (en) Method and apparatus for ultrasound diagnostic imaging
US20190072671A1 (en) 3-D Ultrasound Imaging with Multiple Single-Element Transducers and Ultrasound Signal Propagation Correction
US6322510B1 (en) Ultrasonic imaging method and apparatus
JP3093823B2 (ja) 超音波ドプラ診断装置
US20100010350A1 (en) Ultrasonic diagnostic apparatus and ultrasonic image acquisition method
US20070073152A1 (en) Systems and methods for acquiring images simultaneously
US20060122510A1 (en) Ultrasonic diagnostic apparatus and data processing method therefor
US6425868B1 (en) Ultrasonic imaging system
JP2004329609A (ja) 超音波診断装置
JP2002034985A (ja) 超音波診断装置
JP3281435B2 (ja) 超音波ドプラ診断装置
JPH0767449B2 (ja) 超音波ドプラ診断装置
JPH10118065A (ja) 超音波診断装置
JP2823252B2 (ja) 超音波診断装置
JP2002186615A (ja) 超音波診断装置
JP4537754B2 (ja) 超音波診断装置及びパルスドプラ計測装置
JP2563656B2 (ja) 超音波ドプラ映像装置
US11744552B2 (en) Ultrasound diagnostic apparatus, medical image processing apparatus, and medical image processing method
JP2011240006A (ja) 超音波診断装置
JPH03272751A (ja) 超音波診断装置

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
C56 Change in the name or address of the patentee

Owner name: HITACHI ALOKA MEDICAL LTD.

Free format text: FORMER NAME: AROCKAR K.K.

CP01 Change in the name or title of a patent holder

Address after: Tokyo, Japan

Patentee after: Hitachi Aloka Medical Ltd.

Address before: Tokyo, Japan

Patentee before: Arockar K. K.

CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20110824

Termination date: 20160328