东北大学 计算机科学与工程学院, 辽宁 沈阳 110169
收稿日期:2017-11-10
基金项目:中央高校基本科研业务费专项资金资助项目(N150403002)。
作者简介:张石(1963-),男,辽宁抚顺人,东北大学教授,博士生导师。
摘要:针对低复杂度的最小方差方法虽降低计算复杂度但成像质量不高的问题, 提出了一种改进的低复杂度幅度相位估计波束形成算法.通过抽取协方差矩阵的有效行来计算自适应加权值, 其中, 有效行数的选取由更加准确的高斯相位相干系数确定, 然后再用高斯相位相干因子对加权矢量做进一步修正, 从而达到降低复杂度的同时保证成像质量的目的.通过实验验证, 将改进算法与原算法和另一种降低复杂度的波束域方法进行对比分析, 充分证明了改进算法在降低复杂度和提升成像质量方面的优越性.
关键词:低复杂度相位相干因子幅度相位估计波束形成超声成像
Improved Low Complexity Amplitude and Phase Estimation Beamforming Algorithm
ZHANG Shi, SHE Li-huang, YAN Xin
School of Computer Science & Engineering, Northeastern University, Shenyang 110169, China
Corresponding author: SHE Li-huang, E-mail: shelihuang@ise.neu.edu.cn
Abstract: To improve the imaging quality in the minimum variance method of low complexity, an improved low-complexity amplitude and phase estimation beamforming algorithm was proposed. The adaptive weighted value was calculated by extracting the effective rows of the covariance matrix, where the selection of the effective number of rows was determined by a more accurate Gaussian phase coherence coefficient, which was then used to modify the weighting vector in order to reduce the complexity and ensure the quality of the imaging at the same time. Through experimental verification, the improved algorithm was compared with the original algorithm and another beamspace method to reduce the complexity, which fully proves the superiority of the improved algorithm in reducing the complexity and improving the imaging quality.
Key words: low complexityphase coherence factoramplitude and phase estimationbeamformingultrasound image
目前, 医学超声成像自适应波束形成算法中, 诸多改进算法被提出, 如:前后向空间平滑技术[1]、特征空间法[2-3]、相干因子方法[4]等, 均能够有效提升成像质量, 但却增加了算法的复杂度, 因而对复杂度的研究也引起了学者们的广泛关注.Nilsen等[5]提出波束域方法, 将阵列信号转换为更低维度的波束域信号, 从而降低矩阵维度.Vignon等[6]提出了一种减少发射和接收阵元数目的波束形成方法; Synnevag等[7]提出了一种低复杂度的波束形成方法, 从预先设计好的一些无失真的权值中自适应选择权值; Asl等[8]将Toeplitz结构运用于协方差矩阵的估计中, 大大降低矩阵求逆的计算量.
Deylami等[9]提出一种低复杂度最小方差波束形成算法, 本文为了与改进算法作对比, 称其为GCF-LMV算法, 该方法通过抽取协方差矩阵的某些行来计算加权值, 从而降低计算量, 若保留协方差矩阵的η行, 则计算复杂度为O(η2L), 与O(L3)相比明显降低.
上述方法虽降低了复杂度, 但其成像质量不高, 因而本文提出了一种改进的低复杂度幅度相位估计波束形成算法, 该算法基于稳健性更好的幅度相位估计算法, 并且在重要参数η的选取中采用高斯相位相干因子代替原算法的广义相干系数, 最后再将该因子对权矢量做进一步修正, 从而大幅度提高成像分辨率.为了验证算法的效果, 将本文算法(gPCF-LAPES)与MV, APES, GCF-LMV, 波束域(B-APES)算法进行仿真对比分析, 结果证明了所提算法的优越性.
1 信号模型若超声换能器由M阵元构成, 各阵元间距相等, 那么超声声场内的散射点经过波束形成, 其输出为
(1) |
(2) |
(3) |
2 幅度相位估计算法APES算法[10]采用一种滤波的方法, 将干扰和噪声信号进行抑制, 从而得到准确的期望频率信号.同时, 该算法对期望信号进行去相关处理, 最终将求解权值W(n)的过程转换为约束最小二乘问题进行计算:
(4) |
(5) |
(6) |
(7) |
(8) |
(9) |
(10) |
(11) |
(12) |
(13) |
(14) |
该算法实现的关键在于如何确定Rη, 若选取行数多于源目标点数, 那么计算量将无法达到最优, 若选取行数少于源目标点数, 那么将会影响成像质量.
该算法计算广义相干系数GCF后, 通过式(15)计算选取行数ηmap:
(15) |
4 低复杂度幅度相位估计算法本文将原算法部分公式进行相应改进以满足所提改进算法的需求, 改进算法中, 所求权矢量需满足式(16)才能够满足APES算法的约束条件,
(16) |
(17) |
(18) |
1) ?计算孔径中所有阵元接收的回波信号;
2) ?将延时信号划分子阵列, 经过空间平滑等处理得到采样协方差矩阵;
3) ?根据APES算法计算噪声协方差矩阵;
4) ?计算gPCF[11];
5) ?根据式(18)计算噪声协方差矩阵所需选取的行数;
6) ?计算算法的加权值W;
7) ?将步骤4)计算的gPCF对步骤6)加权值W进行修正, 得到最终的加权矢量.
5 算法仿真本文采用Field II和MATLAB作为仿真平台, 分别对MV, APES, GCF-LMV, B-APES, gPCF-LAPES算法进行仿真, 实验分为点散射子和囊肿模型实验, 并加入信噪比为60 dB的高斯白噪声来模拟真实的超声声场环境.仿真实验设置参数[10]:换能器类型为线阵, 阵元数量为96, 阵元间距0.1 mm, 中心频率为3 MHz, 子阵列阵元为32, 对角加载常数为10-5, 声速为1 540 m·s-1, 动态显示范围为60 dB.
5.1 点散射目标成像测试点散射子成像实验通过对6对等距分布间隔2 mm的点散射子进行测试.声波发射采用定点聚焦, 聚焦点深度50 mm, 接收方式采取动态聚焦.MV, APES, GCF-LMV, B-APES, gPCF-LAPES 5种算法成像结果如图 1所示.
图 1(Fig. 1)
图 1 点散射目标成像结果Fig.1 Imaging results of point scattering targets (a)—MV; (b)—APES; (c)—B-APES; (d)—GCF-LMV; (e)—gPCF-LAPES. |
由图 1a和图 1b可以看出, MV和APES算法成像质量较差, 旁瓣较为明显, 其点目标之间存在粘连状况, 难以区分; 图 1c的B-APES算法成像结果与原始APES算法相比较, 其成像质量略有降低, 在20~40 mm深度出现明显旁瓣及微弱的伪影; 图 1d为GCF-LMV算法在20~30 mm深度旁瓣伪影杂乱, 但40~60 mm深度成像效果较前三种算法较好; 图 1e中成像效果明显提升, 在各深度处点散射子能够明显区分, 不存在粘连情况, 且没有明显旁瓣出现.
图 2通过描绘横向波束幅度响应图对图 1结果进行进一步的论述, 如图 2所示, GCF-LMV算法主瓣宽度和旁瓣幅度与前三种算法相比均有所降低; 图中gPCF-LAPES算法主瓣宽度最窄, 且旁瓣幅度最低, 由图 2b可得, 该算法的旁瓣幅度降低到60 dB以下, 说明gPCF-LAPES算法的成像分辨率相比于GCF-LMV算法都有了很大的提升, 证明了该算法的优越性.
图 2(Fig. 2)
图 2 不同深度横向波束幅度响应图Fig.2 Lateral responses at different depths (a)—30 mm;(b)—50 mm. |
同时, 本文采用PSL和FWHM两个参数来定量评价算法的成像分辨率, PSL为旁瓣幅度, FWHM为半峰宽度, 深度为50 mm.由表 1计算结果显示, gPCF-LAPES算法旁瓣幅度为-66.7 dB, 且其主瓣宽度0.846 mm, 比原算法降低0.24 mm, 充分说明了gPCF-LAPES算法在提升成像分辨率上的优势.
表 1(Table 1)
表 1 点散射目标50 mm处PSL和FWHM值Table 1 Values of PSL and FWHM at 50 mm in point scatter experiment
| 表 1 点散射目标50 mm处PSL和FWHM值 Table 1 Values of PSL and FWHM at 50 mm in point scatter experiment |
5.2 囊肿模型成像测试本文进行了囊肿模型仿真实验, 以此来评价各算法成像的对比分辨率, 将囊肿模型中心设定在深度50 mm, 且设定其半径为4 mm, 同时在囊肿外部随机分布10 000个散射子用以模拟真实的人体组织成分.图 3为不同波束形成算法的成像结果.
图 3(Fig. 3)
图 3 囊肿模型成像对比图Fig.3 Contrast imaging of cyst model (a)—MV; (b)—APES; (c)—B-APES; (d)—GCF-LMV; (e)—gPCF-LAPES. |
如图 3所示, gPCF-LAPES算法所得的囊肿模型形状与边缘清晰可见, 囊肿内部不含白色点, 其余4种算法成像囊肿内部有轻微的白色伪影存在, 但gPCF-LAPES算法外部出现背景区域过抑制的现象.为了定量评价各算法的对比度, 采用CR(囊肿外部平均幅度-包囊内部平均幅度)来作为评价参数, 如表 2所示, gPCF-LAPES内部平均幅度为-67.27 dB, 外部平均幅度为-27.81 dB, 最终CR值显示其对比度最高, 证明了本算法能够有效地提升成像的对比度.
表 2(Table 2)
表 2 囊肿模型成像对比度Table 2 Imaging contrast of cyst model
| 表 2 囊肿模型成像对比度 Table 2 Imaging contrast of cyst model |
6 结论本文基于低复杂度的最小方差算法, 提出了一种改进的低复杂度的幅度相位估计算法, 只选取协方差矩阵的有效行来计算自适应权值, 降低了计算量, 并通过高斯相位相干因子对计算出的权值进行修正, 得到更加准确的自适应权值.通过点散射子和囊肿模型仿真实验对比得出, 本文改进算法能够解决原算法成像质量较低的问题, 且保持其低复杂度的优点, 但其仍存在背景区域过抑制的问题.
参考文献
[1] | Synnevag J F, Austeng A, Holm S. Adaptive beamforming applied to medical ultrasound imaging[J].IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2007, 54(8): 1603–1613. |
[2] | Mehdizadeh S, Austeng A, Johansen T F, et al. Eigenspace based minimum variance beamforming applied to ultrasound imaging of acoustically hard tissues[J].IEEE Transactions on Medical Imaging, 2012, 31(10): 1912–1921.DOI:10.1109/TMI.2012.2208469 |
[3] | Zeng X, Chen C, Wang Y. Eigenspace-based minimum variance beamformer combined with Wiener postfilter for medical ultrasound imaging[J].Ultrasonics, 2012, 52(8): 996–1004.DOI:10.1016/j.ultras.2012.07.012 |
[4] | Wang S L, Li P C. MVDR-based coherence weighting for high-frame-rate adaptive imaging[J].IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2009, 56(10): 2097–2110.DOI:10.1109/TUFFC.2009.1293 |
[5] | Nilsen C I C, Hafizovic I. Beamspace adaptive beamforming for ultrasound imaging[J].IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2009, 56(10): 2187–2197.DOI:10.1109/TUFFC.2009.1301 |
[6] | Vignon F, Burcher M R. Capon beamforming in medical ultrasound imaging with focused beams[J].IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2008, 55(3): 619–628.DOI:10.1109/TUFFC.2008.686 |
[7] | Synnevag J F, Holm S, Austeng A. A low complexity data-dependent beamformer[J].Ultrasonics Symposium, 2008, 58: 1084–1087. |
[8] | Asl B M, Mahloojifar A. A low-complexity adaptive beamformer for ultrasound imaging using structured covariance matrx[J].IEEE Transactions on Ultrasonics, Ferroelectrics and Frequency Control, 2012, 59(4): 660–667.DOI:10.1109/TUFFC.2012.2244 |
[9] | Deylami A M, Asl B M. Low complex subspace minimum variance beamformer for medical ultrasound imaging[J].Ultrasonics, 2016, 66: 43–53.DOI:10.1016/j.ultras.2015.11.012 |
[10] | 佘黎煌, 闫鑫, 张石. 结合子空间投影的幅度相位估计波束形成算法[J].东北大学学报(自然科学版), 2018, 39(7): 927–930. ( She Li-huang, Yan Xin, Zhang Shi. Amplitude and phase estimation beamforming algorithm combined with subspace projection[J].Journal of Northeastern University (Natural Science), 2018, 39(7): 927–930.) |
[11] | Hasegawa H. Enhancing effect of phase coherence factor for improvement of spatial resolution in ultrasonic imaging[J].Journal of Medical Ultrasonics, 2016, 43(1): 1–9. |