CN102770079B - 超声波成像装置和控制延迟的方法 - Google Patents
超声波成像装置和控制延迟的方法 Download PDFInfo
- Publication number
- CN102770079B CN102770079B CN201180010435.9A CN201180010435A CN102770079B CN 102770079 B CN102770079 B CN 102770079B CN 201180010435 A CN201180010435 A CN 201180010435A CN 102770079 B CN102770079 B CN 102770079B
- Authority
- CN
- China
- Prior art keywords
- propagation time
- time value
- depth
- value
- address
- 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
Links
Classifications
-
- G—PHYSICS
- G10—MUSICAL INSTRUMENTS; ACOUSTICS
- G10K—SOUND-PRODUCING DEVICES; METHODS OR DEVICES FOR PROTECTING AGAINST, OR FOR DAMPING, NOISE OR OTHER ACOUSTIC WAVES IN GENERAL; ACOUSTICS NOT OTHERWISE PROVIDED FOR
- G10K11/00—Methods 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/18—Methods or devices for transmitting, conducting or directing sound
- G10K11/26—Sound-focusing or directing, e.g. scanning
- G10K11/34—Sound-focusing or directing, e.g. scanning using electrical steering of transducer arrays, e.g. beam steering
- G10K11/341—Circuits therefor
- G10K11/346—Circuits therefor using phase variation
-
- 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/52034—Data rate converters
-
- 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/52046—Techniques for image enhancement involving transmitter or receiver
- G01S7/52049—Techniques for image enhancement involving transmitter or receiver using correction of medium-induced phase aberration
Abstract
本发明使得能够通过在深度方向上使用递归关系来针对每个接收通道并行计算在校正折射之后的超声波传播时间值。此外,可以通过使用在基准深度处预先获得的精确传播时间值来在每次到达基准深度就校正传播时间值而避免误差的累积。对于这种误差校正,用于计算传播时间值的递归关系可以是近似表达式。例如,可以使用基准深度之间的基准传播时间值的倾度来计算传播时间值。在实际电路中,接收信号被依次存储在存储器中,通过计算与超声波的传播时间值对应的地址位置并且把计算出的地址中存储的接收信号相加来形成接收波束。
Description
技术领域
本发明涉及一种用于接收在被检体中反射或生成的超声波以及对所接收的超声波进行成像的超声波成像装置,更具体地说,涉及一种对从用于形成接收波束的多个接收元件获得的信号根据超声波传播时间而执行延迟校正的数字超声波成像装置。
背景技术
由于超声波不产生暴露并且基本上无创,因此超声波图像诊断被广泛使用。尤其近年来,因为超声波成像装置的数字化的进步,所以可以容易地以高精度来控制用于形成超声波束而生成的延迟时间。据此,可以通过使用用于每个区域的诊断的紧凑装置来获得高质量生物断层图像。除了将超声波辐射到被检体、接收回声并且生成其图像的超声回声诊断装置之外,还已提出了辐射脉冲光并且接收被检体内部生成的光声波(超声波)并且产生其图像的光声成像装置。
在传统超声波成像装置中,通常基于超声波传播通过声音速度为恒定的介质的假设来执行用于形成超声波束的延迟时间控制。另一方面,如果声音速度与生物组织(被检体)不同的材料(比如用于紧固在超声探针的表面上的声学透镜层和生物组织的固定板)存在于生物组织与发送/接收超声波的换能器之间,则超声波在具有不同声音速度的区域的边界处折射。结果,在传统延迟时间控制的情况下,超声波束的焦点未被良好地形成,并且图像质量恶化。当厚脂肪层同样存在于被检体表面上时,由于脂肪层和脂肪层之下的生物组织的声音速度不同,所以超声波以相同方式折射,结果,在脂肪层之下的组织并未清楚地成像。
可以通过使用Snell定律考虑折射路径来解析地确定当具有不同声音速度的层存在于被检体与换能器之间时的延迟量。已经提出了使用Snell定律校正延迟时间的一些方法。
美国专利No.6,607,489(PTL 1)公开了一种具有用于紧固生物组织的固定板的超声装置,其中,针对每个换能器和每个焦点来把归因于在固定板与目标组织之间的声音速度的差导致的基于精确计算的延迟时间预先存储在装置中,并且使用该延迟时间来控制发送/接收波束的延迟。还公开了实时地计算延迟时间,但未公开具体数值计算算法。
WO 01-026555(PTL 2)公开了一种校正关于由组织表面上的脂肪层和声学透镜产生的折射的影响的超声波成像装置,以使得抑制关于在被检体的脂肪层之下的组织区域的图像的质量的恶化。由于对相对厚的并且其厚度取决于目标被检体而不同的脂肪层产生的折射的影响,所以通过使用基于Snell定律的延迟时间校正近似表达式的计算来确定每个换能器中的延迟时间。在WO 01-026555的方法中,基于邻近换能器中的延迟时间使用递归关系来计算下一换能器中的延迟时间。由此减少计算时间。
(PTL 1)美国专利No.6,607,489
(PTL 2)国际公开No.WO 01-026555
发明内容
然而,在美国专利No.6,607,489的情况下,必须预先对于所有换能器存储对于所有焦点的延迟时间。因此,需要用于存储延迟时间值的大容量存储器,以便在接收期间具有大量动态聚焦的焦点,以实现高图像质量。通常,分辨率随着换能器元件的数量增加以及动态接收聚焦的焦点的数量增加而改进。还期望捕获被检体深处的图像。在任一情况下,要存储的延迟时间值的数量增加,并且电路规模变得庞大。
此外,如果所***的固定板的厚度改变,则相应地需要延迟时间值的新集合。美国专利No.6,607,489还公开了一种包括对包含折射的超声波传播路径执行精确计算以确定延迟时间的计算单元的装置,并且可以根据所***的固定板的厚度来计算延迟时间值的集合,但未公开超声波传播路径计算的具体数值算法。因此,使用专用于数值运算的处理器(比如MPU或通用PC的CPU)来计算延迟时间,以用于高效地执行通用数值运算。然而,在接收期间的动态聚焦中,必须根据发送超声波脉冲的传播定时来高速地改变每个换能器中的延迟时间。因此,同样在使用用于执行控制的通用PC来计算延迟时间值并且将这些值读入用于形成超声波束的数字电路的情况下,在数字电路中需要用于存储延迟时间值的大容量存储器。另外,需要高速地将数据从PC传送到数字电路的部件。在两种配置中,电路规模增加,并且装置变得昂贵。专用于数值运算的MPU还导致用于形成超声波束的数字电路的规模和成本的增加。以此方式,根据美国专利No.6,607,489,传统数字电路的规模增加,这使得装置的大小和成本以及处理速度和对可以设置的分辨率的限制变得劣于现有技术。结果,传统超声诊断装置的优点,即图像观测的实时特性以及在由紧凑装置配置实现的宽范围中的灵活诊断,消失。
WO 01-026555中公开的传统技术的优点是:即使具有不同声音速度的介质变得较厚,也可以减少近似的延迟时间的误差;并且可以以相对较高速度来执行处理。然而,如果根据该方法而执行用于实时地改变焦点的动态接收聚焦,则以下问题存在。也就是说,因为通过递归关系依次执行的计算处理使用关于换能器位置的递归关系,所以必须对于所有换能器计算对于各个焦点的延迟量。因此,如果装置在少数焦点上具有很多换能器,则良好地执行实时处理。但在执行动态接收聚焦的情况下,必须对于每个焦点执行使用递归关系的计算,并且处理量增加,因此该方法不适合于实时处理。如果设置很多接收焦点以使得分辨率更高,则该问题尤其明显。
另一方面,在传统超声波成像装置的接收处理的情况下,通过对于每个通道几乎独立地并行处理(比如针对每个换能器的信号处理,特别地,对接收的信号的放大和数字化和延迟时间控制)来实现实时高速图像生成。然而,在根据WO 01-026555的处理的情况下,因为每次接收焦点改变就必须对于所有换能器执行用于重置延迟时间的递归关系处理,如上所述,所以无法有效地使用该并行处理运算。为了实时地执行其中设置很多接收焦点的整个处理,仍然需要用于分离地执行延迟计算的计算处理电路以及用于存储计算结果的大存储器容量。
因此,为了使用根据WO 01-026555的处理来执行动态接收聚焦,必须大大地改变传统数字电路,尤其是必须增加待使用的存储器和计算电路,或必须减少接收焦点的数量,以使得当改变焦点时用于重置延迟时间的时间减少。结果,难以通过增加用于接收聚焦的接收焦点的数量来在接收时改进分辨率。
基于光声成像方法的超声波成像装置同样具有源自上述问题的问题。特别地,在用于生成三维图像的光声成像方法的情况下,必须对于每个深度估计接收超声波的传播时间,恰如超声波脉冲回声方法的动态接收聚焦。优选的是,二维阵列用于换能器阵列(其中,布置很多换能器),以改进分辨率。因此,在基于传统光声成像方法的超声波图像诊断装置的情况下,用于重构三维图像的处理量变得巨大,并且难以实时地生成并且显示图像。因为使用通用PC来依次执行在获得数据之后的处理,所以与使用传统超声波脉冲回声方法的超声波图像诊断装置相比,实时诊断是困难的。
鉴于以上情况,本发明的目的在于提供一种超声波成像装置来在考虑超声波的传播路径的情况下高速地执行延迟量计算处理。
本发明在一个方面中提出一种超声波成像装置,包括:超声波接收机,具有用于接收被检体内部反射或生成的超声波的多个接收元件;接收波束形成器,根据从接收波束的焦点位置到接收元件的传播时间通过对与超声波接收机的接收元件中的至少一部分相对应的接收通道的接收信号执行延迟控制来形成接收波束信号;以及
图像处理器,使用接收波束信号来生成图像,其中,接收波束形成器通过反复进行用于把与焦点深度的改变相对应的传播时间值的改变与已经计算出的传播时间值相加以确定用于下一焦点的传播时间值的处理,来依次计算用于多个焦点的传播时间值,以及反复进行的处理能够对于接收通道中的每一个被单独执行,并且接收波束形成器对于接收通道中的至少一部分并行地计算传播时间值。
根据本发明,对于每个接收通道的深度方向通过递归关系计算与多个焦点对应的传播时间值,因此可以对于多个通道并行地计算传播时间值,并且可以高速地计算延迟时间值。
参照附图根据示例性实施例的以下描述,本发明的其它特征将变得清楚。
附图说明
图1A是描述根据示例1的超声波成像装置的示图。
图1B示出根据示例1的超声波成像装置中使用的超声波束。
图1C示出与图1B所示的超声波束对应的由超声波成像装置生成的图像数据。
图2是描述根据示例1的接收波束形成器的示图。
图3是描述根据示例1的读取地址计算单元的示图。
图4A是描述根据示例1的折射校正地址增量计算块的示图。
图4B是描述根据示例1的用于折射校正地址增量计算的流水线配置的示图。
图5是描述根据示例1的存储器读取控制器的示图。
图6是存储器地址计算处理的流程图。
图7是描述超声波传播时间的示图。
图8A是描述动态接收聚焦的示图。
图8B是描述在存储器中存储接收的信号的状态的示图。
图9A是描述用于根据示例2的光声成像方法的三维图像重构的体素数据结构的示图。
图9B是描述用于在示例2的光声成像方法中在特定体素中创建数据的构思的示图。
图9C是描述根据示例2的处理方法的示图。
图10A是描述根据示例2的基于光声成像方法的超声波图像诊断装置的配置的示图。
图10B是描述示例2中的图像重构的处理的示图。
具体实施方式
本发明是一种可以有效地应用于诊断装置的技术,其通过多个接收元件接收在被检体中反射或生成的超声波,将所接收的超声波转换为接收的信号,并且根据超声波传播时间对于每个接收元件执行延迟控制以生成接收波束信号。该诊断装置的示例是超声回声诊断装置,其允许超声波脉冲从换能器(超声发射机/接收机单元)进入被检体,接收从被检体内部的区域反射的超声波回声,并且生成其图像。另一示例是光声成像装置,其将脉冲光辐射到被检体,接收通过光声效应在内部生成的超声波,将其转换为接收的信号,并且使用该接收的信号来生成图像。在该说明书中,接收来自被检体内部的超声波并且生成图像的装置通常被称为“超声波成像装置”。
<超声波传播时间计算处理>
(推导用于计算超声波传播时间的近似递归关系)
现将描述根据本发明的考虑超声波的折射的情况下近似传播时间的计算方法的示例。图7示出从焦点21传播到换能器203的超声波的折射的状态。当x是换能器203的位置,D是聚焦长度之外的介质26的部分的厚度(即介质26的厚度),z是聚焦长度(焦点深度)之中的介质27的部分,θ是介质26中的折射角,并且θ'是介质27中的折射角时,建立以下关系表达式。如果z是负值,焦点21在介质26中,并且在此情况下,可以基于声音速度是均匀的假设来计算超声波传播时间,恰如现有技术的情况,所以将描述并非负值的z的情况。取决于z是负值还是非负值,可以通过与传统技术进行切换来驱动该方法。z在此以下被称为“深度”。
[数学式1]
Dtan+ztanθ′=x, ···式1(b)
其中,Vtis是介质27中的声音速度,Vm是介质26中的声音速度,T是在焦点21反射或生成的超声波到达换能器203的传播时间。式1(a)表示折射的Snell定律,式1(b)表示焦点21和换能器203的几何关系,式1(c)表示超声波传播路径的长度与声音速度的关系。
考虑处于比焦点21更深的位置中的焦点211。如果介质26的厚度不改变,则焦点的改变仅为dz,其为z的改变分量。因此,当焦点深度对于每个dz依次改变时,通过考虑dz的增量来相对于超声波传播时间Ti建立以下近似递归关系。
[数学式2]
Ti+1=Ti+dTi.···式2
如下表示当焦点深度改变dz时的超声波传播时间的改变dTi,其中,在与各个zi对应的式1中,ξi和ηi是cosθ'和cosθ(ξi=cosθ',ηi=cosθ)。
[数学式3]
ξi+1=ξi+dξi, ···式3(b)
ηi+1=ηi+dηi, ···式3(c)
zi+1=zi+dz. ···式3(d)
这里,使用ri如下表示dξi和dηi。
[数学式4]
dηi=γ2ridξi. ···式4(c)
γ是以下所示的量,并且如果介质26和介质27确定则在递归关系计算中是常数。以此方式,通过根据递归关系确定的参数ξi和ηi来确定超声波传播时间的增量dTi。在该说明书中,用于确定dTi并且由递归关系确定的参数被称为“递归参数”。
[数学式5]
通过使用上述关系表达式(近似递归关系),如果在深度zi处的ξi、ηi和dTi是已知的,则可以计算在深度改变dz的焦点处的ξi+1和ηi+1。使用该表达式,可以计算超声波传播时间的增量dTi,并且可以通过将该增量dTi与计算出的超声波传播时间Ti相加来近似地计算超声波传播时间Ti+1。
通过使用该递归关系在更深的焦点深度(dz>0)的方向上或在更浅的焦点深度(dz<0)的方向上依次进行计算,可以依次近似地计算关于具有不同深度的多个焦点的超声波传播时间值Ti+1的串。例如,使用哪个计算方向无关紧要,但如果该方法应用于动态接收聚焦,则焦点深度根据发送的脉冲的传播而依次增加的方向是优选的。
由于焦点的间隔dz的绝对值较小,所以近似的精度增加。特别地,用于对接收的信号(接收超声波的接收元件生成(转换)的信号)进行采样以生成数字信号的时钟频率的周期是数字接收的信号的时间分辨率的极限,因此优选的是选择与该周期对应的dz。在此情况下,焦点的最小间隔是用于所发送的超声波脉冲在时钟频率的周期期间在介质26中传播的距离,并且是|dz|=Vtis×Tclk,其中,Tclk是时钟频率的周期。
(初始值)
参数γ和焦点的间隔dz对于所有换能器都相同,如果给定每个换能器的位置x和初始值ξ0、η0、T0和z0,则通过以上递归关系来依次计算在焦点的每个间隔dz处的超声波传播时间Ti+1。
对于初始值ξ0、η0和T0,可以使用对于每个换能器预先计算的值。特别地,在较深焦点深度的方向上依次计算的情况下,初始深度z0可以设置为零,也就是说,初始焦点可以设置在介质26和介质27的边界上。此时,超声波仅在介质26中从焦点传播到每个换能器,从而可以对于每个换能器基于相同声音速度而以与传统超声波传播时间相同的方式来计算初始值T0、ξ0和η0,如下所示。
[数学式6]
如果介质26的厚度D不改变,则分离地计算的值可以是固定的,如果D改变,则计算的值可以提供为表。该处理可以被包括在初始化时的计算中。
(误差校正)
通常,如果使用近似递归关系来递归地执行计算,则归因于近似的误差累积,并且误差随着深度增加而增加。为了避免该情况,可以使用在预定基准深度zm处预先高度精确地计算的基准递归参数ξm和ηm以及基准传播时间Tm来校正误差的累积。通过以适当间隔来设置多个基准深度值,可以周期性地校正误差。
为了确定在基准深度zm处的ξm、ηm和Tm,首先求解以下四次方程式以获得h。
[数学式7]
使用在此获得的根中的满足0<h<D的最小实数根,基准值ξm、ηm和Tm可以表示如下。
[数学式8]
以此方式,通过数值计算预先确定在基准深度zm处的高精度值,一旦焦点深度到达基准深度zm,此后就使用ξm、ηm和Tm作为新的初始值来继续递归关系计算。上述基准值的计算,尤其是四次方程式的计算是复杂的,但由于仅对于相对少数量的基准深度值执行计算,因此不花费太多时间。
可以通过实际测量来确定基准传播时间。具体地说,每个换能器使用超声波装置接收来自实际上设置在基准深度zm处的点声源(例如水听器)的超声波脉冲。然后,实际上可以基于接收的脉冲的到达时间而测量基准传播时间Tm,以使得实际测量的基准传播时间Tm用于校正近似误差。在此情况下,通过逆向求解方程式8(c)确定h,并且该h代入方程式8(a)和方程式8(b)来确定基准递归参数ξm和ηm。
基准传播时间Tm是用于校正的量,因此基准传播时间值Tm的数量可以远小于在通过上述近似递归关系确定的在焦点的每个间隔dz处的超声波传播时间值Ti的数量。具体地说,当在接收时在动态聚焦中依次确定在每个接收焦点处的超声波传播时间Ti的值时,与对于所有接收焦点来存储超声波传播时间值相比,存储基准传播时间Tm和基准递归参数值ξm和ηm需要更少容量。
(变型)
在将初始深度z0设置为不为零的值的情况下,可以通过与确定基准值相同的方法来确定初始值。换句话说,将在方程式7中把zm看做z0对于h确定根,并且该h用于计算方程式8,或另外使用实际测量。
如果初始深度z0设置为不为零的值,则可以在更深的焦点深度的方向上或在更浅的焦点深度的方向上使用递归关系来进行计算。如果在更浅的焦点深度的方向上使用递归关系来进行计算,则当深度zi变为负值时,使用递归关系的计算停止。如果深度zi变为负值,则这说明焦点21在介质26中,在此情况下,计算应切换到当声音速度均匀时使用的传播时间的计算,其与传统技术相同。
在以上描述中,焦点的间隔dz在递归关系的每个步骤中是固定的,但如果焦点的间隔对于所有换能器是共同的,则可以在递归关系的每个步骤中改变。具体地说,随着深度变得更深,超声波传播时间的改变量dTi减少,近似精度增加,从而随着深度变得更深,可以通过增加dz的绝对值来省略更多的计算处理。
(该计算的优点)
上述递归关系描述为递归参数的有理多项式,并且其计算仅包括四个算术运算,不包括诸如超越函数和开平方的运算,因此计算处理可以容易地安装为数字电路。特别地,在不安装通用CPU或不使用通用PC的情况下,可以容易地封装该数字电路。可以使用MPU或GPU来封装该数字电路。
可以仅使用取决于每个换能器的参数x和初始值ξ0、η0和T0以及对于所有换能器共同的参数γ和焦点的间隔dz来依次计算上述递归关系。换句话说,无需在递归计算期间使用其它换能器的计算结果,因此每个换能器可以并行地执行该计算。因此,可以通过将该计算应用于由传统超声波图像诊断装置执行的用于每个通道的并行波束形成处理,来高速地执行图像生成执行处理。
在更深的焦点深度的方向上按顺序确定超声波传播时间的情况下,特别地,当在将接收焦点依次移动到与发送的超声波脉冲的传播定时对应的更深方向的同时调整接收波束的焦点时,该配置是合适的。
<动态接收聚焦>
现在,将描述使用超声波脉冲回声方法来执行动态接收聚焦的情况。在超声波脉冲回声方法中,首先发送超声波脉冲。所发送的超声波脉冲行进通过介质。同时,随着超声波脉冲行进而在每个位置处反射的超声波回声传播通过介质并且到达每个换能器,如上所述。在所发送的超声波脉冲从该超声波脉冲的发送到每个反射点的行进时间、以及所反射的超声回声从每个反射点传播到每个换能器的超声波传播时间的逝去之后,超声回声被每个换能器装置接收。因此,每个反射点提供用于每个换能器的每个接收信号的该延迟时间,即所发送的超声波脉冲的行进时间和超声回声的超声波传播时间的总和,并且通过将这些接收信号值相加来获得在每个点上聚焦的延迟求和信号。通过沿着发送脉冲的行进方向布置每个点处的延迟求和信号来获得深度方向上的一维介质反射强度分布。然后,通过改变(扫描)发送脉冲的发送位置和接收换能器的位置来获得与介质的反射强度有关的二维反射强度图像。以此方式,当执行动态接收聚焦时,每个信号值是根据每个换能器的超声波传播时间(该时间是在根据超声波脉冲的行进而改变接收焦点21的深度的同时如上确定的)而被延迟的,并且信号值相加。优选的是,相加的信号值变为通过根据所发送的超声波脉冲的位置按时间顺序排列的时间顺序加法信号。接下来描述该处理。
首先将参照图8描述动态接收聚焦的概述。图8A示出描述接收聚焦的概述的图,图8B示出在存储每个换能器的接收信号的存储器中存储接收信号的状态。为了简化解释,仅考虑三个换能器200、201和202以及两个接收焦点212和213。在此省略具有不同声音速度的层。假设换能器201和202相对于换能器200对称地定位。分别以400、401和402来表示从换能器200、201和202获得(输出)的接收信号。
从每个换能器200、201和202获得的接收信号400、401和402被处理为对应接收通道。下文中对每个换能器的处理又称为“对每个接收通道的处理”。
在根据接收超声波的时间的地址处按时间顺序在存储器中存储从由每个换能器接收的超声波依次转换出的接收信号400、401和402。超声回声深度在地址的序列中增加。来自接收焦点212的超声回声信号在接收信号400的情况下存储在地址位置40中,而在接收信号401和402的情况下存储在地址位置41中。以相同方式,来自接收焦点213的超声回声信号在接收信号400的情况下存储在地址位置42中,而在接收信号401和402的情况下存储在地址位置43中。
地址位置40与所发送的超声波脉冲沿着中心轴线20传播到焦点212的时间和所反射的超声波从焦点212传播到换能器200的时间的总和对应。地址位置41与所发送的超声波脉冲沿着中心轴线20传播到焦点212与所反射的超声波从焦点212传播到换能器201和202的时间的总和对应。通过把通过计算上述递归关系获得的接收的超声波的传播时间与所发送的超声波脉冲的传播时间相加来确定这些地址位置。对于对称位置中的换能器(比如换能器201和202),由于接收的超声波的传播时间相同,所以可以仅对于这些换能器中的一个计算地址位置。
具体地说,每个接收通道的地址位置Paddr可以表示如下。
[数学式9]
在此,i是用于表示动态接收聚焦的焦点的索引。Ti是通过上述近似递归关系确定的值。Paddr与紧接着在发送之后在存储器中存储的接收信号的地址位置的增量相对应。方程式9中的第一项和第二项与直到所发送的超声波脉冲到达焦点的传播时间对应,第三项与接收超声波从焦点到达换能器的传播时间对应。用于与具体i有关的每个接收通道的地址位置的信号值的总和与相应的动态接收聚焦的焦点(即在该点处的所发送的超声波脉冲位置处的反射强度)对应。通过根据所发送的超声波脉冲的行进而依次增加接收焦点i的位置的深度,可以根据所发送的超声波脉冲沿着索引i的行进来获得与反射强度相对应的时间序列信号。特别是在本发明中,Ti是根据深度而递归地确定的,因此可以沿着索引i依次执行计算。
如果焦点的间隔dz设置为小值(dz=Vtis·Tclk)以匹配时钟频率的周期,并且深度的初始值设置为介质26和介质27的边界,则Paddr变为如下。
[数学式10]
通过使用i作为读数根据用于每个接收通道的i依次读取存储器中存储的接收信号值,并且对于所有接收通道将接收信号值相加,获得在边界处的介质27中的延迟求和信号。更具体地说,与用于所发送的超声波脉冲传播到焦点zi的定时相对应地对于每个换能器确定从焦点zi起的超声波传播时间Ti,并且计算地址Paddr。通过从用于每个换能器的接收信号的计算出的地址位置读取接收信号的值并将其相加,获得与焦点zi对应的延迟求和信号。由于焦点的间隔dz匹配时钟频率,所以可以使用粗略地匹配时钟频率的间隔来在根据所发送的超声波脉冲的行进的深度方向上获得具有高分辨率的时间序列信号。
可以通过以下递归关系而不是方程式10来确定地址位置Paddr。
[数学式11]
Paddr(i+1)=Paddr(i)+dPaddr(i)
由于介质26中的延迟求和信号是通过传统延迟控制方法获得的,因此可以在介质26和介质27的边界处切换处理。在此情况下,i可以用作边界下方的区域中的读数。
优选的是,根据发送脉冲的传播定时按时间序列来布置动态接收聚焦的延迟求和时间序列信号。为了减少波束形成处理的时间,优选的是,在较深接收聚焦的序列中计算地址位置。
特别地,延迟求和处理是根据所发送的超声波脉冲的传播定时来执行的,从而可以在将时间序列接收信号存储在存储器中的同时、根据处理在确保预定时间滞后之后开始用于延迟求和处理的读取操作。由此,部分并行处理变为可能,并且通常可以减少波束形成的处理时间。
如果根据从目标焦点接收超声波信号的定时已经完成地址位置的计算,则与读取存储器同时地执行延迟求和处理,这进一步提高速度。用于实现该操作的可能方法是预先计算地址位置,但需要存储器存储地址位置。特别地,如果设置接收焦点的短间隔以获得高图像分辨率,则需要用于存储很多存储器地址的存储器容量。而在本发明的情况下,根据接收焦点的深度通过计算处理来递归地计算接收的超声波的传播时间,从而可以与延迟求和处理的进行同步地并行计算地址位置。因此,除了存储用于误差校正的相对较少数量的基准地址之外,可以在不提供用于存储很多存储器地址值的存储器的情况下实现高速波束形成。此外,消除了对归因于用于存储存储器地址值的存储器容量而导致的对接收焦点的数量的主要限制,从而可以执行用于更高的接收焦点密度的动态接收聚焦,并且与折射校正关联的该优点实现高分辨率。
因此,根据本发明,使用一系列递归操作在增加接收焦点的深度的同时依次计算地址,由此可以高速地用数字电路配置(对传统配置没有主要附加)并且通过执行折射校正(正如现有技术那样)来生成良好质量的图像。
(误差校正)
为了防止上述近似递归关系计算中的误差的累积,还优选的是,包括通过在基准焦点(基准深度)处的基准传播时间的校正处理。图6示出包括动态接收聚焦中的误差校正的存储器地址计算处理的示例。
首先,在步骤S1中,递归关系的初始值ξ0、η0和T0被例如使用方程式6计算并且设置。同时,读取并且设置包括用于近似误差校正的基准深度zm、基准传播时间Tm和基准递归参数ξm和ηm的表。这里,基准传播时间Tm被写为地址值。对与深度对应的读取计数器进行初始化。
然后,确定深度是否为基准深度(步骤S2)。如果不是基准深度(S2:否),则通过计算上述递归关系来计算在该深度处的ξi、ηi和超声波传播时间Ti,并且基于超声波传播时间Ti而确定读取地址值(步骤S3)。如果其为基准深度(S2:是),则基于对应基准传播时间来确定基准地址值,并且该值被设置为地址值(步骤S4)。同时,递归关系的参数被更新为校正基准值ξm和ηm。
实际存储器的地址被根据如此确定的地址值来确定,并且用于接收波束形成的延迟求和处理(步骤S5)。然后,读取计数器递增,并且深度增加达到焦点的间隔量dz(步骤S6)。重复地址计算处理,直到到达读取计数器结束的深度。当到达读取计数器结束的深度时,处理结束。
以此方式,每次读取计数器递增,就在每个接收通道中确定存储接收信号的存储器的地址位置。读取如此确定的地址位置处的接收信号值,并且对于所有接收通道把这些接收信号值相加,由此执行动态接收聚焦的延迟求和处理。
<其它计算方案>
在以上描述中,使用方程式2、方程式3、方程式4和方程式5来计算递归关系。然而,本发明的特征在于,对于每个接收通道独立地计算递归关系,这不一定受限于上述方程式。如果可以随着接收焦点改变而通过针对每个接收通道的一系列递归运算来依次计算地址,则可以实现以上发明的效果。
换句话说,如果在不使用其它接收通道的计算结果的情况下确定超声波传播时间Ti相对于深度间隔dz的增量dTi(如方程式3至方程式5那样),并且使用方程式2将结果递归地相加,则可以使用其它表达式。因此,可以使用除了方程式3、方程式4和方程式5之外的方程式,只要该方程式是如下近似表达式:使用该近似表达式,每个接收通道可以独立地计算超声波传播时间的增量dTi。
(变型1)固定增量dTi
在执行上述近似校正的情况下基准深度值之间的预定增量可以用于每个接收通道。换句话说,如果深度zi满足zm≤zi<zm+i,则dTi可以是dTi=Δm。将再次参照图6描述该处理。
首先,在初始化时(步骤S1),读取初始值ξ0、η0和T0,并且读取包括基准深度zm、基准传播时间Tm和基准传播时间增量Δm的表。该基准数据被设置用于每个接收通道。
为了计算在除了基准深度之外的位置处的地址值(步骤S3),根据方程式2(Ti+1=Ti+dTi)计算传播时间Ti,并且根据该值确定地址值。在此情况下,此时的基准传播时间增量Δm用于增量dTi。另一方面,为了计算在基准深度处的地址值(步骤S4),在基准深度处根据基准传播时间Tm确定地址值。
此后的处理与以上相同。换句话说,从地址值计算实际读取地址(步骤S5),读取计数器值递增(步骤S6),并且重复该处理,直到读取结束。
基准传播时间增量Δm可以是被设置用于每个基准深度zm的每个间隔的任何值,例如对使用方程式7和方程式8对于每个基准深度zm之间的每个间隔计算出的高度精确的超声波传播时间的曲线进行拟合的斜度。可以利用两个邻近基准深度zm与zm+1之间的基准传播时间值Tm和Tm+1的改变。在此情况下,基准传播时间增量Δm可以从相邻的基准传播时间值Tm(Δm=dz·(Tm+1-Tm)/(zm+1-zm))简单地被计算,从而无需在表中存储基准传播时间增量Δm。
在该示例中,可以省略或简化用于确定增量dTi的计算处理,从而计算处理尺寸较小,并且可以极大地省去待实现的数字电路。另一方面,与使用方程式3、方程式4和方程式5的情况相比,近似精度下降,从而需要在更多基准深度zm处的基准传播时间值Tm,这增加用于存储数据的存储器容量。
在以上示例中,由于深度间隔dz是恒定的,所以使用基准传播时间增量Δm,但如果深度间隔dz改变,则应使用基准深度之间的传播时间的增加率。换句话说,通过dTi=dzi·(Tm+1-Tm)/(zm+1-zm)来确定传播时间。在这两种情况下,基于相对于邻近基准深度中的基准传播时间值的改变量来确定基准深度之间的超声波传播时间增量。
(变型2)基于深度的多项式的增量dTi的近似
作为使用方程式3至方程式5的近似递归关系的计算与使用线性插值的计算之间的方法,可以对于每个基准深度zm使用深度的多项式来近似传播时间增量dTi。
[数学式12]
dTi (m)=F(x)(zi-zm)···式12
在此,dTi (m)是在基准深度zm与基准深度zm+1之间使用的超声波传播时间增量。F(x)(zi-zm)是深度zi与基准深度zm之间的差(zi-zm)的多项式,并且其系数是对于每个接收通道和基准深度确定的。例如通过对使用方程式7和方程式8确定的值进行拟合来计算系数。当从诸如水听器之类的点声源接收超声波脉冲时,可以基于信号的实际测量值来执行拟合。
在图6所示的处理流程中,除了基准深度zm和基准传播时间Tm之外,多项式F(x)(zi-zm)的每个系数也必须在初始化(S1)期间以及在基准深度处的更新(S4)期间被初始化或更新。在该方法中,必须对于每个基准深度存储多项式的系数,但与将超声波传播时间增量设置为恒定的情况相比,近似精度增加,从而优点是基准深度的数量可以较低。即使与使用方程式3、方程式4和方程式5的情况相比近似精度下降,但在此的优点是可以简化超声波传播时间增量dTi的计算。
(其它变型)
如果超声波传播时间增量可以通过少量参数来表示并且可以简单地由数值来表示,以及如果可以使用样条函数、小波函数、各种正交函数扩展等而不使用多项式,则用于近似超声波传播时间增量dTi的方法可以是任何方法。在任何情况下,对于使用方程式7和方程式8计算出的值或在邻近的基准深度之间的实际测量值执行参数拟合,以上述相同方式在基准深度处的更新期间更新参数,以及通过参数指定的近似函数可以用于基准深度之间的递归关系表达式2。
在以上描述中,主要示出把本发明应用于存储器地址计算以用于接收期间的波束形成的示例,但本发明可以应用于在发送期间的波束形成以控制延迟。
以上描述主要示出通过将超声波照射到被检体上并且接收其反射的超声波(回声)来生成图像的超声波回声诊断装置。然而,由于需要相似的接收波束形成,所以本发明也可以应用于通过将脉冲化的光辐射到被检体上并且接收通过光声效应生成的光声波(超声波)来生成图像的光声成像装置。
此外,在以上描述中,示出当存在其声音速度不同的一个***板或一个层时的情况,但这种层的数量可以是两个或更多个。本发明也可以应用于当其声音速度不同的这样的层不存在的情况。
[示例1]
现将使用示例详细描述将本发明应用于超声波图像诊断装置的情况。
<总体配置>
现将参照图1描述根据该示例的超声波图像诊断装置(超声波回声诊断装置)。
图1A示出根据该示例的超声波成像装置的配置。图1B示出根据该示例的超声波成像装置中使用的超声波束。图1C示出由与图1B所示的超声波束对应的超声波成像装置生成的图像或图像数据。
在图1A中,标号301表示输入单元,302表示控制器,303表示发射机,304表示孔径选择切换器,305表示换能器阵列,306表示***板,307表示被检体,308表示接收波束形成器,309表示信号处理器,310表示图像处理器,以及311表示显示器。
输入单元301用于基于要捕获被检体的图像的条件而把各种参数输入到控制器302,并且由标准开关和音量调控器或键盘等构成。待输入的参数与标准超声波图像诊断装置相同,但待输入的该示例的特征参数是***板306的厚度和声音速度以及被检体307的声音速度。在把该超声波成像装置构建为超声***造影法的情况下(其中,用于紧固作为被检体307的胸部的压迫板用作***板306),***板306的厚度和声音速度在装置中已经是已知和固定的。在此情况下,例如,可以在控制器302中的非易失性存储器中预先存储这些参数。对于被检体307的声音速度,取决于待成像的被检体307的区域,可以选择预先存储的若干类型的值中的值,或可以使用预先存储的代表性值。
控制器302初始化超声波成像装置,并且基于正输入的参数来执行各种设置,然后控制其它单元中的每一个。控制器302首先把发送控制信号发送到发射机303,生成其中加上了发送延迟的被发送的超声波驱动脉冲,并且同时使用孔径选择切换器304来选择发送/接收孔径并且将其切换到发送侧。然后紧接着在超声波脉冲从换能器阵列305的发送结束之后,控制器302将孔径选择切换器304切换到接收侧。然后,控制器302控制接收波束形成器308和信号处理器309的处理,使用孔径选择切换器304重新选择发送/接收孔径,并且通过重复以上控制操作来扫描超声波发送/接收波束。当扫描的一帧结束时,控制器302控制图像处理器以创建帧图像。
发射机303基于来自控制器302的发送控制信号来生成用于驱动换能器阵列305上的每个换能器并且生成超声波的发送驱动脉冲。同时,发射机303把与换能器阵列305上的每个换能器的位置和发送焦点相对应的延迟添加到每个换能器的发送驱动脉冲,以便生成发送汇聚波束。此外,如下所述,发射机303生成发送驱动脉冲并且仅对于孔径选择切换器304选择的换能器延迟该脉冲,以创建用于发送的孔径。
在此添加的延迟用于发送波束的汇聚。在发送中,脉冲化的超声波即刻发送,从而通常可以仅设置一个焦点。因此,待设置的延迟时间值的群组的数量保持较低,并且即使在控制器302中存储每个延迟时间值,存储器体积也不增加太多。因此,预先计算的折射校正延迟量可以存储在控制器302中,并且用于焦点。
与接收波束的汇聚相比,发送波束的汇聚不影响图像质量,从而可以执行不校正折射的传统延迟控制处理。在此情况下,可以通过简单的计算处理来控制延迟,从而简化装置配置,但使用以下平均声音速度近似的延迟量计算特别有效。
[数学式13]
在此,Vtis是被检体307中的声音速度,Vm是***板306中的声音速度,D是***板306的厚度,f是作为从换能器阵列305到发送焦点的深度的发送焦距。如果这些值对于设备是唯一的,例如使用固定板作为***板306(其为固定板)的情况,则它们从输入单元301被输入作为输入参数,或部分地存储在控制器302中作为预定值。符号x是从发送波束的中心测量的在换能器阵列305上的换能器的位置,并且由换能器阵列305上的换能器元件的间隔以及元件相对于发送波束的中心的位置确定。使用该表达式,可以对于与形成孔径的每个换能器对应的每个发送通道来计算期望的发送焦点位置中的延迟量。
以上表达式用于假设超声波在将换能器和焦点连接的线上传播而考虑声音速度值的差来确定延迟时间。折射的影响未被考虑,但已经近似地校正具有不同声音速度的每个介质中的相应传播时间值。特别地,可以通过将微小改变添加到传统延迟计算数字信号处理电路来实现该表达式,该表达式与在均匀声音速度条件下确定延迟时间的传统方式相似并且不产生计算量的太多增加。执行该计算处理使得不必需预先存储每个延迟时间量。
由发射机303生成并且延迟的发送驱动脉冲经由孔径选择切换器304输入到在换能器阵列305上形成发送孔径的每个换能器,超声波脉冲从每个换能器被发送。
换能器阵列305具有多个换能器(发送/接收元件)。换能器被一维地阵列化并且可以基于发送/接收超声波的延迟控制通过波束形成处理来切换发送/接收波束的焦点位置。换能器阵列用于发送和接收超声波,并且充当超声波发射机和超声波接收机。孔径选择切换器304切换并且连接换能器阵列305上的每个换能器、发射机303的每个发送通道以及接收波束形成器308的每个接收通道。与切换每个换能器的发送/接收同时地,在换能器阵列305上选择孔径部分,并且通过切换待选择的换能器来移位孔径,由此扫描超声波发送/接收波束。例如,形成发送/接收孔径的换能器的数量是几十个单元。
将参照图1B描述超声波发送/接收波束的扫描。为了简化解释,未示出***板306。在换能器阵列305上使用一部分换能器(图1B仅示出三个换能器)形成孔径312。从孔径312的每个换能器发送的超声波被组合,并且形成发送超声波束313。对于接收,从孔径312的每个换能器获得的接收信号被处理,由此以相同方式形成接收超声波束313。超声波束313的中心变为超声波束扫描线314。通过该发送/接收超声波束313,获得沿着超声波束扫描线314的区域中的被检体组织315的超声波回声信号。超声回声信号反映沿着超声波束扫描线314的区域中的被检体的声学结构,随着超声波脉冲传播在根据深度的时间/位置中出现反射信号。
通过切换孔径选择切换器304并且依次选择换能器以构成孔径312,孔径312在换能器阵列的阵列方向上移位。由此平行移动每个超声波束扫描线314,并且可以获得反映被检体的声学结构的超声波的反射信号的截面数据。
在该示例中,通过切换孔径选择切换器304来移位孔径312,从而执行超声波束扫描。应用发射机303中的驱动脉冲的发送通道的数量和处理接收波束形成器308中的接收信号的接收通道的数量可以与构成孔径312的换能器的数量相同。以此方式,该示例具有这样的优点:发送通道和接收通道的数量可以小于换能器阵列305的换能器的总数量,并且可以简化电路配置。然而,发射机303的发送通道的数量和接收波束形成器308的接收通道的数量可以大于孔径312的换能器的数量,或可以与换能器阵列305的换能器的总数量相同。
在以上描述中,孔径312的换能器的数量对于发送和接收是相同的,但孔径312的大小在发送与接收之间可以改变。在图1B中,用于发送和接收的超声波束共同地示出为超声波束313,但波束的形式在发送与接收之间可以不同。
从换能器阵列305发送的超声波脉冲经由***板306发送到被检体307。***板306近似地是平面平行板,并且如果其厚度和声音速度是已知的,则可以使用任何介质,并且由于声音速度与被检体307不同,因此在波束形成时在延迟控制中需要校正。例如,***板可以是在换能器阵列305的阵列方向和垂直于超声波发送方向的方向上会聚超声波的声学透镜,或用于抑制归因于换能器与被检体307之间的声学阻抗的差而引起的超声波的反射的匹配层。以此方式,***板306可以是其声音速度与被检体307不同的薄层。或者,***板306可以是其声音速度与被检体307不同的厚层(例如用于紧固被检体307的固定板或当因为被检体307的形状所以被检体307不直接接触换能器阵列305时使用的支架)。***板306也可以是被检体中的脂肪层。在脂肪层的情况下,通过在没有折射校正的情况下使用正常延迟控制或使用正常脂肪层的平均声音速度发送/接收超声波而获得的图像来首先测量脂肪层的厚度。通过经由输入单元301输入该数据作为***板306的厚度,可以获得具有折射校正的清楚图像。
在穿过***板306之后发送到被检体307的超声波脉冲在被检体307中的每个区域中被反射,反射的超声波经由***板306被换能器阵列305的孔径再次接收。孔径上的每个换能器把接收的超声波转换为电接收信号,通过被切换用于接收的孔径选择切换器304把每个接收信号输入到接收波束形成器308。接收波束形成器308通过相加来生成与接收波束扫描线对应的延迟求和信号,以使得根据输入的每个接收信号来形成接收波束,并且将其输出到信号处理器309。信号处理器309对延迟求和信号执行带通滤波处理、相位检测处理、动态范围转换处理等,以生成具有与超声波的反射强度对应的亮度(强度)的信号,并且将信号输入到图像处理器310。
将参照图1B和图1C描述上述对应关系。图1C中的标号317指定从一个接收通道获得的亮度信号序列。亮度信号序列317是由图1B中的超声波束扫描线314上的接收超声波束313经由接收波束形成器308和信号处理器309生成的。亮度信号串的亮度(强度)的值与接收超声波束313的反射强度相对应,亮度信号串上的位置与接收超声波束313的反射部分的深度对应。作为在被检体307的截面上的超声波反射强度的二维分布的图像数据316的一个帧由发送/接收波束的扫描创建的多个亮度信号串的集合构成。特别地,生成与被检体组织315的每个部分的超声波反射系数对应的图像分布318。
为了执行以上操作,图像处理器310以根据从信号处理器309输入的深度值的顺序把按时间序列排列的亮度信号串临时存储。信号处理器310根据超声波束扫描来存储多个亮度信号串,并且当扫描的一个帧结束时执行扫描转换处理。显示器311通过格栅扫描来显示图像,其中,在图像的横向方向上设置显示扫描线的方向,但时间序列的亮度信号串在深度方向(即被检体207的纵向方向)上排列。因此,图像处理器310执行扫描转换处理,以便使数据阵列与显示器311的显示扫描相匹配。同时,图像处理器310还执行各种插值处理、坐标变换处理和图像滤波处理(不限于上述处理)以使亮度信号串的数据与显示单元的像素相匹配。在此仅示出图像处理器310的基本操作以简化描述,但可以添加用于改进图像质量的各种操作,操作的细节的差异不应干扰应用本发明。
除了由发送/接收超声波生成的图像显示数据之外,如上所述,图像处理器310还可以组合各种数据并且将其输出到显示器311作为图像显示数据。这些数据的示例是:控制器302中存储的每个参数;对于观测被检体307的图像有效的各种信息,如比例;以及用于协助输入的显示数据。
显示器311可以是CRT、液晶显示器等,其显示从图像处理器310输出的图像。
通过重复以上操作,几乎实时地在显示器311上显示被检体307的断层图像。
<接收波束形成器>
现将参照图2描述执行动态接收聚焦的接收波束形成器308。
具有放大器的A/D转换器11放大从孔径选择切换器304输出的接收信号,并且使用时钟频率来把信号数字化。存储器12存储来自具有放大器的A/D转换器11的输出数字信号。根据时钟频率而被依次数字化的接收信号被存储在存储器12中作为与接收时间对应的离散时间序列信号。读取地址计算单元13计算与在每个接收信号的每个接收焦点处的传播时间对应的读取地址,以执行动态接收聚焦。存储器读取控制器14使用来自读取地址计算单元13的读取地址值来读取存储器12中的数字信号值。与发送脉冲传播定时对应地执行读取操作,并且对于各个接收深度聚焦的接收信号幅度值被输出到各个接收通道。加法器15将这些信号值相加,由此获得动态接收聚焦的延迟求和信号,其输出到信号处理器309。对于多个接收通道提供具有放大器的A/D转换器11、存储器12、读取地址计算单元13和存储器读取控制器14,使得并行执行处理并且实时生成图像。
在图2中,对于所有接收通道提供地址计算单元13。然而,如果换能器阵列305的配置使得:换能器以相等间隔布置,并且构成接收孔径的换能器相对于其中心对称地布置,则传播时间对于在对称位置中布置的换能器是相同的。因此,可以仅对于接收孔径中对称位置中布置的接收通道中的一个布置地址计算单元13,使得另一个接收通道共享该输出值,于是可以简化电路配置。
图2示出当接收波束形成器308由数字电路构成时的每个块的概述。可以使用固定数字电路或可编程数字电路(如FPGA)来实现读取地址计算单元13和存储器读取控制器14。可以替代地使用具有标准处理器(CPU)的计算机。在接收波束形成器308中,每个接收通道并行地执行延迟处理,并且形成接收波束。本发明可以被看作这样的延迟控制方法。如果接收波束形成器308由FPGA或CPU构成,则该程序可以被看作本发明的延迟控制程序。
<读取地址计算单元>
现将参照图3详细描述读取地址计算单元(存储器地址计算单元)13。图3示出当读取地址计算单元13由数字电路构成时的每个块的概述。图3中的配置是当迭代地计算读取地址时,如上所述使用基于基准传播时间在多个基准深度(基准焦点)处的校准的处理的示例。在该示例中,与用于在超声波图像诊断装置中同步处理的***时钟频率的周期相对应地设置动态接收波束形成中的接收焦点间隔。***时钟频率的周期与存储器12中存储的数字接收信号的时间分辨率相对应,可以通过将接收动态波束形成的接收焦点的间隔对应于***时钟频率的周期来执行其中在接收信号的所有深度中匹配焦点的延迟求和处理。由此超声波成像装置中生成的图像的分辨率改进。设dz是接收焦点的间隔,并且设Tclk是时钟频率的周期,则dz由dz=Vtis×Tclk给出。时钟信号还用于延迟求和处理的读取信号的基准。
控制器302将时钟信号6000、基准深度时钟值6001、地址校准值6002、用于校准的参数值6003、开始/结束控制信号6006和地址增量操作参数6007输出到读取地址计算单元13。从控制器302接收这些输入,读取地址计算单元13将与超声波传播时间对应的读取地址值6005输出到存储器读取控制器14。
基准深度时钟值6001是在时钟频率的每个周期处测量的基准深度。时钟计数器61通过累计时钟信号6000来创建时钟累计值,并且在时钟频率的每个周期处对逝去的时间进行计数。
地址校准值6002是在与每个基准深度时钟值6001对应的基准传播时间处的地址值。用于校准的参数值6003是用于与每个基准深度时钟值对应的校准的参数值。如果使用上述方程式3、方程式4和方程式5来执行一系列递归运算,则用于校准的参数值6003是基准递归参数ξm和ηm。如果使用基准传播时间增量Δm,则用于校准的参数值6003是转换为时钟计数的基准传播时间增量Δm。如果使用近似多项式,则用于校准的参数值6003是多项式的系数。在控制器302中使用基准深度时钟值作为索引把这些用于校准的参数值6003存储为查找表。从控制器302的查找表读取的值分别存储在临时存储器601、602和603中。当下述基准深度控制信号6011被发送到控制器302时,控制器302更新从查找表读取的值。在查找表的开始,存储扫描开始时钟值、折射校正开始时钟值、扫描初始地址值和折射校正初始地址值。在此,扫描开始时钟值是与用于开始动态接收聚焦的时间对应的时钟值。折射校正开始时钟值是与超声波往复地传播***板306的厚度的时间对应的时钟值。扫描初始地址值是与开始动态聚焦的时间对应的地址值,折射校正初始地址值是与开始折射校正的时间对应的地址值。对于接收焦点存在于***板306中的动态接收聚焦开始点和***板306的边缘,当以均匀声音速度执行正常延迟控制时,这些地址值被确定为接收波束形成的读取时钟定时和地址值。以此方式,可以指定用于开始动态接收聚焦和折射校正的时钟计数值及其初始地址值。
开始/结束控制信号6006是用于控制读取地址计算单元13的操作的开始/结束的信号。
对于地址增量计算参数6007,用于在下述均匀声音速度地址增量计算块65和折射地址增量计算块64中的计算的计算参数被控制器302预先在每个块中设置。具体地说,设置声音速度Vtis和Vm、***板306的厚度D以及动态聚焦的焦点的间隔等。
现将描述用于该示例的地址值。用于读取存储器12中存储的信号的地址是整数。然而,如果与***时钟对应地设置动态接收聚焦的焦点的间隔,如在该示例的情况下那样,则与超声波传播时间对应的时间是实数值。在该示例中,递归地计算地址,从而如果地址写为整数,则在地址值中生成取整舍入误差。如果生成取整舍入误差,则延迟求和时间的精度恶化,并且图像分辨率下降。为了防止该情况,根据该示例,以包括整数部分和小数部分的实数形式来计算地址值。换句话说,地址值的高有效位用作整数部分,其低有效位用作小数部分。地址值的高有效位整数部分被称为“地址值整数部分”,地址值的低有效位小数部分被称为“地址值小数部分”。只要保证所需精度,则地址值小数部分的位的数量就可以是任何值,并且优选地是2至8位,尽管这取决于时钟频率、待捕获的图像的深度以及待使用的超声波的频率。
下述存储器读取控制器14使用在通过对实数地址取整舍入生成的两个地址处的用于接收信号值的地址值小数部分来执行插值。由此,与存储器12中存储的接收信号的***时钟频率的周期对应的时间分辨率***值,可以执行高度精确的延迟求和处理,并且伴随着动态接收聚焦匹配***时钟频率的周期的效果分辨率被进一步改进。以此方式,读取地址计算单元13内部使用的地址值是包括从控制器302的查找表获得的地址校准值6002的实数地址值。这对于以下描述中使用的地址值也是正确的。
现将连同整个(装置)的操作一起描述地址计算单元13的每个块。
在基准深度时钟值6001和地址校准值6002被从控制器302发送并且设置在临时存储器601和602中之后,在开始的定时从控制器302发送开始/结束控制信号6006。当接收到开始/结束控制信号6006时,比较器62输出增量计算控制信号6010,其命令均匀声音速度地址增量计算块65开始操作。比较器62设置切换器66,以使得将均匀声音速度地址增量计算块65输出到临时存储器604。
如果时钟计数器61的时钟累计值匹配临时存储器601中的基准深度时钟值6001(在此情况下,扫描开始时钟值),则比较器62将基准深度控制信号6011输出到切换器67和68以及控制器302。切换器67通过基准深度控制信号6011的输入来将地址校准值6002(此时,扫描初始地址值)存储在临时存储器605中,并且将地址校准值6002设置为初始地址值。接收到基准深度控制信号6011的输入,控制器302更新在临时存储器601、602和603中的基准深度时钟值6001、地址校准值6002和用于校准的参数值6003。此时,折射校正开始时钟值作为基准深度时钟值6001存储在临时存储器601中,折射校正初始地址值作为地址校准值6002存储在临时存储器602中。折射地址增量计算的初始参数作为用于校准的参数值6003存储在临时存储器603中。
均匀声音速度地址增量计算块65与时钟信号6000同步地计算***板306中的地址增量值6004。如上所述,当动态聚焦处理开始时,待从均匀声音速度地址增量计算块65输出的地址增量值从切换器66输出,并且存储在临时存储器604中。加法器69将临时存储器604中的地址增量值和在先前时钟处计算出并且存储在临时存储器605中的输出地址值相加,以使得计算输出地址值6005。以此方式,与时钟信号6000同步,可以递归地计算并且依次输出与接收焦点的深度对应的地址值。
随着时间逝去,接收焦点变得更深,接收焦点最终到达***板306和被检体307的边界。此时,时钟计数器62的时钟累计值是折射校正开始时钟值。如果时钟累计值与存储器601中的基准深度时钟值6001(此时,折射校正开始时钟值)匹配,则比较器62输出增量计算控制信号6010。这里,增量计算控制信号6010用于指示均匀声音速度地址增量计算块65的操作的结束以及折射地址增量计算块64的操作的开始。比较器62还被设置为使得:切换器66把折射地址增量计算块64的输出输出到临时存储器604。下文中,增量计算控制信号6010的输出停止,直到开始/结束控制信号6006再次输入。
如果时钟累计值与折射校正开始时钟值匹配,则比较器62将基准深度控制信号6011输出到切换器67和68以及控制器302。通过基准深度控制信号6011,切换器67将临时存储器602中的折射校正初始地址值存储到临时存储器605,并且将临时存储器603中的折射地址增量计算的初始参数输入到折射地址增量计算块64。控制器302接收基准深度控制信号6011的输入,并且把临时存储器601、602和603中的基准深度时钟信号6001、地址校准值6002和用于校准的参数值6003更新为查找表中存储的下一基准值。
折射地址增量计算块64与时钟信号6000同步地通过考虑归因于***板306的折射来计算地址增量值,并且将值存储在临时存储器604中。然而,恰如***板306中的动态接收聚焦那样,计算出的地址增量值依次与先前时钟处的地址值相加,以随着接收焦点的深度增加而确定输出地址值6005。
以此方式,与时钟信号6000同步地,与接收焦点的深度对应的地址值可以被迭代地计算并且依次输出。如果时钟累计值匹配临时存储器601中的基准深度时钟值6001,则更新用于校正的临时存储器,以使得在下一基准深度中也可以执行相同校正。换句话说,基准深度时钟值6001、地址校准值6002和用于校准的参数值6003被更新为控制器302的查找表中存储的接下来的值。重复以上操作,并且操作在开始/结束控制信号6006中接收到结束信号的点停止。
只要值与时钟信号同步,就无需对于时钟的每个增量执行加法器69与均匀声音速度地址增量计算块65或折射地址增量计算块64的输出的加法以及读取临时存储器605中的输出地址值6005。如果均匀声音速度地址增量计算块65或折射地址增量计算块64中的计算处理花费时间,则可以根据每个块的计算时间和输出来每隔时钟信号的若干时钟计时来执行加法/输出操作。在此情况下,在控制器302中的查找表中写入的基准深度时钟值6001应在匹配该输出操作的定时中写入。
在地址计算单元13的以上示例中,描述了根据在基于多个基准深度的基准传播时间来执行校准的同时进行处理的情况。然而,如果不执行该校正,则可以简化一部分配置,这是由于无需超声波传播时间的计算和对应地址的更新。
均匀声音速度地址增量计算块65中的处理可以与传统技术相同,因此省略其描述。
<折射地址增量计算块>
现将参照图4描述折射地址增量计算块64中的处理的示例。该示例用于使用方程式3、方程式4和方程式5来计算递归关系。图4A示出用于使用方程式3、方程式4和方程式5的递归关系计算处理的块的概述。在此,其中与表达式的对应性相对容易理解的配置被用于简化解释,但本发明的实施例不限于该配置。
该框图通常划分为计算数据的计算子块(由实线表示),用于存储在计算中的值的临时存储器(由虚线表示)以及用于存储在计算中不改变的参数的临时存储器(由双线表示)。
临时存储器701、702、703、704、705、706、707、708和709存储在计算中的中间值。这些临时存储器主要被布置用于在计算中匹配处理的定时。这些临时存储器包括存储器以使得更容易解释方程式3、方程式4和方程式5的每个变量,并且具体地说,例如,临时存储器704、705和706可以省略。
临时存储器710、711、712、713、714和715存储其值在使用递归关系的一系列递归计算中不改变的计算参数。当计算操作开始时,控制器302提供这些计算参数作为地址增量计算参数6007。这些临时存储器尤其在下文中被称为“参数存储器”。
计算子块71、72、73、74、75、76、77和78是表示为子块的其中执行各个计算的部分。
在递归操作的开始,被给出作为校准的参数值6003的值存储在临时存储器701、702和703中。使用存储器701中存储的值和参数存储器711中存储的值,计算机子块71计算方程式4(b)的右手边上的分子部分,并且将结果输出到计算子块73。计算子块72使用临时存储器701和702中存储的值来计算方程式4(a)的右手边,并且将结果存储在临时存储器704中。计算子块73使用来自计算子块71的输出、临时存储器703和704中存储的值以及参数存储器712中存储的值来计算方程式4(b)的右手边,并且将结果存储在临时存储器705中。计算子块74使用临时存储器704和705中存储的值来计算方程式4(c)的右手边,并且将结果存储在临时存储器706中。计算子块75使用临时存储器701和705中存储的值来计算方程式3(b)的右手边,并且将结果存储在临时存储器707中。计算子块76使用临时存储器702和706中存储的值来计算方程式3(c)的右手边,并且将结果存储在临时存储器708中。计算子块77使用临时存储器703中存储的值以及参数存储器710中存储的值来计算方程式3(d)的右手边,并且将结果存储在临时存储器709中。
此时临时存储器707、708和709中存储的值是在使用递归关系在第一步骤中的递归计算已经完成的情况下的值。临时存储器707、708和709中存储的值再次存储在临时存储器701、702和703中,并且执行在第二步骤中的递归计算。同时,临时存储器707中存储的值被输入到计算子块78,其中,使用参数存储器714和715中存储的值来计算方程式3(a),以在折射校正之后确定传播时间增量。结果然后通过方程式11等转换为地址,从而获得地址增量值6004。
在依次执行周期性递归计算的步骤的同时,地址增量值6004当接收焦点在动态接收聚焦的接收焦点间隔处加深时被依次计算出来。所确定的地址增量值6004存储在临时存储器604中。
根据该示例,如果折射地址增量计算块64中的计算量大并且与***时钟同步地执行该计算,则计算处理可能不在一个时钟周期内完成。在此情况下,必须在地址计算单元13中调整每个块中的时钟同步。此外,地址值的计算时间变为大于***时钟周期若干倍,因此处理减慢。为了防止该情况,折射地址增量计算块64可以由流水线(pipeline)构成,如图4B所示。图4B示出由FIFO存储器710和多个块部分79(其中的一个在图4A中示出)构成的流水线。这里,流水线配置具有四个级,但用于流水线的级的数量是根据块79的计算时间来设置的。在流水线处理中,在读取加法器69的输出所需的定时之前处理开始,曾存储在FIFO存储器710中的地址增量值6004在加法器69的处理定时被读取并且存储在临时存储器604中。在此情况下,需要用于递归处理的初始值,因此应首先读取临时存储器603中的用于校准的参数值6003。
折射地址增量计算块64的该示例是执行计算处理算法的示例,并且尤其是用于协助理解根据方程式3、方程式4和方程式5的处理的描述的示例,因此可以通过省略或改变处理配置的一部分来简化电路配置。
不使用方程式3、方程式4和方程式5,可以使用其中地址增量为部分恒定的基准传播时间增量dTi或多项式近似来简化折射地址增量计算块64的处理。
<存储器读取控制单元>
现将参照图5描述存储器读取控制单元14的示例。在该示例中,存储器读取地址被计算为实数,从而接收波束可以比现有技术更精确地聚焦。以下待描述的存储器读取控制用于有效地利用该优点。具体地说,使用时间上连续并且从存储器12读取的两个接收信号值来执行与输出地址值6005的小数地址部分对应的插值。这允许精度比存储器12中存储的接收信号的时间分辨率更精细的计算。以下解释的使用接收信号插值的接收信号读取对于本发明并非是必要的,并且插值方法也仅是示例,实际的插值计算处理不限于以下方法,例如,可以插值两个或更多个信号值。
图5中的标号81、82和83表示临时存储器,其中,来自读取地址计算单元13的输出地址值6005和从存储器12读取的两个连续接收信号值被存储。临时存储器81、82和83不是必要的,但被布置用于同步处理以及用于协助理解在此的描述。标号84表示减法器,85表示乘法器,86表示加法器。
临时存储器81中存储的输出地址值6005包括作为高有效位部分的地址整数部分801和作为低有效位部分的地址小数部分208。地址整数部分801用于指定并读取与存储器12中存储的时间序列数字接收信号当中的具体时间相对应的信号值。为了解释,这里在下文中,临时存储器81中存储的地址整数部分由n表示,地址小数部分由q表示。
在由地址整数部分801指定的地址处的接收信号值以及在时间上连续的与该部分对应的下一地址的接收信号值分别从存储器12存储在临时存储器82和83中。在此,由Sn表示临时存储器82中存储的接收信号值,由Sn+1表示临时存储器83中存储的接收信号值。下标n表示由地址整数部分801指定的存储器地址。
减法器84从临时存储器82中存储的接收信号值Sn减去临时存储器83中存储的接收信号值Sn+1,并且将结果输入到乘法器85。乘法器85将减法器84的输出与临时存储器81中存储的地址值的地址小数部分q相乘,并且将结果输入到加法器86。加法器86将临时存储器83的输出与乘法器85中存储的接收信号值Sn+1相加,并且将结果输出到加法器15。
通过以上计算,计算以下表达式。
[数学式14]
Sn·q+Sn+1·(1-q).···式14
由此,在通过地址小数部分q对时间上邻近的接收信号值Sn和Sn+1插值之后的值被输出到加法器15。在存储器12中,接收信号按由***时钟的间隔确定的时间间隔存储,但通过该插值,可以获得与其间的超声波传播时间值对应的接收信号。结果,可以在不取决于时钟频率的周期对延迟求和处理操作的目标接收信号值的时间精度取整舍入的情况下高精度地执行延迟求和处理。
然后,通过加法器15进行的上述处理以及稍后的处理对于接收通道的接收信号执行延迟求和,并且超声波断层图像经由信号处理器309和图像处理器310而显示在显示器311上。
在该示例中,超声波发送/接收波束的扫描已经被描述为与选择构成换能器阵列305的孔径的换能器的一部分以及通过使换能器偏移进行的线扫描有关。这用于简化描述,并且整个换能器阵列305可以用作孔径。在此情况下,本发明可以应用于基于分区扫描方法的超声波图像诊断装置,其中,使用每个换能器的发送/接收控制来改变超声波发送/接收波束的发送/接收方向。在此情况下,接收焦点的位置可以设置在从换能器阵列305的中心径向地扩展的接收超声波束上,而不是直接在换能器阵列305的中心正下方的位置。以接收焦点的深度变深的顺序、使用近似递归关系以相对于接收超声波束中的一个而确定用于每个换能器或接收通道的超声波传播时间,来递归地依次确定接收存储器的地址位置。如果对于使用该地址位置而被读取的每个换能器或接收通道的接收信号值执行延迟求和,则获得相对于超声波束的延迟求和转换信号。通过改变接收超声波束和扫描波束的方向,可以创建并显示超声波断层图像,如上所述。
[示例2]
本发明的示例2是基于使用光声波的光声成像方法的超声波成像装置。根据光声成像方法,通过将脉冲层波束辐射到被检体上并且检测作为热膨胀而引起的光声波的超声波,对被检体的内部结构成像。
例如,在日本专利申请公开(PCT申请的翻译)No.2001-507952中如下公开了光声成像方法的具体过程。
(1)二维换能器阵列(二维接收元件阵列)位于被检体表面上,并且单脉冲化电磁能量被辐射到被检体上。
(2)在电磁能量的辐射之后采样并且立即存储每个换能器(每个接收元件)的接收光声信号。
(3)对于待成像的被检体中的点r’计算用于光声波从点r’到达位置r的传播时间,并且与传播时间对应的每个换能器的信号被加在一起以获得点r’的图像值。
(4)对于待成像的每个点r’重复步骤(3)。
日本专利申请公开No.2005-21380公开了使用公共换能器交替地重构基于光声成像方法的图像和正常超声波回声图像的装置。
以此方式,同样根据光声成像方法,必须计算从光声波生成点到每个换能器的光声波的传播时间,并且在使用该计算的传播时间来匹配相位的情况下按时间序列把接收光声信号加在一起。因此,如果在光声波的传播路径上存在其中声音速度是不同的层并且生成了折射,则接收光声信号的相位不匹配,并且创建的图像分辨率下降,恰如基于超声波脉冲回声方法的接收的情况。在该示例中,本发明的超声波传播时间的校正被应用于光声成像方法,因此抑制待创建的图像的分辨率的下降。
将参照图9详细描述基于光声成像方法的三维图像重构方法。
图9A是描述用于根据光声成像方法的三维图像重构的体素数据结构的概念示图。图9B是描述用于在光声成像方法中在特定体素中创建数据的构思的图。图9C是描述根据该示例的处理方法的图。
在图9中,标号305表示换能器阵列。与上述超声波回声方法不同,二维阵列用于重构三维图像。此外,在光声成像方法中,换能器阵列的每个换能器仅接收作为由激光辐射生成的光声波的超声波。
在图9A中,标号91表示三维体素阵列。其中被检体中的成像区域被虚拟地分段为格栅的数据结构被称为“体素阵列”。与该数据结构对应的格栅结构被称为体素阵列,以及构成体素阵列的单元格栅被称为“体素”。体素数据阵列中存储的数据组被称为“体素阵列数据”。在基于光声成像方法的三维图像重构中,每个体素中生成的光声波的强度被评估,被转换为图像亮度作为辐射的激光吸收量,并且存储在体素数据阵列中。基于如此获得的图像亮度数据(亮度体素阵列数据),显示三维图像,或重构并显示断层图像。
在图9B中,标号91a表示三维体素阵列91中的特定体素。标号92表示从体素91a到换能器阵列305的每个换能器的光声波传播路径。基于由换能器阵列305的每个换能器接收的以时间序列的接收光声信号来评估体素91a中生成的光声波的强度。换句话说,与从体素91a到换能器的传播时间对应的部分中的接收信号的幅度值被针对每个换能器确定并相加,由此确定体素91a中生成的光声波的强度。
如果此时被检体的声音速度是均匀的,则可以通过体素91a和每个换能器的几何位置关系以及被检体的声音速度来简单地确定传播时间。另一方面,如果具有不同声音速度的***板存在于被检体与换能器阵列305之间,或被检体的表面区域是层并且其声音速度与较低层部分不同,则通过基于本发明的处理来确定包括折射效应的传播时间。
在此,省略***板的说明以简化附图。为了校准取决于每个换能器的位置以及孔径的大小和方向而评估体素91a的立体角的改变,并且为了增加体素91a中生成的光声波强度的评估中的精度,可以在把每个换能器信号相加时执行加权。可以执行考虑取决于换能器自身的灵敏度分布的角方向性的加权。这些问题的细节对于本发明是次要的,因此在此省略详细描述,本发明关注当其声音速度与被检体不同的***层存在时所需的对传播时间执行的折射校正。
将参照图9C描述当本发明应用于光声成像方法时的处理的概述。标号91b表示体素阵列91当中的在被检体深度方向上排列的体素列。标号93表示体素列91b的中心轴。体素列的中心轴93在近似垂直于二维换能器阵列305的阵列平面的被检体深度方向上伸展。
标号94是换能器阵列305中的任意换能器。在此,假设体素列的中心轴93与换能器94之间的距离(即从换能器94到体素列的中心轴93垂直绘制的长度)是x。下文中x被称为“换能器-体素列距离”。
如果深度方向上的体素列91b的每个体素之间的距离被看作深度方向上的焦点的间隔dz,并且换能器位置x被看作垂线的长度,则光声波的传播时间T变为与方程式1相同。因此,如果确定声音速度不同的***板和被检体的界面上的体素之间的传播时间,则通过按顺序递归地确定光声波传播时间增量,可以对于其深度沿着体素列的中心轴93每次变深dz的各个体素按顺序递归地确定传播时间。具体地说,使用方程式2、方程式3、方程式4和方程式5按顺序递归地计算光声波传播时间值。恰如使用超声回声方法的上述装置,初始光声波传播时间可以确定为均匀声音速度的情况。对于在任意深度处的体素使用通过方程式7和方程式8计算的值作为初始值,可以在更浅的深度的方向上按顺序递归地计算光声波传播时间值。或者,在递归计算中,可以执行用于抑制近似误差的累计的校准。对于递归关系的近似,可以使用多项式近似和正交函数***。可以以与使用超声回声方法的装置相同的方式来执行这些操作。
以此方式,使用根据深度递归地确定的光声波的传播时间,可以从由换能器94接收的按时间序列的接收光声信号中提取从体素列91b上的每个体素生成的光声波强度的分量。如此从每个体素生成的光声信号的贡献与对应体素数据阵列的数据相加,结果被再次存储。然后,体素串91b偏移,并且重复相同操作。通过重复该操作,换能器94接收的接收光声信号的贡献可以与和三维体素阵列91对应的体素数据阵列的每个数据相加。通过对于换能器阵列305的每个换能器执行该操作,可以对于与三维体素阵列91对应的体素数据阵列的每个数据获得其相位对准的接收光声信号的相加值。
在以上描述中,对于体素和换能器的所有组合计算光声波传播时间,以简化解释。但取决于体素阵列91的设置以及换能器阵列305的换能器的间隔的设置,可以使用公共传播时间计算。换句话说,传播时间根据光声波的传播路径被确定,并且根据每个体素、每个换能器和***板的几何位置被确定。例如,针对换能器和体素的集合的传播时间与在位置方面平行偏移达到相同量的换能器和体素的集合的传播时间相同。考虑到该情况,优选的是,在与换能器的二维阵列表面平行的平面上布置的体素的间隔与换能器的间隔匹配或处于整数的比率。然后,通过针对特定换能器和体素列91b的组合在深度方向上递归地计算而确定的光声波传播时间值的序列可以直接应用于处于相同位置关系的另一组合。
现将参照图10描述根据该示例的基于光声成像方法的超声波图像诊断装置。图10A示出根据该示例的基于光声成像方法的超声波图像诊断装置的配置。与示例1相同的构成元件以相同标号表示,其中省略多余描述。
输入单元301用于输入操作设置。控制器302控制整个装置。激光发射控制器101控制并且驱动激光的脉冲辐射。激光束辐射单元102将脉冲激光束辐射到被检体上。透明固定板103用于在该透明固定板103与***板306之间紧固被检体307,并且以允许激光束透射通过的材料制成。换能器阵列305包括二维排列的换能器,用于接收光声波。其中,接收光声信号处理单元104放大由换能器阵列305接收的接收光声信号以增加SNR(信噪比),将放大的接收光声信号转换为数字信号,并且由模拟放大器和A/D转换器等构成。为了进一步改进SNR,可以包括频率滤波处理。图像重构计算单元105基于数字化的接收光声信号使用光声成像方法来重构三维图像,并且创建其中评估了每个体素中生成的光声波的强度的体素阵列数据。图像处理器310使用由图像重构计算单元105创建的体素阵列数据根据来自控制器302的指令来创建待显示在显示器311上的图像数据。图像处理器310显示三维图像数据,并且特别地,创建二维投影数据以及各种断层数据(比如最大强度亮度投影图像(MIP))。显示器311显示图像处理器310的输出。
该示例与示例1的差异在于可以从输入单元301向控制器302指定待显示在显示器311上的数据的类型,以便表示三维图像数据的配置。通过该信息,图像处理器310执行坐标转换、插值、图像滤波等,以基于三维体素阵列数据来创建各种投影图像和断层图像。可以使用评估在每个体素中生成的光声波的强度的数据来创建亮度体素阵列数据。该亮度转换可由图像重构计算单元105执行。
换能器阵列305是用于获得三维数据的二维阵列,如上所述。L表示换能器阵列305的换能器的数量。
***板306是具有与被检体307不同的声音速度的平面平行板,并且允许超声波透射通过,恰如示例1,但在该示例中是特别为了连同透明固定板103一起紧固被检体307的目的而安装的。透明固定板103在该激光束的波长范围中应是透明的,以便辐射脉冲激光,但形状可以仅是适合于保持并且紧固被检体307的形状。
该示例与示例1之间的最大差异在于,由于使用光声成像方法,所以把脉冲激光束辐射到被检体307上而不发射超声波。激光辐射单元102包括可以以高功率来发射脉冲的激光光源。激光束的波长是在被检体307中的透射率相对较高的范围中的波长。特别地,当生物组织是被检体307时,优选的是,发射其波长近似是800nm到1200nm的近红外光,然后可以安装SHG转换器,以用于YAG激光器装置或钛-蓝宝石激光器装置。也可以使用高功率LD阵列光源。激光辐射单元102具有各种光学装置,诸如波束扩展器和扩散板等,以便把发射的激光束均匀地辐射到被检体307上。
激光发射控制单元101通过基于来自控制器302的发射控制信号控制并且驱动激光辐射单元102来辐射脉冲激光。然后,来自控制器302的控制信号开始光声波接收操作。通过激光脉冲辐射在被检体307中生成的光声波由换能器阵列305接收,并且转换为接收光声信号(其由接收光声信号处理单元104放大),然后输入到图像重构计算单元105作为数字采样数据。图像重构计算单元105执行图像重构处理,并且将体素阵列数据输出到图像处理单元310,图像处理器310执行上述图像处理并且将指定的图像显示在显示单元311上。
现将参照图10B详细描述图像重构计算单元105的处理。在图10B中,标号104指示接收光声信号处理器,标号12指示存储器,标号13指示读取地址计算单元,标号14指示存储器读取控制器,标号106指示用于存储体素阵列数据的体素阵列存储器。体素阵列存储器106可以由随机存取DRAM构成或由FIFO存储器与随机存取DRAM的组合构成。在图10B中,体素阵列存储器106概念地示出为具有N行乘M列的矩阵阵列存储器,以简化对操作的描述。在此,N是体素阵列的与换能器阵列表面平行的平面上的划分的总数,M是深度方向上体素阵列的划分的总数。换句话说,作为N乘M阵列存储器的体素阵列存储器106的每个元件与体素阵列中的每个体素对应。关于与图9中的概念示图的对应性,体素阵列91具有数量为N的体素列91b,每个体素列91b包括数量为M的体素。在此情况下,N×M体素阵列存储器106的每个行与体素列91b对应,该行中的每个存储器元件与体素列91b的每个体素对应。特别优选的是,与体素阵列存储器106的每个行的存储器元件对应的每个体素沿着体素列91b的深度按顺序设置。
根据该示例的优选配置,分别针对换能器阵列305的每个换能器,模拟放大器和A/D转换器的集合被设置在接收光声信号处理器104中。由此,每个换能器接收的接收光声信号的时间序列数据串被数字化并且通过并行处理而存储在存储器12中。在存储器12中,对于每个换能器存储数字化的接收光声信号的时间序列数据串。例如,存储器12可以由对于每个换能器单独提供的一维存储器构成。存储器被构建为使得可以对于各个换能器访问二维阵列。
还优选的是,布置数量为N的存储器读取控制器14,以便并行处理特定换能器。数量为N的存储器读取控制器14中的每一个被布置分别与阵列存储器106中的N行对应。恰如示例1,存储器读取控制器14使用由读取地址计算单元13输出的存储器读取地址来在存储器12中与光声波传播时间对应的定时读取接收光声信号的值。读取的信号值被加到阵列存储器106的预定存储器元件(存储器元件中的值和读取的信号值被加一起,并且其结果再次存储在存储器元件中)。在此,在阵列存储器106中待引用的存储器元件在与存储器读取控制器14对应的行上。
读取地址计算单元13使用换能器-体素列距离x和体素的深度的间隔dz来计算与深度依次改变的体素的光声波传播时间相对应的地址,恰如示例1那样。使用该地址,存储器读取控制器14对于阵列存储器106的对应行中的存储器元件依次执行上述加法处理。
可以布置与存储器读取控制器14相同数量的读取地址计算单元13,但可以取决于该换能器和对应体素列91b的位置关系而由公共读取地址计算单元13计算地址。例如,公共读取地址计算单元13可以计算多个存储器读取控制器14的地址,以用于具有相同换能器-体素列距离x的换能器。如果体素阵列91在其下表面上垂直地并且水平地相等地划分,则因为几何对称性,所以同心地围绕作为中心的换能器的体素列91b具有相同换能器-体素列距离x。以此方式,读取地址计算单元13的数量可以小于N。
每个存储器读取控制器14使用由每个读取地址计算单元13递归地计算的存储器读取地址来读取接收光声信号值的值,由此阵列存储器106的每个行中的存储器元件中的值被依次并行地更新。换句话说,阵列存储器106的列由每个读取地址计算单元13根据递归计算处理而依次更新。取决于如何为递归计算提供初始值,更新序列可以在更深的深度或更浅的深度的方向上。
为了更新阵列存储器106中的数据,阵列存储器106可以仅由随机存取DRAM构成,以便直接引用、加上或存储值。或者,值可以一次存储在例如对于每个存储器读取控制器14设置的行存储器或FIFO缓冲存储器中,以便在处理对所有深度结束之后执行随机存取DRAM中的相加和更新处理。
通过对于不同换能器重复上述处理,可以获得其中基于光声成像方法而重构三维图像的体素阵列数据。具体地说,在接收光声信号处理器104处理来自每个换能器的输出并且在存储器12中存储接收光声信号之后,与特定换能器对应的数量为N的存储器读取控制器14更新阵列存储器106中的数据。或者可以提供数量为L×N的存储器读取控制器14,以并行地对于所有换能器处理数据。具体地说,在后一情况下,公共读取地址计算单元13可以对于具有相同换能器-体素列距离x的换能器和体素列的集合来计算存储器读取地址。近来,可以使用GPU中的很多内核进行体积数据的并行处理,因此包括GPU的数字电路可以执行该示例的并行处理。
该示例是以下装置:该装置使用近似递归关系依次对于在不同深度处的目标递归地确定超声波(作为包括归因于具有不同声音速度的***板而引起的折射效应的光声波)的传播时间值,并且使用结果基于光声成像方法来重构三维图像。通过按顺序对于具有不同深度的体素递归地确定传播时间值并且基于所确定的传播时间值来并行地对于每个换能器和每个体素列执行处理,使整个处理时间减少,并且可以实时重构三维图像。换句话说,通常可以对于具有不同换能器-体素列距离x的换能器和体素列的组合来单独递归地计算传播时间值,从而可以对于这些换能器和体素列容易地执行并行处理,并且可以减少整个处理时间。如果读取地址计算单元13对于其换能器-体素列距离x相同的换能器和体素列的组合是公共的,则可以抑制电路大小的增加。
与传统公开技术相比,本发明具有以下优点。
在基于光声成像方法的三维图像重构的情况下,换能器阵列是二维阵列,并且体素布置中的体素数量越大,分辨率越好。这说明,处理所需的光声波传播时间值的数量变多。由于使用很多类型的路径,所以需要长处理时间的美国专利No.6,607,489中所示的用于折射路径的普通计算不适用于实时处理。如果预先计算并且存储延迟时间值,则必须存储很多延迟时间值。特别地,取决于图像捕获条件,每次改变体素阵列就需要新的延迟时间值的集合,并且假设所有可能的情况来存储所有这些值是不实际的。因此,每次体素阵列改变就必须计算延迟时间值,结果,不可避免的是增加装置规模。
在递归地计算延迟时间的WO 01-026555中公开的方法的情况下,恰如本发明,解决了上述问题,但在改变该示例中提到的换能器-体素列距离x的方向上执行递归计算。递归计算基本上需要顺序计算,因此为了通过该方法容易地执行并行处理,必须在深度方向上执行递归计算。这说明,如果体素阵列中的体素的数量在垂直方向上是Nx,在水平方向上为Ny并且在深度方向上为Nz,则可以并行地执行数量为Nz的计算。而在本发明的情况下,可以对于换能器和体素列的集合执行并行处理,也就是说,可以并行地执行数量的Nx×Ny的计算。如果使用适合于生成三维图像的均匀体素阵列,则区段的数量在每个方向上对于Nx、Ny和Nz大约相同,从而该示例的并行处理的程度(可被并行执行的处理的数量)可以高于WO 01-026555中公开的方法。并行处理程度的这种差异随着为了改进分辨率而增加区段数量并更多增加换能器的数量而增大。
此外,在标准格栅类型体素阵列或二维换能器阵列布置的情况下,每个换能器-体素列距离x变为不均匀的。因此,在x方向上的递归关系计算中,近似准确度变为不一致的,并且处理参数增多。在该示例中,可通过将焦点的深度的间隔设置为体素阵列的深度的间隔来执行用于均匀间隔的近似递归关系处理,因此可以简化处理。
因此,如果本发明被用于光声成像方法,则可以执行适合于并行处理的用于延迟时间计算的递归处理,并且作为该示例的情况,可以提供可以实时执行整个处理的基于光声成像方法的超声波图像诊断装置。
在此,描述了基于光声成像方法的超声波图像诊断装置,但可以通过集成示例1的超声发送单元和孔径选择切换器来构成联合地使用超声回声方法的图像诊断装置,因此可以同时显示光声图像和超声波图像。
虽然已经参照示例性实施例描述了本发明,但应理解,本发明不限于公开的示例性实施例。所附权利要求的范围将要被赋予最宽泛的解释,从而包括所有这样的修改以及等效结构和功能。
本发明的方面也可以通过读出并执行存储器设备上记录的程序以执行上述实施例的功能的***或装置的计算机(或诸如CPU或MPU之类的设备)以及其步骤由通过例如读出并且执行存储器设备上记录的程序以执行上述实施例的功能的***或装置的计算机执行的方法来实现。
本申请要求于2010年2月23日提交的日本专利申请No.2010-036890的权益,其全部内容通过引用而合并到此。
Claims (13)
1.一种超声波成像装置,包括:
超声波接收机,具有用于经过介质接收被检体的区域中反射或生成的超声波的多个接收元件,所述介质具有与该区域的声音速度不同的声音速度;
接收波束形成器,根据从接收波束的焦点位置到接收元件的传播时间通过对与所述超声波接收机的接收元件中的至少一部分相对应的接收通道的接收信号执行延迟控制来形成接收波束信号;以及
图像处理器,使用所述接收波束信号来生成图像,其中,
所述接收波束形成器通过反复进行用于把与焦点深度的改变量相对应的传播时间值的改变与已经计算出的传播时间值相加以确定用于下一焦点的传播时间值的处理,来依次计算用于多个焦点的传播时间值,
反复进行的处理能够对于所述接收通道中的每一个被独立执行,并且所述接收波束形成器对于所述接收通道中的至少一部分并行地计算传播时间值,以及
所述接收波束形成器使用包括ξi=cosθ′的递归参数并根据以下等式或以下等式的近似来计算传播时间值的改变:
其中dTi是传播时间值的改变,dz是焦点深度的改变量,Vtis是区域中的声音速度,θ′是被检体的区域中的折射率,并且递归参数ξi是通过使用焦点深度的改变量、该区域中的声音速度以及介质中的声音速度的递归关系而被确定的。
2.如权利要求1所述的超声波成像装置,其中,
递归参数是根据在焦点的相邻深度处的递归参数值的多项式中描述的递归关系而被确定的。
3.如权利要求1所述的超声波成像装置,其中,
所述传播时间值的改变被确定为焦点深度的多项式。
4.如权利要求1所述的超声波成像装置,其中,传播时间值的改变是常数。
5.如权利要求1至权利要求4中任一项所述的超声波成像装置,其中,
所述接收波束形成器存储基准传播时间值,基准传播时间值是针对特定焦点深度预先确定的传播时间值,以及
当在所述反复进行的处理中焦点深度到达所述特定焦点深度时,在该特定焦点深度处的传播时间值被视为所存储的基准传播时间值,并且基于该基准传播时间值来计算在下一焦点深度处的传播时间值。
6.如权利要求2所述的超声波成像装置,其中,
所述接收波束形成器存储基准传播时间值和基准递归参数,基准传播时间值和基准递归参数是针对特定焦点深度预先确定的传播时间值和递归参数,以及
当在反复进行的处理中焦点深度到达所述特定焦点深度时,在该特定焦点深度处的传播时间值和递归参数被视为所存储的基准传播时间值和基准递归参数,并且基于基准传播时间值和基准递归参数来计算在下一焦点深度处的传播时间值和递归参数。
7.如权利要求3所述的超声波成像装置,其中,
所述接收波束形成器存储针对特定焦点深度预先确定的传播时间值以及要被用于所述特定焦点深度的后续焦点的多项式的系数,以及
当在反复进行的处理中焦点深度到达所述特定焦点深度时,在该特定焦点深度处的传播时间值被视为所存储的基准传播时间值,并且在随后处理中使用所存储的系数通过多项式来确定传播时间值的改变。
8.如权利要求4所述的超声波成像装置,其中,
所述接收波束形成器存储针对多个基准焦点而被预先确定的基准传播时间值,以及
当在反复进行的处理中焦点到达基准焦点中的任一个时,在所到达的该焦点处的传播时间值被视为所存储的基准传播时间值,并且通过在所到达的该焦点的深度处的基准传播时间值和在下一焦点深度处的基准传播时间值来确定在随后处理中的传播时间值的改变。
9.如权利要求1至权利要求4中任一项所述的超声波成像装置,其中,
所述接收波束形成器针对接收通道中的每一个都包括:
存储器,其按时间顺序存储从与接收通道相对应的接收元件依次接收的接收信号;
存储器读取控制器,控制从存储器进行读取;以及
存储器地址计算单元,计算存储器的读取地址,
接收波束信号是通过把根据由存储器地址计算单元计算出的读取地址由存储器读取控制器从存储器读取的接收通道中的每一个的接收信号相加来形成的,以及
存储器地址计算单元计算与针对多个焦点依次计算出的传播时间值相对应的读取地址位置,以便依次计算与多个焦点相对应的读取地址。
10.如权利要求9所述的超声波成像装置,其中,
所述存储器地址计算单元以实数形式计算读取地址,以及
接收波束形成器根据计算出的实数形式的读取地址的小数部分来对与计算出的实数形式的读取地址的整数部分相对应的地址和下一地址中存储的两个接收信号进行插值,并且使用所插值的接收信号来形成接收波束信号。
11.如权利要求1至权利要求4中任一项所述的超声波成像装置,其中,
所述超声波接收机还充当用于将超声波发送到被检体的超声波发射机,以及
所述超声波成像装置通过把超声波辐射到被检体上并且接收在被检体中反射了的反射超声波来形成被检体内部的图像。
12.一种超声波成像装置,包括:
超声波接收机,具有用于经过介质接收从被检体的区域中生成的声学波的多个接收元件,所述介质具有与该区域的声音速度不同的声音速度;
图像重构计算单元,基于由所述超声接收机接收的超声波来获得评估每个体素中生成的声学波的强度的体素阵列数据;以及
图像处理器,使用由所述图像重构计算单元获得的所述体素阵列数据来获得图像数据,其中,
所述图像重构计算单元通过反复进行用于把与体素深度的改变量相对应的传播时间值的改变与已经计算出的传播时间值相加以确定用于下一体素的传播时间值的处理,来依次计算用于多个体素的传播时间值,
反复进行的处理能够对于每个体素列独立执行,并且所述图像重构计算单元对于所述体素列的至少一部分并行地计算传播时间值,以及
所述图像重构计算单元使用包括ξi=cosθ′的递归参数并根据以下等式或以下等式的近似来计算传播时间值的改变:
其中,dTi是传播时间值的改变,dz是焦点深度的改变量,Vtis是区域中的声音速度,θ′是被检体的区域中的折射率,并且递归参数ξi是通过使用焦点深度的改变量、该区域中的声音速度以及介质中的声音速度的递归关系而被确定的。
13.一种控制超声波成像装置中的延迟的方法,所述超声波成像装置经过介质接收被检体的区域中反射或生成的超声波,并且通过根据在接收波束的焦点位置与接收元件之间的传播时间执行延迟处理来形成接收波束,所述介质具有与该区域的声音速度不同的声音速度,所述方法包括:
通过反复进行用于计算与焦点深度的改变量相对应的传播时间值的改变的处理以及用于把计算出的传播时间值的改变与已经计算出的传播时间值相加以确定用于下一焦点的传播时间值的处理,来依次计算用于多个焦点的传播时间值,
使得反复进行的处理能够对于每一个接收通道被独立执行,并且对于所述接收通道中的至少一部分并行地计算传播时间值,以及
使用包括ξi=cosθ′的递归参数并根据以下等式或以下等式的近似来计算传播时间值的改变:
其中dTi是传播时间值的改变,dz是焦点深度的改变量,Vtis是区域中的声音速度,θ′是被检体的区域中的折射率,并且递归参数ξi是通过使用焦点深度的改变量、该区域中的声音速度以及介质中的声音速度的递归关系而被确定的。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010036890A JP5645421B2 (ja) | 2010-02-23 | 2010-02-23 | 超音波画像装置および遅延制御方法 |
JP2010-036890 | 2010-02-23 | ||
PCT/JP2011/054422 WO2011105593A1 (en) | 2010-02-23 | 2011-02-22 | Ultrasonic imaging apparatus and method of controlling delay |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102770079A CN102770079A (zh) | 2012-11-07 |
CN102770079B true CN102770079B (zh) | 2015-03-11 |
Family
ID=44063956
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201180010435.9A Expired - Fee Related CN102770079B (zh) | 2010-02-23 | 2011-02-22 | 超声波成像装置和控制延迟的方法 |
Country Status (5)
Country | Link |
---|---|
US (1) | US9165552B2 (zh) |
EP (1) | EP2538846B1 (zh) |
JP (1) | JP5645421B2 (zh) |
CN (1) | CN102770079B (zh) |
WO (1) | WO2011105593A1 (zh) |
Families Citing this family (35)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2009052293A1 (en) * | 2007-10-16 | 2009-04-23 | Zymogenetics, Inc. | Combination of blys inhibition and anti-cd 20 agents for treatment of autoimmune disease |
JP5637725B2 (ja) | 2010-04-12 | 2014-12-10 | キヤノン株式会社 | 音響波イメージング装置 |
JP5721462B2 (ja) | 2011-02-09 | 2015-05-20 | キヤノン株式会社 | 被検体情報取得装置 |
JP5655021B2 (ja) | 2011-03-29 | 2015-01-14 | 富士フイルム株式会社 | 光音響画像化方法および装置 |
JP2012235832A (ja) | 2011-05-10 | 2012-12-06 | Canon Inc | 被検体情報取得装置 |
JP2013102960A (ja) * | 2011-11-14 | 2013-05-30 | Fujifilm Corp | 超音波診断装置及び方法 |
JP5865050B2 (ja) | 2011-12-15 | 2016-02-17 | キヤノン株式会社 | 被検体情報取得装置 |
JP5905856B2 (ja) * | 2012-07-13 | 2016-04-20 | 富士フイルム株式会社 | 超音波検査装置 |
JP5840181B2 (ja) * | 2012-08-17 | 2016-01-06 | 富士フイルム株式会社 | 光音響画像生成装置及び方法 |
JP6039305B2 (ja) | 2012-08-23 | 2016-12-07 | キヤノン株式会社 | 被検体情報取得装置、情報処理装置および被検体情報取得方法 |
JP2014068756A (ja) * | 2012-09-28 | 2014-04-21 | Fujifilm Corp | 超音波検査装置、超音波画像データ生成方法およびプログラム |
EP2749209A1 (en) * | 2012-12-28 | 2014-07-02 | Canon Kabushiki Kaisha | Object information acquisition apparatus, display method, and program |
US20140182383A1 (en) | 2012-12-28 | 2014-07-03 | Canon Kabushiki Kaisha | Object information obtaining device, display method, and non-transitory computer-readable storage medium |
CA2903479C (en) | 2013-03-15 | 2023-10-10 | Butterfly Network, Inc. | Monolithic ultrasonic imaging devices, systems and methods |
JP6508867B2 (ja) | 2013-03-29 | 2019-05-08 | キヤノン株式会社 | 被検体情報取得装置およびその制御方法 |
US20140350406A1 (en) * | 2013-05-24 | 2014-11-27 | Siemens Medical Solutions Usa, Inc. | Dynamic Operation for Subarrays in Medical Diagnostic Ultrasound Imaging |
JP6111870B2 (ja) * | 2013-05-31 | 2017-04-12 | セイコーエプソン株式会社 | 超音波測定装置、超音波画像装置及び超音波画像の処理方法 |
EP2818117A1 (en) * | 2013-06-26 | 2014-12-31 | Canon Kabushiki Kaisha | Object information obtaining system, signal processing method, and program |
US10090627B2 (en) * | 2013-08-21 | 2018-10-02 | Te Connectivity Corporation | Filters for terminal crimping devices using ultrasonic signals |
JP6180850B2 (ja) * | 2013-08-30 | 2017-08-16 | 富士フイルム株式会社 | 超音波診断装置および超音波画像生成方法 |
JP2015112326A (ja) * | 2013-12-12 | 2015-06-22 | キヤノン株式会社 | プローブ、被検体情報取得装置 |
US10405829B2 (en) | 2014-12-01 | 2019-09-10 | Clarius Mobile Health Corp. | Ultrasound machine having scalable receive beamformer architecture comprising multiple beamformers with common coefficient generator and related methods |
JP2017000547A (ja) * | 2015-06-12 | 2017-01-05 | 株式会社日立製作所 | 超音波診断装置 |
US10079633B2 (en) | 2015-09-29 | 2018-09-18 | The United States Of America, As Represented By The Secretary Of The Army | Time-based and frequency-based radio beamforming waveform transmission |
US10320467B2 (en) | 2015-09-29 | 2019-06-11 | The United States Of America, As Represented By The Secretary Of The Army | Frequency-based radio beamforming waveform transmission |
US10193612B2 (en) * | 2015-09-29 | 2019-01-29 | The United States Of America, As Represented By The Secretary Of The Army | Time-based radio beamforming waveform transmission |
JP6545089B2 (ja) * | 2015-11-27 | 2019-07-17 | エイブリック株式会社 | 超音波診断用送信回路および超音波送信方法 |
CN108267745A (zh) * | 2016-06-22 | 2018-07-10 | 安溪县景宏技术咨询有限公司 | 图像形成装置及其图像形成方法 |
US20180055369A1 (en) * | 2016-08-31 | 2018-03-01 | Qualcomm Incorporated | Layered sensing including rf-acoustic imaging |
US11439369B1 (en) * | 2017-10-26 | 2022-09-13 | United States Of America As Represented By The Secretary Of The Air Force | Method and apparatus for imaging with reduced level of off-axis artifacts |
JP7213491B2 (ja) * | 2019-05-21 | 2023-01-27 | 国立大学法人豊橋技術科学大学 | 超音波画像構築方法、超音波画像構築装置、超音波画像構築プログラム |
CN110693524B (zh) * | 2019-10-16 | 2022-08-05 | 深圳蓝影医学科技股份有限公司 | 一种超声医学成像聚焦校正方法和装置 |
WO2021216723A1 (en) * | 2020-04-22 | 2021-10-28 | Bfly Operations, Inc. | Methods and apparatuses for beamforming in ultrasound systems |
US11927565B2 (en) * | 2021-05-27 | 2024-03-12 | Evident Canada, Inc. | Techniques for adapting time delays of ultrasound inspection system during data acquisition |
CN113598719A (zh) * | 2021-09-01 | 2021-11-05 | 上海科技大学 | 光声重建算法的硬件加速方法、装置及*** |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0696791A2 (en) * | 1994-08-09 | 1996-02-14 | Hewlett-Packard Company | Delay generator for phased array ultrasound beamformer |
US5724972A (en) * | 1996-05-02 | 1998-03-10 | Acuson Corporation | Method and apparatus for distributed focus control with slope tracking |
CN101371789A (zh) * | 2007-08-22 | 2009-02-25 | 深圳迈瑞生物医疗电子股份有限公司 | 实时计算波束合成聚焦延迟参数的方法及装置 |
Family Cites Families (29)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS58144767A (ja) * | 1982-02-24 | 1983-08-29 | Hitachi Medical Corp | 超音波受信装置 |
EP0256481B1 (de) * | 1986-08-20 | 1994-01-12 | Siemens Aktiengesellschaft | Verfahren und Einrichtung zur adaptiven Fokussierung bei einem medizinischen Ultraschall-Bildgabegerät |
US5928152A (en) * | 1994-08-05 | 1999-07-27 | Acuson Corporation | Method and apparatus for a baseband processor of a receive beamformer system |
US5817024A (en) | 1996-06-28 | 1998-10-06 | Sonosight, Inc. | Hand held ultrasonic diagnostic instrument with digital beamformer |
US5713356A (en) | 1996-10-04 | 1998-02-03 | Optosonics, Inc. | Photoacoustic breast scanner |
JP4316700B2 (ja) * | 1998-08-18 | 2009-08-19 | 株式会社日立メディコ | 超音波診断装置 |
US6123671A (en) * | 1998-12-31 | 2000-09-26 | General Electric Company | Method and apparatus for distributed, agile calculation of beamforming time delays and apodization values |
US6685645B1 (en) * | 2001-10-20 | 2004-02-03 | Zonare Medical Systems, Inc. | Broad-beam imaging |
WO2001026555A1 (fr) * | 1999-10-15 | 2001-04-19 | Hitachi Medical Corporation | Dispositif d'imagerie ultrasonore |
JP2003521341A (ja) * | 2000-01-31 | 2003-07-15 | アー.ヤー. アンゲルセン、ビョルン | 医療用超音波イメージングにおける位相面収差およびパルス残響の補正 |
JP4416256B2 (ja) * | 2000-03-10 | 2010-02-17 | 株式会社日立メディコ | 超音波撮像装置 |
US6443897B1 (en) * | 2000-12-28 | 2002-09-03 | Ge Medical Systems Global Technology Company, Llc | Refraction delay error correction using agile beamformer |
US6607489B2 (en) * | 2001-04-05 | 2003-08-19 | General Electric Company | Focus correction for ultrasound imaging through mammography compression plate |
US6705994B2 (en) * | 2002-07-08 | 2004-03-16 | Insightec - Image Guided Treatment Ltd | Tissue inhomogeneity correction in ultrasound imaging |
JP3884370B2 (ja) * | 2002-11-15 | 2007-02-21 | アロカ株式会社 | 超音波診断装置 |
US8088067B2 (en) * | 2002-12-23 | 2012-01-03 | Insightec Ltd. | Tissue aberration corrections in ultrasound therapy |
JP4406226B2 (ja) | 2003-07-02 | 2010-01-27 | 株式会社東芝 | 生体情報映像装置 |
US20050107694A1 (en) * | 2003-11-17 | 2005-05-19 | Jansen Floribertus H. | Method and system for ultrasonic tagging of fluorescence |
CN101112320B (zh) * | 2006-07-28 | 2010-07-21 | 深圳迈瑞生物医疗电子股份有限公司 | 波束合成的接收聚焦参数的实时计算方法及其装置 |
JP5016326B2 (ja) * | 2007-03-06 | 2012-09-05 | 株式会社日立メディコ | 超音波診断装置 |
CA2697343C (en) * | 2007-08-23 | 2017-06-20 | Verasonics, Inc. | Adaptive ultrasound image reconstruction based on sensing of local media motion |
JP5294687B2 (ja) | 2008-05-02 | 2013-09-18 | キヤノン株式会社 | 超音波測定装置およびその制御方法 |
JP5495607B2 (ja) | 2008-05-27 | 2014-05-21 | キヤノン株式会社 | 超音波診断装置 |
US8240597B2 (en) | 2008-08-06 | 2012-08-14 | Honeywell International Inc. | UAV ducted fan lip shaping |
JP5419404B2 (ja) * | 2008-09-04 | 2014-02-19 | キヤノン株式会社 | 光音響装置 |
JP2010068885A (ja) * | 2008-09-17 | 2010-04-02 | Canon Inc | 測定装置 |
JP5400466B2 (ja) * | 2009-05-01 | 2014-01-29 | キヤノン株式会社 | 画像診断装置、画像診断方法 |
JP5393256B2 (ja) | 2009-05-25 | 2014-01-22 | キヤノン株式会社 | 超音波装置 |
CN101858972B (zh) * | 2010-03-23 | 2013-01-30 | 深圳市蓝韵实业有限公司 | 基于延时参数实时计算和流水线的多波束合成方法和装置 |
-
2010
- 2010-02-23 JP JP2010036890A patent/JP5645421B2/ja not_active Expired - Fee Related
-
2011
- 2011-02-22 WO PCT/JP2011/054422 patent/WO2011105593A1/en active Application Filing
- 2011-02-22 US US13/520,410 patent/US9165552B2/en not_active Expired - Fee Related
- 2011-02-22 CN CN201180010435.9A patent/CN102770079B/zh not_active Expired - Fee Related
- 2011-02-22 EP EP11708114.1A patent/EP2538846B1/en not_active Not-in-force
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0696791A2 (en) * | 1994-08-09 | 1996-02-14 | Hewlett-Packard Company | Delay generator for phased array ultrasound beamformer |
US5724972A (en) * | 1996-05-02 | 1998-03-10 | Acuson Corporation | Method and apparatus for distributed focus control with slope tracking |
CN101371789A (zh) * | 2007-08-22 | 2009-02-25 | 深圳迈瑞生物医疗电子股份有限公司 | 实时计算波束合成聚焦延迟参数的方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN102770079A (zh) | 2012-11-07 |
WO2011105593A1 (en) | 2011-09-01 |
US9165552B2 (en) | 2015-10-20 |
EP2538846B1 (en) | 2017-10-11 |
JP2011172611A (ja) | 2011-09-08 |
EP2538846A1 (en) | 2013-01-02 |
JP5645421B2 (ja) | 2014-12-24 |
US20120281902A1 (en) | 2012-11-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102770079B (zh) | 超声波成像装置和控制延迟的方法 | |
Jakovljevic et al. | Local speed of sound estimation in tissue using pulse-echo ultrasound: Model-based approach | |
EP3581961A1 (en) | Method and apparatus for ultrasound imaging with improved beamforming | |
JP6053860B2 (ja) | 医療用超音波イメージングシステムにおいて剪断波情報を求めるためのコンピュータプログラム | |
EP1281074B1 (en) | Correction of phasefront aberrations and pulse reverberations in medical ultrasound imaging | |
US9572549B2 (en) | Calibration of multiple aperture ultrasound probes | |
US5531117A (en) | Closed loop maximum likelihood phase aberration correction in phased-array imaging systems | |
US6705994B2 (en) | Tissue inhomogeneity correction in ultrasound imaging | |
EP2561380B1 (en) | Imaging method and apparatus using shear waves | |
CN111819467B (zh) | 波传播和散射参数的估计方法和仪器 | |
RU2458342C1 (ru) | Способ ультразвуковой томографии и устройство для его осуществления | |
US20180143165A1 (en) | Local Speed Of Sound Estimation Method For Medical Ultrasound | |
JP2014512923A (ja) | オーバーラップする送信ビームにおける適格と評価された領域を使用する増強された超音波画像形成 | |
US11029289B2 (en) | Ultrasonic TFM with calculated angle beams | |
AU2021209270A1 (en) | Method and system for ultrasonic characterization of a medium | |
AU2021209268A1 (en) | Method and system for ultrasonic characterization of a medium | |
RU2232547C2 (ru) | Способ и устройство для получения ультразвуковых изображений структур и сосудов головного мозга | |
US11628470B2 (en) | Method for locating elements on a flexible ultrasound phased array | |
Måsøy et al. | Estimation of ultrasound wave aberration with signals from random scatterers | |
Dziewierz et al. | A design methodology for 2D sparse NDE arrays using an efficient implementation of refracted-ray TFM | |
AU2021209265A1 (en) | Method and system for ultrasonic characterization of a medium | |
Freiburger et al. | Parallel processing techniques for the speckle brightness phase aberration correction algorithm | |
RU2788389C1 (ru) | Способ обнаружения аберраций при ультразвуковом исследовании | |
Lacefield et al. | Spatial coherence analysis applied to aberration correction using a two-dimensional array system | |
JP2669204B2 (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150311 |
|
CF01 | Termination of patent right due to non-payment of annual fee |