CN102028499A - 超声探头自适应成像系统 - Google Patents

超声探头自适应成像系统 Download PDF

Info

Publication number
CN102028499A
CN102028499A CN 201010603305 CN201010603305A CN102028499A CN 102028499 A CN102028499 A CN 102028499A CN 201010603305 CN201010603305 CN 201010603305 CN 201010603305 A CN201010603305 A CN 201010603305A CN 102028499 A CN102028499 A CN 102028499A
Authority
CN
China
Prior art keywords
probe
primitive
frequency
signal
demodulation
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
CN 201010603305
Other languages
English (en)
Other versions
CN102028499B (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.)
Vinno Technology Suzhou Co Ltd
Original Assignee
Vinno Technology Suzhou 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 Vinno Technology Suzhou Co Ltd filed Critical Vinno Technology Suzhou Co Ltd
Priority to CN2010106033053A priority Critical patent/CN102028499B/zh
Publication of CN102028499A publication Critical patent/CN102028499A/zh
Application granted granted Critical
Publication of CN102028499B publication Critical patent/CN102028499B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Ultra Sonic Daignosis Equipment (AREA)

Abstract

本发明公开了超声探头自适应成像系统,该成像系统根据探头中每一个基元的几何参数和声学参数进行自适应成像,具体来说,在发射波束合成模块或接收波束合成模块中,根据每个基元的几何参数和声学参数,来计算每个基元的发射延时值或接收延时值和变迹系数,通过针对每个基元的发射延时或接收延时的控制,使得发射波束合成或接收波束合成的聚焦的焦点更集中,焦点处的信号更强,提高图像质量;通过针对每个基元的变迹系数的控制,抑制发射波束或接收波束的旁瓣,优化图像效果;解调滤波模块中解调中心频率与解调带宽与探头中存储的中心频率、带宽相匹配,进一步提高图像质量,综上,通过成像系统中发射波束合成模块、接收波束合成模块、解调滤波模块与探头的较好的匹配,提高了超声图像质量,使用该超声探头自适应成像系统,也减少了探头的控制成本和次品率。

Description

超声探头自适应成像系统
技术领域
本发明涉及超声探头自适应成像系统。
背景技术
医学超声成像集成了超声物理学、现代电子学、信息学和生物医学等多学科的原理和技术手段,是继X线成像技术后,在医学中发展最迅速,应用最广泛的技术之一。影响超声图像质量的因素很多,包括换能器,主机系统以及显示器等等。其中,换能器(或称为探头)是超声成像关键因素之一,是超声图像之源。再好的主机,如果换能器的性能不好,也不会得到好的超声图像。探头生产厂商将探头的设计参数控制在一定的误差范围之内就算合格。越好的探头厂商,这个误差范围控制的越好。超声成像系统和超声探头参数性能如果没有很好的匹配,就得不到最优的图像,甚至得不到好的图像。超声成像系统根据探头的统计均值来进行设计,而不会根据每一个探头的实际测量参数进行设计。这样就会导致用户接收到的主机系统和探头工作在一起的时候不能获得最优的图像质量。
影响超声图像质量的关键因素包括探头的几何参数、声学参数。这些参数直接和探头生产的工艺水平相关。好的工艺水平和严格的筛选标准可以将探头的几何参数和声学参数控制在小的误差范围之类,但是这样会大大增加探头的生产成本,从而提高用户的购买成本。其中,几何参数会影响到成像链路上的发射和接收的波束合成。声学参数则影响整个成像链路上的信号质量。如果探头的实际几何参数和声学参数与主机系统设计所依据的均值参数有较大偏离,就会导致图像质量下降。
目前业内的现状就是,好的探头生产厂商对于几何参数和声学参数的控制筛选较为严格,图像质量较好,但生产成本高,卖给用户的价格也相对较高;而相对弱一些的探头生产厂商对于几何参数和声学参数的控制筛选要求较低,图像质量会相对差一些,生产成本较低,卖给用户的价格较低。不论是好的探头生产厂商还是弱一些的探头生产厂商,都没有给出一个好的主机系统与探头最优匹配问题的解决方案,因为生产厂商不可能生产出与目标值零偏差的探头。
目前,探头出厂前,会测试每一个探头的几何参数和声学参数,在误差范围之内的就算合格,否则就是次品。探头生产厂商能够提供给超声成像系统生产厂商的是探头的设计规格书,里面会给出探头的几何参数和声学参数的统计平均值和统计方差。因此,超声成像设备厂商无法得到每一个探头的实际几何参数和声学参数。即便探头生产厂商能够提供每一个探头以及每一个探头所有基元的几何参数和声学参数的测试报告,超声成像设备厂商目前也无法把每一个探头的测试参数录入成像系统。即便能够录入,也得保证所有的探头出厂在成像系统设计之前完成。这在实际的系统设计中是不现实的。
当前,有些探头生产厂商将有些探头的一些信息如基元的传播时间TOF,基元的灵敏度,存入到探头中,系统根据探头来做一些设计补偿来提高成像质量,但是这些是远远不够的。出于提高图像质量的考虑,申请人对超声探头成像进行了相关研究,诞生了本发明。
发明内容
本发明的目的是提供一种基于探头的几何参数和声学参数来进行整个成像链路的自适应优化以达到较优的图像质量的成像系统。
为达到上述目的,本发明采用的技术方案是:超声探头自适应成像系统,所述成像系统包括
探头,所述探头具有多个基元,所述探头中存储有所述多个基元的几何参数和声学参数,所述几何参数包括每个基元的几何位置,所述声学参数包括每个基元的脉冲响应或频率响应、中心频率、带宽;
可与所述探头相通信的发射波束合成模块,所述发射波束合成模块用于从所述探头中调用所述多个基元的几何参数和声学参数后,计算每个基元的发射延时值和变迹系数,然后对发射信号进行波束合成,并将合成好的发射波束施加到所述探头;
可与所述探头相通信的接收波束合成模块,所述接收波束合成模块用于从所述探头中调用所述多个基元的几何参数和声学参数后,计算每个基元的接收延时值和变迹系数,然后对所述探头接收的回波信号进行接收波束合成;
设置在所述接收波束合成模块下游的解调滤波模块,所述解调滤波模块用于从所述探头中调用所述中心频率、带宽来对所述接收波束合成模块输出的接收波束信号进行解调滤波。
超声探头自适应成像系统,所述成像系统包括
探头,所述探头具有多个基元,所述探头中存储有所述多个基元的几何参数和声学参数,所述几何参数包括每个基元的几何位置,所述声学参数包括每个基元的脉冲响应或频率响应、中心频率、带宽;
可与所述探头相通信的发射波束合成模块,所述发射波束合成模块用于从所述探头中调用所述多个基元的几何参数和声学参数后,计算每个基元的发射延时值和变迹系数,然后对发射信号进行波束合成,并将合成好的发射波束施加到所述探头;
可与所述探头相通信的解调滤波模块,所述解调滤波模块用于从所述探头中调用所述中心频率、带宽来对所述探头接收的每一基元的回波信号进行解调滤波;
设置在所述解调滤波模块下游的接收波束合成模块,所述接收波束合成模块用于从所述探头中调用所述多个基元的几何参数和声学参数后,计算每个基元的接收延时值和变迹系数,然后对所述解调滤波模块输出的信号进行接收波束合成。
设置笛卡尔坐标系,坐标原点(0,0,0)定于所述探头表面的中心位置,x轴定义为所述探头表面的切线方向,y轴定义为所述探头的高度方向,z轴定义为所述探头的扫描深度方向,
基元的几何位置设定为(xj,yj,zj),发射或接收焦点位置(focx,focy,focz),扫描线顶点到焦点的距离R,则
D ( j ) = R + offset - ( x j - foc x ) 2 + ( y j - foc y ) 2 + ( z j - foc z ) 2
其中,D(j)为所述发射延时值或所述接收延时值,offset为设定的发射延时偏移值或所述接收延时偏移值。
对基元j施加一个带宽远大于探头带宽的脉冲,利用示波器记录水槽中水听器获取的回波信号即为该基元的脉冲响应hj(n),所述示波器的采样频率为探头带宽的10倍以上,对该基元的脉冲响应hj(n)进行频域表达即可获得该基元的频率响应wj(f),
设定所述探头中多个基元的平均脉冲响应h(n)或平均频率响应w(f),该平均脉冲响应h(n)或平均频率响应w(f)是对多个基元的脉冲响应hj(n)或频率响应wj(f)通过统计取平均获取的,设定施加到所述探头的发射波束信号频谱Hxmit(f),该发射波束信号在被检测组织中传播时的衰减和散射在频域表达为Hatt(f,r)和Hscatter(f,r),则所述探头接收的回波信号频谱Hrx(f,r)
Hrx(f,r)=Hxmit(f)w(f)Hatt(f,r)w(f)Hscatter(f,r)
其中r为深度,f为频率,
设定所述发射波束或所述接收波束的孔径为A,则
A = R F num ( f )
其中R为扫描线顶点到焦点的距离,Fnum(f)是所述回波信号的中心频率f的函数,
F num ( f ) = 1 2 tan α ( f ) + 1 2 R sin α ( f ) · F
其中α(f)为基元的接收角,R为探头的半径,F为焦距,
设定s为基元在探头表面的坐标,孔径A中工作的基元为第L号基元到第M号基元,那么
A = Σ j = L + 1 M | s j - s j - 1 |
对于线阵探头或相控阵探头,sj=xj;对于弧阵探头,则
s j = x j 2 + z j 2 2 z j tan - 1 ( z j 2 x j )
设定变迹窗W(n),窗长为T,变迹窗可以采用任何的窗函数,则基元的变迹系数w(j)
w ( j ) = W [ round ( s j A * T ) ] s j > = 0 W [ round ( s j A * T ) + T ] s j < 0 .
对基元j施加一个带宽远大于探头带宽的脉冲,利用示波器记录水槽中水听器获取的回波信号即为该基元的脉冲响应hj(n),所述示波器的采样频率为探头带宽的10倍以上,对该基元的脉冲响应hj(n)进行频域表达即可获得该基元的频率响应wj(f),设定施加到所述探头的发射波束信号频谱Hxmit(f),该发射波束信号在被检测组织中传播时的衰减和散射在频域表达为Hatt(f,r)和Hscatter(f,r),则所述基元j接收的回波信号频谱Hrx(j)(f,r)
Hrx(j)(f,r)=Hxmit(f)wj(f)Hatt(f,r)wj(f)Hscatter(f,r)
其中r为深度,f为频率,
根据基元j的回波信号的频谱,对每一基元的回波进行解调滤波,之后再对每一基元解调后的回波进行接收波束合成;
设定所述发射波束或所述接收波束的孔径为A,则
A = R F num ( f )
其中R为扫描线顶点到焦点的距离,Fnum(f)是所述回波信号的中心频率f的函数,
F num ( f ) = 1 2 tan &alpha; ( f ) + 1 2 R sin &alpha; ( f ) &CenterDot; F
其中α(f)为基元的接收角,R为探头的半径,F为焦距,
设定s为基元在探头表面的坐标,孔径A中工作的基元为第L号基元到第M号基元,那么
A = &Sigma; j = L + 1 M | s j - s j - 1 |
对于线阵探头或相控阵探头,sj=xj;对于弧阵探头,则
s j = x j 2 + z j 2 2 z j tan - 1 ( z j 2 x j )
设定变迹窗W(n),窗长为T,变迹窗可以采用任何的窗函数,则基元的变迹系数w(j)
w ( j ) = W [ round ( s j A * T ) ] s j > = 0 W [ round ( s j A * T ) + T ] s j < 0 .
由于上述技术方案的运用,本发明与现有技术相比具有下列优点:在发射波束合成模块或接收波束合成模块中,根据每个基元的几何参数和声学参数,来计算每个基元的发射延时值或接收延时值和变迹系数,通过针对每个基元的发射延时或接收延时的控制,使得发射波束合成或接收波束合成的聚焦的焦点更集中,焦点处的信号更强,提高图像质量;通过针对每个基元的变迹系数的控制,抑制发射波束或接收波束的旁瓣,优化图像效果;解调滤波模块中解调中心频率与解调带宽与探头中存储的中心频率、带宽相匹配,进一步提高图像质量,综上,该成像系统根据探头中每一个基元的几何参数和声学参数进行自适应成像,具体来说,成像系统中发射波束合成模块、接收波束合成模块均引入探头中每一个基元的几何参数和声学参数来进行聚焦控制和变迹控制,解调滤波模块也引入探头的声学参数进行解调滤波,实现成像系统中发射波束合成模块、接收波束合成模块、解调滤波模块与探头的较好的匹配,提高了超声图像质量,使用该超声探头自适应成像系统,也减少了探头的控制成本和次品率。
附图说明
附图1为采用笛卡尔坐标系来表示本发明中探头的几何参数的示意图;
附图2为本发明的流程示意图,其中解调滤波模块设置在接收波束合成模块的下游;
附图3为本发明的流程示意图,其中接收波束合成模块设置在解调滤波模块的下游。
具体实施方式
下面结合附图来进一步阐述本发明的结构。
在图2中,超声探头自适应成像系统,包括探头、可与探头相通信的发射波束合成模块和接收波束合成模块、设置在接收波束合成模块下游的解调滤波模块,探头具有多个基元,探头中存储有多个基元的几何参数和声学参数,几何参数包括每个基元的几何位置,声学参数包括每个基元的脉冲响应或频率响应、中心频率、带宽;发射波束合成模块用于从探头中调用多个基元的几何参数和声学参数后,计算每个基元的发射延时值和变迹系数,然后对发射信号进行波束合成,并将合成好的发射波束施加到探头;接收波束合成模块用于从探头中调用多个基元的几何参数和声学参数后,计算每个基元的接收延时值和变迹系数,然后对探头接收的回波信号进行接收波束合成;解调滤波模块用于从探头中调用中心频率、带宽来对接收波束合成模块输出的接收波束信号进行解调滤波,经解调滤波的信号再经过信号及图像处理,得到最终的超声图像。
在图3中的超声探头自适应成像系统中,接收波束合成模块设置在解调滤波模块的下游,其余与图2相同,这样在图3中,解调滤波模块用于从探头中调用中心频率、带宽来对探头接收的每一基元或通道的回波信号进行解调滤波;接收波束合成模块用于从探头中调用多个基元的几何参数和声学参数后,计算每个基元的接收延时值和变迹系数,然后对解调滤波模块输出的信号进行接收波束合成,再经过信号及图像处理,得到最终的超声图像。
在图2和图3中,在发射波束合成模块或接收波束合成模块中,根据每个基元的几何参数和声学参数,来计算每个基元的发射延时值或接收延时值和变迹系数,通过针对每个基元的发射延时或接收延时的控制,使得发射波束合成或接收波束合成的聚焦的焦点更集中,焦点处的信号更强,提高图像质量;通过针对每个基元的变迹系数的控制,抑制发射波束或接收波束的旁瓣,优化图像效果;解调滤波模块中解调中心频率与解调带宽与探头中存储的中心频率、带宽相匹配,进一步提高图像质量,综上,该成像系统根据探头中每一个基元的几何参数和声学参数进行自适应成像,具体来说,成像系统中发射波束合成模块、接收波束合成模块均引入探头中每一个基元的几何参数和声学参数来进行聚焦控制和变迹控制,解调滤波模块也引入探头的声学参数进行解调滤波,实现成像系统中发射波束合成模块、接收波束合成模块、解调滤波模块与探头的较好的匹配,提高了超声图像质量,使用该超声探头自适应成像系统,也减少了探头的控制成本和次品率。
对于探头中基元的几何参数,建立空间坐标系并给出探头的每个基元在坐标系中的空间坐标。坐标系可以定义为任何的坐标系,如笛卡尔坐标系(x,y,z)或者是极坐标系(R,Ф,θ)。探头可以是任何形状的探头,如弧阵探头、线阵探头、相控阵探头。坐标系的原点可以任意选取。在图1中,坐标系采用笛卡尔坐标系,坐标原点(0,0,0)定于换能器表面的中心位置,x轴定义为探头表面的切线方向,y轴定义为探头的高度方向,z轴定义为探头的扫描深度方向。设探头共有N个基元,这样将探头的每个基元的坐标(xj,yj,zj)(1<=j<=N)存入到探头的存储器中。超声成像系统识别到探头时,可以读取每一基元j的几何坐标(xj,yj,zj),这样系统发射和接收部分可以针对不同的探头有不同波束合成控制,而不像现在的系统设计中假设基元的排布是均匀并且对所有探头的几何参数都是相同的。对于声学参数,将每个基元j的脉冲响应hj(n)或者频率响应wj(f)存入到探头的存储器中。基元的脉冲响应按如下方法测得:依次对探头的每一个基元j施加一个窄脉冲(脉冲带宽远大于探头带宽),利用示波器依次记录水槽中水听器获取的回波信号即为基元j的脉冲响应hj(n),示波器的采样频率须在探头带宽的10倍以上;对获取的基元脉冲响应hj(n)进行频域表达即可获得基元j的频率响应wj(f)。当成像系统识别到探头时,可以读入探头存储器中每一基元j的脉冲响应hj(n)或频率响应wj(f),参与到整个信号链路上的频率控制和信号处理。
基元的几何位置设定为(xj,yj,zj),发射或接收焦点位置(focx,focy,focz),扫描线顶点到焦点的距离R,则
D ( j ) = R + offset - ( x j - foc x ) 2 + ( y j - foc y ) 2 + ( z j - foc z ) 2
其中,D(j)为发射延时值或接收延时值,offset为设定的发射延时偏移值或接收延时偏移值。
设定探头中多个基元的平均脉冲响应h(n)或平均频率响应w(f),该平均脉冲响应h(n)或平均频率响应w(f)是对多个基元j的脉冲响应hj(n)或频率响应wj(f)通过统计取平均获取的,设定施加到探头的发射波束信号频谱Hxmit(f),该发射波束信号在被检测组织中传播时的衰减和散射在频域表达为Hatt(f,r)和Hscatter(f,r),则在图2中,探头接收的回波信号频谱Hrx(f,r)
Hrx(f,r)=Hxmit(f)w(f)Hatt(f,r)w(f)Hscatter(f,r)
在图3中,探头接收的每一基元j的回波信号频谱Hrx(j)(f,r)
Hrx(j)(f,r)=Hxmit(f)wj(f)Hatt(f,r)wj(f)Hscatter(f,r)
其中r为深度,f为频率,
根据基元j的回波信号的频谱,对每一基元的回波进行解调滤波,之后再对每一基元解调后的回波进行接收波束合成;
设定所述发射波束或所述接收波束的孔径为A,则
A = R F num ( f )
其中R为扫描线顶点到焦点的距离,Fnum(f)是所述回波信号的中心频率f的函数,
F num ( f ) = 1 2 tan &alpha; ( f ) + 1 2 R sin &alpha; ( f ) &CenterDot; F
其中α(f)为基元的接收角,R为探头的半径,F为焦距,
设定s为基元在探头表面的坐标,孔径A中工作的基元为第L号基元到第M号基元,那么
A = &Sigma; j = L + 1 M | s j - s j - 1 |
对于线阵探头或相控阵探头,sj=xj;对于弧阵探头,则
s j = x j 2 + z j 2 2 z j tan - 1 ( z j 2 x j )
设定变迹窗W(n),窗长为T,变迹窗可以采用任何的窗函数,如汉明窗或高斯窗等,则基元的变迹系数w(j)
w ( j ) = W [ round ( s j A * T ) ] s j > = 0 W [ round ( s j A * T ) + T ] s j < 0 .

Claims (5)

1.超声探头自适应成像系统,其特征在于:所述成像系统包括
探头,所述探头具有多个基元,所述探头中存储有所述多个基元的几何参数和声学参数,所述几何参数包括每个基元的几何位置,所述声学参数包括每个基元的脉冲响应或频率响应、中心频率、带宽;
可与所述探头相通信的发射波束合成模块,所述发射波束合成模块用于从所述探头中调用所述多个基元的几何参数和声学参数后,计算每个基元的发射延时值和变迹系数,然后对发射信号进行波束合成,并将合成好的发射波束施加到所述探头;
可与所述探头相通信的接收波束合成模块,所述接收波束合成模块用于从所述探头中调用所述多个基元的几何参数和声学参数后,计算每个基元的接收延时值和变迹系数,然后对所述探头接收的回波信号进行接收波束合成;
设置在所述接收波束合成模块下游的解调滤波模块,所述解调滤波模块用于从所述探头中调用所述中心频率、带宽来对所述接收波束合成模块输出的接收波束信号进行解调滤波。
2.超声探头自适应成像系统,其特征在于:所述成像系统包括
探头,所述探头具有多个基元,所述探头中存储有所述多个基元的几何参数和声学参数,所述几何参数包括每个基元的几何位置,所述声学参数包括每个基元的脉冲响应或频率响应、中心频率、带宽;
可与所述探头相通信的发射波束合成模块,所述发射波束合成模块用于从所述探头中调用所述多个基元的几何参数和声学参数后,计算每个基元的发射延时值和变迹系数,然后对发射信号进行波束合成,并将合成好的发射波束施加到所述探头;
可与所述探头相通信的解调滤波模块,所述解调滤波模块用于从所述探头中调用所述中心频率、带宽来对所述探头接收的每一基元的回波信号进行解调滤波;
设置在所述解调滤波模块下游的接收波束合成模块,所述接收波束合成模块用于从所述探头中调用所述多个基元的几何参数和声学参数后,计算每个基元的接收延时值和变迹系数,然后对所述解调滤波模块输出的信号进行接收波束合成。
3.根据权利要求1或2所述的超声探头自适应成像系统,其特征在于:
设置笛卡尔坐标系,坐标原点(0,0,0)定于所述探头表面的中心位置,x轴定义为所述探头表面的切线方向,y轴定义为所述探头的高度方向,z轴定义为所述探头的扫描深度方向,
基元的几何位置设定为(xj,yj,zj),发射或接收焦点位置(focx,focy,focz),扫描线顶点到焦点的距离R,则
D ( j ) = R + offset - ( x j - foc x ) 2 + ( y j - foc y ) 2 + ( z j - foc z ) 2
其中,D(j)为所述发射延时值或所述接收延时值,offset为设定的发射延时偏移值或接收延时偏移值。
4.根据权利要求1所述的超声探头自适应成像系统,其特征在于:对基元j施加一个带宽远大于探头带宽的脉冲,利用示波器记录水槽中水听器获取的回波信号即为该基元的脉冲响应hj(n),所述示波器的采样频率为探头带宽的10倍以上,对该基元的脉冲响应hj(n)进行频域表达即可获得该基元的频率响应wj(f),
设定所述探头中多个基元的平均脉冲响应h(n)或平均频率响应w(f),该平均脉冲响应h(n)或平均频率响应w(f)是对多个基元的脉冲响应hj(n)或频率响应wj(f)通过统计取平均获取的,设定施加到所述探头的发射波束信号频谱Hxmit(f),该发射波束信号在被检测组织中传播时的衰减和散射在频域表达为Hatt(f,r)和Hscatter(f,r),则所述探头接收的回波信号频谱Hrx(f,r)
Hrx(f,r)=Hxmit(f)w(f)Hatt(f,r)w(f)Hscatter(f,r)
其中r为深度,f为频率
设定所述发射波束或所述接收波束的孔径为A,则
A = R F num ( f )
其中R为扫描线顶点到焦点的距离,Fnum(f)是所述回波信号的中心频率f的函数,
F num ( f ) = 1 2 tan &alpha; ( f ) + 1 2 R sin &alpha; ( f ) &CenterDot; F
其中α(f)为基元的接收角,R为探头的半径,F为焦距,
设定s为基元在探头表面的坐标,孔径A中工作的基元为第L号基元到第M号基元,那么
A = &Sigma; j = L + 1 M | s j - s j - 1 |
对于线阵探头或相控阵探头,sj=xj;对于弧阵探头,则
s j = x j 2 + z j 2 2 z j tan - 1 ( z j 2 x j )
设定变迹窗W(n),窗长为T,变迹窗可以采用任何的窗函数,则基元的变迹系数w(j)
w ( j ) = W [ round ( s j A * T ) ] s j > = 0 W [ round ( s j A * T ) + T ] s j < 0 .
5.根据权利要求2所述的超声探头自适应成像系统,其特征在于:对基元j施加一个带宽远大于探头带宽的脉冲,利用示波器记录水槽中水听器获取的回波信号即为该基元的脉冲响应hj(n),所述示波器的采样频率为探头带宽的10倍以上,对该基元的脉冲响应hj(n)进行频域表达即可获得该基元的频率响应wj(f),设定施加到所述探头的发射波束信号频谱Hxmit(f),该发射波束信号在被检测组织中传播时的衰减和散射在频域表达为Hatt(f,r)和Hscatter(f,r),则所述基元j接收的回波信号频谱Hrx(j)(f,r)
Hrx(j)(f,r)=Hxmit(f)wj(f)Hatt(f,r)wj(f)Hscatter(f,r)
其中r为深度,f为频率,
根据基元j的回波信号的频谱,对每一基元的回波进行解调滤波,之后再对每一基元解调后的回波进行接收波束合成;
设定所述发射波束或所述接收波束的孔径为A,则
A = R F num ( f )
其中R为扫描线顶点到焦点的距离,Fnum(f)是所述回波信号的中心频率f的函数,
F num ( f ) = 1 2 tan &alpha; ( f ) + 1 2 R sin &alpha; ( f ) &CenterDot; F
其中α(f)为基元的接收角,R为探头的半径,F为焦距,
设定s为基元在探头表面的坐标,孔径A中工作的基元为第L号基元到第M号基元,那么
A = &Sigma; j = L + 1 M | s j - s j - 1 |
对于线阵探头或相控阵探头,sj=xj;对于弧阵探头,则
s j = x j 2 + z j 2 2 z j tan - 1 ( z j 2 x j )
设定变迹窗W(n),窗长为T,变迹窗可以采用任何的窗函数,则基元的变迹系数w(j)
w ( j ) = W [ round ( s j A * T ) ] s j > = 0 W [ round ( s j A * T ) + T ] s j < 0 .
CN2010106033053A 2010-12-24 2010-12-24 超声探头自适应成像系统 Expired - Fee Related CN102028499B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010106033053A CN102028499B (zh) 2010-12-24 2010-12-24 超声探头自适应成像系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010106033053A CN102028499B (zh) 2010-12-24 2010-12-24 超声探头自适应成像系统

Publications (2)

Publication Number Publication Date
CN102028499A true CN102028499A (zh) 2011-04-27
CN102028499B CN102028499B (zh) 2012-10-31

Family

ID=43882339

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010106033053A Expired - Fee Related CN102028499B (zh) 2010-12-24 2010-12-24 超声探头自适应成像系统

Country Status (1)

Country Link
CN (1) CN102028499B (zh)

Cited By (16)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102579074A (zh) * 2011-12-23 2012-07-18 深圳市理邦精密仪器股份有限公司 一种应用于波束合成的边缘补偿方法和装置
CN102860843A (zh) * 2012-09-29 2013-01-09 深圳市理邦精密仪器股份有限公司 一种胎心信号采集的方法及装置
CN102871685A (zh) * 2012-09-14 2013-01-16 飞依诺科技(苏州)有限公司 超声探头几何参数的校正方法和装置及系统
CN103034979A (zh) * 2012-11-30 2013-04-10 声泰特(成都)科技有限公司 一种超声图像清晰度提升方法
CN103405247A (zh) * 2013-08-12 2013-11-27 郑力鹏 超声检查距离控制器
CN103584886A (zh) * 2013-11-20 2014-02-19 无锡祥生医学影像有限责任公司 一种基于相位相干信息的自适应变迹方法
CN103784166A (zh) * 2014-03-03 2014-05-14 哈尔滨工业大学 多功能一体化数字超声诊断系统
CN104688271A (zh) * 2015-03-27 2015-06-10 清华大学 合成聚焦超声成像方法和装置
CN104997533A (zh) * 2015-06-23 2015-10-28 武汉超信电子工程有限公司 超声探头几何参数自动校正方法和装置
CN105640587A (zh) * 2014-11-12 2016-06-08 Ge医疗系统环球技术有限公司 一种用于增强超声影像中的介入式器械的方法及装置
CN108652660A (zh) * 2017-03-28 2018-10-16 美国西门子医疗解决公司 用于医学诊断超声中的衰减估计的衍射校正
WO2018195873A1 (zh) * 2017-04-27 2018-11-01 深圳迈瑞生物医疗电子股份有限公司 一种超声连续波多普勒成像方法及装置、存储介质
CN109974843A (zh) * 2017-12-27 2019-07-05 声博科技股份有限公司 测量声学换能器的宽带环路灵敏度的方法及系统
CN109974844A (zh) * 2017-12-27 2019-07-05 声博科技股份有限公司 测量声学换能器的特征环路灵敏度的方法与系统
CN111544038A (zh) * 2020-05-12 2020-08-18 上海深至信息科技有限公司 一种云平台超声成像系统
CN114098809A (zh) * 2021-11-29 2022-03-01 无锡海鹰电子医疗系统有限公司 一种超声诊断仪发射系统及其快速优化方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101209211A (zh) * 2006-12-30 2008-07-02 深圳迈瑞生物医疗电子股份有限公司 接收孔径可调整的数字化超声波束合成方法及装置
US20090054780A1 (en) * 2007-08-24 2009-02-26 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Method and device for real-time computation of point-by-point apodization coefficients
CN101858972A (zh) * 2010-03-23 2010-10-13 深圳市蓝韵实业有限公司 基于延时参数实时计算和流水线的多波束合成方法和装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101209211A (zh) * 2006-12-30 2008-07-02 深圳迈瑞生物医疗电子股份有限公司 接收孔径可调整的数字化超声波束合成方法及装置
US20090054780A1 (en) * 2007-08-24 2009-02-26 Shenzhen Mindray Bio-Medical Electronics Co., Ltd. Method and device for real-time computation of point-by-point apodization coefficients
CN101858972A (zh) * 2010-03-23 2010-10-13 深圳市蓝韵实业有限公司 基于延时参数实时计算和流水线的多波束合成方法和装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
《电子与信息学报》 20100430 郑驰超等 基于编码发射与自适应波束形成的超声成像 第959-962页 1-5 第32卷, 第4期 2 *

Cited By (28)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102579074B (zh) * 2011-12-23 2014-04-30 深圳市理邦精密仪器股份有限公司 一种应用于波束合成的边缘补偿方法和装置
CN102579074A (zh) * 2011-12-23 2012-07-18 深圳市理邦精密仪器股份有限公司 一种应用于波束合成的边缘补偿方法和装置
CN102871685A (zh) * 2012-09-14 2013-01-16 飞依诺科技(苏州)有限公司 超声探头几何参数的校正方法和装置及系统
CN102871685B (zh) * 2012-09-14 2015-04-08 飞依诺科技(苏州)有限公司 超声探头几何参数的校正方法和装置及系统
CN102860843A (zh) * 2012-09-29 2013-01-09 深圳市理邦精密仪器股份有限公司 一种胎心信号采集的方法及装置
CN103034979A (zh) * 2012-11-30 2013-04-10 声泰特(成都)科技有限公司 一种超声图像清晰度提升方法
CN103034979B (zh) * 2012-11-30 2015-03-25 声泰特(成都)科技有限公司 一种超声图像清晰度提升方法
CN103405247A (zh) * 2013-08-12 2013-11-27 郑力鹏 超声检查距离控制器
CN103584886A (zh) * 2013-11-20 2014-02-19 无锡祥生医学影像有限责任公司 一种基于相位相干信息的自适应变迹方法
CN103584886B (zh) * 2013-11-20 2015-07-15 无锡祥生医学影像有限责任公司 一种基于相位相干信息的自适应变迹方法
CN103784166B (zh) * 2014-03-03 2015-08-19 哈尔滨工业大学 多功能一体化数字超声诊断系统
CN103784166A (zh) * 2014-03-03 2014-05-14 哈尔滨工业大学 多功能一体化数字超声诊断系统
CN105640587A (zh) * 2014-11-12 2016-06-08 Ge医疗系统环球技术有限公司 一种用于增强超声影像中的介入式器械的方法及装置
CN104688271A (zh) * 2015-03-27 2015-06-10 清华大学 合成聚焦超声成像方法和装置
US11307297B2 (en) 2015-03-27 2022-04-19 Tsinghua University Method and device for ultrasonic imaging by synthetic focusing
CN104997533B (zh) * 2015-06-23 2017-11-14 武汉超信电子工程有限公司 超声探头几何参数自动校正方法和装置
CN104997533A (zh) * 2015-06-23 2015-10-28 武汉超信电子工程有限公司 超声探头几何参数自动校正方法和装置
US11166699B2 (en) 2017-03-28 2021-11-09 Siemens Medical Solutions Usa, Inc. Diffraction correction for attenuation estimation in medical diagnostic ultrasound
CN108652660A (zh) * 2017-03-28 2018-10-16 美国西门子医疗解决公司 用于医学诊断超声中的衰减估计的衍射校正
WO2018195873A1 (zh) * 2017-04-27 2018-11-01 深圳迈瑞生物医疗电子股份有限公司 一种超声连续波多普勒成像方法及装置、存储介质
CN110087554A (zh) * 2017-04-27 2019-08-02 深圳迈瑞生物医疗电子股份有限公司 一种超声连续波多普勒成像方法及装置、存储介质
CN109974844A (zh) * 2017-12-27 2019-07-05 声博科技股份有限公司 测量声学换能器的特征环路灵敏度的方法与系统
CN109974843B (zh) * 2017-12-27 2021-10-22 声博科技股份有限公司 测量声学换能器的宽带环路灵敏度的方法及系统
CN109974843A (zh) * 2017-12-27 2019-07-05 声博科技股份有限公司 测量声学换能器的宽带环路灵敏度的方法及系统
CN111544038A (zh) * 2020-05-12 2020-08-18 上海深至信息科技有限公司 一种云平台超声成像系统
CN111544038B (zh) * 2020-05-12 2024-02-02 上海深至信息科技有限公司 一种云平台超声成像系统
CN114098809A (zh) * 2021-11-29 2022-03-01 无锡海鹰电子医疗系统有限公司 一种超声诊断仪发射系统及其快速优化方法
CN114098809B (zh) * 2021-11-29 2024-07-12 无锡海鹰电子医疗系统有限公司 一种超声诊断仪发射系统及其快速优化方法

Also Published As

Publication number Publication date
CN102028499B (zh) 2012-10-31

Similar Documents

Publication Publication Date Title
CN102028499B (zh) 超声探头自适应成像系统
US8167806B2 (en) Ultrasonograph
CN111819467B (zh) 波传播和散射参数的估计方法和仪器
Prada et al. Time reversal techniques in ultrasonic nondestructive testing of scattering media
CN101191786B (zh) 超声波检查方法
US20100142781A1 (en) Systems and Method for Adaptive Beamforming for Image Reconstruction and/or Target/Source Localization
IL281888B1 (en) An image reconstruction method based on non-linear mapping
KR20050058365A (ko) 개선된 고조파 영상 처리를 위한 시스템 및 방법
EP1046056A1 (en) Method and apparatus for improving sidelobe performance of sparse array using harmonic imaging
Wilcox Ultrasonic arrays in NDE: Beyond the B-scan
CN101966088A (zh) 基于柔性相控超声阵列的口腔综合检测方法与装置
Allevato et al. Air-coupled ultrasonic spiral phased array for high-precision beamforming and imaging
JP4643591B2 (ja) 超音波撮像装置
Fink Acoustic time-reversal mirrors
CN110392553B (zh) 用于定位声学传感器的定位设备和系统
CN108415020A (zh) 一种改进的时域超声信号合成孔径算法
US20240264296A1 (en) Method and system for processing beamformed data
JPH03500454A (ja) 人為構造を除外した超音波反射伝送映像化方法および装置
US11885917B2 (en) Methods and instrumentation for estimation of wave propagation and scattering parameters
EP4154031A1 (en) Methods and instrumentation for estimation of wave propagation and scattering parameters
Salari et al. Increasing Penetration Depth and Field-of-view in 3D Ultrasound Imaging with a Lensed Row-column Array
CN114047256A (zh) 基于动态阵元合成孔径聚焦的平板陶瓷膜缺陷超声成像方法
CN105548363A (zh) 基于多途识别的超声检测成像方法
JP3492095B2 (ja) 超音波撮像装置
Kumru et al. Ultrasonic array characterization in multiscattering and attenuating media using pin targets

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C53 Correction of patent for invention or patent application
CB03 Change of inventor or designer information

Inventor after: Tian Yuan

Inventor before: Tian Lin

COR Change of bibliographic data

Free format text: CORRECT: INVENTOR; FROM: TIAN LIN TO: TIAN YUAN

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: 20121031

Termination date: 20171224

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