Low-dose CT reconstruction via L1 dictionary learning regularization using iteratively reweighted least-squares

被引:32
作者
Zhang, Cheng [1 ,2 ,3 ]
Zhang, Tao [2 ]
Li, Ming [1 ]
Peng, Chengtao [4 ]
Liu, Zhaobang [1 ]
Zheng, Jian [1 ]
机构
[1] Chinese Acad Sci, Suzhou Inst Biomed Engn & Technol, Suzhou 215163, Peoples R China
[2] Chinese Acad Sci, Changchun Inst Opt Fine Mech & Phys, Changchun 130033, Peoples R China
[3] Univ Chinese Acad Sci, Beijing 100049, Peoples R China
[4] Univ Sci & Technol China, Dept Elect Sci & Technol, Hefei 230061, Peoples R China
基金
中国国家自然科学基金;
关键词
Dictionary learning; Image reconstruction; L-1-norm; Iteratively reweighted least squares; IMAGE-RECONSTRUCTION; SIGNAL RECOVERY;
D O I
10.1186/s12938-016-0193-y
中图分类号
R318 [生物医学工程];
学科分类号
0831 ;
摘要
Background: In order to reduce the radiation dose of CT (computed tomography), compressed sensing theory has been a hot topic since it provides the possibility of a high quality recovery from the sparse sampling data. Recently, the algorithm based on DL (dictionary learning) was developed to deal with the sparse CT reconstruction problem. However, the existing DL algorithm focuses on the minimization problem with the L-2-norm regularization term, which leads to reconstruction quality deteriorating while the sampling rate declines further. Therefore, it is essential to improve the DL method to meet the demand of more dose reduction. Methods: In this paper, we replaced the L-2-norm regularization term with the L-1-norm one. It is expected that the proposed L-1-DL method could alleviate the over-smoothing effect of the L-2-minimization and reserve more image details. The proposed algorithm solves the L-1-minimization problem by a weighting strategy, solving the new weighted L-2-minimization problem based on IRLS (iteratively reweighted least squares). Results: Through the numerical simulation, the proposed algorithm is compared with the existing DL method (adaptive dictionary based statistical iterative reconstruction, ADSIR) and other two typical compressed sensing algorithms. It is revealed that the proposed algorithm is more accurate than the other algorithms especially when further reducing the sampling rate or increasing the noise. Conclusion: The proposed L-1-DL algorithm can utilize more prior information of image sparsity than ADSIR. By transforming the L-2-norm regularization term of ADSIR with the L-1-norm one and solving the L-1-minimization problem by IRLS strategy, L-1-DL could reconstruct the image more exactly.
引用
收藏
页数:21
相关论文
共 19 条
[1]   K-SVD: An algorithm for designing overcomplete dictionaries for sparse representation [J].
Aharon, Michal ;
Elad, Michael ;
Bruckstein, Alfred .
IEEE TRANSACTIONS ON SIGNAL PROCESSING, 2006, 54 (11) :4311-4322
[2]   Robust uncertainty principles:: Exact signal reconstruction from highly incomplete frequency information [J].
Candès, EJ ;
Romberg, J ;
Tao, T .
IEEE TRANSACTIONS ON INFORMATION THEORY, 2006, 52 (02) :489-509
[3]   Near-optimal signal recovery from random projections: Universal encoding strategies? [J].
Candes, Emmanuel J. ;
Tao, Terence .
IEEE TRANSACTIONS ON INFORMATION THEORY, 2006, 52 (12) :5406-5425
[4]   Projected Cancer Risks From Computed Tomographic Scans Performed in the United States in 2007 [J].
de Gonzalez, Amy Berrington ;
Mahesh, Mahadevappa ;
Kim, Kwang-Pyo ;
Bhargavan, Mythreyi ;
Lewis, Rebecca ;
Mettler, Fred ;
Land, Charles .
ARCHIVES OF INTERNAL MEDICINE, 2009, 169 (22) :2071-2077
[5]   Statistical image reconstruction for polyenergetic X-ray computed tomography [J].
Elbakri, IA ;
Fessler, JA .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 2002, 21 (02) :89-99
[6]  
Gunawan Dadang, 2010, ADV COMP CONTR TEL T, DOI [10.1109/ACT.2010.31, DOI 10.1109/ACT.2010.31]
[7]   SHANNON SAMPLING THEOREM - ITS VARIOUS EXTENSIONS AND APPLICATIONS - TUTORIAL REVIEW [J].
JERRI, AJ .
PROCEEDINGS OF THE IEEE, 1977, 65 (11) :1565-1596
[8]   Accelerated iterative transmission CT reconstruction using an ordered subsets convex algorithm [J].
Kamphuis, C ;
Beekman, FJ .
IEEE TRANSACTIONS ON MEDICAL IMAGING, 1998, 17 (06) :1101-1105
[9]   Sparse representations for limited data tomography [J].
Liao, Hstau Y. ;
Sapiro, Guillermo .
2008 IEEE INTERNATIONAL SYMPOSIUM ON BIOMEDICAL IMAGING: FROM NANO TO MACRO, VOLS 1-4, 2008, :1375-+
[10]   Learning multiscale sparse representations for image and video restoration [J].
Mairal, Julien ;
Sapiro, Guillermo ;
Elad, Michael .
MULTISCALE MODELING & SIMULATION, 2008, 7 (01) :214-241