CN103237500B - 用于波束形成中变迹的动态孔径控制和归一化 - Google Patents

用于波束形成中变迹的动态孔径控制和归一化 Download PDF

Info

Publication number
CN103237500B
CN103237500B CN201180058809.4A CN201180058809A CN103237500B CN 103237500 B CN103237500 B CN 103237500B CN 201180058809 A CN201180058809 A CN 201180058809A CN 103237500 B CN103237500 B CN 103237500B
Authority
CN
China
Prior art keywords
data
apodization
channel
block
dynamic
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.)
Active
Application number
CN201180058809.4A
Other languages
English (en)
Other versions
CN103237500A (zh
Inventor
D·P·马吉
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.)
Texas Instruments Inc
Original Assignee
Texas Instruments Inc
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 Texas Instruments Inc filed Critical Texas Instruments Inc
Publication of CN103237500A publication Critical patent/CN103237500A/zh
Application granted granted Critical
Publication of CN103237500B publication Critical patent/CN103237500B/zh
Active 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/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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/54Control of the diagnostic device
    • 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
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B8/00Diagnosis using ultrasonic, sonic or infrasonic waves
    • A61B8/48Diagnostic techniques
    • A61B8/486Diagnostic techniques involving arbitrary m-mode
    • 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

Abstract

一种用于超声波束形成的变迹方法(500)包括提供(501)在多个(k)数据通道中表示来自靶组织的回声数据的滤波数字通道数据和预定数量的有效通道(Nact)。基于软件的集成变迹算法利用选择的变迹函数h[n,k]动态变迹数字通道数据,其中n是样本数。该集成变迹算法通过产生参数numApertureChannels(Nap[n])其中Nap[n]<Nact且基于用于波束形成的动态波束聚焦位置选择Nap[n]个特定数字通道,来应用动态孔径控制创建有效的孔径。应用的动态数据缩放通过利用h[n,k]和比例因子之间的向量内积产生归一化的变迹因子hnorm[n,k],来提供数据归一化。归一化的变迹因子被应用(503)于Nap[n]个特定数字通道中的数字通道数据以产生经变迹和归一化的数字通道数据。

Description

用于波束形成中变迹的动态孔径控制和归一化
技术领域
公开的实施例涉及用于超声波和相关波束形成算法的波束形成器、实现公开的波束形成算法的集成电路(IC)及其超声系统。
背景技术
波束形成(beamforming)是在用于指向性信号发射或接收的传感器阵列中使用的信号处理技术。通过使用自适应的或固定的接收/发射波束图形获得空间选择性。
使用波束形成的一种医学应用是超声诊断学。超声波能量由发射波束形成器聚焦在靶组织处,并且靶组织调制和返回的超声波能量由接收波束形成器聚焦。该接收波束形成器可以提供信号,以产生B超图像、表示靶组织的彩色多普勒或频谱多普勒信息或其组合。这种波束形成系统可以提供人体组织或其他对象的组织的实时横断面(层析成像)二维(2D)图像。
在传统的超声成像系统中,定制的专用集成电路(ASIC)器件执行波束形成器计算。作为软件波束形成解决方案的一部分,变迹(apodization)涉及对每个数据通道(每个数据通道包括专用的换能器元件)的适当增益应用,并且需要使每个数据通道内的输入信号成形。这种成形减少波束形成信号中的栅旁瓣效应,其中该栅旁瓣效应源于侧压幅度变化和换能器元件间隔。
在信号处理数据路径中期望的其他特性包括动态孔径控制和数据归一化(normalization)。动态孔径控制的功能是通过动态调整有效换能器元件的数量来保持遍及整体超声扫描线的期望孔径。关于归一化,由于在每个样本数n(每个样本数n对应于不同的时刻)可以使用不同数量的换能器元件,归一化可以考虑这种缩放效应从而为波束形成信号维持恒定信号电平。
在大部分ASIC波束形成器系统设计中,经由使数据通道特定组合的选择作为预期孔径的函数的模拟开关,处理孔径控制。图1示意性示出用于成像靶组织的常规超声波束形成器系统100的简化结构图绘图,其包括来自从靶返回的回声的接收回声信号14的基于切换的时间延迟转向和聚焦。系统100包括包含多个换能器元件11的接收阵列10,如图所示其配置成线性阵列。
在图1中动态孔径控制示意性示出为多对接收器通道开关16a-16e,这些通道开关在回声接收时段期间被孔径控制电路17按顺序关闭。孔径控制电路17在任意给定的采样时间(时刻)动态选择换能器元件11中的哪些是有效换能器元件。对于处于接收模式的每个有效换能器元件11,建立数据通道以将靶数据从换能器元件11转移到波束形成器中的信号处理元件。接收通道开关16a-16e是模拟电子开关。数据通道与开关16a、16b、16c、16d和16e中的每一个关联。
在回声接收时段期间,来自每个数据通道内的相应有效换能器元件11的数据被延迟相应的处理通道时间延迟,包括图1所示的转向延迟和聚焦延迟,其起到动态聚焦接收的信号回声的作用。尽管没有在图1中示出,在通道时间延迟的应用之后,相应的延迟通道信号被变迹以便采用变迹因子(apodization factor)使每个接收的通道信号缩放期望的值,从而减少如上所述的在随后形成的波束形成信号中的栅旁瓣效应。然后通过加法放大器15对变迹信号求和以产生图1中示为“聚焦回声信号”的波束形成信号。多个聚焦回声信号可以用于形成扫描线,且多个扫描线可以被组合以在适合的显示装置上形成靶组织的图像。
尽管没有示出,加法放大器15输出的在每个采样时间n的聚焦回声信号然后可以被最优归一化。归一化可以用于使平均信号水平保持在采样时刻n中在整个扫描线上基本相同(恒定)。
由常规超声波束形成器系统100实现的动态孔径控制类型可能难以根据有效通道的数量和采样时间进行同步,尤其如果该实现是模拟的。由于这种同步困难,作为采样时间的函数确定对于使平均信号水平保持沿整个扫描线相同(恒定)所需的数据归一化通常变得存在更多问题。
发明内容
公开的方法包括用于超声波束形成的变迹方法、由此形成的超声成像系统和诸如用于实现超声成像块和功能的数字信号处理器(DSP)的硬件,其包含滤波块、变迹块、动态孔径控制块、控制器和加法器。提供在多个(k)数据通道中表示来自靶组织的回声数据的数字通道数据和从该多个数据通道选出的预定数量(≤k)的有效通道(Nact)。
基于软件的集成变迹算法利用选择的变迹函数h[n,k]动态变迹数字通道数据。通过产生参数numApertureChannels(Nap[n])其中Nap[n]≤Nact且基于用于超声波束形成的动态波束聚焦位置选择Nap[n]个特定数据通道,该集成变迹算法应用动态孔径控制以创建有效的孔径。通过利用h[n,k]和比例因子之间的向量内积产生归一化的变迹因子hnorm[n,k],应用的动态数据缩放提供数据归一化。
所述归一化变迹因子可以仅应用于Nap[n]个特定数据通道中数字通道数据,以产生经变迹和归一化的数字通道数据。在动态变迹之前或者在其之后,滤波也应用于Nap[n]个特定数据通道中的通道数据。因为变迹和滤波可以仅应用于Nap[n]个特定数据通道中的数据,其中该Nap[n]个特定数据通道是用于在每个样本数n(时刻)形成孔径,所以与变迹和滤波来自该系统中所有数据通道的数据这种布置相反,数据处理效率显著提高(例如更少的循环)。
公开的实施例通过提供基于软件的方法解决对集成变迹方法的未满足的需要,其中该基于软件的方法执行动态孔径控制和归一化两方面以便为超声波束形成器系统产生缩放的波束形成数据。公开的实施例能够在无数据通道切换的情况下变迹,这是因为使用向量内积处理所述归一化/缩放,由此避免了在数据通道之间进行切换的需要。
因为在每个样本数(n)的孔径通道和数据归一化之间的关系是已知的,所以部分提供更准确计时的集成变迹算法。以软件为基础,公开的集成变迹算法可以用于任意期望的变迹函数。由于公开的集成变迹方法涉及较少循环,他们也足以在硬件上实现,以产生等价的数值结果(相比于更多乘法密集型的传统算法)。公开的集成变迹算法可以用于为整个超声图像中给定(单个)扫描线或者所有扫描线执行数据归一化。
附图说明
图1示意性示出用于成像靶组织的常规超声波束形成器系统的简化结构图绘图,其包括对来自换能器元件接收阵列的接收回声信号的基于切换的时间延迟转向和聚焦。
图2根据示例实施例示意性示出用于成像靶组织的超声波束形成器系统的简化结构图绘图,其包括公开的动态孔径控制和归一化。
图3根据本发明的实施例示出可以实现图2中虚线范围内的所有系统元件的DSP IC的简化结构图。
图4是根据示例实施例的示例性超声系统的结构图,该系统可以实现公开的动态孔径控制和归一化。
图5是根据示例实施例的变迹超声信号的方法的流程图,该方法包括动态孔径控制和归一化。
图6提供根据示例实施例在德克萨斯仪器股份有限公司的C64+体系结构DSP上运行的基于汉宁窗的集成变迹算法的软件流水信息。
具体实施方式
如上面所述,常规的超声系统例如图1中所示常规超声波束形成器系统100使用模拟通道开关执行动态孔径控制、对通道数据执行变迹(如果实现)以及对求和通道数据执行归一化(如果实现),对求和通道数据的归一化用于缩放平均信号水平以使其在波束形成器系统内分开位置(例如,作为独立的数据处理块)中在整个扫描线上相同(恒定)。信号处理的某些部分可以以模拟方式执行(例如,孔径控制17)且某些部分可以以数字方式执行(变迹,如果实现的话)。相比之下,公开的集成变迹方法提供基于软件的方法,该方法执行动态孔径控制和归一化,由此将前述独立的数据处理操作结合成单个操作以确保数据缩放准确度,也消除了开关的需要。
公开的基于软件的集成变迹算法可实现在诸如ASIC、现场可编程门阵列(FPGA)和DSP的IC上。由于如下所述用向量内积处理缩放,因此消除对开关的需要。同样,由于在每个样本数n(时刻)上用于动态孔径归一化的相应数据通道之间的关系是已知的,因此该数据可以在每个样本数n上准确缩放。
为了更加准确和确定地实现动态孔径控制和数据归一化两方面,下面公开一种示例性的基于软件的集成变迹方法。通常,变迹处理的离散时间实现方式可以在数学上表示为:
y[n,k]=h[n,k]*x[n,k]                 (1)
其中h[n,k]是关于第k个数据通道和样本数n的变迹函数,x[n,k]是以样本数n从第k个数据通道接收的信号,而y[n,k]是以样本数n来自第k个数据通道的经滤波的接收信号。对于简单的单位增益变迹函数h[n,k],该变迹函数h[n,k]可以定义为:
h [ n , k ] = 0 k &NotElement; k ap [ n ] , n &GreaterEqual; 0 1 k &Element; k ap [ n ] , n &GreaterEqual; 0 - - - ( 2 )
其中Kap[n]将孔径通道指数集合定义为样本数n的函数,其中每个n对应于具体的时刻。对个更复杂的变迹函数h[n,k]而言,这种特定的高斯窗函数在本领域中称为汉宁窗,该变迹函数变成下面的表示:
h [ n , k ] = 0 k &NotElement; k ap [ n ] , n &GreaterEqual; 0 1 2 ( 1 - cos ( 2 &pi; ( k + 1 ) N ap [ n ] + 1 ) ) k &Element; k ap [ n ] , n &GreaterEqual; 0 - - - ( 3 )
其中cos()是余弦函数且Nap[n]将孔径通道(有效数据通道)的数量(N)定义为样本数n的函数。
为了确定给定样本数n(时刻)的有效孔径通道的数量,此处称为Nact,按照换能器参数定义孔径的f#(即f数量)。Nact通常是用户指定的且涉及选择用于波束形成的数据通道(每个数据通道与换能器关联)的最大数量,其中Nact小于或等于波束形成系统内的数据通道/换能器元件的数量。在光学中,f#被定义为透镜的焦点距离(Rfp)与透镜直径(D)的比率,可以数学表示为:
f # = R fp D - - - ( 4 )
在大多数的应用中,当f#等于1时获得良好的聚焦。然而,公开的实施例包括f#不等于1。
f#可以按照超声换能器参数表示为:
f # = R fp N ap [ n ] * elemPitch - - - ( 5 )
其中Rfp是到透镜焦点的径向距离,而elemPitch是各换能器元件之间的间距。等式(5)可以写成考虑复杂度增加的非恒定元件间距。由于波束形成器的目标是动态调整离散时间系统在每个样本数n的焦点,因此到透镜的焦点距离Rfp可以写成:
R fp = c 2 &CenterDot; f s &CenterDot; n - - - ( 6 )
其中c是在感兴趣组织中的声速,且fs是离散时间信号处理系统的采样率。把等式(6)代入等式(5),推导出以下描述数据通道的数量Nap[n]的关系,其定义了孔径通道的数量并且选自有效通道Nact,该有效通道Nact是用于创建在给定样本数n的期望孔径以获得良好的聚焦所需要的。
N ap [ n ] = c 2 &CenterDot; f s &CenterDot; f # &CenterDot; elemPitch &CenterDot; n - - - ( 7 )
关于Nap[n]的等式(7)能够确定用于创建在每个样本数n的孔径所需要的指数集合(即Kap[n]),以便选择的变迹函数(例如高斯窗函数诸如汉宁或汉明变迹窗)仅可以应用于数据通道的Nap[n]个特定数据通道,由此节省数据处理资源。Kap[n]在上述等式2中示出。
例如,相控阵列换能器的波束聚焦位置被设置在构成换能器的换能器元件的中间(中点)。因此,通过使用以换能器中点为中心的Nap[n]个通道来创建期望的孔径。然而,对于线性阵列换能器而言,线性阵列换能器的波束聚焦位置沿换能器的表面移动。因此,通过使用围绕动态波束聚焦位置的Nap[n]个通道创建期望的孔径。
一旦变迹已应用于通道数据,经变迹的通道数据可以被归一化以保持在样本数n范围沿整个扫描线的平均信号水平相同(有效恒定)。用于实现归一化的比例因子Krat[n]在本文称为归一化比例因子Krat[n],其可以写成:
K rat [ n ] = N act N ap [ n ] - - - ( 8 )
其中,上述的Nact是有效换能器通道的数量。应注意,关于比例因子Krat[n]最大值的值应被指定以处理1除以0的情况且处理给定超声图像中的所有扫描线的归一化。由此,归一化比例因子Krat[n]可能值的范围可以表示为:
1≤Krat[n]≤Kmax                   (9)
使用(8,9)中的归一化比例因子Krat[n]和(7)中关于孔径通道数量Nap[n]的关系,产生如下(10)显示的广义内积关系并使用该内积关系计算关于任意选择的变迹函数h[n,k]的归一化变迹实现hnorm[n,k]:
hnorm[n,k]=Krat[n]·h[n,k]其中0≤k≤N              (10)
其中h[n,k]是期望的变迹函数,且hnorm[n,k]是期望的变迹函数的归一化形式。
图2根据示例实施例示意性示出用于成像靶组织的超声波束形成器系统200的简化结构图绘图,其包括公开的动态孔径控制和归一化。在典型的实施例中,波束形成器系统200是延迟和求和超声波束形成器系统。
尽管图2示出以线性方式布置的换能器元件,正如在某些相控阵扫描系统中一样,但是公开的实施例也可以用于以其他配置安置的换能器元件,例如弯曲布置、柱状布置以及2D(即二维)和3D(即三维)布置。另外,所公开的实施例也可应用于其他阵列类型的换能器,包括曲线或弯曲的换能器、柱形换能器、环形阵列换能器和矩形阵列换能器。
波束形成器系统200包括换能器阵列112,其包含如所示元件1121-1128的多个(k)换能器元件,每个元件定义的数据通道包含将靶组织105反射的声波转换成电子传感信号的压电换能器。换能器元件1121-1128由此形成8个独立的数据通道1到8(每个数据通道一个换能器)中的第一元件,该8个独立的数据通道从换能器元件1121-1128延伸到下面所述的numActiveChannels(Nact)块120a。
在换能器阵列112中的各换能器均各自耦合到电压控制放大器(VCA)116,然后耦合到用于将放大的换能器信号数字转换的模数转换器(ADC)117。图2中的ADC117耦合到示为动态孔径控制块120的Nact的块120a上,该动态孔径控制块120显示为由所示控制器241控制。Nact设置波束形成器系统200使用的k数据通道的最大数量。该最大数量通常是预定的且保持恒定,诸如由用户输入提供并存储在存储器中。Nact具有以下的限制:
0≤Nact≤numChannels(如图2所示numChannels=8)
图2中显示的Nact块120a在该图2的示例中从图2所示的可能的8个数据通道中选择出6个。动态孔径控制块120也包括numApertureChannels(Nap)计算块120b。上述的动态Nap[n]计算如所示在控制器241的控制下由Nap计算块120b实现。Nap[n]计算块120b可操作以动态确定用于创建有效孔径的特定数据通道。
Nap[n]具有以下限制:
0≤Nap≤Nact
由Nap计算块120b提供的输出被耦合到包含多个滤波器的滤波器块118上,这些滤波器应用适当的时间延迟以补偿由于靶组织105和各换能器1121-1128之间路径长度差异导致的不同回声到达时间。滤波器块118中的滤波器耦合到变迹增益块119中的变迹块,这些变迹块为基于用于波束形成的动态波束聚焦位置从多个数据通道选择的Nap[n]个特定数据通道中的数据产生经变迹和归一化的通道数据。变迹增益块119由此提供对提供的各数据通道信号的缩放(即加权)。由变迹增益块119提供到各自Nap[n]个特定数据通道的变迹增益是基于上述期望的变迹函数hnorm[n,k]的归一化形式,由此在逐样本(n个接n个)的基础上被改变。
各变迹增益块119的输出耦合到加法器121上,该加法器对变迹增益块119输出的经变迹和归一化的信号求和,以产生示为z[n]的表示聚焦和归一化回声信号的期望波束形成信号。多个聚焦和归一化回声信号z[n]可以用于形成扫描线,且多个扫描线可以被组合以在合适的显示装置上形成靶组织105的图像。
系统200还包括包含计算结构246的控制器241。控制器241可以由包括一个或多个DSP、FPGA或ASIC的装置提供。例如,DPS可以提供图2中所示的虚线范围内的系统200的所有部件。系统200包括耦合到控制器241的存储器248,该存储器用于存储包含所公开的集成变迹算法和波束形成期间样本计数(数量)数据的数据。控制器241可操作以从存储器248加载数据用于处理。
可以看到控制器241被耦合以控制滤波器块118和变迹块119两者。提供给变迹增益块119中各数据通道的加权因子基于上面如(10)提到的关于选择变迹函数h[n,k]的归一化变迹实现hnorm[n,k],下面再重复一遍:
hnorm[n,k]=Krat[n]·h[n,k]其中0≤k≤N               (10)
其中h[n,k]是期望的变迹函数,而hnorm[n,k]是期望的变迹函数的归一化形式。
图3根据示例实施例示出可以实现图2中所示虚线范围内所有系统元件的DSP IC300的简化结构图。所示DSP IC300形成在具有顶部半导体表面的衬底310(例如硅衬底)上并且包括乘法累加(MAC)单元320,该乘法累加单元可操作产生控制信号以实现公开的动态孔径控制和归一化。DSP IC300通常包括易失性存储器(例如RAM)325和非易失性存储器(例如ROM)330。根据本发明实施例的算法可以存储在非易失性存储器330内。所示的DSP IC300也包括用于输入和输出的(多个)接口端口340、计数器/计时器345、存储器控制器350和总线355。
如同常规DSP一样,该DSP IC300可以执行指令以实现一个或者多个数字信号处理算法或处理。举例来说,该指令数据可以包括当载入和初始化到DSP IC300中时可以提示DSP IC300实现不同数字信号处理算法或处理的多种系数以及指令。该DSP IC300可以接收来自图2中所示ADC117的数据,并且然后根据现有配置对该数据应用算法。
MAC单元320通常包括延迟与变迹电路,用于处理数字化超声传感信号以形成经延迟、经变迹和经归一化的数字超声传感信号。MAC单元320也可以提供控制器241和计算结构246。易失性存储器325可以提供需要的任何存储器。
图4是根据示例实施例可以实现公开的动态孔径控制和归一化的示例性超声系统400的结构图。系统400包括包含发射(Tx)波束形成器525的发射部分520和包含接收(Rx)波束形成器545的接收部分540,它们共享换能器550的共同阵列。
系统400包括耦合到Tx波束形成器525和Rx波束形成器545的波束形成器中央控制单元510。波束形成器中央控制单元510可以具体为DSP,诸如上述参考图3的DSP IC300,用于实现上述由图2中所示超声波束形成器系统200提供的所公开动态孔径控制和归一化。接收部分540的Rx波束形成器545耦合到后端成像DSP560。后端成像DSP560耦合到显示器570。
图5是根据示例实施例用于超声波束形成的变迹从多个数据通道(k)接收的超声信号的方法500的流程图,每个所述数据通道与多个换能器元件的不同换能器元件关联。步骤501包括提供在多个数据通道的每一个数据通道中表示从靶组织接收的回声数据的数字通道数据以及来自多个数据通道的预定数量的有效数据通道(Nact)。
步骤502包括应用基于软件的集成变迹算法,以便通过使用选择的变迹函数h[n,k]来动态变迹所述数字通道数据,其中n是样本数,其通过使用计算系统(诸如图2中所示包含计算结构246的控制器241)来实现公开的集成变迹算法,其中所示计算系统包括非暂时性机器可读存储装置(诸如图2中所示的存储器248)。该集成变迹算法包括(i)应用动态孔径控制创建有效孔径,其包括产生参数numApertureChannels(Nap[n])其中Nap[n]≤Nact,和基于用于波束形成的动态波束聚焦位置从多个数据通道中选择Nap[n]个特定数据通道,和(ii)应用动态数据缩放来提供数据归一化,其通过使用变迹函数h[n,k]和比例因子之间的向量内积关系产生归一化变迹因子hnorm[n,k]。
步骤503包括将归一化变迹因子hnorm[n,k]应用到数字通道数据,从而在至少一个样本数[n]时为由Nap[n]个特定数据通道定义的数据通道产生经变迹和归一化的数字通道数据。方法500还包括应用适当的时间延迟到相应数据通道以补偿不同回声到达时间的滤波。滤波可以在数据处理流程中的变迹之前或之后实现。尽管数学上的滤波可以跟在变迹之后,事实上相比于在变迹之前滤波的数据处理流程(例如图2中所示的系统200),在变迹之后的滤波会导致变迹更多的数据值(例如要求更多的乘法)。来自Nap[n]个特定数据通道的经变迹和归一化的滤波数字通道数据然后可以用于产生表示聚焦与归一化回声信号的期望波束形成信号(图2中所示z[n]),其中该Nap[n]个特定数据通道是基于用于波束形成的动态波束聚焦位置选择的。
公开的包括动态孔径控制和归一化的集成变迹算法已经用C代码实现并且已经在由德克萨斯仪器股份有限公司制造的TMS320C6474DSP上示范。汉宁变迹在目标为TI DSP C64+体系结构的软件上实现。图6中示出在C64+体系结构上运行的关于集成变迹算法的软件流水信息。
发现该特定集成变迹算法花费15个循环来计算每个变迹值,用以实现动态孔径控制和归一化两方面。应注意到,实现公开的集成变迹算法所涉及的体系结构不是有限的乘法器(M单元),而作为替代是有限的移位/装载(L,S,D单元)。形成有限的L,S,D单元的优势在于数据处理是较少数据密集型。公开的算法可以由离线计算的变迹表支持并且存储在表中,并且被参考以将流水线调度减少到所示的11个循环的理想范围。
公开的动态孔径控制(因为对于发射波束形成而言不需要归一化)可以用于发射超声波束形成。另外,所公开的波束形成可以被扩展至超出超声应用范围。本领域的技术人员应认识到在要求保护的发明的范围内,可以对示例实施例和实现的许多其他实施例作出其他改进。

Claims (13)

1.一种用于超声波束形成的变迹从多个数据通道接收的超声信号的方法,其中所述多个数据通道即k个数据通道,所述多个数据通道中的每个与多个换能器元件的不同换能器元件关联,所述方法包括:
在所述多个数据通道的每一个数据通道中提供表示从靶组织接收的回声数据的数字通道数据和提供从所述多个数据通道中选择的预定数量的有效数据通道,即Nact
应用基于软件的集成变迹算法以使用选择的变迹函数h[n,k]来动态变迹所述数字通道数据,其中n是样本数,使用包含非暂时性机器可读存储装置的计算系统来实现所述集成变迹算法,所述集成变迹算法包括:
(i)应用动态孔径控制创建有效孔径包括:产生参数numApertureChannels,即Nap[n],其中Nap[n]≤Nact,和基于用于所述超声波束形成的动态波束聚焦位置从所述有效数据通道中选择Nap[n]个特定数据通道,
(ii)通过使用所述变迹函数h[n,k]和比例因子之间的向量内积关系产生归一化变迹因子hnorm[n,k],应用动态数据缩放提供数据归一化,和
将所述归一化变迹因子hnorm[n,k]应用到所述数字通道数据,从而在至少一个所述n时为由所述Nap[n]个特定数据通道定义的所述数据通道产生经变迹和归一化的数字通道数据,
其中用Krat[n]表示所述比例因子,且所述向量内积关系包括:hnorm[n,k]=Krat[n]·h[n,k]其中0<k≤Nact
其中 K rat [ n ] = N act N ap [ n ] .
2.根据权利要求1所述的方法,其中所述归一化变迹因子hnorm[n,k]仅专门提供给所述多个数据通道的所述Nap[n]个特定数据通道。
3.根据权利要求1所述的方法,其中所述选择的变迹函数h[n,k]包括高斯窗函数。
4.根据权利要求2所述的方法,进一步包括:
仅为所述多个数据通道的所述Nap[n]个特定数据通道专门滤波所述数字通道数据,以便所述经变迹和归一化的数字通道数据包括经变迹和归一化的滤波数字通道数据;
对所述经变迹和归一化的滤波数字通道数据求和以产生波束形成信号;
沿直线改变所述动态波束聚焦位置;
多次重复包括所述求和以及所述改变的所述方法以产生多个所述波束形成信号,和
从所述多个波束形成信号产生扫描线。
5.根据权利要求4所述的方法,进一步包括产生多个所述扫描线,和结合所述多个扫描线以在显示装置上形成所述靶组织的图像。
6.根据权利要求1所述的方法,进一步包括使用数字信号处理器即DSP,用于仅为所述多个数据通道的所述Nap[n]个特定数据通道专门滤波所述数字通道数据,以便所述经变迹和归一化的数字通道数据包括经变迹和归一化的滤波数字通道数据,并且用于应用所述基于软件的集成变迹算法以动态变迹所述滤波数字通道数据。
7.一种用于超声波束形成的超声成像系统,其包括:
多个换能器元件,用于向靶组织发射超声发射脉冲,并且接收来自所述靶组织的回声信号作为响应,其中所述多个换能器元件即k个换能器元件;
发射部分,其用于为所述超声发射脉冲的所述发射驱动所述多个换能器元件,以及
接收部分,其用于处理由所述多个换能器响应于所述回声信号产生的多个传感信号,所述接收部分定义多个数据通道,所述多个数据通道中的每个与所述多个换能器元件的不同换能器元件关联,所述接收部分包括:
数字化块,其用于在所述多个数据通道的每个数据通道中产生表示来自所述回声信号的回声数据的数字通道数据,
动态孔径控制块,其耦合到所述数字化块的输出端,且包含用于设置所述系统采用的所述多个数据通道的数量的有效数据通道Nact块和用于动态确定用于创建有效孔径的所述多个数据通道的数量Nap[n]的动态numApertureChannels计算块,即Nap[n]计算块;
滤波块,其耦合到所述Nap[n]计算块的输出端,用以滤波所述数字通道数据从而产生滤波数字通道数据,
耦合到所述滤波块的变迹块;
控制器和相关联的存储器,它们经耦合以控制所述滤波块、所述变迹块和所述动态孔径控制块;
其中所述控制器执行基于软件的集成变迹算法,以利用选择的变迹函数h[n,k]动态变迹所述滤波数字通道数据,其中n是样本数,其中包含非暂时性机器可读存储装置的计算系统用于实现所述集成变迹算法,所述集成变迹算法包括:
(i)动态孔径控制,用于控制所述动态孔径控制块以创建所述有效孔径,且包括基于用于所述波束形成的动态波束聚焦位置从所述有效数据通道中选择Nap[n]个特定数据通道,以及
(ii)通过使用所述变迹函数h[n,k]与比例因子Krat[n]之间的向量内积关系产生归一化变迹因子hnorm[n,k],动态数据缩放以提供数据归一化;
其中所述控制器将所述归一化变迹因子hnorm[n,k]提供给所述变迹块以处理所述滤波数字通道数据从而在至少一个所述n时产生经变迹和归一化的滤波数字通道数据;
用于求和的加法器,其耦合到所述变迹块的输出端以形成波束形成信号;
后端成像显示处理器,其经耦合以接收和处理所述波束形成信号从而产生显示信号,所述显示信号适用于使显示装置产生图像,以及
显示装置,其用于接收所述显示信号并产生所述图像,
其中由Krat[n]表示所述比例因子且所述向量内积关系包括:
hnorm[n,k]=Krat[n]·h[n,k]其中0<k≤Nact,其中
8.根据权利要求7所述的超声成像系统,其中所述归一化变迹因子hnorm[n,k]仅专门提供给所述多个数据通道的所述Nap[n]个特定数据通道。
9.根据权利要求7所述的超声成像系统,其中所述系统可操作用于沿直线改变所述动态波束聚焦位置,且多次重复所述求和与所述改变以产生多个所述波束形成信号。
10.根据权利要求7所述的超声成像系统,进一步包括使用数字信号处理器(DSP)实现所述滤波块、所述变迹块、所述动态孔径控制块、所述控制器和所述加法器。
11.一种用于超声信号处理的数字信号处理器集成电路即DSP IC,其包括:
衬底,其具有半导体表面;
耦合到数字化块的输出端的动态孔径控制块,所述数字化块在多个数据通道的每个数据通道中提供表示来自回声信号的回声数据的数字通道数据,所述多个数据通道的每个与多个换能器元件的不同换能器元件相关联,其中所述多个数据通道即k个数据通道;
所述动态孔径控制块包括用于设置采用的所述多个数据通道的最大数量的有效数据通道Nact块和用于动态确定所述多个数据通道的数量Nap[n]从而创建有效孔径的动态numApertureChannels计算块,即Nap[n]计算块;
滤波块,其耦合到所述Nap[n]计算块的输出端以滤波所述数字通道数据从而产生滤波数字通道数据;
耦合到所述滤波块的变迹块;
控制器和相关联的存储器,它们经耦合以控制所述滤波块、所述变迹块和所述动态孔径控制与归一化块;
其中所述控制器执行基于软件的集成变迹算法以利用选择的变迹函数h[n,k]动态变迹所述滤波数字通道数据,其中n是样本数,实现所述集成变迹算法包括:
(i)动态孔径控制,其用于控制所述动态孔径控制与归一化块以创建所述有效孔径,且包括基于用于波束形成的动态波束聚焦位置从所述有效数据通道中选择Nap[n]个特定数据通道,以及
(ii)通过使用所述变迹函数h[n,k]与比例因子Krat[n]之间的向量内积关系产生归一化变迹因子hnorm[n,k],动态数据缩放以提供数据归一化;
其中所述控制器为所述变迹块提供所述归一化变迹因子hnorm[n,k]以处理所述滤波数字通道数据从而在至少一个所述n时产生经变迹和归一化的滤波数字通道数据,以及
用于求和的加法器,其耦合到所述变迹块的输出端以形成波束形成信号,
其中用Krat[n]表示所述比例因子,且所述向量内积关系包括:hnorm[n,k]=Krat[n]·h[n,k]其中0<k≤Nact,其中
12.根据权利要求11所述的DSP IC,其中所述归一化变迹因子hnorm[n,k]仅专门提供给所述多个数据通道的所述Nap[n]个特定数据通道。
13.根据权利要求11所述的DSP IC,其中所述选择的变迹函数h[n,k]包括高斯窗函数。
CN201180058809.4A 2010-12-06 2011-12-06 用于波束形成中变迹的动态孔径控制和归一化 Active CN103237500B (zh)

Applications Claiming Priority (5)

Application Number Priority Date Filing Date Title
US42006810P 2010-12-06 2010-12-06
US61/420,068 2010-12-06
US13/305,378 2011-11-28
US13/305,378 US8545406B2 (en) 2010-12-06 2011-11-28 Dynamic aperture control and normalization for apodization in beamforming
PCT/US2011/063504 WO2012078610A2 (en) 2010-12-06 2011-12-06 Dynamic aperture control and normalization for apodization in beamforming

Publications (2)

Publication Number Publication Date
CN103237500A CN103237500A (zh) 2013-08-07
CN103237500B true CN103237500B (zh) 2015-06-10

Family

ID=46162874

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201180058809.4A Active CN103237500B (zh) 2010-12-06 2011-12-06 用于波束形成中变迹的动态孔径控制和归一化

Country Status (4)

Country Link
US (1) US8545406B2 (zh)
JP (1) JP5911506B2 (zh)
CN (1) CN103237500B (zh)
WO (1) WO2012078610A2 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2011095896A1 (en) 2010-02-08 2011-08-11 Dalhousie University Ultrasound imaging system using beamforming techniques for phase coherence grating lobe suppression
US9739875B2 (en) 2013-02-12 2017-08-22 Urs-Us Medical Technology Inc. Analog store digital read ultrasound beamforming system and method
KR20160046669A (ko) * 2014-10-21 2016-04-29 알피니언메디칼시스템 주식회사 빔포밍 장치, 초음파 이미징 장치 및 빔포밍 방법
CN104914440B (zh) * 2015-06-23 2017-06-16 重庆大学 一种融合动态孔径的压缩感知超声波束合成方法
EP3389500A4 (en) * 2015-12-18 2019-01-02 Ursus Medical LLC Ultrasound beamforming system and method with reconfigurable aperture
US11547389B2 (en) * 2017-09-22 2023-01-10 Philips Image Guided Therapy Corporation Methods and systems for ultrasound contrast enhancement
US20190361102A1 (en) * 2018-05-22 2019-11-28 Analog Devices, Inc. Delay and apodization control interface for ultrasound beamformer
CN109223035B (zh) * 2018-08-21 2021-09-28 青岛海信医疗设备股份有限公司 超声信号处理方法、装置、设备及存储介质
CN109558543A (zh) * 2018-12-11 2019-04-02 拉扎斯网络科技(上海)有限公司 一种样本采样方法、样本采样装置、服务器和存储介质
EP4228515A4 (en) * 2020-11-18 2024-04-03 Wuhan United Imaging Healthcare Co Ltd ULTRASONIC IMAGING SYSTEM AND METHOD
WO2022105843A1 (zh) * 2020-11-18 2022-05-27 武汉联影医疗科技有限公司 一种超声波成像方法、系统和存储介质
US11504093B2 (en) 2021-01-22 2022-11-22 Exo Imaging, Inc. Equalization for matrix based line imagers for ultrasound imaging systems
WO2023097538A1 (en) * 2021-12-01 2023-06-08 Qualcomm Incorporated Techniques for performing uniform beamforming
WO2024054589A1 (en) 2022-09-09 2024-03-14 Exo Imaging, Inc. Coherent matrix of digital imaging systems on chip

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5647365A (en) * 1995-11-13 1997-07-15 Siemens Medical Systems, Inc. Apodization parameter generator for ultrasound medical imaging system
US6436044B1 (en) * 2000-02-14 2002-08-20 Siemens Medical Solutions Usa, Inc. System and method for adaptive beamformer apodization
CN1431624A (zh) * 2003-02-28 2003-07-23 清华大学 基于△-σ变换的超声动态接收变迹方法
CN1568230A (zh) * 2001-09-17 2005-01-19 Ge帕拉莱尔设计公司 换能器的频率和振幅旁瓣缩减
CN101373181A (zh) * 2007-08-24 2009-02-25 深圳迈瑞生物医疗电子股份有限公司 实时计算逐点变迹系数的方法及装置

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4180790A (en) 1977-12-27 1979-12-25 General Electric Company Dynamic array aperture and focus control for ultrasonic imaging systems
FR2596269A1 (fr) * 1986-03-28 1987-10-02 Labo Electronique Physique Echographe ultrasonore apodise a barrette lineaire de transducteurs piezoelectriques et procede de realisation d'une telle barrette
US5685308A (en) * 1994-08-05 1997-11-11 Acuson Corporation Method and apparatus for receive beamformer system
US6629929B1 (en) * 2002-11-08 2003-10-07 Koninklijke Philips Electronics N.V. Method and apparatus for automatically setting the transmit aperture and apodization of an ultrasound transducer array
US20060173313A1 (en) 2005-01-27 2006-08-03 Siemens Medical Solutions Usa, Inc. Coherence factor adaptive ultrasound imaging
JP5028041B2 (ja) * 2006-07-14 2012-09-19 親良 炭 ビームフォーミング装置及び方法
CN101116622B (zh) * 2006-08-02 2010-12-01 深圳迈瑞生物医疗电子股份有限公司 波束合成的接收变迹参数的实时计算方法及其装置
KR101055577B1 (ko) * 2007-11-23 2011-08-23 삼성메디슨 주식회사 초음파 시스템
US8834369B2 (en) 2008-06-27 2014-09-16 Texas Instruments Incorporated Receive beamformer for ultrasound

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5647365A (en) * 1995-11-13 1997-07-15 Siemens Medical Systems, Inc. Apodization parameter generator for ultrasound medical imaging system
US6436044B1 (en) * 2000-02-14 2002-08-20 Siemens Medical Solutions Usa, Inc. System and method for adaptive beamformer apodization
CN1568230A (zh) * 2001-09-17 2005-01-19 Ge帕拉莱尔设计公司 换能器的频率和振幅旁瓣缩减
CN1431624A (zh) * 2003-02-28 2003-07-23 清华大学 基于△-σ变换的超声动态接收变迹方法
CN101373181A (zh) * 2007-08-24 2009-02-25 深圳迈瑞生物医疗电子股份有限公司 实时计算逐点变迹系数的方法及装置

Also Published As

Publication number Publication date
US20120143059A1 (en) 2012-06-07
CN103237500A (zh) 2013-08-07
WO2012078610A3 (en) 2013-01-10
US8545406B2 (en) 2013-10-01
JP5911506B2 (ja) 2016-04-27
WO2012078610A2 (en) 2012-06-14
JP2013544622A (ja) 2013-12-19

Similar Documents

Publication Publication Date Title
CN103237500B (zh) 用于波束形成中变迹的动态孔径控制和归一化
US8416643B2 (en) Receive beamformer for ultrasound having delay value sorting
US8834369B2 (en) Receive beamformer for ultrasound
JP5813776B2 (ja) 超音波撮像装置
CN110063749B (zh) 超声波测定装置、超声波图像装置及超声波测定方法
US5844139A (en) Method and apparatus for providing dynamically variable time delays for ultrasound beamformer
EP0467690B1 (en) Ultrasonic imaging by radial scan of trapezoidal sector
US6910380B2 (en) Ultrasonic transmitting and receiving apparatus
CN109363714B (zh) 一种超声成像设备及其超声成像方法
JP5575907B2 (ja) 超音波プローブ及び超音波診断装置
JP3803374B2 (ja) 位相偏差補正のための2dアレーの作動方法及び接続装置
JP6416884B2 (ja) ダイナミックフォーカスを有する超音波プローブならびに関連するシステムおよび方法
JP2009536853A (ja) マルチラインビーム生成器による超音波合成送信フォーカシング
JP2000157548A (ja) 超音波散乱体をイメ―ジングするための方法及びシステム
US20180299537A1 (en) Beamforming apparatus, beamforming method, and ultrasonic imaging apparatus
US10098613B2 (en) Image processing module, ultrasound imaging apparatus, image processing method, and control method of ultrasound imaging apparatus
US8926514B2 (en) Iterative time delay values for ultrasound beamforming
EP2944976B1 (en) Beam forming apparatus, method for forming beams, ultrasonic imaging apparatus, and ultrasonic probe
US20070083109A1 (en) Adaptive line synthesis for ultrasound
JP2002058671A (ja) 超音波診断装置
JP4130604B2 (ja) 超音波診断装置
JP2001104307A (ja) 柔軟性のある受信ビームフォーミング用集積回路アーキテクチャ
CN117084717A (zh) 一种血流成像方法和超声成像装置
PAUL SIMULATION AND ANALYSIS OF A DIGITAL FOCUSED BEAMFORMER FOR SONAR (SIGNAL PROCESSING)

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