Robust Low-Dose CT Sinogram Preprocessing via Exploiting Noise-Generating Mechanism

被引:63
作者
Xie, Qi [1 ,2 ]
Zeng, Dong [3 ,4 ]
Zhao, Qian [1 ,2 ]
Meng, Deyu [1 ,2 ]
Xu, Zongben [1 ,2 ]
Liang, Zhengrong [5 ,6 ]
Ma, Jianhua [3 ,4 ]
机构
[1] Xi An Jiao Tong Univ, Sch Math & Stat, Xian 710049, Shaanxi, Peoples R China
[2] Xi An Jiao Tong Univ, Minist Educ, Key Lab Intelligent Networks & Network Secur, Xian 710049, Shaanxi, Peoples R China
[3] Southern Med Univ, Sch Biomed Engn, Guangzhou 510515, Guangdong, Peoples R China
[4] Southern Med Univ, Guangzhou Key Lab Med Radiat Imaging & Detect Tec, Guangzhou 510515, Guangdong, Peoples R China
[5] SUNY Stony Brook, Dept Radiol, Stony Brook, NY 11794 USA
[6] SUNY Stony Brook, Dept Biomed Engn, Stony Brook, NY 11794 USA
基金
美国国家卫生研究院; 中国国家自然科学基金;
关键词
Computed tomography; noise modeling; maximum a posteriori (MAP); statistical model; regularization; RAY COMPUTED-TOMOGRAPHY; TOTAL-VARIATION MINIMIZATION; IMAGE-RECONSTRUCTION; RESTORATION; REGULARIZATION; REDUCTION; SHRINKAGE; ALGORITHM;
D O I
10.1109/TMI.2017.2767290
中图分类号
TP39 [计算机的应用];
学科分类号
081203 ; 0835 ;
摘要
Computed tomography (CT) image recovery from low-mAs acquisitions without adequate treatment is always severely degraded due to a number of physical factors. In this paper, we formulate the low-dose CT sinogram preprocessing as a standard maximum a posteriori (MAP) estimation, which takes full consideration of the statistical properties of the two intrinsic noise sources in low-dose CT, i.e., the X-ray photon statistics and the electronic noise background. In addition, instead of using a general image prior as found in the traditional sinogram recovery models, we design a new prior formulation to more rationally encode the piecewise-linear configurations underlying a sinogram than previously used ones, like the TV prior term. As compared with the previous methods, especially the MAP-based ones, both the likelihood/loss and prior/regularization terms in the proposed model are ameliorated in a more accurate manner and better comply with the statistical essence of the generation mechanism of a practical sinogram. We further construct an efficient alternating direction method of multipliers algorithm to solve the proposed MAP framework. Experiments on simulated and real low-dose CT data demonstrate the superiority of the proposed method according to both visual inspection and comprehensive quantitative performance evaluation.
引用
收藏
页码:2487 / 2498
页数:12
相关论文
共 45 条
[1]  
[Anonymous], FOUND TRENDS MACH LE
[2]  
[Anonymous], 2010, 100920105055 ARXIV
[3]  
[Anonymous], 2004, CONVEX OPTIMIZATION
[4]  
[Anonymous], P IEEE NSS MIC
[5]  
[Anonymous], 2009, Computed Tomography: Principles, Design Artifacts and Recent Advances
[6]   A Fast Iterative Shrinkage-Thresholding Algorithm for Linear Inverse Problems [J].
Beck, Amir ;
Teboulle, Marc .
SIAM JOURNAL ON IMAGING SCIENCES, 2009, 2 (01) :183-202
[7]   Image recovery via total variation minimization and related problems [J].
Chambolle, A ;
Lions, PL .
NUMERISCHE MATHEMATIK, 1997, 76 (02) :167-188
[8]   The direct extension of ADMM for multi-block convex minimization problems is not necessarily convergent [J].
Chen, Caihua ;
He, Bingsheng ;
Ye, Yinyu ;
Yuan, Xiaoming .
MATHEMATICAL PROGRAMMING, 2016, 155 (1-2) :57-79
[9]   The statistical sinogram smoothing via adaptive-weighted total variation regularization for low-dose X-ray CT [J].
Cui, Xueying ;
Gui, Zhiguo ;
Zhang, Quan ;
Liu, Yi ;
Ma, Ruifen .
OPTIK, 2014, 125 (18) :5352-5356
[10]  
Fessler JA., 2000, Handbook of Medical Imaging, V2, P1, DOI [DOI 10.1117/3.831079.CH1, DOI 10.1117/3.831079]