CN102100567A - 彩色多普勒超声波诊断装置 - Google Patents

彩色多普勒超声波诊断装置 Download PDF

Info

Publication number
CN102100567A
CN102100567A CN2010105985463A CN201010598546A CN102100567A CN 102100567 A CN102100567 A CN 102100567A CN 2010105985463 A CN2010105985463 A CN 2010105985463A CN 201010598546 A CN201010598546 A CN 201010598546A CN 102100567 A CN102100567 A CN 102100567A
Authority
CN
China
Prior art keywords
mentioned
blood flow
orthogonal
orthogonal transformation
output
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
Application number
CN2010105985463A
Other languages
English (en)
Other versions
CN102100567B (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.)
Canon Medical Systems Corp
Original Assignee
Toshiba Corp
Toshiba Medical Systems Corp
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 Toshiba Corp, Toshiba Medical Systems Corp filed Critical Toshiba Corp
Publication of CN102100567A publication Critical patent/CN102100567A/zh
Application granted granted Critical
Publication of CN102100567B publication Critical patent/CN102100567B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/52Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/5215Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data
    • A61B8/5238Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of medical diagnostic data for combining image data of patient, e.g. merging several images from different acquisition modes into one image
    • 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/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/488Diagnostic techniques involving Doppler signals
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/08Detecting organic movements or changes, e.g. tumours, cysts, swellings
    • A61B8/0891Detecting organic movements or changes, e.g. tumours, cysts, swellings for diagnosis of blood vessels

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Surgery (AREA)
  • Pathology (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Animal Behavior & Ethology (AREA)
  • General Health & Medical Sciences (AREA)
  • Public Health (AREA)
  • Veterinary Medicine (AREA)
  • Hematology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

一种彩色多普勒超声波诊断装置。要解决的问题:收发宽频带超声波时杂波成分和血流成分的多普勒频带扩展而重叠,难以从血流成分中分离杂波成分,成为杂波成分多的血流图像。彩色多普勒超声波诊断装置具备:超声波探头;收发部,经由上述超声波探头在与被检体之间收发超声波;正交检波部,对收发部的输出进行正交检波;滤波单元,实施二维滤波,该二维滤波用于在深度方向上对正交检波部的输出进行第一正交变换,在整体方向上对第一正交变换的结果进行第二正交变换,从所得到的频率数据中提取血流成分,对二维滤波的输出进行第二正交变换的逆变换和第一正交变换的逆变换;及血流成分信息计算部,根据上述滤波单元的输出计算血流成分信息。

Description

彩色多普勒超声波诊断装置
本申请基于2009年12月21日年在日本提交的申请号为2009-289634号的专利申请,要求该在先申请的优先权,且将其全部内容均引入本申请。
技术领域
本发明涉及一种彩色多普勒超声波诊断装置。
背景技术
超声波诊断中的彩色多普勒法为,在同一方向上对生物体多次照射超声波,通过基于多普勒效应的偏移频率来提取血流成分信息。血流成分信息包含血流速度、功率和方差。例如,如图11所示,在同一方向(同一光栅)上以1/PRF(PRF:重复频率)的周期照射N次超声波。在同一方向上将与同一深度相对应的N个回波成分沿时间轴排列的数据列被称作整体(ensemble)方向或多普勒方向。深度方向也被称作距离方向。整体数据被进行傅立叶变换。由此求出多普勒频率。
组织信号(杂波成分信号)的强度比血流成分的信号强度大40~100dB左右。为了描绘在组织中流动的血流成分,必须抑制组织信号。被称作壁滤波器的高通滤波器抑制杂波成分。装置要求100dB以上的S/N(信噪比)。
生物体存在频率依存衰减。即使通过宽频带来发送超声波,回波信号的中心也成为低波频带。即使通过宽频带进行发送,在接收时也难以S/N(信噪比)良好地利用其能量。因此,在要求较高的S/N的彩色多普勒中,通过窄频带来发送超声波。
近年的装置实现了较高的S/N。即使是彩色多普勒,即使通过相当于B模式的宽频带来发送超声波,也能够微细地表现血流成分。宽频带化的距离分辨率优良。
但是,在彩色多普勒模式中通过宽频带发送超声波时,发生如下问题。
在设超声波的频率为fRF、音速为c,反射体以一定速度v在朝向探头的方向上移动的情况下,所得到的多普勒位移的频率fDop由公式(1)表示。
f Dop = 2 v c f RF - - - ( 1 )
即,多普勒频率与超声波的发送频率成正比例,例如通过2MHz的超声波所得到的多普勒频率与通过4MHz的超声波所得到的多普勒频率差两倍。图12上的图表示超声波频率和多普勒频率的关系。横轴是多普勒频率,纵轴是超声波发送频率。在某个多普勒频率、超声波频率下的多普勒信号的振幅,用未图示的与上述2个轴正交的第三轴来表现。血流成分的多普勒频率和杂波成分的多普勒频率均与超声波频率成正比例。
在现有技术中,仅对多普勒频率进行信号处理。因此,成为对于如下信号的处理,该信号为图12上方的频率轴被积分为一个而血流成分和杂波成分部分地重叠的信号。图12下方的图表的横轴是多普勒频率、纵轴是多普勒信号的振幅。即,即使血流成分和杂波成分被分离,在频率轴上血流成分和杂波成分的多普勒频率也重叠,并成为几乎不能分离。
在一直使用的窄频带的收发的情况下,仅存在于发送超声波的中心频率f0的频率附近,所以血流成分和杂波成分的多普勒频率不重叠。但是,当进行宽频带的收发时,有时血流成分和多普勒频率重叠。于是,产生的问题为,通过壁滤波器不能够分离杂波成分和血流成分,成为杂波成分较多的图像。
并且,在通过中心频率f0宽频带地进行了发送的情况下,实际的接收信号受到生物体中的频率依存衰减的影响,而如图13所示那样以比f0低的频率f1接收,但是一般不能够得知f1。因此,将公式(1)变形,
v = 2 c f RF f Dop - - - ( 2 )
在公式(2)中,由于fRF不确定,所以不能够求出正确的速度v。在一直使用的窄频带的收发的情况下,发送频率f0和接收频率f1几乎相同,所以没有这样的问题。
另一方面,专利文献1公开了利用两种超声波收发频率的方法。在该专利文献中公开的方法为:使用不同的超声波收发频率,利用同一速度的血流成分的多普勒频率不同的情况,推断超过折返速度的速度。但是,没有提及从血流成分中分离杂波成分的情况。
专利文献1:日本特许第2953083号公报
发明内容
本发明的目的在于解决如下技术问题:在对宽频带的超声波进行了收发的情况下,杂波成分的多普勒波频带和血流成分的多普勒波频带重叠,难以从血流成分中分离杂波成分,由此成为杂波成分较多的血流图像。
本发明的彩色多普勒超声波诊断装置具备:超声波探头;收发部,经由上述超声波探头在与被检体之间收发超声波;正交检波部,对上述收发部的输出进行正交检波;滤波单元,实施二维滤波,该二维滤波用于在深度方向上对上述正交检波部的输出进行第一正交变换,在整体方向上对上述第一正交变换的结果进行第二正交变换,从所得到的频率数据中提取血流成分,并且该滤波单元对上述二维滤波的输出依次进行上述第二正交变换的逆变换和上述第一正交变换的逆变换;以及血流成分信息计算部,根据上述滤波单元的输出计算血流成分信息。
即使在对宽频带的超声波进行了收发的情况下,也能够从血流成分中有效地分离杂波成分而得到抑制了杂波成分的血流图像。
附图说明
图1是示出第一实施方式的彩色多普勒超声波诊断装置的结构的图。
图2是示出图1的二维滤波单元的处理步骤的图。
图3是示出接收了从图1的接收电路输出的正交检波的复信号(IQ信号)的二维频率特性的一个例子的图。
图4是示出第一实施方式的第一变形例的处理步骤的流程图。
图5是示出第一实施方式的第二变形例的处理步骤的流程图。
图6是示出图5的BPF的三种通过频带的关系的图。
图7是示出第二实施方式的二维滤波单元的处理步骤的流程图。
图8是示出在图7中基于第一正交变换的傅立叶频率fn和速度Vn的关系的图。
图9是示出图8的折返解除后的速度Vn的图。
图10是示出第一实施方式的第二变形例的二维滤波单元的结构的图。
图11是示出包的图。
图12是示出在通过宽频带超声波进行了收发时的以往的带通滤波中不能够从血流成分中分离杂波成分的问题的图。
图13是将生物体中的频率依存衰减的接收频率与发送频率一起示出的图。
符号的说明:
1超声波探头;2发送电路;3控制单元;4接收电路;5B模式处理单元;6高通滤波器;7二维滤波单元;8彩色多普勒处理单元;9坐标变换单元;10图像合成单元;11监视器。
具体实施方式
根据一个实施方式,彩色多普勒超声波诊断装置具备:超声波探头;收发部,经由超声波探头在与被检体之间收发超声波;正交检波部,对收发部的输出进行正交检波;滤波单元,实施二维滤波,该二维滤波用于在深度方向上对正交检波部的输出进行第一正交变换,在整体方向上对第一正交变换的结果进行第二正交变换,从所得到的频率数据中提取血流成分,对二维滤波的输出依次进行第二正交变换的逆变换和第一正交变换的逆变换;以及血流成分信息计算部,根据滤波单元的输出计算血流成分信息。
下面,参照附图对本实施方式进行说明。
图1示出第一实施方式的彩色多普勒超声波诊断装置的结构。超声波探头1具有多个压电振子。在压电振子的前侧配置有音响匹配层。在压电振子的背侧配置有背衬部件等。多个压电振子排列成二维或三维状,以便能够进行电气的二维或三维扫描。发送电路2对探头1的压电振子施加脉冲信号,由此产生超声波。超声波探头1的压电振子将来自被检体的超声波回波变换为电信号。
发送电路2构成为,能够根据操作者的模式选择,在控制单元3的控制下,将发送超声波的频带在比较宽频带和比较窄频带之间切换。在B模式下选择比较宽频带的超声波。在第一彩色多普勒模式下选择比较窄频带的超声波。在第二彩色多普勒模式下选择比较宽频带的超声波。在第二彩色多普勒模式下使用的频带也可以比在B模式下使用的频带为宽频带。另外,在第一彩色多普勒模式下选择的频带与现有的彩色多普勒模式的频带等价。在第一彩色多普勒模式下,通过与现有相同的带通滤波器从血流成分中分离杂波成分。在第二彩色多普勒模式下,通过本实施方式中特征性的二维滤波从血流成分中分离杂波成分。详细情况将后述。
在此,对超声波的宽频带化进行说明。发送频带f0的倒数为发送周期T0,将具有超声波1个周期T0的时间长度的超声波脉冲称作突发1波脉冲,将具有2个周期“2·T0”的时间长度的超声波脉冲称作突发2波脉冲,将具有M个周期“M·T0”的时间长度的超声波脉冲称作突发M波脉冲,并将1波、2波、…M波称作突发波数。对发送周期乘以突发波数而得到的值相当于脉冲长度。典型地说,宽频带通过发送频率为f0、突发波数小于3的突发波数M的超声波来产生。窄频带通过发送频率为f0、突发波数为3以上的突发波数M的超声波来产生。另外,宽频带不一定仅由突发波数小于3的脉冲来定义。例如,也可以由突发波的包络线上的半值宽度ΔW和频率T之比“ΔW/T”来定义。在使用了该比的情况下,宽频带成为“ΔW/T”=2.2以下。此外,也可以由超声波的频谱的比频带=ΔW/fc(ΔW:半值宽度;fc:频普中心频率)来定义宽频带的范围。在该情况下,比频带=ΔW/fc=0.3以上相当于“宽频带”。
从超声波探头1向被检体P发送的超声波,不断地被体内组织的声阻抗的非连续面反射。回波信号被超声波探头1接收。该回波的振幅依存于进行反射的非连续面的声阻抗之差。此外,在所发送的超声波脉冲被移动中的血流以及心脏壁等的表面反射时的回波,由于多普勒效应而依存于移动体的超声波发送方向的速度成分,而受到频率偏移。
发送电路2具有脉冲产生器、发送延迟部和脉冲发生器。脉冲产生器为,以规定的速率频率fr Hz(周期:1/fr秒)重复发生用于生成发送超声波的速率脉冲。发送延迟部为,在每个波道将超声波聚焦为波束状,且将为了决定发送指向性所需的延迟时间赋予各波道的速率脉冲。脉冲发生器在基于速率脉冲的定时,在每个波道对探头1施加驱动脉冲。
接收电路4具有前置放大器、A/D变换器、接收延迟部、加法器以及正交检波电路。前置放大器在每个波道对经由探头1取入的回波信号进行放大。接收延迟部对所放大的回波信号赋予为了决定接收指向性所需的延迟时间,之后在加法器中进行加法处理。通过该加法,来自与回波信号的接收指向性相对应的方向的反射成分被强调,根据接收指向性和发送指向性而形成超声波收发的综合性波束。正交检波电路在第一、第二彩色多普勒模式下起作用。正交检波电路对回波信号混合与发送超声波同相的参照信号,并且对回波信号混合与发送超声波相差π/2相位的参照信号。以下,将该正交检波信号称作IQ信号。
B模式处理单元5通过对数放大、包络线检波处理等处理,从接收电路4的回波信号中产生B模式图像的数据。
在接收电路4的输出上连接高通滤波器6和二维滤波单元7。通过控制单元3的控制,在第一彩色多普勒模式下高通滤波器6起作用,在第二彩色多普勒模式下二维滤波单元7起作用。高通滤波器6通过与现有相同的被称作壁滤波器的时间轴上的一维高通滤波器来抑制杂波成分,提取血流成分。关于二维滤波单元7将后述。
彩色多普勒处理单元8取出从高通滤波器6或二维滤波单元7接受的回波信号的基于多普勒效应的偏移频率,作为移动体主要提取血流成分,并对于多个点逐个地求出平均速度、方差、功率等血流数据。被称作数字扫描转换器(DSC)的坐标变换单元9,将平均速度图像、方差图像、功率图像等血流图像的与超声波扫描顺序相对应的坐标系,变换为一般的正交坐标系。血流图像在图像合成单元10中被与B模式图像合成,并显示在监视器11上。
以下,对二维滤波单元7进行详细说明。本实施方式的二维滤波单元7为,即使在选择了通过高通滤波器6不能够分离血流成分和杂波成分的程度的宽频带时,也实现血流成分和杂波成分的分离。高通滤波器6仅通过多普勒频率的一维进行观测,但当着眼于超声波频率和多普勒频率的二维滤波时,能够将通过高通滤波器6不能分离的血流成分和杂波成分也进行分离。
在本实施方式中,实施如下的滤波:在深度方向上进行第一正交变换,在用于多普勒检测的包方向上进行第二正交变换,由此分割为二维正交空间,并在该二维空间上区分血流成分和杂波成分。作为第一正交变换,在乘以窗口函数之后进行离散傅立叶变换。作为第二正交变换,进行离散傅立叶变换、向正交多项式空间的变换、离散余弦变换、卡洛(Karhunen-Loeve)变换等变换。在二维空间中实施了滤波之后,进行第二正交变换的逆变换,并且进行上述第一正交变换的逆变换,并在全部返回时间轴上之后,通过彩色多普勒处理单元8推断速度、功率、方差。彩色多普勒处理单元8的处理在第一、第二彩色模式下共通。还能够不进行第一正交变换的逆变换,而对超声波频率的每个频带推断速度、功率、方差,尤其关于速度利用与超声波频率成正比例的情况,而成为更高精度的值。
在进行了宽频带的收发的情况下的彩色多普勒超声波诊断装置中,在本实施方式中,通过增加称为超声波频率的又一个维,能够对从在一维方向上观察的多普勒频率中不能分离的血液成分和杂波成分进行分离,能够抑制在现有技术的方法中含有较多杂波成分的图像的杂波成分。
并且,在进行了宽频带的收发的情况下,由于存在生物体的频率依存衰减,所以一般难以得知接收中心频率,因此难以推断正确的速度值,但是在本实施方式中,由于分割为每个频率来求出速度值,所以能够推断正确的速度值。并且,能够推断比折返速度快的速度。
如上所述,如图3所示那样,在接收电路4的正交检波电路中,对接收信号被以频率f0混合而产生的IQ信号的杂波成分和血流成分进行分离。
如图2所示,首先,在深度方向上进行第一正交变换(S1)。接着,对于进行了第一正交变换的结果在整体方向上进行第二正交变换(S2)。在第一及第二正交变换均为离散傅立叶变换的情况下,其结果如图3所示。在RF频率f0即IQ频率为0下的壁滤波器的截止频率在多普勒频率上为fc0的情况下,构成IQ频率f中的截止频率fc(f)成为
fc ( f ) = fc 0 f 0 ( | f | + f 0 ) - - - ( 3 )
的高通滤波器(S3)。即,在图3中,为将杂波部分切除的滤波。在此所说的滤波是在频率轴上的处理,所以对二维的各频率成分乘以适当的系数。接着,在返回时间区域时,首先进行第二正交变换的逆变换(S4),接着进行第一正交变换的逆变换(S5)。由此,信号成为除去了杂波成分的IQ信号,所以之后能够通过与现有技术相同的方式,通过自相关法推断速度、方差、功率(S6)。
在第一及第二正交变换为离散傅立叶变换的情况下,也可以代替分别进行一维变换,而进行二维离散傅立叶变换。逆变换也能够使用二维离散逆傅立叶变换。在离散傅立叶变换的情况下,在变换之前需要窗口函数,所以在乘以规定的窗口函数之后进行变换。
在整体方向上通常仅存在6~20左右的数据。当针对这样较小的数据尺寸乘以窗口函数并进行离散傅立叶变换时,血流成分的灵敏度降低。此外,在离散傅立叶变换中,原理上不能够分离观测时间以上的低频成分,所以不能够充分除去由低频成分构成的杂波成分。
因此,作为第二正交变换进行基于正交多项式空间的变换。所谓正交多项式是对勒让德(Le Gendre)多项式乘以适当的系数之后的多项式作为正交基底的多项式。利用正交多项式空间构成壁滤波器的方法,记载在以下的文献中。
Steinar Bjaerun,Hans Torp,“Clutter filter design for ultrasound color flow imaging”,IEEE Transaction on ultrasonics,ferroelectrics,and frequency control,Vol.49,No.2,pp.204-216,Feb.,2002
正交多项式空间的优点是,不需要为频率函数,所以不需要窗口函数。并且,所谓一次系数是直线成分,所以在观测时间内能够检测出与具有某个斜度的直线能够近似的、较迟运动的杂波成分,所以能够除去这种较迟的杂波成分。该正交空间与傅立叶频率空间不同,所以不能够如图6那样表现,但正交多项式空间的频率(在此将正交变换后的结果称作频率。将通常意义上的频率记为傅立叶频率)在傅立叶空间具有怎样的频率,如果使用Z变换则能够简单地赋予对应。根据该对应,如果在“深度方向傅立叶频率+整体方向正交多项式频率”的二维空间中,实施在傅立叶频率轴上具有成为fc(f)的截止频率的滤波器,则能够实施有效地除去杂波成分的滤波。
除此之外,作为第二正交变换,能够使用离散余弦变换。离散余弦变换也未假设周期函数,所以不需要窗口函数。此外,作为第二正交变换,能够使用卡洛(Karhunen-Loeve)变换。卡洛(Karhunen-Loeve)变换是根据信号的统计性质来分解为包含最大的主成分的正交成分的方法。由于正交变换一定存在逆变换,所以能够在正交空间的频率轴上实施了滤波之后进行逆变换,而返回到原来的时间轴。
图4示出与第一实施方式相关的第一变形例的步骤。该第一变形例与第一实施方式的不同点在于,不进行第二正交变换(S2),而是对第一正交变换输出的各频率按照每个距离进行滤波,该滤波是在时间轴上进行与第二正交空间上的频率轴上的滤波等价的处理(S7、S8、S9)。在该方式中,第一正交变换是进行深度方向的离散傅立叶变换,并对各个频率,加上改变截止频率而实施了现有技术的壁滤波器(时间轴上的壁滤波器)后的输出(S10),之后进行离散逆傅立叶变换。第一实施方式和该第一变形例是数学上等价的处理。
图5示出与第一实施方式相关的第二变形例的步骤。该第二变形例与第一实施方式的不同点在于,不进行第一、第二正交变换(S1、S2),而全部是时间轴的处理。如图6、图10所示,由通过频带不同的多个带通滤波器(BPF),将IQ信号按照每个距离分割为第一正交变换的频带的组(S11、S12、S13)。对该分割后的信号,之后进行与第一实施方式的第一变形例相同的处理。图6在频率轴上示出了该情况。
(第二实施方式)
本实施方式的结构与第一实施方式相同。图7示出本实施方式的二维滤波处理的步骤。到进行第二正交变换的逆变换(S1~S4)为止与第一实施方式相同。
在进行了第二正交变换的逆变换之后,不实施第一正交变换的逆变换。保持在第一频率轴上的状态,在各个频率成分中推断速度、方差、功率(S21、S22、S23)。速度、方差、功率是多普勒频率轴上的参数,所以能够通过与第一实施方式相同的处理来推断。
在此,使第一正交变换为离散傅立叶变换。另外,不限定第二正交变换的方法,可以是任意的正交变换方法。在第一正交变换的傅立叶频率为fn的情况下,假设能够得到速度Vn、方差Tn、功率Pn。其中,Vn为将自相关函数的偏角标准化到-0.5~+0.5的范围内的速度。当将fn和Vn描绘在图表上时,如图8所示产生折返。当将该折返解除(展开)时,成为如图9所示。当通过最小二乘法求图9的直线的斜度而得到了a这个值时,由于用PRF除公式(1)两边得到的值为标准化速度,且fRF=f0+fn,所以
vn = 2 v c · PRF ( f 0 + fn ) - - - ( 4 )
由此,图10(b)的直线的斜度a为
a = 2 v c · PRF - - - ( 5 )
由此,通过以下公式求出速度v。
v = a · c · PRF 2 - - - ( 6 )
该速度v是通过最小二乘法根据将宽频带的接收回波分割为窄频带而得到的各多普勒频率求得的速度,与根据宽频带的多普勒信号直接求得的速度相比精度更高,同时能够得到折返速度以上的速度。即,关于最终的速度、方差、功率为,通过最小二乘法根据每个傅立叶频率的速度值来求得直线的斜度,并根据该斜度求得最终的速度。在图8中,在正确得到展开的程度为止能够正确得到较快的流速。
作为与第二实施方式相关的第一、第二变形例,可以考虑分别在第一实施方式的第一、第二变形例中,将“实施利用了第二正交变换的时间轴上的滤波”的块之后置换为第二实施方式的“通过自相关法推断速度、方差、功率”之后的块的方法。
虽然上面对特定的实施例进行了说明,但这些实施例仅为例子,而并不试图限定本发明的范围。以上说明的新实施方式能够通过各种其他的方式来实施,进而在不脱离本发明的精神的范围内,能够对上面说明的实施方式中的方式进行各种省略、替代和变更。所附的申请范围和与其等同的技术方案覆盖这样的方式或变形例,其落入本发明的范围中。

Claims (12)

1.一种彩色多普勒超声波诊断装置,其特征在于,具备:
超声波探头;
收发部,经由上述超声波探头在与被检体之间收发超声波;
正交检波部,对上述收发部的输出进行正交检波;
滤波单元,实施二维滤波,该二维滤波用于在深度方向上对上述正交检波部的输出进行第一正交变换,在整体方向上对上述第一正交变换的结果进行第二正交变换,从所得到的频率数据中提取血流成分,并且该滤波单元对上述二维滤波的输出依次进行上述第二正交变换的逆变换和上述第一正交变换的逆变换;以及
血流成分信息计算部,根据上述滤波单元的输出计算血流成分信息。
2.根据权利要求1所述的彩色多普勒超声波诊断装置,其特征在于,
上述第一正交变换和上述第二正交变换是离散傅立叶变换。
3.根据权利要求1所述的彩色多普勒超声波诊断装置,其特征在于,
通过离散二维傅立叶变换和离散二维逆傅立叶变换,一次性地进行上述第一、第二正交变换和逆变换。
4.根据权利要求1所述的彩色多普勒超声波诊断装置,其特征在于,
上述第一正交变换采用离散傅立叶变换,上述第二正交变换采用正交多项式。
5.根据权利要求1所述的彩色多普勒超声波诊断装置,其特征在于,
上述第一正交变换采用离散傅立叶变换,上述第二正交变换采用离散余弦变换。
6.根据权利要求1所述的彩色多普勒超声波诊断装置,其特征在于,
上述第一正交变换采用离散傅立叶变换,上述第二正交变换采用卡洛变换。
7.一种彩色多普勒超声波诊断装置,其特征在于,具备:
超声波探头;
收发部,经由上述超声波探头在与被检体之间收发超声波;
正交检波部,对上述收发部的输出进行正交检波;
滤波单元,在深度方向上对上述正交检波部的输出进行正交变换,为了提取血流成分,在整体方向上按照每个距离对上述正交变换的结果实施滤波特性不同的时间轴上的滤波,对上述滤波器的输出进行上述正交变换的逆变换;以及
血流成分信息计算部,根据上述滤波单元的输出计算血流成分信息。
8.一种彩色多普勒超声波诊断装置,其特征在于,具备:
超声波探头;
收发部,经由上述超声波探头在与被检体之间收发超声波;
正交检波部,对上述收发部的输出进行正交检波;
滤波单元,在深度方向上对上述正交检波部的输出实施带通滤波,为了从上述带通滤波的结果提取血流成分,在整体方向上按照每个距离实施滤波特性不同的时间轴上的滤波;以及
血流成分信息计算部,根据上述滤波单元的输出计算血流成分信息。
9.一种彩色多普勒超声波诊断装置,其特征在于,具备:
超声波探头;
收发部,经由上述超声波探头在与被检体之间收发超声波;
正交检波部,对上述收发部的输出进行正交检波;
滤波单元,实施二维滤波,该二维滤波用于在深度方向上对上述正交检波部的输出进行第一正交变换,在整体方向上对上述第一正交变换的结果进行第二正交变换,从所得到的频率数据中提取血流成分,并且该滤波单元对上述二维滤波的输出进行上述第二正交变换的逆变换;以及
血流成分信息计算部,根据从上述滤波单元输出的基于上述第一正交变换的频率成分计算血流成分信息。
10.一种彩色多普勒超声波诊断装置,其特征在于,具备:
超声波探头;
收发部,经由上述超声波探头在与被检体之间收发超声波;
正交检波部,对上述收发部的输出进行正交检波;
滤波单元,实施二维滤波,该二维滤波用于在深度方向上对上述正交检波部的输出进行第一正交变换,在整体方向上对上述第一正交变换的结果进行第二正交变换,从所得到的频率数据中提取血流成分,并且该滤波单元对上述二维滤波的输出进行上述第二正交变换的逆变换;以及
血流成分信息计算部,根据上述滤波单元的输出,对上述第一正交变换的成分推断包含速度、方差、功率在内的血流成分信息,并根据这些信息来推断最终的速度、方差、功率。
11.根据权利要求10所述的彩色多普勒超声波诊断装置,其特征在于,
上述血流成分信息包含血流速度,上述血流速度为,通过最小二乘法根据每个傅立叶频率的速度值求出直线的斜度,并根据该斜度求出最终的速度。
12.一种彩色多普勒超声波诊断装置,其特征在于,具备:
超声波探头;
收发部,经由上述超声波探头在与被检体之间收发超声波;
正交检波部,对上述收发部的输出进行正交检波;
滤波单元,在深度方向上对上述正交检波部的输出实施带通滤波,为了从上述带通滤波的结果提取血流成分,在整体方向上按照每个距离实施滤波特性不同的时间轴上的滤波;以及
血流成分信息计算部,根据上述滤波单元的输出,对各频带推断速度、方差、功率等的血流成分信息,并根据这些信息来推断最终的速度、方差、功率。
CN 201010598546 2009-12-21 2010-12-21 彩色多普勒超声波诊断装置 Active CN102100567B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP289634/2009 2009-12-21
JP2009289634A JP5398514B2 (ja) 2009-12-21 2009-12-21 カラードプラ超音波診断装置

Publications (2)

Publication Number Publication Date
CN102100567A true CN102100567A (zh) 2011-06-22
CN102100567B CN102100567B (zh) 2013-12-25

Family

ID=44152056

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201010598546 Active CN102100567B (zh) 2009-12-21 2010-12-21 彩色多普勒超声波诊断装置

Country Status (3)

Country Link
US (1) US8858446B2 (zh)
JP (1) JP5398514B2 (zh)
CN (1) CN102100567B (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106466192A (zh) * 2015-08-17 2017-03-01 富士胶片株式会社 超声波诊断装置和多普勒波形图像生成方法
JP2017153592A (ja) * 2016-02-29 2017-09-07 東芝メディカルシステムズ株式会社 超音波診断装置、信号処理装置および信号処理プログラム
CN108354629A (zh) * 2018-03-30 2018-08-03 苏州佳世达电通有限公司 超音波成像方法
CN109963513A (zh) * 2016-11-17 2019-07-02 皇家飞利浦有限公司 用于使用闪烁伪影来检测肾结石的超声系统和方法
RU220404U1 (ru) * 2023-07-13 2023-09-12 Федеральное государственное бюджетное учреждение "Национальный медицинский исследовательский центр радиологии" Министерства здравоохранения Российской Федерации (ФГБУ "НМИЦ радиологии" Минздрава России) Устройство для определения возможности УЗИ аппарата использовать эффект артефакта мерцания для обнаружения мочевых камней

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20140350406A1 (en) * 2013-05-24 2014-11-27 Siemens Medical Solutions Usa, Inc. Dynamic Operation for Subarrays in Medical Diagnostic Ultrasound Imaging
CN107106142B (zh) * 2014-12-26 2020-08-04 古野电气株式会社 超声波体组织检测装置与方法、以及存储介质
CN111352433B (zh) * 2018-12-20 2021-04-06 中国科学院沈阳自动化研究所 一种无人机水平姿态角的故障诊断方法
JP7191789B2 (ja) * 2019-08-08 2022-12-19 富士フイルム株式会社 超音波観察ユニット、超音波診断装置、超音波振動子の駆動方法、超音波振動子の駆動プログラム
US20220054104A1 (en) * 2020-07-30 2022-02-24 Canon Medical Systems Corporation Ultrasonic diagnostic device and image processing device

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1039902A (zh) * 1988-07-26 1990-02-21 菲利浦光灯制造公司 傅里叶处理脉冲多普勒信号的方法和装置
CN1133166A (zh) * 1995-04-10 1996-10-16 深圳安科高技术有限公司 一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法及仪器
US5910118A (en) * 1996-10-01 1999-06-08 Kabushiki Kaisha Toshiba Removal of clutter components in ultrasound color-doppler imaging
CN101022765A (zh) * 2004-08-31 2007-08-22 华盛顿州大学 用于评估狭窄血管的壁振动的超声技术
CN101347343A (zh) * 2007-07-17 2009-01-21 株式会社东芝 超声波图像取得装置
US20090209861A1 (en) * 2005-02-17 2009-08-20 Panasonic Corporation Ultrasonic Doppler Blood Flow Meter

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5477858A (en) * 1986-07-30 1995-12-26 Siemens Medical Systems, Inc. Ultrasound blood flow/tissue imaging system
JP3029315B2 (ja) * 1990-03-30 2000-04-04 株式会社東芝 超音波血流イメージング装置
JP4801912B2 (ja) * 2004-03-24 2011-10-26 株式会社東芝 超音波診断装置
JP4634738B2 (ja) * 2004-04-30 2011-02-16 株式会社東芝 超音波診断装置
JP4080490B2 (ja) * 2005-03-04 2008-04-23 アロカ株式会社 超音波診断装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1039902A (zh) * 1988-07-26 1990-02-21 菲利浦光灯制造公司 傅里叶处理脉冲多普勒信号的方法和装置
CN1133166A (zh) * 1995-04-10 1996-10-16 深圳安科高技术有限公司 一种使用经颅多普勒超声技术检测颅内脑血流状态分布的方法及仪器
US5910118A (en) * 1996-10-01 1999-06-08 Kabushiki Kaisha Toshiba Removal of clutter components in ultrasound color-doppler imaging
CN101022765A (zh) * 2004-08-31 2007-08-22 华盛顿州大学 用于评估狭窄血管的壁振动的超声技术
US20090209861A1 (en) * 2005-02-17 2009-08-20 Panasonic Corporation Ultrasonic Doppler Blood Flow Meter
CN101347343A (zh) * 2007-07-17 2009-01-21 株式会社东芝 超声波图像取得装置

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106466192A (zh) * 2015-08-17 2017-03-01 富士胶片株式会社 超声波诊断装置和多普勒波形图像生成方法
JP2017153592A (ja) * 2016-02-29 2017-09-07 東芝メディカルシステムズ株式会社 超音波診断装置、信号処理装置および信号処理プログラム
CN109963513A (zh) * 2016-11-17 2019-07-02 皇家飞利浦有限公司 用于使用闪烁伪影来检测肾结石的超声系统和方法
CN109963513B (zh) * 2016-11-17 2023-03-21 皇家飞利浦有限公司 用于使用闪烁伪影来检测肾结石的超声系统和方法
CN108354629A (zh) * 2018-03-30 2018-08-03 苏州佳世达电通有限公司 超音波成像方法
RU220404U1 (ru) * 2023-07-13 2023-09-12 Федеральное государственное бюджетное учреждение "Национальный медицинский исследовательский центр радиологии" Министерства здравоохранения Российской Федерации (ФГБУ "НМИЦ радиологии" Минздрава России) Устройство для определения возможности УЗИ аппарата использовать эффект артефакта мерцания для обнаружения мочевых камней

Also Published As

Publication number Publication date
JP2011125635A (ja) 2011-06-30
US20110152689A1 (en) 2011-06-23
US8858446B2 (en) 2014-10-14
JP5398514B2 (ja) 2014-01-29
CN102100567B (zh) 2013-12-25

Similar Documents

Publication Publication Date Title
CN102100567B (zh) 彩色多普勒超声波诊断装置
JP2002539877A (ja) コード化送信パルスを用いた医療診断用超音波イメージングシステム
JP4931910B2 (ja) 超音波画像装置
EP2154548A1 (en) Ultrasound diagnostic apparatus for nonlinear contrast imaging
JP2003038490A (ja) 位相反転超音波画像処理システムおよび方法
JP2001286472A (ja) 超音波診断装置
EP2169426A1 (en) Ultrasound signal processing apparatus and method for determining attenuation
JP3093823B2 (ja) 超音波ドプラ診断装置
US7803114B2 (en) Ultrasonic diagnostic apparatus and data processing method therefor
JP4698003B2 (ja) 超音波診断装置
JP2007507271A (ja) 超音波イメージングにおけるスモールアンサンブル長によるクラッタフィルタリング
JP5355327B2 (ja) 超音波診断装置
JP4627675B2 (ja) 超音波診断装置
JP2916219B2 (ja) 超音波診断装置
JP6780509B2 (ja) 超音波診断装置の制御装置、超音波診断装置およびクラッター成分低減方法
WO2017047328A1 (ja) 超音波診断装置、及び超音波撮像方法
JP4568080B2 (ja) 超音波診断装置
JP3465197B2 (ja) 超音波診断装置
JP2008148794A (ja) 画像生成方法および超音波診断装置
JP2004041324A (ja) 超音波診断装置
JP4219626B2 (ja) 超音波断層画像生成装置
JPH05237103A (ja) 超音波診断装置
JPH05200024A (ja) 超音波ドプラ診断装置
JPS63311938A (ja) 超音波血流計測装置
Udesena et al. A frequency splitting method for blood flow estimation.

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
C41 Transfer of patent application or patent right or utility model
TR01 Transfer of patent right

Effective date of registration: 20160718

Address after: Japan Tochigi

Patentee after: Toshiba Medical System Co., Ltd.

Address before: Tokyo, Japan, Japan

Patentee before: Toshiba Corp

Patentee before: Toshiba Medical System Co., Ltd.