CN101116622A - 波束合成的接收变迹参数的实时计算方法及其装置 - Google Patents

波束合成的接收变迹参数的实时计算方法及其装置 Download PDF

Info

Publication number
CN101116622A
CN101116622A CNA2006100619695A CN200610061969A CN101116622A CN 101116622 A CN101116622 A CN 101116622A CN A2006100619695 A CNA2006100619695 A CN A2006100619695A CN 200610061969 A CN200610061969 A CN 200610061969A CN 101116622 A CN101116622 A CN 101116622A
Authority
CN
China
Prior art keywords
apodization
memory
depth
channel
aperture
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
CNA2006100619695A
Other languages
English (en)
Other versions
CN101116622B (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.)
Shenzhen Mindray Bio Medical Electronics Co Ltd
Shenzhen Mindray Scientific Co Ltd
Original Assignee
Shenzhen Mindray Bio Medical Electronics Co 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 Shenzhen Mindray Bio Medical Electronics Co Ltd filed Critical Shenzhen Mindray Bio Medical Electronics Co Ltd
Priority to CN2006100619695A priority Critical patent/CN101116622B/zh
Priority to US11/607,573 priority patent/US7800979B2/en
Publication of CN101116622A publication Critical patent/CN101116622A/zh
Application granted granted Critical
Publication of CN101116622B publication Critical patent/CN101116622B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • G10K11/348Circuits therefor using amplitude variation

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Multimedia (AREA)
  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

一种波束合成的接收变迹参数的实时计算方法和装置,所述装置作为一个独立附件接入超声成像***的接收波束合成模块的信号输入端,包括:第一至第三乘法器、第一至第四加法器、减法器、除法器、抽取率存储器、接收孔径存储器、抽取起点存储器、预置变迹曲线存储器、变迹深度计数器、通道计数器、系数选择器、数据选择器、地址处理器和定时控制模块。变迹参数的实时计算包括:在***内预置一条长度为N的预置的变迹参照曲线,对所述预置的变迹参照曲线按深度进行不同起点,不同抽取率的抽取,抽取结果作为各信号接收处理通道不同深度的变迹曲线。采用本发明的技术方案能节省***的存储器资源,***切换探头时装载参数的速度较快。

Description

波束合成的接收变迹参数的实时计算方法及其装置
技术领域
本发明涉及本发明涉及医用超声波诊断***中数字波束的合成,尤其涉及旨在节省***存储器资源的在接收进程中波束合成变迹参数的实时计算方法及其装置。
背景技术
超声成像设备中,变迹系数是波束合成不可缺少的一个参数。对于绝大多数超声成像***,都采取把变迹参数提前算好并存储在硬件存储器的方式。由于变迹系数是一个随接收波束深度变化而变化的量,对于线阵梯形扫描及相控阵扫描,变迹系数还和接收线位置相关,所以这些变迹系数的总量可能有几十Mbit,甚至更多。这样必然导致***要增加额外的存储器,造价变高,另外在切换探头时载入新参数需要的时间也较长,不方便使用。
美国6,123,671号专利公开了不同阵元如何从预存的变迹参数读取出对应变迹值的方法,但没有涉及变迹参数的实时计算,由于存储的变迹参数和深度相关,美国6,123,671号专利方案仍然要耗费较大的存储空间,从该专利的内容来看,仅适用于凸阵和线阵垂直发射,因此美国6,123,671号专利方法的适用面受到一定的限制。
硬件实时参数计算也许是一种解决机制。超声成像以扫描线为基本单位,只需要预存少量参数,并在每条扫描线发射接收前写入几个控制参数,即可以在接收期间由硬件自动生成所需的变迹参数,这将大大节省***的存储器资源。
发明内容
本发明要解决的技术问题是提供一种在波束合成中实时计算变迹参数的方法及其实现该方法的硬件装置。
本发明采用如下技术方案:
设计一种波束合成的接收变迹参数的实时计算方法,该方法执行如下步骤:
a.探头各阵元接收来自接收线的反射回波信号;
b.所述探头各阵元把接收到的反射回波信号送入各自的信号接收处理通道进行放大处理和模数转换;
c.经模数转换后的各通道数字化回波数据分别被送入FIFO存储器;
d.变迹参数实时计算装置实时计算生成各通道FIFO存储器的数字化回波信号的变迹参数;
e.接收波束合成模块调用变迹参数实时计算装置产生的变迹参数进行波束合成。
在上述步骤d中实时计算所述变迹参数包括如下步骤:
a.在***内预置一条长度为N的函数,如高斯窗或汉明窗作为变迹参照曲线;
b.根据预设的参数对所述预置的变迹参照曲线按深度进行不同起点,不同抽取率的抽取,抽取结果作为各信号接收处理通道不同深度的变迹曲线。
利用所述变迹参照曲线的对称特性,在***内只保存左半部分或右半部分的N/2点。为保证所需的计算精度,所述N应大于32,N的一个典型值为1024。
利用所述变迹参照曲线抽取不同深度的变迹参数时,首先判断对应阵元是否在孔径以内,孔径以外的变迹系数为0,孔径以内的阵元从StartPj开始抽取,抽取率为Δnj,j为代表不同的接收深度,孔径以内的变迹系数Rapoj,n用下式计算:
其中Pj,n=StartPj+(n-ChanNum/2+Aperj/2)*Δnj+Offsetj
if(Pj,n>N/2-1)Pj,n=N-1-Pj,n
Offsetj=-Fxj/d*Δnj d为阵元间距
Fxj=Fxj-1+ΔFx=Fx0+j*ΔFx,
上式中Win(Pjn)为具有正态分布的函数,Aperj为接收孔径大小,Starpj为不同变迹深度下接收孔径内的变迹曲线抽取的起点值,ΔFx是指两个变迹曲线发生变化的接收焦点之间的横坐标的间距,j为深度,n为阵元序号,对接收通道为ChanNum的***,n的取值范围为0~ChanNum-1。
上述实时计算方法的硬件实现方式包括如下过程:
预先把各有关参数存入各存储器,即把变迹曲线抽取率值Δnj存入抽取率存储器,接收孔径大小Aperj存入接收孔径存储器,不同变迹深度下接收孔径内的变迹曲线抽取的起点值Starpj存入抽取起点存储器,归一化后变迹参照曲线的左半部分存入变迹曲线存储器,则实时计算所述接收变迹参数包括如下步骤:
i、在定时控制模块的控制下,变迹深度计数器按照变迹变化的时间间隔进行加1计数;
ii、所述变迹深度计数器把其计数值分别赋予抽取率存储器、接收孔径存储器和抽取起点存储器作为它们的读地址并分别从所述抽取率存储器读取该变迹深度下的抽取率参数Δnj,从所述接收孔径存储器中读取该变迹深度下的接收孔径参数Aperj,从所述抽取起点存储器中读取该变迹深度下的抽取起点参数StartPj;
iii、所述第二孔径偏转参数ΔFx/d与变迹深度计数值经过第一乘法器相乘,得到的结果再与第一孔径偏转参数Fx0/d经过第一加法器相加得到该变迹深度下的因子Fxj/d,然后与所述抽取率参数Δnj在第二乘法器处相乘,相乘的结果就是在此模式下的偏移量Offsetj
iv、如果与所述通道计数器的计数值对应的通道在接收孔径Apefj之内,则系数选择信号选择第二乘法器计算输出偏移量;如果与所述通道计数器的计数值对应的通道在接收孔径Aperj外,则系数选择信号选择第二乘法器计算输出0,该选择控制在系数选择器处实现;
v、该深度下的接收孔径值从接收孔径存储器中读出后除以2,相除后的结果与通道计数器的计数值减去ChanNum/2的结果在第二加法器处相加,相加后的结果在第三乘法器处乘以抽取率Δnj,相乘的结果即为因子(n-ChanNum/2+Aperj)*Δnj
vi、经过第三加法器与第四加法器将偏移量Offsetj、抽取起点StartPj和因子(n-ChanNum/2+Aperj)*Δnj相加,得到在变迹深度j通道n的抽取坐标,该抽取坐标对应因子Pj,n,该抽取坐标在地址处理模块中进行地址处理,如果抽取坐标Pj,n大于变迹曲线长度的一半,即N/2时,地址取坐标取N-1-Pj,n;如果抽取坐标小于等于变迹曲线长度的一半,地址取Pj,n;
vii、地址处理模块输出的地址用做变迹曲线存储器的读地址,读出来的数据经过一个数据选择模块,如果该通道在接收孔径内,数据选择信号选择存储器数据输出,该地址读出来的数据就是该变迹深度j通道n的变迹参数值;如果该通道在接收孔径外,数据选择信号选择0数据,则在该变迹深度j下该通道n的变迹参数值为0;
viii、定时控制模块在变迹深度j时控制通道计数器以1的时间间隔从通道0计数到通道ChanNum-1,循环执行步骤i至步骤viii从而获得ChanNum个通道的变迹参数实时计算。
所述接收通道数ChanNum取值范围为1至512,其典型值为32、64或128。
本发明孙采用的技术方案还包括一种波束合成的接收变迹参数的实时计算装置,该装置作为一个独立附件接入超声成像***的接收波束合成模块的信号输入端,其特征在于所述计算装置包括:
第一至第三乘法器、第一至第四加法器、减法器、除法器、抽取率存储器、接收孔径存储器、抽取起点存储器、变迹曲线存储器、变迹深度计数器、通道计数器、系数选择器、数据选择器、地址处理器和定时控制模块:
所述定时控制模块的信号输出端分别连接所述通道计数器和变迹深度计数器,所述变迹深度计数器把其计数值分别送给抽取率存储器、接收孔径存储器和抽取起点存储器,所述第一乘法器的输入端连接第二孔径偏转参数ΔFx/d与变迹深度计数器,其输出值与第一孔径偏转参数Fx0/d一起被送入所述第一加法器;
所述第二乘法器的输入端分别连接第一加法器和抽取率存储器的输出端,其输出端连接系数选择器;所述除法器的输入端接所述接收孔径存储器,其输出端接所述第二加法器;所述减法器把所述通道计数器的计数值减去ChanNum/2,其输出数值被送入所述第二加法器;
所述第三乘法器的输入端分别连接第二加法器和抽取率存储器的输出端,其输出端与系数选择器的输出端一起连接第三加法器,所述第三加法器的输出与抽取起点存储器的输出在第四加法器处相加,被送入所述地址处理模块中进行地址处理,地址处理模块输出的地址用做所述变迹曲线存储器的读地址,读出来的数据经过所述数据选择模块,该数据选择模块输出变迹深度j通道n的变迹值。
与现有技术相比,本发明波束合成的接收变迹参数的实时计算方法及其装置具有如下优点:在保证波束合成质量的前提下,***需要存储的参数量更少,因此节省了***的存储器资源,***切换探头时装载参数的速度也较快。
附图说明
图1是采用本发明装置的一个超声波成象***的原理方框图;
图2是本发明所采用的变迹参照曲线图;
图3是本发明是超声波成象***接收变迹系数获取方式示意图;
图4是是本发明考虑波束校正的波束合成的接收变迹参数的实时计算方法的硬件构成框图。
具体实施方式
下面将结合附图和实施例对本发明方法及其装置进行详细的说明。
图1是B型超声成像***的原理方框图,本发明的接收变迹参数的实时计算装置在图1中作为作为一个独立附件接入超声成像***的接收波束合成模块的信号输入端。
下面主要以64通道单波束的B型超声成像***为例进行说明,很显然,本发明的技术方案很容易推广到其它通道数(如24、48或更多通道)以及多波束的情况。
一、波束合成原理可以用公式表示为:
BFecho ( j ) = Σ i = 1 n [ rapo ( i ) ( j ) * A ( j ) * Gecho ( i ) ( j ) + rapo ( i ) ( j ) * B ( j ) * Gecho ( i ) ( j + 1 ) ] - - - ( 1 )
式中:BFecho(j)为波束合成结果,j代表时间,i为接收通道号,n为接收孔径,Gecho为经接收粗延时调整后各通道的回波信号,A和B为插值系数。rapo即为变迹参数,从(1)式可以看出,它的作用在于给不同通道的回波以不同的权重。一般情况下,各个通道的变迹参数各不相同,并且随深度变化而变化。
在实际***设置中,不同深度的变迹曲线(即不同通道的变迹值组成的曲线)尽管各不相同,但是都可以认为是同一条曲线(大多数情况下为高斯窗或者汉明窗)经过拉伸、压缩后的一段或者全部。基于此,本发明提出了一种基于对同一条预置曲线进行不同抽取,得到不同深度变迹曲线的方法。
变迹系数由系数预存的一条变迹曲线中抽取得到(参见图2)。预存的变迹曲线的长度为N,考虑到一般情况下该曲线肯定为对称曲线,因此只保存N/2点(如图2曲线的左半部分),每点8bit。在硬件存储量许可的情况下,N值取大一点比较好,在本发明的一个实现中,N值取1024点。
硬件由变迹曲线抽取变迹系数时,首先判断对应阵元是否在孔径以内,孔径以外的变迹系数为0;孔径以内的阵元从StartPj开始抽取,抽取率为Δnj,此处j是代表不同的接收深度。StartPj和Δnj由软件计算好并预存在RAM中。假定变迹系数定义做Rapoj,n,j为深度,n为阵元序号(对接收通道为ChanNum的***,n取值范围是0~ChanNum-1)。Rapoj,n计算方式如下:
其中Pj,n=StartPj+(n-ChanNum/2+Aperj/2)*Δnj    (2)
if(Pj,n>N/2-1)Pj,n=N-1-Pj,n
上式中win为预存的变迹曲线,Aperj代表深度j处的接收孔径。设置参数时,通过Aperj,StartPj和Δnj保证按照式(1)抽取出来的变迹曲线中心点始终为最大,并且基本对称。
式(1)满足凸阵和线阵垂直发射和接收的情况,并且要求接收线起点位于接收孔径的中心。如图3中接收线L,O为孔径中心,接收线L的起点和O重合。
二、考虑孔径偏移
实际情况下,接收线往往并不位于接收孔径的中心,或者说孔径发生了偏移(图3中接收线M,O为孔径中心,O’为接收线M起点),这时候就不能直接使用式(1)计算接收线M的变迹。
接收线M和接收线L主要区别在于与接收线M最接近的阵元不一样,一般来说我们希望离接收焦点最近的阵元有最大的变迹,只需要将式(1)中Pj,n的计算方式变化一下即可,式(1)变为:
其中Pj,n=StartPj+(n-ChanNum/2+Aperj/2)*Δnj+Offsetj    (2)
if(Pj,n>N/2-1)Pj,n=N-1-Pj,n
Offsetj=O′O/d*Δnj  d为阵元间距
式(2)相对于式(1),在计算Pj,n时多加了一个偏移量Offsetj,该偏移量和深度相关,在实现时可以由软件提前计算出O’O/d,在扫描开始前写入硬件,然后实时计算Offsetj。式(2)适用于凸阵和线阵的垂直扫描。
三、偏转接收
对于相控阵和线阵梯形扫描,接收扫描线并非垂直于探头表面,(图3中接收线N,O为孔径中心,O”为接收线N起点,α为接收线偏转的角度)。由于接收线有一个偏转角度,所以对每一接收焦点,与其最邻近的探头阵元均不相同。所以式(2)不能适用于接收线N。对于接收焦点坐标位于(Fxj,Fzj)的接收焦点,变迹曲线Rapoj,n计算方式如下:
Figure A20061006196900131
其中Pj,n=StartPj+(n-ChanNum/2+Aperj/2)*Δnj+Offsetj    (3)
if(Pj,n>N/2-1)Pj,n=N-1-Pj,n
Offsetj=-Fxj/d*Δnj   d为阵元间距
式(1)和式(2)都可以统一到式(3)中,对于式(1),Fxj始终为0,对于式(2),Fxj深度不变,始终为00’,而更广泛的情况,Fxj随深度发生变化。
Fxj的值可以由硬件实时计算。计算方法为:
Fxj=Fxj-1+ΔFx=Fx0+j*ΔFx   (4)
ΔFx是指两个变迹曲线发生变化的接收焦点之间的横坐标间距(参见图3),因为变迹曲线产生变化的时间间隔是固定的,所以ΔFx也是个定值。
式(3)中,计算Offsetj需要用到除法运算,一般硬件计算电路均希望尽量避免除法,所以可以按如下方式计算Fxj/d:
Fxj/d=(Fx0+j*ΔFx)/ d=Fx0/d+j*ΔFx/d    (5)
Fx0/d和ΔFx/d的值由软件计算好,在扫描开始前写入硬件的寄存器,硬件通过累加的方式计算Fxj/d。
四、硬件实现方式
以一个64通道单波束的超声成像***为例,在进行波束合成时需要进行最多64通道的变迹参数计算,实现装置如图4所示。
在该实现装置中包括变迹曲线存储器、抽取起点存储器、抽取率存储器以及接收孔径存储器,包括第一偏转参数以及第二偏转参数,还包括深度计数器、通道计数器以及定时控制模块。第一偏转参数对应于算法中的Fx0/d因子,第二偏转参数对应于算法中的ΔFx/d因子。
抽取率存储器存储了不同变迹深度下的变迹曲线抽取率值,这个存储器的地址对应变迹深度,地址中的值对应该深度下的抽取率值。接收孔径存储器存储了不同变迹深度下的接收孔径大小,这个存储器的地址对应变迹深度,地址中的值对应该变迹深度下的接收孔径大小。抽取起点存储器存储了不同变迹深度下接收孔径内的变迹曲线抽取的起点值,这个存储器的地址对应变迹深度,地址中的值对应抽取起点值。变迹曲线存储器存储了一条左右对称归一化曲线的左半部分,地址对应曲线的x轴坐标,地址中的数值对应与曲线的幅度值。
在定时控制模块的控制下,变迹深度计数器按照变迹变化的时间间隔进行加1计数,在某个变迹深度下,控制通道计数器按一定的时间间隔从通道0计数到通道63。
变迹深度计数器的计数值是抽取率存储器、接收孔径存储器以及抽取起点存储器的读地址,根据变迹深度计数器的计数值,可以分别从抽取率存储器读取该变迹深度下的抽取率参数Δnj,从接收孔径存储器中读取该变迹深度下的接收孔径参数Aperj,从抽取起点存储器中读取该变迹深度下的抽取起点参数StartPj。
第二孔径偏转系数与变迹深度计数值经过第一乘法器相乘,得到的结果再与第一孔径偏转系数经过第一加法器相加可以得到该变迹深度下的因子Fxj/d。然后与抽取率参数Δnj在第二乘法器处相乘,相乘的结果就是在此模式下的偏移量。
如果与通道计数器的计数值对应的通道在接收孔径内,则系数选择信号选择第二乘法器计算的偏移量输出;如果与通道计数器的计数值对应的通道在接收孔径外,则系数选择信号选择输出0。该选择控制在系数选择器处实现。
该深度下的接收孔径值从接收孔径存储器中读出后除以2(该除法器可以用右移一位简单实现)。相除后的结果与通道计数器的计数值减去ChanNum/2的结果在第二加法器处相加,相加后的结果在第三乘法器处乘以抽取率Δnj,相乘的结果即为公式(3)中的因子(n-ChanNum/2+Aperj)*Δnj。
经过第三加法器与第四加法器将偏移量、抽取起点StartPj以及因子(n-ChanNum/2+Aperj)*Δnj相加,得到在变迹深度j通道n的抽取坐标,该抽取坐标对应公式(3)中的因子Pj,n。该抽取坐标在地址处理模块中进行地址处理。如果抽取坐标大于变迹曲线长度的一半(即N/2)时,地址取坐标取N-1-Pj,n;如果抽取坐标小于等于变迹曲线长度的一半,地址取Pj,n。地址处理模块输出的地址用做于变迹曲线存储器的读地址,读出来的数据经过一个数据选择模块,如果该通道在接收孔径内,数据选择信号选择存储器数据输出,该地址读出来的数据就是该变迹深度j通道n的变迹值;如果该通道在接收孔径外,数据选择信号选择0数据,则在该变迹深度下该通道的变迹值为0。
定时控制模块在变迹深度j时控制通道计数器从0计数到63便可以完成64通道的变迹参数计算。当定时控制模块控制变迹深度在整个波束合成期从变迹深度0计数到最大的扫描深度时就可以完成64通道单波束***全部变迹参数的实时计算。

Claims (8)

1.一种波束合成的接收变迹参数的实时计算方法,其特征在于该方法执行如下步骤:
a.探头各阵元接收来自接收线的反射回波信号;
b.所述探头各阵元把接收到的反射回波信号送入各自的信号接收处理通道进行放大处理和模数转换;
c.经模数转换后的各通道数字化回波数据分别被送入FIFO存储器;
d.变迹参数实时计算装置实时计算生成各通道FIFO存储器的数字化回波信号的变迹参数;
e.接收波束合成模块调用变迹参数实时计算装置产生的变迹参数进行波束合成。
2.根据权利要求1所述的波束合成的接收变迹参数的实时计算方法,其特征在于:在步骤d中实时计算所述变迹参数包括如下步骤:
a.在***内预置一条长度为N的函数,如高斯窗或汉明窗作为变迹参照曲线;
b.根据预设的参数对所述预置的变迹参照曲线按深度进行不同起点,不同抽取率的抽取,抽取结果作为各信号接收处理通道不同深度的变迹曲线。
3.根据权利要求2所述的波束合成的接收变迹参数的实时计算方法,其特征在于:利用所述变迹参照曲线的对称特性,在***内只保存左半部分或右半部分的N/2点。
4.根据权利要求2或3所述的波束合成的接收变迹参数的实时计算方法,其特征在于:为保证所需的计算精度,所述N应大于32,所述N的一个典型值为1024。
5.根据权利要求3所述的波束合成的接收变迹参数的实时计算方法,其特征在于:利用所述变迹参照曲线抽取不同深度的变迹参数时,首先判断对应阵元是否在孔径以内,孔径以外的变迹系数为0,孔径以内的阵元从StartPj开始抽取,抽取率为Δnj,j为代表不同的接收深度,孔径以内的变迹系数Rapoj,n用下式计算:
Figure A2006100619690003C1
其中Pj,n=StartPj+(n-ChanNum/2+Aperj/2)*Δnj+Offsetj
if(Pj,n>N/2-1)Pj,n=N-1-Pj,n
Offsetj=-Fxj/d*Δnj    d为阵元间距
Fxj=Fxj-1+ΔFx=Fx0+j*ΔFx,
上式中Win(Pjn)为预置的变迹参数曲线,Aperj为接收孔径大小,Starpj为不同变迹深度下接收孔径内的变迹曲线抽取的起点值,ΔFx是指两个变迹曲线发生变化的接收焦点之间的横坐标的间距,j为深度,n为阵元序号,对接收通道为ChanNum的***,n的取值范围为0~ChanNum-1。
6.根据权利要求5所述的波束合成的接收变迹参数的实时计算方法,其特征在于该实时计算方法的硬件实现方式包括如下过程:
预先把各有关参数存入各存储器,即把变迹曲线抽取率值Δnj存入抽取率存储器,接收孔径大小Aperj存入接收孔径存储器,不同变迹深度下接收孔径内的变迹曲线抽取的起点值Starpj存入抽取起点存储器,归一化后变迹参照曲线的左半部分存入变迹曲线存储器,则实时计算所述接收变迹参数包括如下步骤:
i、在定时控制模块的控制下,变迹深度计数器按照变迹变化的时间间隔进行加1计数;
ii、所述变迹深度计数器把其计数值分别赋予抽取率存储器、接收孔径存储器和抽取起点存储器作为它们的读地址并分别从所述抽取率存储器读取该变迹深度下的抽取率参数Δnj,从所述接收孔径存储器中读取该变迹深度下的接收孔径参数Aperj,从所述抽取起点存储器中读取该变迹深度下的抽取起点参数StartPj
iii、所述第二孔径偏转参数ΔFx/d与变迹深度计数值经过第一乘法器相乘,得到的结果再与第一孔径偏转参数Fx0/d经过第一加法器相加得到该变迹深度下的因子Fxj/d,然后与所述抽取率参数Δnj在第二乘法器处相乘,相乘的结果就是在此模式下的偏移量Offsetj
iv、如果与所述通道计数器的计数值对应的通道在接收孔径Aperj之内,则系数选择信号选择第二乘法器计算输出偏移量;如果与所述通道计数器的计数值对应的通道在接收孔径Aperj之外,则系数选择信号选择第二乘法器计算输出0,该选择控制在系数选择器处实现;
v、该深度下的接收孔径值从接收孔径存储器中读出后除以2,相除后的结果与通道计数器的计数值减去ChanNum/2的结果在第二加法器处相加,相加后的结果在第三乘法器处乘以抽取率Δnj,相乘的结果即为因子(n-ChanNum/2+Aperj)*Δnj
vi、经过第三加法器与第四加法器将偏移量Offsetj、抽取起点StartPj和因子(n-ChanNum/2+Aperj)*Δnj相加,得到在变迹深度j通道n的抽取坐标,该抽取坐标对应因子Pj,n,该抽取坐标在地址处理模块中进行地址处理,如果抽取坐标Pj,n大于变迹曲线长度的一半,即N/2时,地址取坐标取N-1-Pj,n;如果抽取坐标小于等于变迹曲线长度的一半,地址取Pj,n
vii、地址处理模块输出的地址用做变迹曲线存储器的读地址,读出来的数据经过一个数据选择模块,如果该通道在接收孔径内,数据选择信号选择存储器数据输出,该地址读出来的数据就是该变迹深度j通道n的变迹参数值;如果该通道在接收孔径外,数据选择信号选择0数据,则在该变迹深度j下该通道n的变迹参数值为0;
viii、定时控制模块在变迹深度j时控制通道计数器以1的时间间隔从通道0计数到通道ChanNum-1,循环执行步骤i至步骤viii从而获得ChanNum个通道的变迹参数实时计算。
7.根据权利要求6所述的波束合成的接收变迹参数的实时计算方法,其特征在于:所述接收通道数ChanNum取值范围为1至512,其典型值为32、64或128。
8.一种波束合成的接收变迹参数的实时计算装置,该装置作为一个独立附件接入超声成像***的接收波束合成模块的信号输入端,其特征在于所述计算装置包括:
第一至第三乘法器、第一至第四加法器、减法器、除法器、抽取率存储器、接收孔径存储器、抽取起点存储器、变迹曲线存储器、变迹深度计数器、通道计数器、系数选择器、数据选择器、地址处理器和定时控制模块:
所述定时控制模块的信号输出端分别连接所述通道计数器和变迹深度计数器,所述变迹深度计数器把其计数值分别送给抽取率存储器、接收孔径存储器和抽取起点存储器,所述第一乘法器的输入端连接第二孔径偏转参数ΔFx/d与变迹深度计数器,其输出值与第一孔径偏转参数Fx0/d一起被送入所述第一加法器;
所述第二乘法器的输入端分别连接第一加法器和抽取率存储器的输出端,其输出端连接系数选择器;所述除法器的输入端接所述接收孔径存储器,其输出端接所述第二加法器;所述减法器把所述通道计数器的计数值减去ChanNum/2,其输出数值被送入所述第二加法器;
所述第三乘法器的输入端分别连接第二加法器和抽取率存储器的输出端,其输出端与系数选择器的输出端一起连接第三加法器,所述第三加法器的输出与抽取起点存储器的输出在第四加法器处相加,被送入所述地址处理模块中进行地址处理,地址处理模块输出的地址用做所述变迹曲线存储器的读地址,读出来的数据经过所述数据选择模块,该数据选择模块输出变迹深度j通道n的变迹值。
CN2006100619695A 2006-08-02 2006-08-02 波束合成的接收变迹参数的实时计算方法及其装置 Active CN101116622B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN2006100619695A CN101116622B (zh) 2006-08-02 2006-08-02 波束合成的接收变迹参数的实时计算方法及其装置
US11/607,573 US7800979B2 (en) 2006-08-02 2006-12-01 Beamforming method and apparatus used in ultrasonic imaging system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2006100619695A CN101116622B (zh) 2006-08-02 2006-08-02 波束合成的接收变迹参数的实时计算方法及其装置

Publications (2)

Publication Number Publication Date
CN101116622A true CN101116622A (zh) 2008-02-06
CN101116622B CN101116622B (zh) 2010-12-01

Family

ID=39030128

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2006100619695A Active CN101116622B (zh) 2006-08-02 2006-08-02 波束合成的接收变迹参数的实时计算方法及其装置

Country Status (2)

Country Link
US (1) US7800979B2 (zh)
CN (1) CN101116622B (zh)

Cited By (26)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101940479A (zh) * 2010-07-15 2011-01-12 重庆大学 基于分段动态变迹系数的超声成像波束合成方法及其装置
CN101961251A (zh) * 2010-09-29 2011-02-02 深圳市蓝韵实业有限公司 一种医学超声诊断***中实时计算变迹曲线的方法及装置
CN103181782A (zh) * 2011-12-29 2013-07-03 三星麦迪森株式会社 超声***和提供多普勒频谱图像的方法
CN104582582A (zh) * 2012-08-21 2015-04-29 毛伊图像公司 超声成像***存储器架构
CN105520752A (zh) * 2014-10-21 2016-04-27 爱飞纽医疗机械贸易有限公司 波束成形设备、超声成像设备以及波束成形方法
US9526475B2 (en) 2006-09-14 2016-12-27 Maui Imaging, Inc. Point source transmission and speed-of-sound correction using multi-aperture ultrasound imaging
US9572549B2 (en) 2012-08-10 2017-02-21 Maui Imaging, Inc. Calibration of multiple aperture ultrasound probes
US9582876B2 (en) 2006-02-06 2017-02-28 Maui Imaging, Inc. Method and apparatus to visualize the coronary arteries using ultrasound
US9668714B2 (en) 2010-04-14 2017-06-06 Maui Imaging, Inc. Systems and methods for improving ultrasound image quality by applying weighting factors
US9788813B2 (en) 2010-10-13 2017-10-17 Maui Imaging, Inc. Multiple aperture probe internal apparatus and cable assemblies
US9883848B2 (en) 2013-09-13 2018-02-06 Maui Imaging, Inc. Ultrasound imaging using apparent point-source transmit transducer
CN107669294A (zh) * 2017-09-22 2018-02-09 青岛海信医疗设备股份有限公司 波束合成中的变迹系数的实时计算方法及装置
CN108113703A (zh) * 2017-12-21 2018-06-05 飞依诺科技(苏州)有限公司 一种用于超声扫描的生成变迹值的方法和装置
US10130333B2 (en) 2006-10-25 2018-11-20 Maui Imaging, Inc. Method and apparatus to produce ultrasonic images using multiple apertures
US10206662B2 (en) 2009-04-14 2019-02-19 Maui Imaging, Inc. Calibration of ultrasound probes
CN109394268A (zh) * 2018-12-07 2019-03-01 余姚市华耀工具科技有限公司 息肉危害程度映射平台
CN109419532A (zh) * 2017-08-30 2019-03-05 佳能株式会社 超声探头和包括超声探头的光声装置
US10226234B2 (en) 2011-12-01 2019-03-12 Maui Imaging, Inc. Motion detection using ping-based and multiple aperture doppler ultrasound
US10267913B2 (en) 2013-03-13 2019-04-23 Maui Imaging, Inc. Alignment of ultrasound transducer arrays and multiple aperture probe assembly
US10401493B2 (en) 2014-08-18 2019-09-03 Maui Imaging, Inc. Network-based ultrasound imaging system
US10617384B2 (en) 2011-12-29 2020-04-14 Maui Imaging, Inc. M-mode ultrasound imaging of arbitrary paths
US10675000B2 (en) 2007-10-01 2020-06-09 Maui Imaging, Inc. Determining material stiffness using multiple aperture ultrasound
CN111772679A (zh) * 2020-08-04 2020-10-16 深圳蓝韵医学影像有限公司 一种超声成像的方法、装置、计算机设备及计算机可读介质
US10835208B2 (en) 2010-04-14 2020-11-17 Maui Imaging, Inc. Concave ultrasound transducers and 3D arrays
US10856846B2 (en) 2016-01-27 2020-12-08 Maui Imaging, Inc. Ultrasound imaging with sparse array probes
US12048587B2 (en) 2020-12-07 2024-07-30 Maui Imaging, Inc. Ultrasound imaging with sparse array probes

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101116622B (zh) * 2006-08-02 2010-12-01 深圳迈瑞生物医疗电子股份有限公司 波束合成的接收变迹参数的实时计算方法及其装置
KR100891056B1 (ko) 2007-05-31 2009-03-31 이재성 개방형 pet-mri통합 유닛
KR100891057B1 (ko) 2007-05-31 2009-03-31 이재성 탈장착형 pet-mri통합 유닛
CN101373181B (zh) * 2007-08-24 2012-03-21 深圳迈瑞生物医疗电子股份有限公司 实时计算逐点变迹系数的方法及装置
US8545406B2 (en) * 2010-12-06 2013-10-01 Texas Instruments Incorporated Dynamic aperture control and normalization for apodization in beamforming
US8929177B2 (en) * 2013-03-14 2015-01-06 Fujifilm Sonosite, Inc. System and method for performing progressive beamforming
US11933892B2 (en) 2016-02-04 2024-03-19 Koninklijke Philips N.V. Ultrasound imaging system and method

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6252975B1 (en) 1998-12-17 2001-06-26 Xerox Corporation Method and system for real time feature based motion analysis for key frame selection from a video
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
CN2387863Y (zh) * 1999-07-19 2000-07-19 开曼迈瑞医疗电子(深圳)有限公司 全数字b超实时接收动态变迹波束合成装置
EP1683266B1 (en) 2003-06-18 2009-04-15 Telefonaktiebolaget LM Ericsson (publ) Power amplifier pre-distortion
CN100390564C (zh) 2003-10-30 2008-05-28 吴鸿明 一种大量程、高精度无线电测量高度的方法
CN101116622B (zh) * 2006-08-02 2010-12-01 深圳迈瑞生物医疗电子股份有限公司 波束合成的接收变迹参数的实时计算方法及其装置
CN101373181B (zh) * 2007-08-24 2012-03-21 深圳迈瑞生物医疗电子股份有限公司 实时计算逐点变迹系数的方法及装置

Cited By (38)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9582876B2 (en) 2006-02-06 2017-02-28 Maui Imaging, Inc. Method and apparatus to visualize the coronary arteries using ultrasound
US9986975B2 (en) 2006-09-14 2018-06-05 Maui Imaging, Inc. Point source transmission and speed-of-sound correction using multi-aperture ultrasound imaging
US9526475B2 (en) 2006-09-14 2016-12-27 Maui Imaging, Inc. Point source transmission and speed-of-sound correction using multi-aperture ultrasound imaging
US10130333B2 (en) 2006-10-25 2018-11-20 Maui Imaging, Inc. Method and apparatus to produce ultrasonic images using multiple apertures
US10675000B2 (en) 2007-10-01 2020-06-09 Maui Imaging, Inc. Determining material stiffness using multiple aperture ultrasound
US10206662B2 (en) 2009-04-14 2019-02-19 Maui Imaging, Inc. Calibration of ultrasound probes
US11998395B2 (en) 2010-02-18 2024-06-04 Maui Imaging, Inc. Point source transmission and speed-of-sound correction using multi-aperture ultrasound imaging
US9668714B2 (en) 2010-04-14 2017-06-06 Maui Imaging, Inc. Systems and methods for improving ultrasound image quality by applying weighting factors
US11172911B2 (en) 2010-04-14 2021-11-16 Maui Imaging, Inc. Systems and methods for improving ultrasound image quality by applying weighting factors
US10835208B2 (en) 2010-04-14 2020-11-17 Maui Imaging, Inc. Concave ultrasound transducers and 3D arrays
CN101940479A (zh) * 2010-07-15 2011-01-12 重庆大学 基于分段动态变迹系数的超声成像波束合成方法及其装置
CN101940479B (zh) * 2010-07-15 2012-10-03 重庆大学 基于分段动态变迹系数的超声成像波束合成方法及其装置
CN101961251A (zh) * 2010-09-29 2011-02-02 深圳市蓝韵实业有限公司 一种医学超声诊断***中实时计算变迹曲线的方法及装置
US9788813B2 (en) 2010-10-13 2017-10-17 Maui Imaging, Inc. Multiple aperture probe internal apparatus and cable assemblies
US10226234B2 (en) 2011-12-01 2019-03-12 Maui Imaging, Inc. Motion detection using ping-based and multiple aperture doppler ultrasound
CN103181782A (zh) * 2011-12-29 2013-07-03 三星麦迪森株式会社 超声***和提供多普勒频谱图像的方法
US10617384B2 (en) 2011-12-29 2020-04-14 Maui Imaging, Inc. M-mode ultrasound imaging of arbitrary paths
US10064605B2 (en) 2012-08-10 2018-09-04 Maui Imaging, Inc. Calibration of multiple aperture ultrasound probes
US9572549B2 (en) 2012-08-10 2017-02-21 Maui Imaging, Inc. Calibration of multiple aperture ultrasound probes
US11253233B2 (en) 2012-08-10 2022-02-22 Maui Imaging, Inc. Calibration of multiple aperture ultrasound probes
CN104582582A (zh) * 2012-08-21 2015-04-29 毛伊图像公司 超声成像***存储器架构
US9986969B2 (en) 2012-08-21 2018-06-05 Maui Imaging, Inc. Ultrasound imaging system memory architecture
US10267913B2 (en) 2013-03-13 2019-04-23 Maui Imaging, Inc. Alignment of ultrasound transducer arrays and multiple aperture probe assembly
US9883848B2 (en) 2013-09-13 2018-02-06 Maui Imaging, Inc. Ultrasound imaging using apparent point-source transmit transducer
US10653392B2 (en) 2013-09-13 2020-05-19 Maui Imaging, Inc. Ultrasound imaging using apparent point-source transmit transducer
US10401493B2 (en) 2014-08-18 2019-09-03 Maui Imaging, Inc. Network-based ultrasound imaging system
CN105520752A (zh) * 2014-10-21 2016-04-27 爱飞纽医疗机械贸易有限公司 波束成形设备、超声成像设备以及波束成形方法
CN105520752B (zh) * 2014-10-21 2019-11-05 爱飞纽医疗机械贸易有限公司 波束成形设备、超声成像设备以及波束成形方法
US10856846B2 (en) 2016-01-27 2020-12-08 Maui Imaging, Inc. Ultrasound imaging with sparse array probes
CN109419532A (zh) * 2017-08-30 2019-03-05 佳能株式会社 超声探头和包括超声探头的光声装置
CN107669294B (zh) * 2017-09-22 2020-03-20 青岛海信医疗设备股份有限公司 波束合成中的变迹系数的实时计算方法及装置
CN107669294A (zh) * 2017-09-22 2018-02-09 青岛海信医疗设备股份有限公司 波束合成中的变迹系数的实时计算方法及装置
CN108113703B (zh) * 2017-12-21 2021-02-09 飞依诺科技(苏州)有限公司 一种用于超声扫描的生成变迹值的方法和装置
CN108113703A (zh) * 2017-12-21 2018-06-05 飞依诺科技(苏州)有限公司 一种用于超声扫描的生成变迹值的方法和装置
CN109394268B (zh) * 2018-12-07 2021-05-11 刘志红 息肉危害程度映射平台
CN109394268A (zh) * 2018-12-07 2019-03-01 余姚市华耀工具科技有限公司 息肉危害程度映射平台
CN111772679A (zh) * 2020-08-04 2020-10-16 深圳蓝韵医学影像有限公司 一种超声成像的方法、装置、计算机设备及计算机可读介质
US12048587B2 (en) 2020-12-07 2024-07-30 Maui Imaging, Inc. Ultrasound imaging with sparse array probes

Also Published As

Publication number Publication date
CN101116622B (zh) 2010-12-01
US20080033299A1 (en) 2008-02-07
US7800979B2 (en) 2010-09-21

Similar Documents

Publication Publication Date Title
CN101116622A (zh) 波束合成的接收变迹参数的实时计算方法及其装置
CN101112320B (zh) 波束合成的接收聚焦参数的实时计算方法及其装置
EP1998676B1 (en) Method and system for controlling delay in an ultrasound receive beamformer
KR100355719B1 (ko) 다단계 지연 소자를 이용한 초음파 수신 비임 형성 장치
US5844139A (en) Method and apparatus for providing dynamically variable time delays for ultrasound beamformer
US4989143A (en) Adaptive coherent energy beam formation using iterative phase conjugation
US4835689A (en) Adaptive coherent energy beam formation using phase conjugation
US8545406B2 (en) Dynamic aperture control and normalization for apodization in beamforming
EP0696791A2 (en) Delay generator for phased array ultrasound beamformer
US9022937B2 (en) Ultrasound device and method for real-time computation of point-by-point apodization coefficients with trigonometric functions
CN101900808B (zh) 一种多通道扫描延时参数计算方法
CN108024798B (zh) 超声波诊断装置以及延迟数据生成方法
CN108784737B (zh) 一种超声成像的波束合成方法和装置
Nikolov et al. Fast parametric beamformer for synthetic aperture imaging
US5647365A (en) Apodization parameter generator for ultrasound medical imaging system
US8926514B2 (en) Iterative time delay values for ultrasound beamforming
KR20070113084A (ko) 다중 수신 스캔라인 형성을 위한 초음파 진단 장치 및 방법
US20070239013A1 (en) Delay controller for ultrasound receive beamformer
US5779640A (en) Ultrasound multi-beam distortion correction system and method
Walker et al. Aberrator integration error in adaptive imaging
JP2017000547A (ja) 超音波診断装置
US5891041A (en) Ultrasonic imaging system adapted for use with ultrasonic probes having different center frequencies
US11633172B2 (en) Synthetic transmit focusing ultrasound system with speed of sound aberration correction
JP2002336249A (ja) 多段構造の遅延素子を用いる超音波受信ビーム成形装置
KR100886932B1 (ko) 단일 보간기를 이용한 다채널 빔집속장치 및 방법

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
TR01 Transfer of patent right

Effective date of registration: 20180528

Address after: 518057 MINDRAY science and technology south twelve road MINDRAY high tech Industrial Park, Shenzhen, Guangdong

Co-patentee after: SHENZHEN MINDRAY SCIENTIFIC Co.,Ltd.

Patentee after: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

Address before: 518057 MINDRAY science and technology south twelve road MINDRAY high tech Industrial Park, Shenzhen, Guangdong

Patentee before: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

TR01 Transfer of patent right
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20080206

Assignee: Shenzhen Mindray Animal Medical Technology Co.,Ltd.

Assignor: SHENZHEN MINDRAY BIO-MEDICAL ELECTRONICS Co.,Ltd.

Contract record no.: X2022440020009

Denomination of invention: Real-time calculation method and device for receiving apodization parameters of beamforming

Granted publication date: 20101201

License type: Common License

Record date: 20220804

EE01 Entry into force of recordation of patent licensing contract