计算机科学与应用    Vol.4 No.10 (October 2014)

一种基于各向异性总变分最小化的CT图像重建算法
CT Image Reconstruction Algorithm Based on Anisotropic Total Variation Minimization

作者: 张 莹 , 王 丹 ：北京工业大学计算机学院，北京;

关键词: 计算机断层成像； 不完全角度重建； 总变分最小化； 交替方向法； 各向异性总变分最小化； CT Image Reconstruction； Few-View Reconstruction； Total Variation Minimization； Alternating Direction Method； Anisotropic Total Variation Minimization

摘要:
由于受数据采集时间、照射剂量、成像系统扫描的几何位置等因素的约束，计算机断层成像(CT)技术目前只能在有限角度范围或在较少的投影角度得到数据，这些都属于不完全角度重建问题。图像重建问题中的总变分(Total-Variation, TV)最小化模型使用基于交替方向法(alternating direction method, ADM)的稀疏优化算法能够在不完全角度的图像重建中获得较优的重建结果。然而，在极稀疏的角度数量下，各向同性TV最小化算法的重建精度不是很理想，存在进一步改善空间。本文针对该问题，通过基于稀疏优化的交替方向方法推导基于各向异性TV最小化的CT图像重建算法。实验结果表明，在稀疏角度重建中，本文提出的基于各向异性TV最小化重建算法与各向同性TV最小化重建算法相比，在稀疏性保持良好的基础上，重建精度上存在优势，综合性能方面表现更优异。

Abstract: In many practical applications, due to the data acquisition time, dose, and geometric constraint scanning, only in a limited Angle range, various data is available to acquire. It is the so-called few-view problem. In recent years, the Total Variation (TV) minimization model, using alternating direction method (ADM) in sparse optimization algorithm shows better reconstruction results among these TV-based algorithms. However, Isotropic TV minimization based algorithms for few- view reconstruction has not so good accuracy and there is further improvement to achieve. Aiming at this problem, Anisotropic TV minimization algorithm is proposed in this paper. The algorithm is based on ADM and uses sparse optimization theory. Experimental results demonstrate that the proposed method compared with Isotropic TV minimization algorithm, has higher reconstruction accuracy and a more excellent comprehensive performance.

文章引用: 张 莹 , 王 丹 (2014) 一种基于各向异性总变分最小化的CT图像重建算法。 计算机科学与应用， 4， 240-247. doi: 10.12677/CSA.2014.410033

参考文献

[1] Andersen, A.H. (1989) Algebraic reconstruction in CT from limited views. IEEE Transactions on Medical Imaging, 8, 50-55.

[2] Candes, E., Romberg, J. and Tao, T. (2006) Robust uncertainty principles: Exact signal reconstruction from highly incomplete frequency Information. IEEE Transactions on Information Theory, 52, 489-509.

[3] Donoho, D. (2006) Compressed sensing. IEEE Transactions on Information Theory, 52, 1289-1306.

[4] Candes, E. and Tao, T. (2005) Decoding by linear programming. IEEE Transactions on Information Theory, 59, 1207-1223.

[5] Goldstein, T. and Osher, S. (2009) The split Bregman method for l1 regularized problems. SIAM Journal on Imaging Sciences, 2, 323-343.

[6] Li, C. (2009) An efficient algorithm for total variation regularization with applications to the single pixel camera and compressive sensing. Rice University, Houston.

[7] Vandeghinste, B., Goossens, B., Beenhouwer, J., et al. (2011) Split-Bregman-based sparse-view CT reconstruction. 11th International Meeting on Fully Three-Dimensional Reconstruction in Radiology and Nuclear Medicine, 431-434.

[8] Deng, W. and Yin, W. (2012) On the global and linear convergence of the generalized alternating direction method of multipliers. Rice CAAM Report, 12-14.

[9] Chen, C. and Xu, G. (2013) The linearized split Bregman iterative algorithm and its convergence analysis for robust tomographic image reconstruction. UCLA CAM Report, 13-66.

[10] 喻玲娟, 谢晓春 (2008) 压缩感知理论简介. 数字视频, 12, 16-18．

[11] 傅迎华 (2008) 可压缩传感重构算法与近似QR分解. 计算机应用, 9, 2300-2302.

[12] Chen, Z., Jin, X., Li, L., et al. (2013) A limited-angle CT reconstruction method based on anisotropic TV minimization. Physics in Medicine and Biology, 58, 2119-2141.

分享
Top