CN102727256B - 基于虚拟阵元的双聚焦波束合成方法及装置 - Google Patents
基于虚拟阵元的双聚焦波束合成方法及装置 Download PDFInfo
- Publication number
- CN102727256B CN102727256B CN201210255474.1A CN201210255474A CN102727256B CN 102727256 B CN102727256 B CN 102727256B CN 201210255474 A CN201210255474 A CN 201210255474A CN 102727256 B CN102727256 B CN 102727256B
- Authority
- CN
- China
- Prior art keywords
- virtual array
- focusing
- sample point
- data
- following formula
- 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
Links
Images
Landscapes
- Ultra Sonic Daignosis Equipment (AREA)
Abstract
本发明公开了一种基于虚拟阵元的双聚焦波束合成方法及装置,首先利用单一固定焦点进行第一次聚焦;然后利用虚拟阵元的概念计算第二次聚焦的延时参数;最后再根据延时叠加波束合成的原理及动态聚焦技术进行第二次聚焦,得到最终成像的扫描线数据。本发明采用将虚拟阵元技术引入到超声成像系统,结合延时叠加波束合成方法,利用2次延时叠加实现双聚焦波束合成超声成像。该方法本质就是利用两次延时叠加进行多个重叠声场的相关有效信息的提取,从而增加声波信号的信息量和能量,使得成像效果获得大幅度的提高。该方法利用两个波束合成器BF1和BF2以及一个FIFO缓存器通过两次延时叠加,无需再存储大量的低分辨率回波数据。
Description
技术领域
本发明涉及医学超声成像领域,特别涉及一种基于虚拟阵元的双聚焦波束合成方法及装置。
背景技术
在医学超声成像领域,基于声场叠加或图像叠加的合成孔径聚焦(synthetic aperturefocusing,SAF)方法得到了广泛应用。与传统的延时叠加波束合成方法相比,合成孔径聚焦成像是一种比较有潜力的超声成像方法,能够明显提高图像的分辨率,但是该方法要求采样和存储每一个孔径点的整个回波信号,即全波采样,其重建理论较复杂,且成像的分辨率对于成像深度有很大的依赖性,无论对软件和硬件要求都比较高;其次单一利用合成孔径聚焦技术,超声成像的回波信号信噪比SNR很低,具有较高的旁瓣等级,图像的对比度、分辨率较差。因此,为进一步提高超声图像的分辨率和增加探测深度,研究学者开展了一系列新方法、新技术的研究。
近年来,虚拟阵元的概念已广泛应用于军事、民用通信以及雷达天线探测系统等阵列信号处理领域,基本思想是在保持真实阵元数目不变的情况下,应用虚拟阵元技术使得阵元数目在虚拟上得到增加,从而减小波束宽度,提高分辨率。然而,目前在国内超声成像领域,虚拟阵元技术应用甚少。随着临床诊断需求的不断增加,对超声成像系统的成像质量提出了更高的要求。
因此急需一种提高超声成像的分辨率的聚焦波束合成方法及装置。
发明内容
有鉴于此,本发明所要解决的技术问题是提供一种提高超声成像的分辨率的聚焦波束合成方法及装置。
本发明的目的之一是提出一种基于虚拟阵元的双聚焦波束合成方法;本发明的目的之二是提出一种基于虚拟阵元的双聚焦波束合成装置。
本发明的目的之一是通过以下技术方案来实现的:
本发明提供的一种基于虚拟阵元的双聚焦波束合成方法,包括以下步骤:
S1:利用单一固定焦点进行发射和接收超声波;
S2:将回波数据进行第一次延时叠加波束合成聚焦;
S3:利用虚拟阵元方法根据第一次聚焦形成的数据来计算聚焦延时参数;
S4:根据聚焦延时参数进行第二次延时叠加波束合成聚焦;
S5:得到最终成像的扫描线回波数据。
进一步,所述第一次波束合成形成的数据存储于缓存器FIFO中。
进一步,所述第一次聚焦包括以下步骤:
S21:通过以下公式来计算滑动子孔径中第i个虚拟阵元的位置xi为:
其中,d为虚拟阵元间距,Zv为虚拟阵元深度,L=Zv/F#/d为滑动子孔径线性传感器个数,D=L*d为滑动子孔径大小,F#为聚焦系数;
S22:通过以下公式来计算虚拟阵元i的延时τi:
其中,c为声速,在虚拟阵元处聚焦以虚拟子孔径中心作为延迟时间参考点;
S23:通过以下公式来计算各个通道的回波数据形成扫描线数据:
其中,Sn(t)为N条扫描线数据,ω(i)为幅度变迹系数,si(t)为子孔径中阵元i接收到的回波信号,r/c为声波从像点到子孔径中心的传播时间,τi为第i个虚拟阵元施加的延时。
进一步,所述第二次聚焦包括以下步骤:
S41:读取缓存器FIFO中的数据并通过以下公式判断数据样本点是否为有效样本点:
其中,d_x为样本点到相应虚拟阵元的侧向距离,d_z为轴向距离,ii为数据样本点编号,jj为虚拟阵元位置点;若Kjj,ii=1,则样本点位于该虚拟阵元的声场内,该样本点即为有效样本点:
S42:通过以下公式计算半孔径角θa:
其中,D=L*d为滑动子孔径的大小,Zv表示虚拟阵元的深度,Zq、Zs分别为样本点Q、S的深度,真实阵元间距d以及虚拟阵元间距Δ,聚焦系数F#=Zv/D;
S43:通过以下公式计算样本点ii到虚拟阵元jj的延时参数:
其中,“±”表示样本点在相应虚拟阵元的下方或者上方,c为声速;
S44:采用逐点接收聚焦波束合成第n条扫描线数据为:
其中,ω(jj)为相应的幅度变迹系数,M为样本点总数,N为BF1所形成的扫描线总数,Sjj,ii(t-τjj,ii)为波束合成器BF1所形成的第jj条扫描线上的第ii个样本点,jj,ii为样本点ii相对于编号是jj的虚拟阵元的延时参数。
本发明的目的之二是通过以下技术方案来实现的:
本发明提供的基于虚拟阵元的双聚焦波束合成装置,包括波束发射接收器、第一波束形成器BF1、第二波束形成器BF2和缓存器FIFO;
所述波束发射接收器,用于利用滑动子孔径采用固定单一焦点进行发射和接收聚焦;
所述第一波束形成器BF1,用于将各个通道的回波数据通过波束形成器BF1进行DAS波束合成;
所述缓存器FIFO,用于存储BF1输出数据;
所述第二波束形成器BF2,用于接收缓存器FIFO中的数据并结合动态聚焦方法进行延时叠加合成波束,得到最终成像的回波数据;
所述成像显示器,用于显示最终成像的回波数据。
进一步,所述第一波束形成器BF1包括以下模块:虚拟阵元位置确定模块、虚拟阵元延时确定模块和回波数据扫描模块,
所述虚拟阵元位置确定模块,用于通过以下公式来计算滑动子孔径中第i个虚拟阵元的位置xi为:
所述虚拟阵元延时确定模块,用于通过以下公式来计算虚拟阵元i的延时τi:
所述回波数据扫描模块,用于通过以下公式来计算各个通道的回波数据形成扫描线数据:
进一步,所述第二波束形成器BF2包括以下模块:有效样本点判断模块、半孔径角θa计算模块,虚拟阵元延时参数模块和第二回报数据扫描模块,
所述有效样本点判断模块,用于读取缓存器FIFO中的数据并通过以下公式判断数据样本点是否为有效样本点:
所述半孔径角θa计算模块,用于通过以下公式计算半孔径角θa:
所述虚拟阵元延时参数模块,用于通过以下公式计算样本点ii到虚拟阵元j j的延时参数:
所述第二回报数据扫描模块,用于采用逐点接收聚焦波束合成第n条扫描线数据为:
本发明的优点在于:本发明采用将虚拟阵元技术引入到超声成像系统,结合延时叠加(delay-and-sum,DAS)波束合成方法,利用2次延时叠加实现双聚焦波束合成(dual focusingbeamforming,DFB)超声成像。该方法本质就是利用两次延时叠加进行多个重叠声场的相关有效信息的提取,从而增加声波信号的信息量和能量,使得成像效果获得大幅度的提高。该方法通过引入虚拟阵元的概念,结合动态孔径聚焦技术,利用两个波束合成器BF1和BF2以及一个FIFO缓存器通过两次延时叠加,无需再存储大量的低分辨率回波数据。
与传统DRF及SF算法相比,该方法的分辨率在一定范围内不会随着探测深度的增加而急剧下降,在一定程度上解决了分辨率与探测深度之间的矛盾,进一步提高超声成像的分辨率,同时不降低成像的深度,并且可以针对不同应用场合来调节虚拟阵元的位置Zv以及聚焦系数F#,实现高质量的超声成像。
本发明的其它优点、目标和特征在某种程度上将在随后的说明书中进行阐述,并且在某种程度上,基于对下文的考察研究对本领域技术人员而言将是显而易见的,或者可以从本发明的实践中得到教导。本发明的目标和其它优点可以通过下面的说明书,权利要求书,以及附图中所特别指出的结构来实现和获得。
附图说明
为了使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明作进一步的详细描述,其中:
图1为传统聚焦发射声场与虚拟阵元声场对比图,(1a)传统发射聚焦,(1b)一个VE的声场,(1c)2个VE的声场;
图2为基于虚拟阵元的双聚焦波束合成DFB原理框图;
图3为BF1形成扫描线及延时参数计算示意图,3(3a)扫描线示意图,图3(3b)给出了形成第n条扫描线时延时参数τ的计算示意图;
图4为虚拟阵元延时参数计算及声场叠加示意图;
图5为动态范围40dB时的不同深度散射点仿真成像示意图;
图6为DRF、SF与DFB的分辨率及能量衰减对比图,(6a)DRF、SF与DFB的侧向分辨率对比,(6b)DRF、SF与DFB在-25dB处主瓣宽度对比。
具体实施方式
以下将结合附图,对本发明的优选实施例进行详细的描述;应当理解,优选实施例仅为了说明本发明,而不是为了限制本发明的保护范围。
图1为传统聚焦发射声场与虚拟阵元声场对比图,(1a)传统发射聚焦,(1b)一个VE的声场,(1c)2个VE的声场,VE表示虚拟阵元;图2为基于虚拟阵元的双聚焦波束合成DFB原理框图;图3为BF1形成扫描线及延时参数计算示意图,3(3a)扫描线示意图,图3(3b)给出了形成第n条扫描线时延时参数τ的计算示意图,如图所示:本发明提供的一种基于虚拟阵元的双聚焦波束合成方法,包括以下步骤:
S1:利用单一固定焦点进行发射和接收超声波;
S2:将回波数据进行第一次延时叠加波束合成聚焦,并将第一次波束合成形成的数据存储于缓存器FIFO中;所述第一次聚焦包括以下步骤:
S21:通过以下公式来计算滑动子孔径中第i个虚拟阵元的位置xi为:
其中,d为虚拟阵元间距,Zv为虚拟阵元深度,F#为聚焦系数;L=Zv/F#/d为滑动子孔径线性传感器个数,L=Zv/F#/d表达式的含义:当聚焦系数一定的情况下,发射子孔径中阵元的个数。D=L*d为滑动子孔径大小。
S22:通过以下公式来计算虚拟阵元i的延时τi:
其中,c为声速,在虚拟阵元处聚焦以虚拟子孔径中心作为延迟时间参考点;
S23:通过以下公式来计算各个通道的回波数据形成扫描线数据:
其中,Sn(t)为N条扫描线数据,ω(i)为幅度变迹系数,si(t)为子孔径中阵元i接收到的回波信号,r/c为声波从像点到子孔径中心的传播时间,τi为第i个虚拟阵元施加的延时。
S3:利用虚拟阵元方法根据第一次聚焦形成的数据来计算聚焦延时参数;
S4:根据聚焦延时参数进行第二次延时叠加波束合成聚焦;所述第二次聚焦包括以下步骤:
S41:读取缓存器FIFO中的数据并通过以下公式判断数据样本点是否为有效样本点:
其中,d_x为样本点到相应虚拟阵元的侧向距离,d_z为轴向距离,ii为数据样本点编号,jj为虚拟阵元位置点;若Kjj,ii=1,则样本点位于该虚拟阵元的声场内,该样本点即为有效样本点:本发明提供的实施例中判断样本的有效性是根据采样点的位置是否处于虚拟阵元的声场范围内来完成的,即:位于虚拟阵元声场范围之内的为有效样本,进而在下一步将所有有效样本求和,增大回波信息量。
S42:通过以下公式计算半孔径角θa:
其中,D=L*d为滑动子孔径的大小,Zv表示虚拟阵元的深度,真实阵元间距d以及虚拟阵元间距Δ,聚焦系数F#=Zv/D;
其中,计算半孔径角的作用是为了确定虚拟阵元的声场范围。图4中Zq、Zs是样本点Q、S的延时时间的计算。真实阵元指的是超声探头中实际存在的超声传感器阵列,而虚拟阵元就是把发射焦点假想成一个超声传感器阵元。
S43:通过以下公式计算样本点ii到虚拟阵元jj的延时参数:
其中,“±”表示样本点在相应虚拟阵元的下方或者上方,c为声速;
S44:采用逐点接收聚焦波束合成第n条扫描线数据为:
其中,ω(jj)为相应的幅度变迹系数,M为样本点总数,N为BF1所形成的扫描线总数,Sjj,ii(t-τjj,ii)为波束合成器BF1所形成的第jj条扫描线上的第ii个样本点,jj,ii为样本点ii相对于编号是jj的虚拟阵元的延时参数;
S5:得到最终成像的扫描线回波数据。
本发明实施例还提供了一种基于虚拟阵元的双聚焦波束合成装置,包括波束发射接收器、第一波束形成器BF1、第二波束形成器BF2和缓存器FIFO;
所述波束发射接收器,用于利用滑动子孔径采用固定单一焦点进行发射和接收聚焦;
所述第一波束形成器BF1,用于将各个通道的回波数据通过波束形成器BF1进行DAS波束合成;
所述缓存器FIFO,用于存储BF1输出数据;
所述第二波束形成器BF2,用于接收缓存器FIFO中的数据并结合动态聚焦方法进行延时叠加合成波束,得到最终成像的回波数据;
所述成像显示器,用于显示最终成像的回波数据。
所述第一波束形成器BF1包括以下模块:虚拟阵元位置确定模块、虚拟阵元延时确定模块和回波数据扫描模块,
所述虚拟阵元位置确定模块,用于通过以下公式来计算滑动子孔径中第i个虚拟阵元的位置xi为:
所述虚拟阵元延时确定模块,用于通过以下公式来计算虚拟阵元i的延时τi:
所述回波数据扫描模块,用于通过以下公式来计算各个通道的回波数据形成扫描线数据:
所述第二波束形成器BF2包括以下模块:有效样本点判断模块、半孔径角θa计算模块,虚拟阵元延时参数模块和第二回报数据扫描模块,
所述有效样本点判断模块,用于读取缓存器FIFO中的数据并通过以下公式判断数据样本点是否为有效样本点:
所述半孔径角θa计算模块,用于通过以下公式计算半孔径角θa:
所述虚拟阵元延时参数模块,用于通过以下公式计算样本点ii到虚拟阵元j j的延时参数:
所述第二回报数据扫描模块,用于采用逐点接收聚焦波束合成第n条扫描线数据为:
对于有限传感器阵列,随着样本点深度的增加,有效虚拟阵元个数将会超出实际虚拟阵元的个数。因此,越接近传感器阵列边缘的地方,有效样本个数不会再继续增加,此时侧向分辨率会略微下降。
图4为虚拟阵元延时参数计算及声场叠加示意图;图5为动态范围40dB时的不同深度散射点仿真成像示意图;图6为DRF、SF与DFB的分辨率及能量衰减对比图,(6a)DRF、SF与DFB的侧向分辨率对比,(6b)DRF、SF与DFB在-25dB处主瓣宽度对比,如图所示,本发明利用Field II进行点散射目标仿真实验,对比了传统动态接收聚焦DRF、合成聚焦SF以及基于虚拟阵元的双聚焦波束合成算法,并重点分析了这三种算法在分辨率和探测深度方面的差异。
DFB表示双聚焦波束合成、BF1表示波束合成器、DRF表示动态接收聚焦、SF表示合成孔径聚焦。
仿真过程中均采用定点发射和分段动态聚焦接收模式,成像采用线性阵列,阵元总数128,发射信号中心频率3.5MHz,采样频率50MHz,阵元中心间距为一个波长,声速1540m/s。信号仿真时加入了60dB的高斯白噪声,成像的动态范围为40dB[16]。设置虚拟源深度Zv=20mm,聚焦系数F#=2,BF1的幅度变迹系数采用hamming窗,BF2的幅度变迹系数采用boxcar窗。为提高成像效果,采用了滑动子孔径技术,子孔径阵元数目为48,
目标散射点共14个:其中深度70mm处设定5个散射点,两相邻散射点横向间隔为2mm,其余散射点分布在深度为10~100mm的区域内,轴向间隔为10mm。图5为不同方法对不同深度散射点的成像结果。从图5可以看出,传统的DRF波束合成算法侧向分辨率差,不能将深度70mm处的5个散射点分辨开来,且随着深度增加,散射点的旁瓣逐渐增加;相比于DRF成像,SF算法虽然提高了侧向分辨率,但是,图像的对比度随着探测深度的增加而急剧下降,且系统复杂度方面远高于DRF;本发明提出的DFB算法,在提高图像的侧向分辨率和对比度方面都远远优于SF以及DRF,系统复杂度略高于DRF,但远低于SF。
为了进一步深入研究不同方法对图像质量的影响情况,图6分别给出了散射点回波数据的归一化幅值在-3dB和-25dB处主瓣宽度随探测深度变化的曲线。从图6中可以看出,虽然DRF及SF算法在近场区域成像的侧向分辨率略优于DFB算法,但随着深度的增加,DFB算法在-3dB处的主瓣宽度基本保持不变,成像效果明显优于DRF和SF。由于信号的主瓣宽度在-25dB处仍然基本保持均匀,从而能量比较集中,因此,DFB算法能够有效地增加探测深度。
以上所述仅为本发明的优选实施例,并不用于限制本发明,显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
Claims (5)
1.一种基于虚拟阵元的双聚焦波束合成方法,其特征在于:包括以下步骤:
S1:利用单一固定焦点进行发射和接收超声波;
S2:将回波数据进行第一次延时叠加波束合成聚焦;
S3:利用虚拟阵元方法根据第一次聚焦形成的数据来计算聚焦延时参数;
S4:第二次聚焦,即根据聚焦延时参数进行第二次延时叠加波束合成聚焦;
S5:得到最终成像的扫描线回波数据;
所述第一次聚焦形成的波束合成形成的数据存储于缓存器FIFO中。
2.根据权利要求1所述的基于虚拟阵元的双聚焦波束合成方法,其特征在于:所述第一次聚焦包括以下步骤:
S21:通过以下公式来计算滑动子孔径中第i个虚拟阵元的位置xi为:
其中,L=Zv/F#/d为滑动子孔径阵元线性传感器个数,d为虚拟阵元间距,Zv为虚拟阵元深度,D=L*d为滑动子孔径大小,F#为聚焦系数;
S22:通过以下公式来计算虚拟阵元i的延时τi:
其中,c为声速,在虚拟阵元处聚焦以虚拟子孔径中心作为延迟时间参考点,Zv表示虚拟阵元的深度;
S23:通过以下公式来计算各个通道的回波数据形成扫描线数据:
其中,Sn(t)为N条扫描线数据,ω(i)为幅度变迹系数,si(t)为子孔径中阵元i接收到的回波信号,r/c为声波从像点到子孔径中心的传播时间,τi为第i个虚拟阵元施加的延时,t表示时间。
3.根据权利要求1所述的基于虚拟阵元的双聚焦波束合成方法,其特征在于:所述第二次延时叠加波束合成聚焦包括以下步骤:
S41:读取缓存器FIFO中的数据并通过以下公式判断数据样本点是否为有效样本点:
其中,d_x为样本点到相应虚拟阵元的侧向距离,d_z为轴向距离,ii为数据样本点编号,jj为虚拟阵元位置点;若Kjj,ii=1,则样本点位于该虚拟阵元的声场内,该样本点即为有效样本点,tanθa表示半孔径角θa的正切值;
S42:通过以下公式计算半孔径角θa:
其中,D=L*d为滑动子孔径的大小,Zv表示虚拟阵元的深度,真实阵元的间距d以及虚拟阵元间距Δ,聚焦系数F#=Zv/D;
S43:通过以下公式计算样本点ii到虚拟阵元jj的延时参数:
其中,“±”表示样本点在相应虚拟阵元的下方或者上方,c为声速;
S44:采用逐点接收聚焦波束合成第n条扫描线数据为:
其中,ω(jj)为相应的幅度变迹系数,M为样本点总数,N为BF1所形成的扫描线总数,Sjj,ii(t-τjj,ii)为波束合成器BF1所形成的第jj条扫描线上的第ii个样本点,jj,ii为样本点ii相对于编号是jj的虚拟阵元的延时参数。
4.基于虚拟阵元的双聚焦波束合成装置,其特征在于:包括波束发射接收器、第一波束形成器BF1、第二波束形成器BF2、缓存器FIFO和成像显示器;
所述波束发射接收器,用于利用滑动子孔径采用固定单一焦点进行发射和接收聚焦;
所述第一波束形成器BF1,用于将各个通道的回波数据通过波束形成器BF1进行DAS波束合成;
所述缓存器FIFO,用于存储BF1输出数据;
所述第二波束形成器BF2,用于接收缓存器FIFO中的数据并结合动态聚焦方法进行延时叠加合成波束,得到最终成像的回波数据;
所述成像显示器,用于显示最终成像的回波数据;
所述第一波束形成器BF1包括以下模块:虚拟阵元位置确定模块、虚拟阵元延时确定模块和回波数据扫描模块,
所述虚拟阵元位置确定模块,用于通过以下公式来计算滑动子孔径中第i个虚拟阵元的位置xi为:
所述虚拟阵元延时确定模块,用于通过以下公式来计算虚拟阵元i的延时τi:
所述回波数据扫描模块,用于通过以下公式来计算各个通道的回波数据形成扫描线数据:
其中,L=Zv/F#/d为滑动子孔径阵元线性传感器个数,d为虚拟阵元间距,Zv为虚拟阵元深度,c为声速,ω(i)为幅度变迹系数,si(t)为子孔径中阵元i接收到的回波信号,r/c为声波从像点到子孔径中心的传播时间。
5.根据权利要求4所述的基于虚拟阵元的双聚焦波束合成装置,其特征在于:所述第二波束形成器BF2包括以下模块:有效样本点判断模块、半孔径角θa计算模块,虚拟阵元延时参数模块和第二回波数据扫描模块,
所述有效样本点判断模块,用于读取缓存器FIFO中的数据并通过以下公式判断数据样本点是否为有效样本点:
所述半孔径角θa计算模块,用于通过以下公式计算半孔径角θa:
所述虚拟阵元延时参数模块,用于通过以下公式计算样本点ii到虚拟阵元jj的延时参数:
所述第二回波数据扫描模块,用于采用逐点接收聚焦波束合成第n条扫描线数据为:
其中,d_x为样本点到相应虚拟阵元的侧向距离,d_z为轴向距离,D=L*d为滑动子孔径大小,F#为聚焦系数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210255474.1A CN102727256B (zh) | 2012-07-23 | 2012-07-23 | 基于虚拟阵元的双聚焦波束合成方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210255474.1A CN102727256B (zh) | 2012-07-23 | 2012-07-23 | 基于虚拟阵元的双聚焦波束合成方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102727256A CN102727256A (zh) | 2012-10-17 |
CN102727256B true CN102727256B (zh) | 2014-06-18 |
Family
ID=46984059
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210255474.1A Active CN102727256B (zh) | 2012-07-23 | 2012-07-23 | 基于虚拟阵元的双聚焦波束合成方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102727256B (zh) |
Families Citing this family (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102895000B (zh) * | 2012-11-06 | 2014-06-18 | 重庆大学 | 一种基于自适应加权的双聚焦波束合成方法 |
CN103033816B (zh) * | 2012-12-07 | 2014-06-04 | 清华大学 | 基于圆弧扫描转换的合成孔径聚焦超声成像实现方法 |
CN103425894B (zh) * | 2013-08-28 | 2016-08-10 | 重庆大学 | 基于误差约束的动态聚焦最优分段方法及装置 |
CN105629248B (zh) * | 2015-12-30 | 2018-06-05 | 中国科学院声学研究所 | 自适应复合成像方法 |
CN108618799B (zh) * | 2018-04-24 | 2020-06-02 | 华中科技大学 | 一种基于空间相干性的超声ct成像方法 |
CN109350115A (zh) * | 2018-12-04 | 2019-02-19 | 珠海医凯电子科技有限公司 | 交错扫描双波束方法 |
CN114129185B (zh) * | 2019-03-18 | 2023-11-07 | 深圳蓝影医学科技股份有限公司 | 波束合成方法、超声成像方法、装置及设备 |
CN110279430B (zh) * | 2019-06-26 | 2021-01-22 | 北京交通大学 | 基于接收端虚拟多角度复合的平面波成像处理方法及装置 |
CN111772679B (zh) * | 2020-08-04 | 2023-03-31 | 深圳蓝影医学科技股份有限公司 | 一种超声成像的方法、装置、计算机设备及计算机可读介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1736964A1 (en) * | 2005-06-24 | 2006-12-27 | Nederlandse Organisatie voor toegepast-natuurwetenschappelijk Onderzoek TNO | System and method for extracting acoustic signals from signals emitted by a plurality of sources |
CN101697010A (zh) * | 2009-10-27 | 2010-04-21 | 电子科技大学 | 一种多功能线型阵列三维合成孔径雷达系统构建方法 |
CN101893705A (zh) * | 2010-06-30 | 2010-11-24 | 重庆大学 | 基于超声成像系统的动态孔径控制方法 |
CN102247168A (zh) * | 2011-08-08 | 2011-11-23 | 重庆博恩富克医疗设备有限公司 | 一种多波束逐点聚焦延时参数处理方法及装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050195103A1 (en) * | 2004-01-13 | 2005-09-08 | Davis Dennis W. | Phased arrays exploiting geometry phase and methods of creating such arrays |
-
2012
- 2012-07-23 CN CN201210255474.1A patent/CN102727256B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1736964A1 (en) * | 2005-06-24 | 2006-12-27 | Nederlandse Organisatie voor toegepast-natuurwetenschappelijk Onderzoek TNO | System and method for extracting acoustic signals from signals emitted by a plurality of sources |
CN101697010A (zh) * | 2009-10-27 | 2010-04-21 | 电子科技大学 | 一种多功能线型阵列三维合成孔径雷达系统构建方法 |
CN101893705A (zh) * | 2010-06-30 | 2010-11-24 | 重庆大学 | 基于超声成像系统的动态孔径控制方法 |
CN102247168A (zh) * | 2011-08-08 | 2011-11-23 | 重庆博恩富克医疗设备有限公司 | 一种多波束逐点聚焦延时参数处理方法及装置 |
Non-Patent Citations (3)
Title |
---|
Catherine H.Frazier等.Synthetic Aperture Techniques with a Virtual Source Element.《IEEE TRANSACTION ON ULTRASONICS,FERROELECTRICS, AND FREQUENCY CONTROL》.1998,第15卷(第4期),第1924-1927页. * |
Synthetic Aperture Flow Imaging Using a Dual Stage Beamformer Approach;Ye Li 等;《2010IEEE International Ultrasonics Symposium Proceedings》;20101231;第196-207页 * |
Ye Li 等.Synthetic Aperture Flow Imaging Using a Dual Stage Beamformer Approach.《2010IEEE International Ultrasonics Symposium Proceedings》.2010,第196-207页. |
Also Published As
Publication number | Publication date |
---|---|
CN102727256A (zh) | 2012-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102727256B (zh) | 基于虚拟阵元的双聚焦波束合成方法及装置 | |
EP3548920B1 (en) | Methods and systems for filtering ultrasound image clutter | |
US8550998B2 (en) | Nonlinear elastic wave measurement and imaging with two-frequency elastic wave pulse complexes | |
CN102895000B (zh) | 一种基于自适应加权的双聚焦波束合成方法 | |
EP1004894B1 (en) | Method and apparatus for high-frame-rate high-resolution ultrasonic image data acquisition | |
CN106680825B (zh) | 一种声学阵列成像系统与方法 | |
JP5411699B2 (ja) | 超音波撮像装置 | |
CN104546003B (zh) | 一种基于平面波发射的可变帧率彩色超高速超声成像方法 | |
CN104414687A (zh) | 超声波测定装置、超声波图像装置及超声波测定方法 | |
CN103969651A (zh) | 自适应声学成像方法 | |
CN113543721B (zh) | 用于采集复合3d超声图像的方法和系统 | |
CN105167802A (zh) | 多普勒成像方法及装置 | |
US20130165783A1 (en) | Providing motion mode image in ultrasound system | |
WO2006057092A1 (ja) | 超音波撮像装置 | |
Hemmsen et al. | Tissue harmonic synthetic aperture ultrasound imaging | |
KR101652727B1 (ko) | 초음파 도플러 영상 장치 및 그 제어 방법 | |
US8235906B2 (en) | System and method for accelerated focused ultrasound imaging | |
US20190046162A1 (en) | Ultrasonic signal processor, ultrasonic diagnostic device, and ultrasonic signal processing method | |
Liu et al. | A multiscale residual U-net architecture for super-resolution ultrasonic phased array imaging from full matrix capture data | |
TWI252929B (en) | Ultrasonic imaging method and device thereof | |
KR101551469B1 (ko) | 평면파를 이용한 초음파 송수신집속 방법 및 시스템 | |
Han et al. | 3D ultrasound imaging in frequency domain based on concepts of array beam and synthetic aperture | |
CN105548363A (zh) | 基于多途识别的超声检测成像方法 | |
EP4155767A1 (en) | Method and system for processing a set of signals received by a transducer element | |
CN113436109B (zh) | 一种基于深度学习的超快高质量的平面波超声成像方法 |
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 |