CN104902824A - 超声波诊断装置 - Google Patents

超声波诊断装置 Download PDF

Info

Publication number
CN104902824A
CN104902824A CN201380070054.9A CN201380070054A CN104902824A CN 104902824 A CN104902824 A CN 104902824A CN 201380070054 A CN201380070054 A CN 201380070054A CN 104902824 A CN104902824 A CN 104902824A
Authority
CN
China
Prior art keywords
assessed value
dimension
orientation
transducer
computing unit
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
CN201380070054.9A
Other languages
English (en)
Other versions
CN104902824B (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
Hitachi Aloka Medical Ltd
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 Hitachi Aloka Medical Ltd filed Critical Hitachi Aloka Medical Ltd
Publication of CN104902824A publication Critical patent/CN104902824A/zh
Application granted granted Critical
Publication of CN104902824B publication Critical patent/CN104902824B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/5207Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving processing of raw data to produce diagnostic data, e.g. for generating an image
    • 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/5269Devices using data or image processing specially adapted for diagnosis using ultrasonic, sonic or infrasonic waves involving detection or reduction of artifacts
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4488Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer the transducer being a phased array
    • 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/8909Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration
    • G01S15/8915Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array
    • G01S15/8925Short-range imaging systems; Acoustic microscope systems using pulse-echo techniques using a static transducer configuration using a transducer array the array being a two-dimensional transducer configuration, i.e. matrix or orthogonal linear arrays
    • 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/52046Techniques for image enhancement involving transmitter or receiver
    • G01S7/52047Techniques for image enhancement involving transmitter or receiver for elimination of side lobes or of grating lobes; for increasing resolving power
    • 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
    • GPHYSICS
    • G10MUSICAL INSTRUMENTS; ACOUSTICS
    • G10KSOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
    • G10K11/00Methods or devices for transmitting, conducting or directing sound in general; Methods or devices for protecting against, or for damping, noise or other acoustic waves in general
    • G10K11/18Methods or devices for transmitting, conducting or directing sound
    • G10K11/26Sound-focusing or directing, e.g. scanning
    • G10K11/34Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
    • G10K11/341Circuits therefor
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/13Tomography
    • A61B8/14Echo-tomography
    • A61B8/145Echo-tomography characterised by scanning multiple planes
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/44Constructional features of the ultrasonic, sonic or infrasonic diagnostic device
    • A61B8/4483Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer
    • A61B8/4494Constructional features of the ultrasonic, sonic or infrasonic diagnostic device characterised by features of the ultrasound transducer characterised by the arrangement of the transducer elements

Abstract

评估值计算单元(50)基于从延迟处理器(30)获得的多个接收的波信号,来评估整相的程度,从而计算与多个振荡元件(12)的二维阵列相关的二维评估值。评估值计算单元(50)根据在x方向上的一维评估值和在y方向上的一维评估值,来获得xy平面上的二维评估值。乘法单元(60)将二维评估值和求和处理器(40)输出的接收波束相乘,以调整接收的波束的增益。非必要的信号分量因而被减少。

Description

超声波诊断装置
技术领域
本发明涉及一种超声波诊断装置,特别涉及一种减少非必要的信号分量的技术。
背景技术
为了改善超声波图像的图像质量,减少比如包含在接收的信号中的旁瓣(side lobe)和栅瓣(grating lobe)的非必要的信号分量是可取的。在超声波的接收波束的形成中,从多个换能器元件(transducer element)处获得的多个接收的波信号根据焦点位置被进行延迟处理,从而匹配到来自焦点位置的反射波的相位。换言之,对多个接收的波信号进行整相。施加延迟处理的多个接收的波信号进行求和处理,以形成接收波束信号。
然而,从每个换能器元件处获得的接收的波信号还包含非必要的反射波分量,这些非必要的反射波分量来自焦点位置之外的位置。因此,由对施加了延迟处理的多个接收的波信号进行求和处理而获得的接收波束信号,由于非必要的反射波分量,而包含非必要的信号分量。在相关技术中提出了减少非必要的信号分量的技术。
例如,专利文献1公开了使用表示整相(phasing)的程度的CF(相干因子)来减少非必要的信号分量的技术。另外,作为表示整相的程度的值,非专利文献1公开了PCF(相位相干因子)和SCF(符号相干因子),专利文献2公开了STF(符号渡越因子,Sign Transit Factor),以及非专利文献2公开了GCF(广义相干因子)。
相关技术参考:
专利文献:
专利文献1:USP 5,910,115
专利文献2:JP 2012-81114A
非专利文献:
非专利文献1:J.Camecho,et al.,“Phase Coherence Imaging”,IEEEtrans.UFFC,vol.56,No.5,2009
非专利文献2:Pai-Chi Li,et al.,“Adaptive Imaging Using theGeneralized Coherence Factor”,IEEE trans.UFFC,vol.50,No.2,2003
发明内容
技术问题
表示整相的程度的值,例如CF,在减少与其中的多个换能器元件以一维方式排列的一维阵列换能器相关的非必要的信号分量中特别地优选。
在这种情况下,本发明的发明人进行了关于减少在其中的多个换能器元件以二维方式排列的二维阵列换能器中的非必要的信号分量的研发。
本发明是在研发过程中构思的,并且本发明的优势在于提供一种减少在二维阵列换能器中的非必要的信号分量的技术。
技术问题的解决方案
根据本发明的一个方面,提供了一种超声波诊断装置,包括:多个换能器元件,其以二维方式排列;延迟处理器,其对从所述多个换能器元件获得的多个接收的波信号施加延迟处理,以对接收的波信号进行整相;评估值计算单元,其基于施加了延迟处理的所述多个接收的波信号来评估整相的程度,以获得与以二维方式排列的所述多个换能器元件相关的二维评估值;求和处理器(summation processor),其对施加了延迟处理的多个接收的波信号施加求和处理,以获得接收信号;以及信号调整单元,其基于所述二维评估值来调整所述接收信号,以减少非必要的信号分量。根据本发明的另一个方面,优选地,所述评估值计算单元,对于互不相同的多个排列方向中的每个排列方向,针对以二维方式排列的多个换能器元件,评估整相的程度,以计算一维评估值,以及基于从多个排列方向获得的多个所述一维评估值,来计算表示二维排列整体的整相的程度的二维评估值。
在上述的装置中,所述多个换能器元件以晶格形状排列;例如,沿着横向(x方向)和纵向(y方向),并形成二维阵列换能器。多个换能器元件形成的换能器平面可以具有矩形形状或圆形形状。另外,相邻的换能器元件可以互相偏移地排列,或者换能器元件中的一些换能器元件可以被设置为无效的元件,以形成二维的稀疏型阵列换能器。
所述评估值计算单元对于多个换能器元件的每个排列方向评估整相的程度,以计算一维评估值。整相的程度是指,与施加了延迟处理的多个接收的波信号的相位的状态有关的程度,以及,例如,整相的程度的评估包括评估相位的匹配程度或者相位的偏移程度。因为评估值计算单元以一维的方式计算整相的程度是足够的,例如专利文献1的CF(相干因子)、非专利文献1的PCF(相位相干因子)和SCF(符号相干因子)、以及专利文献2的STF(符号渡越因子)等,可以作为一维评估值来使用。可替换地,除了这些值以外的一维评估值也可以使用。
评估值计算单元,基于从多个排列方向获得的多个一维评估值,来计算表示二维排列的整体的整相程度的二维评估值。例如,可以通过相对简单的计算,比如对多个一维评估值的求和、相乘、平方和等,来计算二维评估值。
基于二维评估值来调整通过将施加了延迟处理的多个接收的波信号进行求和处理所得到的接收信号。例如,信号调整单元,根据二维评估值的大小来改变接收信号的增益(振幅)。可替换地,必要时也可以基于二维评估值来调整接收信号的相位等。通过信号调整单元基于二维评估值来调整接收信号,从而减少了非必要的信号分量。
根据本发明的另一个方面,优选地,所述评估值计算单元,对于每个排列方向,基于从沿着该排列方向的至少一行换能器元件获得的多个接收的波信号来计算该排列方向的一维评估值。
根据本发明的另一个方面,优选地,所述评估值计算单元对于相互正交的两个排列方向中的每个排列方向计算所述一维评估值,以及基于从所述两个排列方向获得的两个一维评估值来计算所述二维评估值。
本发明的有益效果
根据本发明的各个方面,提供了一种在二维阵列换能器中减少非必要的信号分量的技术。例如,根据本发明的优选的方式,基于从多个排列方向获得的多个一维评估值来计算表示多个换能器元件的二维排列的整体的整相程度的二维评估值,并且基于二维评估值来调整接收信号,从而减少非必要的信号分量。
附图说明
图1为示出了根据本发明的优选的实施例的超声波诊断装置的整体结构的示意图。
图2为示出了以二维方式排列的多个换能器元件12的具体示例的示意图。
图3为示出了参照二维阵列换能器10的坐标系的示意图。
图4为示出了延迟处理之后的接收的波信号的波面的示意图。
具体实施方式
图1为示出了根据本发明的优选的实施例的超声波诊断装置的整体结构的框图。多个换能器元件12中的每个均为发射和接收超声波的元件,多个换能器元件12以二维方式排列,以形成二维阵列换能器10。二维阵列换能器10是在三维诊断区以三维方式扫描超声波波束的三维图像的超声波探测器。二维阵列换能器10可以以三维方式,电性地或者以电性扫描和机械扫描结合地,扫描超声波波束。
发射单元20输出发射信号到二维阵列换能器10的多个换能器元件12中的每个,对多个换能器元件12进行发射控制以形成发射波束,并在诊断区扫描发射波束。换言之,发射单元20具有发射波束形成器的功能。
每个换能器元件12由发射单元20进行发射控制以发射超声波,并接收来自诊断区的响应发射的波而获得的超声波。由每个换能器元件12接收超声波所获得的接收的波信号,从每个换能器元件12通过前置放大器14发送到延迟处理器30。
延迟处理器30包括多个延迟电路32。每个延迟电路32对从对应的换能器元件12通过对应的前置放大器14获得的接收的波信号施加延迟处理。在这种构造下,根据焦点位置对从多个换能器元件12处获得的多个接收的波信号进行延迟处理,并对来自焦点位置的反射波的相位进行匹配。换言之,对多个接收的波信号进行整相。施加了延迟处理的多个接收的波信号在求和处理器40中进行求和处理,以形成接收波束信号。
在这种方式中,整相求和处理由延迟处理器30和求和处理器40来执行,以实现接收波束形成器的功能。在诊断区的整个区域内扫描接收信号,以跟随发射波束的扫描,并沿着接收波束来收集接收波束信号。
评估值计算单元50基于从延迟处理器30获得的多个接收的波信号来评估整相的程度,从而计算与多个换能器元件12的二维排列相关的二维评估值。乘法单元60将从求和处理器40输出的接收波束信号和二维评估值相乘,来调整接收波束信号的增益。通过该处理,减少非必要的信号分量。换言之,乘法单元60作为一个基于二维评估值来调整接收信号的信号调整单元起作用。
图像形成单元70基于被乘法单元60调整过的接收波束信号来形成三维超声波图像。图像形成单元70使用如体绘制方法(volume renderingmethod),来形成以三维方式表征诊断目标的超声波图像。可替换地,图像形成单元70可以形成诊断目标的层析成像的图像(B模式图像)或多普勒图像。图像形成单元70所形成的超声波图像被显示在液晶显示器或类似物实现的显示器72上。控制器80综合地控制图1的超声波诊断装置的整体。
在如图1所示的结构(功能块)中,前置放大器14、发射单元20、延迟处理器30、求和处理器40、评估值计算单元50、乘法单元60、和图像形成单元70中的每一个均可以由硬件实现,例如,电子/电气电路、处理器或类似物,以及例如存储器的设备可以在实现这些结构时根据需要来使用。另外,控制器80可以通过例如硬件(例如,CPU、处理器、和存储器)和定义了CPU和处理器的运行的软件(程序)之间的配合来实现。
已经概述了图1的超声波诊断装置。接下来,将会详细描述图1的超声波诊断装置中使用的二维评估值。图1所示的结构(部件)在图1的附图标记和下文的说明书中被提及。
图2为示出了以二维方式排列的多个换能器元件12的具体示例的示意图。图2示出了二维阵列换能器10的多个换能器元件12的一种排列状态。在图2的具体示例中,多个换能器元件12沿着相互正交的x方向和y方向中的每个方向排列,以晶格形状排列,从而形成二维阵列换能器10的换能器平面。基于x方向和y方向的一维评估值来计算二维评估值。
图3为示出了参照二维阵列换能器10的坐标系的示意图。图3示出了以二维阵列换能器10的换能器平面的中心作为原点的xyz直角坐标系和极坐标系。
超声波从接收焦点到换能器元件12(x,y)的传播距离通过等式1来计算。在等式1中,c表示超声波的速度,tf表示超声波的传播时间。
[等式1]
当x和y的绝对值均足够地小于r时,等式1可以近似为等式2。
[等式2]
具有与从接收焦点到换能器元件12(x,y)相同距离的从原点到点P(r,α,β)的超声波的传播距离,通过等式3来表示,其与等式2具有相似的推导。
[等式3]
c · t p ≅ r + x 2 + y 2 2 r - s i n α · ( x c o s β + y s i n β )
因此,在换能器元件12(x,y)处,来自接收焦点的反射波(接收的波信号)与来自点P(r,α,β)的反射波(接收的波信号)之间的相位差,用等式4来表示,其中,超声波的频率为f。
[等式4]
在二维阵列换能器10中,当沿着x方向排列的换能器元件12的元件数和沿着y方向排列的换能器元件12的元件数均为N时,相邻的换能器元件12之间的元件间隔为λ/2,当沿着x方向排列的换能器元件的元件编号为m、且沿着y方向排列的换能器元件12的元件编号为n时,可以得到等式5。
[等式5]
x = ( m - N - 1 2 ) · λ 2 , y = ( n - N - 1 2 ) · λ 2 , 0 ≤ m ≤ N - 1 , 0 ≤ n ≤ N - 1 ,
将等式5代入等式4,以得到等式6。
[等式6]
进一步地,等式6可以简化,使得元件编号为(m,n)的换能器元件处的来自接收焦点的反射波(接收的波信号)与来自点P(r,α,β)的反射波(接收的波信号)之间的相位差,可以用等式7来表示。
[等式7]
Δψ(m,n)=Am+Bn+D
D = - N - 1 2 ( A + B )
等式7为(m,n,ΔΨ)坐标系的平面等式。在该平面中,在x=y=0时,ΔΨ等于0(ΔΨ=0)。因此,通过分别估计x方向和y方向中的每个方向的相位改变,可以确定由等式7表示的平面。
当接收焦点位于的情况下执行延迟处理时,因为在延迟处理后来自点P(r,α,β)的反射波(接收的波信号)会被等式7所示的相位差移相,反射波用等式8表示。等式8中的G表示复振幅。
[等式8]
S(m,n)=G·exp(-j(Am+Bn))
使用等式8的结果,来检查二维评估值的计算。二维评估值,是表示以二维方式排列的多个换能器元件12的整体的接收波信号的整相程度的值,并且其基于在x方向和y方向上的一维评估值来计算出。作为一维评估值,可以使用CF(相干因子)、PCF(相位相干因子)、SCF(符号相干因子)、STF(符号渡越因子)、或其他诸如此类的值。
<CF(相干因子)>
当以一维方式排列的第i个换能器元件12的接收的波信号为s(i)时,一维CF通过等式9来计算。
[等式9]
CF 1 D = | &Sigma; i = 0 N - 1 s ( i ) | 2 N &CenterDot; &Sigma; i = 0 N - 1 | s ( i ) | 2
当以二维方式排列的第(m,n)个换能器元件12处的接收的波信号为s(m,n)时,将等式9扩展成二维,可得到表示二维CF的等式10。
[等式10]
CF 2 D = | &Sigma; m = 0 N - 1 &Sigma; n = 0 N - 1 s ( m , n ) | 2 N 2 &CenterDot; &Sigma; m = 0 N - 1 &Sigma; n = 0 N - 1 | s ( m , n ) | 2
当将等式8代入等式10时,消去复振幅G,并且由于m和n是相互独立的值,因此得到等式11。
[等式11]
CF 2 D = | &Sigma; m = 0 N - 1 &Sigma; n = 0 N - 1 exp ( - j ( A m + B n ) ) | 2 N 2 &CenterDot; | &Sigma; m = 0 N - 1 &Sigma; n = 0 N - 1 exp ( - j ( A m + B n ) ) | 2 = | &Sigma; m = 0 N - 1 exp ( - j A m ) | 2 &CenterDot; | &Sigma; n = 0 N - 1 exp ( - j B n ) | 2 N 2 &CenterDot; &Sigma; m = 0 N - 1 | exp ( - j A m ) | 2 &CenterDot; &Sigma; n = 0 N - 1 | exp ( - j B n ) | 2 = CF 1 D x &CenterDot; CF 1 D y
评估值计算单元50,基于从延迟处理器30输出的施加了延迟处理(参考等式8)的多个接收的波信号,根据等式11来应用计算,以得到xy平面中的二维评估值CF2D,其基于在x方向上的一维评估值CF1Dx和在y方向上的一维评估值CF1Dy
<PCF(相位相干因子)>
使用换能器元件12的排列方向上的接收的波信号的相位的标准偏差σ(ΔΨ(i)),通过等式12来计算一维PCF。
[等式12]
PCF 1 D = 1 - &gamma; &sigma; 0 &sigma; ( &Delta; &psi; ( i ) )
&sigma; ( &Delta; &psi; ( i ) ) = 1 N &CenterDot; &Sigma; i = 0 N - 1 ( &Delta; &psi; ( i ) ) 2 - { 1 N &CenterDot; &Sigma; i = 0 N - 1 &Delta; &psi; ( i ) } 2 , i = 0 , 1 , ... , N - 1
其中,σ0=π/31/2为归一化标准偏差的常数,γ为调整参数。
当换能器元件12为二维方式排列时的标准偏差σ(ΔΨ(m,n)),基于等式7的结果,用等式13来表示。
[等式13]
σ(Δψ(m,n))=σ(Am+Bn+D)
等式13可以利用等式14示出的方差的特性,变换为等式15。
[等式14]
σ2(Am+Bn+D)=σ2(Am)+σ2(Bn)
[等式15]
&sigma; ( &Delta; &psi; ( m , n ) ) = &sigma; 2 ( A m ) + &sigma; 2 ( B n )
因此,当表示一维PCF的等式12扩展成二维时,就得到表示二维PCF的等式16。在等式16中,为了将扩展成二维的标准偏差归一化,在平方根中引入了因子1/2。
[等式16]
PCF 2 D = 1 - &gamma; &sigma; 0 &sigma; x 2 + &sigma; y 2 2 &sigma; x = &sigma; ( A m ) , &sigma; y = &sigma; ( B n )
利用在x方向上的一维PCF1Dx和在y方向上的一维PCF1Dy,可以从等式16得到等式17。
[等式17]
PCF 2 D = 1 - &gamma; &sigma; 0 &sigma; x 2 + &sigma; y 2 2 = 1 - 1 2 &CenterDot; { ( &gamma; &sigma; 0 &sigma; x ) 2 + ( &gamma; &sigma; 0 &sigma; y ) 2 } = 1 - ( 1 - PCF 1 D x ) 2 + ( 1 - PCF 1 D y ) 2 2
评估值计算单元50,基于从延迟处理器30输出的施加了延迟处理的多个接收的波信号,根据等式17来应用计算,以得到在xy平面中的二维评估值PCF2D,其基于在x方向上的一维评估值PCF1Dx和在y方向上的一维评估值PCF1Dy
<SCF(符号相干因子)>
SCF与PCF基于相同的原则,接收的波信号被二进制化,以使得二进制化的信号b(i)的标准偏差用作索引(index),不用计算相位。在这种情况下,由于二进制化的信号b(i)的标准偏差的值为0~1,用于PCF的等式12中的γ和σ0可以省略,一维SCF由等式18表示。
[等式18]
SCF1D=1-σ(b(i))
等式16和17可以类似地计算,从而得到如等式19的二维SCF。
[等式19]
SCF 2 D = 1 - &sigma; x 2 + &sigma; y 2 2 = 1 - ( 1 - SCF 1 D x ) 2 + ( 1 - SCF 1 D y ) 2 2
评估值计算单元50,通过将从延迟处理器30输出的施加了延迟处理的多个接收的波信号二进制化,根据等式19来应用计算,以得到xy平面中的二维评估值SCF2D,其基于在x方向上的一维评估值SCF1Dx和在y方向是的一维评估值SCF1Dy。可替换地,评估值计算单元50可以利用调整系数p来得到的SCF2D的p次幂,从而得到二维评估值。
<STF(符号渡越因子)>
一维STF通过等式20基于换能器元件12的排列方向上的接收的波信号的过零密度(对应于平均频率)来计算。
[等式20]
STF 1 D = 1 N - 1 &Sigma; i = 0 N - 2 c ( i )
c ( i ) = 1 i f s i g n ( s ( i ) ) &NotEqual; s i g n ( s ( i + 1 ) ) 0 i f s i g n ( s ( i ) ) = s i g n ( s ( i + 1 ) )
在二维阵列换能器10的情况下,如等式7所示,在元件编号为(m,n)的换能器元件12处,来自接收焦点F的反射波(接收的波信号)与来自点P(r,α,β)的反射波(接收的波信号)之间的相位差,是平面等式。此外,来自点P(r,α,β)的经过延迟处理之后的接收的波信号,如等式8所示。
图4为示出施加了延迟处理的接收的波信号的波平面的示意图。在图4中,xy坐标系对应于二维阵列换能器10的换能器平面,xy坐标系中的实线示出了相位为2nπ[rad](其中n为整数)的部分。此外,xy坐标系中的虚线示出了相位为(2n+1)π[rad](其中n为整数)的部分。
作为表示整相程度的指标,在图4所示的波平面中,使用了在最大频率的方向上的平均频率f2D。平均频率f2D可以基于x方向上的频率fx和y方向上的频率fy,通过等式21来计算。
[等式21]
f 2 D = f x 2 + f y 2 f x &CenterDot; f y
由于过零密度对应于平均频率,所以,如果基于等式21将STF扩展成二维,则得到如等式22所示的二维STF。
[等式22]
STF 2 D = ( STF 1 D x ) 2 + ( STF 1 D y ) 2 STF 1 D x &CenterDot; STF 1 D y
评估值计算单元50,基于从延迟处理器30输出的施加了延迟处理的多个接收的波信号,根据等式22来应用计算,以得到xy平面中的二维评估值STF2D,其基于在x方向上的一维评估值STF1Dx和在y方向上的一维评估值STF1Dy。由于越高的过零密度表示信号距离主瓣越远(即,该信号是应当被减少的信号),例如,可以通过利用调整系数p来得到的(1-STF2D)的p次幂,从而得到二维评估值。
利用CF、PCF、SCF和STF计算二维评估值如下。由于等式7示出的相位差是平面等式,因此当在x方向上计算一维评估值时,例如,理论上,在沿着x方向的任何一行都能得到相同的结果。因此,例如,图2中沿着x方向排列的一行可以设为代表行,x方向上的一维评估值可以基于从包含在该代表行中的多个换能器元件12获得的施加了延迟处理的多个接收的波信号来计算。作为代表行,例如,穿过换能器平面的中心或靠近换能器平面的中心的行是优选的。类似地,当为y方向计算一维评估值时,沿着y方向排列的一行设为代表行。可替换地,可以将多个行设定为代表行,可以将该多个行的均值设定为一维评估值。
由于在二维阵列换能器10中的换能器元件12的元件数目大,一些情况下会在探测器中进行信道缩减(channel reduction)。此时,评估值计算单元50可以基于信道缩减之后的接收的波信号来计算一维评估值,并且可以基于一维评估值来得到二维评估值。
已经描述了本发明的优选的实施例。然而,在上述实施例中,仅在各个方面进行了示例,并不限制本发明的范围。本发明在发明的范围和主旨内,包括各种变形的构造。
[附图标记说明]
10二维阵列换能器;12换能器元件;20发射单元;30延迟处理器;32延迟电路;40求和处理器;50评估值计算单元;60乘法单元;70图像形成单元;72显示器;80控制器。

Claims (12)

1.一种超声波诊断装置,包括:
多个换能器元件,其以二维方式排列;
延迟处理器,其对从所述多个换能器元件获得的多个接收的波信号施加延迟处理,以对接收的波信号进行整相;
评估值计算单元,其基于施加了延迟处理的多个接收的波信号来评估整相的程度,以获得与以二维方式排列的所述多个换能器元件相关的二维评估值;
求和处理器,其对施加了延迟处理的多个接收的波信号施加求和处理,以获得接收信号;以及
信号调整单元,其基于所述二维评估值来调整所述接收信号,以减少非必要的信号分量。
2.根据权利要求1所述的超声波诊断装置,其中:
所述评估值计算单元基于施加了所述延迟处理的多个接收的波信号,来计算多个一维评估值,并基于所述多个一维评估值来计算所述二维评估值。
3.根据权利要求1所述的超声波诊断装置,其中:
所述评估值计算单元,对于多个互不相同的排列方向中的每个排列方向,针对以二维方式排列的所述多个换能器元件来评估整相的程度,以计算一维评估值,以及基于从多个排列方向获得的多个所述一维评估值,来计算表示二维排列整体的整相的程度的所述二维评估值。
4.根据权利要求3所述的超声波诊断装置,其中:
所述评估值计算单元,对于每个排列方向,基于从沿着该排列方向的至少一行换能器元件获得的多个接收的波信号来计算该排列方向的一维评估值。
5.根据权利要求3所述的超声波诊断装置,其中:
所述评估值计算单元对于互不相同的两个排列方向中的每个排列方向计算一维评估值,以及基于从所述两个排列方向获得的两个一维评估值来计算所述二维评估值。
6.根据权利要求5所述的超声波诊断装置,其中:
所述评估值计算单元,对于每个排列方向,基于从沿着该排列方向的至少一行换能器元件获得的多个接收的波信号来计算该排列方向的一维评估值。
7.根据权利要求3所述的超声波诊断装置,其中:
所述评估值计算单元对于相互正交的两个排列方向中的每个排列方向计算一维评估值,以及基于从所述两个排列方向获得的两个一维评估值来计算所述二维评估值。
8.根据权利要求7所述的超声波诊断装置,其中:
所述评估值计算单元,对于每个排列方向,基于从沿着该排列方向的至少一行换能器元件获得的多个接收的波信号来计算该排列方向的一维评估值。
9.根据权利要求3所述的超声波诊断装置,其中:
所述评估值计算单元计算作为所述一维评估值的CF(相干因子)、PCF(相位相干因子)、SCF(符号相干因子)、或STF(符号渡越因子)。
10.根据权利要求5所述的超声波诊断装置,其中:
所述评估值计算单元计算作为所述一维评估值的CF(相干因子)、PCF(相位相干因子)、SCF(符号相干因子)、或STF(符号渡越因子)。
11.根据权利要求7所述的超声波诊断装置,其中:
所述评估值计算单元计算作为所述一维评估值的CF(相干因子)、PCF(相位相干因子)、SCF(符号相干因子)、或STF(符号渡越因子)。
12.根据权利要求1所述的超声波诊断装置,其中:
所述信号调整单元将从所述求和处理器获得的施加了求和处理的接收信号与从所述评估值计算单元获得的所述二维评估值相乘,来调整施加了求和处理的接收信号,以及减少所述非必要的信号分量。
CN201380070054.9A 2013-01-09 2013-12-25 超声波诊断装置 Expired - Fee Related CN104902824B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2013001824A JP5504357B1 (ja) 2013-01-09 2013-01-09 超音波診断装置
JP2013-001824 2013-01-09
PCT/JP2013/084685 WO2014109234A1 (ja) 2013-01-09 2013-12-25 超音波診断装置

Publications (2)

Publication Number Publication Date
CN104902824A true CN104902824A (zh) 2015-09-09
CN104902824B CN104902824B (zh) 2017-05-31

Family

ID=50941835

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201380070054.9A Expired - Fee Related CN104902824B (zh) 2013-01-09 2013-12-25 超声波诊断装置

Country Status (5)

Country Link
US (1) US20150342570A1 (zh)
EP (1) EP2944265A4 (zh)
JP (1) JP5504357B1 (zh)
CN (1) CN104902824B (zh)
WO (1) WO2014109234A1 (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110599405A (zh) * 2018-06-13 2019-12-20 青岛海信医疗设备股份有限公司 一种超声图像增强方法、装置和计算机设备
WO2020155348A1 (zh) * 2019-02-01 2020-08-06 深圳先进技术研究院 一种超声波信号调整方法、装置及超声阵列

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117503203B (zh) * 2024-01-03 2024-03-22 之江实验室 一种用于超声环阵成像的相位畸变校正方法和系统

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0765146A (ja) * 1993-08-24 1995-03-10 Toshiba Corp 超音波診断装置
CN1575772A (zh) * 2003-06-25 2005-02-09 阿洛卡株式会社 超声波诊断装置
JP2010119481A (ja) * 2008-11-18 2010-06-03 Aloka Co Ltd 超音波診断装置
CN102440807A (zh) * 2010-10-13 2012-05-09 日立阿洛卡医疗株式会社 超声波诊断装置
CN102599929A (zh) * 2011-01-25 2012-07-25 日立阿洛卡医疗株式会社 超声波诊断装置
US20120277589A1 (en) * 2011-04-28 2012-11-01 Konica Minolta Medical & Graphic, Inc. Ultrasound diagnostic device

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0256481B1 (de) * 1986-08-20 1994-01-12 Siemens Aktiengesellschaft Verfahren und Einrichtung zur adaptiven Fokussierung bei einem medizinischen Ultraschall-Bildgabegerät
US6111816A (en) * 1997-02-03 2000-08-29 Teratech Corporation Multi-dimensional beamforming device
US5910115A (en) 1997-09-22 1999-06-08 General Electric Company Method and apparatus for coherence filtering of ultrasound images
US6524254B2 (en) * 2001-06-20 2003-02-25 Bae Systems Information And Electronic Systems Integration, Inc. Orthogonally reconfigurable integrated matrix acoustical array
US7285094B2 (en) * 2002-01-30 2007-10-23 Nohara Timothy J 3D ultrasonic imaging apparatus and method
WO2005008280A1 (en) * 2003-07-17 2005-01-27 Angelsen Bjoern A J Corrections for wavefront aberrations in ultrasound imaging
JPWO2006006460A1 (ja) * 2004-07-08 2008-04-24 株式会社日立メディコ 超音波撮像装置
JP5634818B2 (ja) 2010-10-13 2014-12-03 日立アロカメディカル株式会社 超音波診断装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH0765146A (ja) * 1993-08-24 1995-03-10 Toshiba Corp 超音波診断装置
CN1575772A (zh) * 2003-06-25 2005-02-09 阿洛卡株式会社 超声波诊断装置
JP2010119481A (ja) * 2008-11-18 2010-06-03 Aloka Co Ltd 超音波診断装置
CN102440807A (zh) * 2010-10-13 2012-05-09 日立阿洛卡医疗株式会社 超声波诊断装置
CN102599929A (zh) * 2011-01-25 2012-07-25 日立阿洛卡医疗株式会社 超声波诊断装置
US20120277589A1 (en) * 2011-04-28 2012-11-01 Konica Minolta Medical & Graphic, Inc. Ultrasound diagnostic device

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110599405A (zh) * 2018-06-13 2019-12-20 青岛海信医疗设备股份有限公司 一种超声图像增强方法、装置和计算机设备
CN110599405B (zh) * 2018-06-13 2023-06-16 青岛海信医疗设备股份有限公司 一种超声图像增强方法、装置和计算机设备
WO2020155348A1 (zh) * 2019-02-01 2020-08-06 深圳先进技术研究院 一种超声波信号调整方法、装置及超声阵列

Also Published As

Publication number Publication date
WO2014109234A1 (ja) 2014-07-17
JP2014132959A (ja) 2014-07-24
US20150342570A1 (en) 2015-12-03
CN104902824B (zh) 2017-05-31
JP5504357B1 (ja) 2014-05-28
EP2944265A4 (en) 2016-09-28
EP2944265A1 (en) 2015-11-18

Similar Documents

Publication Publication Date Title
Flesch et al. 4D in vivo ultrafast ultrasound imaging using a row-column addressed matrix and coherently-compounded orthogonal plane waves
US8793079B2 (en) Method for imaging of nonlinear interaction scattering
Peralta et al. Coherent multi-transducer ultrasound imaging
JP2016042036A (ja) ビームフォーミング方法、計測イメージング装置、及び、通信装置
US11452498B2 (en) Ultrasound imaging
CN104414687A (zh) 超声波测定装置、超声波图像装置及超声波测定方法
US20140245836A1 (en) Object information acquiring apparatus
Varslot et al. Forward propagation of acoustic pressure pulses in 3D soft biological tissue
EP2766738B1 (en) Three dimensional (3d) transverse oscillation vector velocity ultrasound imaging
CN112513675A (zh) 超声方法和装置
Peralta et al. Impact of aperture, depth, and acoustic clutter on the performance of coherent multi-transducer ultrasound imaging
CN104902824A (zh) 超声波诊断装置
CN101677808B (zh) 超声波诊断装置、超声波诊断装置的数据处理装置以及信号处理方法
Peralta et al. Coherent multi-transducer ultrasound imaging in the presence of aberration
Yen Beamforming of sound from two-dimensional arrays using spatial matched filters
CN111580112B (zh) 基于平面波的水下声呐传感器阵列成像方法
CN103622721A (zh) 被检体信息获取装置和方法与信息处理装置
CN102188261B (zh) 基于非衍射阵列波的合成孔径三维超声成像方法
Hazard et al. Effects of motion on a synthetic aperture beamformer for real-time 3D ultrasound
Han et al. 3D ultrasound imaging in frequency domain with 1D array transducer
US11883240B2 (en) Sparse convolutional beamforming for ultrasound imaging
Sampson et al. FPGA implementation of low-power 3D ultrasound beamformer
Li et al. Fast computation of far-field pulse-echo PSF of arbitrary arrays for large sparse 2-D ultrasound array design
Chen et al. Low-complexity motion-compensated beamforming algorithm and architecture for synthetic transmit aperture in ultrasound imaging
Diarra et al. Comparison of different optimized irregular sparse 2D ultrasound arrays

Legal Events

Date Code Title Description
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C41 Transfer of patent application or patent right or utility model
TA01 Transfer of patent application right

Effective date of registration: 20161111

Address after: Tokyo, Japan, Japan

Applicant after: Hitachi Ltd.

Address before: Tokyo, Japan

Applicant before: Hitachi Aloka Medical Ltd.

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

Granted publication date: 20170531

Termination date: 20181225

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