Statistical iterative algorithm based on adaptive weighted total variation for low-dose CT

doi:10.11772/j.issn.1001-9081.2016.10.2916

Journal of Computer Applications ›› 2016, Vol. 36 ›› Issue (10): 2916-2921.DOI: 10.11772/j.issn.1001-9081.2016.10.2916

Previous Articles Next Articles

Statistical iterative algorithm based on adaptive weighted total variation for low-dose CT

HE Lin¹, ZHANG Quan¹, SHANGGUAN Hong¹, ZHANG Wen¹, ZHANG Pengcheng¹, LIU Yi¹, GUI Zhiguo^1,2

1. National Key Laboratory for Electronic Measurement Technology (North University of China), Taiyuan Shanxi 030051, China;
2. Key Laboratory of Instrumentation Science and Dynamic Measurement, Ministry of Education (North University of China), Taiyuan Shanxi 030051, China

Received:2016-03-16 Revised:2016-06-14 Published:2016-10-10
Supported by:
BackgroundThis work is partially supported by the National Natural Science Foundation of China (61271357), the National High Technology Instruments and Equipment Development Program of China (2014YQ24044508), the Natural Science Foundation of Shanxi Province (2015011046), the Science Foundation Plan of North University of China in 2013.

自适应加权全变分的低剂量CT统计迭代算法

何琳¹, 张权¹, 上官宏¹, 张文¹, 张鹏程¹, 刘祎¹, 桂志国^1,2

1. 电子测试技术国家重点实验室(中北大学), 太原 030051;
2. 仪器科学与动态测试教育部重点实验室(中北大学), 太原 030051

通讯作者: 桂志国,E-mail:gzgtg@163.com
作者简介:何琳(1991—),女,山西运城人,硕士研究生,主要研究方向:图像处理与重建;张权(1974—),男,山西大同人,副教授,博士,主要研究方向:图像处理、科学可视化;上官宏(1988—),女,山西临汾人,博士研究生,主要研究方向:图像处理、医学图像重建;张文(1992—),男,湖北孝感人,硕士研究生,主要研究方向:图像处理与重建;张鹏程(1984—),男,内蒙古巴彦淖尔人,讲师,博士,主要研究方向:剂量计算、方案优化;刘祎(1987—),女,河南睢县人,讲师,博士,主要研究方向:图像处理、医学图像重建;桂志国(1972—),男,天津蓟县人,教授,博士,主要研究方向:信号与信息处理、图像处理和识别、图像重建。
基金资助:
国家自然科学基金资助项目（61271357）；国家重大科学仪器设备开发专项（2014YQ24044508）；山西省自然科学基金资助项目（2015011046）；中北大学2013年校科学基金资助项目。

Abstract

Abstract: Concerning the streak artifacts and impulse noise of the Low-Dose Computed Tomography (LDCT) reconstructed images, a statistical iterative reconstruction method based on adaptive weighted Total Variation (TV) for LDCT was presented. Considering the shortage that traditional TV may bring staircase effect while suppressing streak artifacts, an adaptive weighted TV model that combined the weighting factor based on weighted variation and TV model was proposed. Then, the new model was applied to the Penalized Weighted Least Square (PWLS). Different areas of the image were processed with different de-noising intensities, so as to achieve a good effect of noise suppression and edge preservation. The Shepp-Logan model and the digital pelvis phantom were used to test the effectiveness of the proposed algorithm. Experimental results show that the proposed method has smaller Normalized Mean Square Distance (NMSD) and Normal Average Absolute Distance (NAAD) in the two experiment images, compared with the Filtered Back Projection (FBP), PWLS, PWLS-Median Prior (PWLS-MP) and PWLS-TV algorithms. Meanwhile, the proposed method get Peak Signal-To-Noise Ratio (PSNR) of 40.91 dB and 42.25 dB respectively. Experimental results show that the proposed algorithm can well preserve image details and edges, while eliminating streak artifacts effectively.

Key words: Low-Dose Computed Tomography (LDCT), statistical iterative reconstruction, Penalized Weighted Least Square (PWLS), Total Variation (TV), weighted variation

摘要： 针对低剂量计算机断层扫描（LDCT）重建图像时出现条形伪影和脉冲噪声的现象，提出一种自适应加权全变分的LDCT统计迭代重建算法。该算法克服了传统全变分（TV）算法在去除条形伪影的同时引入阶梯效应的缺点，把基于加权方差的加权因子与TV模型相结合提出自适应加权全变分模型，然后再把新模型应用到惩罚加权最小二乘（PWLS）重建算法中，这样就可以对图像的不同区域进行不同强度的去噪，从而取得噪声抑制和边缘保持的良好效果。采用Shepp-Logan模型和数字骨盆体模来验证算法的有效性，实验结果表明，所提算法的归一化均方距离和归一化平均绝对距离均比滤波反投影（FBP）、PWLS、惩罚加权最小二乘的中值先验（PWLS-MP）以及惩罚加权最小二乘的全变分（PWLS-TV）算法的值小，且可分别获得40.91 dB和42.25 dB的峰值信噪比。实验结果表明，该算法重建出的图像在有效去除条形伪影的同时对图像的边缘和细节起到很好的保护作用。

关键词: 低剂量计算机断层扫描, 统计迭代重建, 惩罚加权最小二乘, 全变分, 加权方差

CLC Number:

HE Lin, ZHANG Quan, SHANGGUAN Hong, ZHANG Wen, ZHANG Pengcheng, LIU Yi, GUI Zhiguo. Statistical iterative algorithm based on adaptive weighted total variation for low-dose CT[J]. Journal of Computer Applications, 2016, 36(10): 2916-2921.

何琳, 张权, 上官宏, 张文, 张鹏程, 刘祎, 桂志国. 自适应加权全变分的低剂量CT统计迭代算法[J]. 计算机应用, 2016, 36(10): 2916-2921.

References

[1] LUI D, CAMERON A, MODHAFAR A, et al. Low-dose computed tomography via spatially adaptive Monte-Carlo reconstruction [J]. Computerized Medical Imaging and Graphics, 2013, 37(7/8): 438-449.
[2] RUST G F, AURICH V, REISER M. Noise dose reduction and image improvements in screening virtual colonoscopy with tube currents of 20 mAs with nonlinear Gaussian filter chains [C]//Proceedings of SPIE 4683. Bellingham, WA: SPIE Press, 2002: 186-197.
[3] CHEN Y, LI Y, YU W, et al. Joint-map tomographic reconstruction with patch similarity based mixture prior model [J]. Multiscale Modeling and Simulation, 2011, 9(4): 1399-1419.
[4] ZHANG H, MA J H, WANG J, et al. Statistical image reconstruction for low-dose CT using nonlocal means-based regularization [J]. Computerized Medical Imaging and Graphics, 2014, 38: 423-435.
[5] ZHANG H, MA J H, WANG J, et al. Statistical image reconstruction for low-dose CT using nonlocal means-based regularization, Part II: an adaptive approach [J]. Computerized Medical Imaging and Graphics, 2015, 43: 26-35.
[6] WANG J, LI T F, LU H B, et al. Penalized weighted least-squares approach to sinogram noise reduction and image reconstruction for low-dose X-ray computed tomography [J]. IEEE Transactions on Medical Imaging, 2006, 25(10): 1272-1283.
[7] XU Q, YU H Y, MOU X Q, et al. Low-dose X-ray CT reconstruction via dictionary learning [J]. IEEE Transactions on Medical Imaging, 2012, 31(9): 1682-1697.
[8] TIAN Z, JIA X, YUAN K H, et al. Low-dose CT reconstruction via edge-preserving total variation regularization [J]. Physics in Medicine and Biology, 2011, 56(18): 5949-5967.
[9] ZHU Y, ZHAO M L, ZHAO Y S, et al. Noise reduction with low-dose CT data based on a modified ROF model [J]. Optics Express, 2012, 20(16): 17987-18004.
[10] CHAO S M, TSAI D M. An improved anisotropic diffusion model for detail- and edge-preserving smoothing [J]. Pattern Recognition Letters, 2010, 31(13): 2012-2023.
[11] LOU Y F, ZENG T Y, OSHER S, et al. A weighted difference of anisotropic and isotropic total variation model for image processing [J]. SIAM Journal on Imaging Sciences, 2015, 8(3): 1798-1823.
[12] TANG J, NETT B E, CHEN G H. Performance comparison between total variation-based compressed sensing and statistical iterative reconstruction algorithms [J]. Physics in Medicine & Biology, 2009, 54(19): 5781-5804.
[13] XU Q, MOU X, WANG G, et al. Statistical interior tomography [J]. IEEE Transactions on Medical Imaging, 2011, 30(5): 1116-1128.
[14] LI T F, LI X, WANG J, et al. Nonlinear sinogram smoothing for low-dose X-ray CT [J]. IEEE Transactions on Nuclear Science, 2004, 51(5): 2505-2513.
[15] FESSLER J A. Penalized weighted least-squares image reconstruction for positron emission tomography [J]. IEEE Transactions on Medical Imaging, 1994, 13(2): 290-300.
[16] RUDIN L, OSHER S, FATEMI E. Nonlinear total variation based noise removal algorithms [J]. Physica D — Nonlinear Phenomena, 1992, 60(1/2/3/4): 259-268.
[17] KIANI M, SEID H P. State estimation of nonlinear dynamic systems using weighted variance-based adaptive particle swarm optimization [J]. Applied Soft Computing, 2015, 34(C): 1-17.
[18] LU X Q, SUN Y, YUAN Y, et al. Image reconstruction by an alternating minimization [J]. Neuro Computing, 2011, 74(5): 661-670.
[19] ELBKRI I, FESSLER J. Statistical image reconstruction for polyene-rgetic X-ray computed tomography [J]. IEEE Transactions on Medical Imaging, 2002, 21: 89-99.
[20] HSIAO T, RANGARAJAN A, GINDI G. A new convex edge-preserving median prior with applications to tomography [J]. IEEE Transactions on Medical Imaging, 2003, 22(5): 580-585.
[21] 钱姗姗, 黄静, 马建华, 等. 基于投影数据非单调性全变分恢复的低剂量CT重建 [J]. 电子学报, 2011, 39(7): 1702-1707.(QIAN S S, HUANG J, MA J H, et al. Nonmonotone total variation minimization based projection restoration for low-dose CT reconstruction [J]. Acta Electronica Sinica, 2011, 39(7): 1702-1707.)
[22] 王丽艳, 韦志辉, 低剂量CT 的线性Bregman迭代重建算法[J]. 电子与信息学报, 2013, 35(10): 2418-2424.(WANG L Y, WEI Z H. Linearized bregman iterations for low-dose CT reconstruction [J]. Journal of Electronics and Information, 2013, 35(10): 2418-2424.)

Statistical iterative algorithm based on adaptive weighted total variation for low-dose CT

自适应加权全变分的低剂量CT统计迭代算法

PDF

Knowledge

Abstract

Cite this article

share this article

References

Related Articles 15

Recommended Articles

Metrics

[1]	Chenyu GE, Liang DONG, Yikun XU, Yi CHANG, Hongming ZHANG. Global-scale radar data restoration algorithm based on total variation and low-rank group sparsity [J]. Journal of Computer Applications, 2021, 41(11): 3353-3361.
[2]	YANG Wenxia, WANG Meng, ZHANG Liang. Semantic face image inpainting based on U-Net with dense blocks [J]. Journal of Computer Applications, 2020, 40(12): 3651-3657.
[3]	XIONG Xinglong, YANG Lixiang, MA Yuzhao, ZHUANG Zibo. Alerting algorithm of low-level wind shear based on fuzzy C-means [J]. Journal of Computer Applications, 2018, 38(3): 655-660.
[4]	GAO Jingzhi, LIU Yi, BAI Xu, ZHANG Quan, GUI Zhiguo. Stationary wavelet domain deep residual convolutional neural network for low-dose computed tomography image estimation [J]. Journal of Computer Applications, 2018, 38(12): 3584-3590.
[5]	ZHAN Yi, LI Meng. Variable exponent variational model for image interpolation [J]. Journal of Computer Applications, 2017, 37(7): 2067-2070.
[6]	YAN Honghai, PU Fangling, XU Xin. Video super resolution method based on structure tensor [J]. Journal of Computer Applications, 2016, 36(7): 1944-1948.
[7]	ZHAO Mingchao, CHEN Zhibin, WEN Youwei. Parallel computation for image denoising via total variation dual model on GPU [J]. Journal of Computer Applications, 2016, 36(5): 1228-1231.
[8]	LIU Jinlong, XIONG Chengyi, GAO Zhirong, ZHOU Cheng, WANG Shuxian. Image compressive sensing reconstruction via total variation and adaptive low-rank regularization [J]. Journal of Computer Applications, 2016, 36(1): 233-237.
[9]	HE Lin, ZHANG Quan, SHANGGUAN Hong, ZHANG Fang, ZHANG Pengcheng, LIU Yi, SUN Weiya, GUI Zhiguo. Adaptive total generalized variation denoising algorithm for low-dose CT images [J]. Journal of Computer Applications, 2016, 36(1): 243-247.
[10]	WU Yulian FENG Xiangchu. Multiplicative noise removal via adaptive kernel regression and total variation [J]. Journal of Computer Applications, 2013, 33(09): 2592-259.
[11]	HU Xuegang LI Yu. Improved image poisson denoising model based on fractional variation [J]. Journal of Computer Applications, 2013, 33(04): 1100-1102.
[12]	LIU Xiaoguang GAO Xingbao ZHOU Dongmei. l2-total variation image restoration based on subspace optimization [J]. Journal of Computer Applications, 2013, 33(04): 1112-1114.
[13]	FAN Bo YANG Xiaomei HU Xuezhu. Super-resolution image reconstruction algorithms based on compressive sensing [J]. Journal of Computer Applications, 2013, 33(02): 480-483.
[14]	HU Wenjin LI Zhanming. Numerical method based on total variation wavelet inpainting [J]. Journal of Computer Applications, 2013, 33(01): 28-30.
[15]	JDUAN Ji-zhong ZHANG Li-yi LIU Yu SUN Yun-shan. Improved MPEG-2 video coding scheme based on compressed sensing [J]. Journal of Computer Applications, 2012, 32(12): 3411-3414.