Low-dose CT reconstruction using spatially encoded nonlocal penalty

被引:32
作者
Kim, Kyungsang
El Fakhri, Georges
Li, Quanzheng [1 ]
机构
[1] Massachusetts Gen Hosp, Gordon Ctr Med Imaging, 125 Nashua St 6th Floor,Suite 660, Boston, MA 02114 USA
关键词
Grand Challenge; Low-dose CT reconstruction; spatially encoded nonlocal penalty; IMAGE-RECONSTRUCTION; ORDERED SUBSETS; REGULARIZATION; ALGORITHMS; QUALITY; ART;
D O I
10.1002/mp.12523
中图分类号
R8 [特种医学]; R445 [影像诊断学];
学科分类号
1002 ; 100207 ; 1009 ;
摘要
Purpose: Computed tomography (CT) is one of the most used imaging modalities for imaging both symptomatic and asymptomatic patients. However, because of the high demand for lower radiation dose during CT scans, the reconstructed image can suffer from noise and artifacts due to the trade-off between the image quality and the radiation dose. The purpose of this paper is to improve the image quality of quarter dose images and to select the best hyperparameters using the regular dose image as ground truth. Methods: We first generated the axially stacked two-dimensional sinograms from the multislice raw projections with flying focal spots using a single slice rebinning method, which is an axially approximate method to provide simple implementation and efficient memory usage. To improve the image quality, a cost function containing the Poisson log-likelihood and spatially encoded nonlocal penalty is proposed. Specifically, an ordered subsets separable quadratic surrogates (OS-SQS) method for the log-likelihood is exploited and the patch-based similarity constraint with a spatially variant factor is developed to reduce the noise significantly while preserving features. Furthermore, we applied the Nesterov's momentum method for acceleration and the diminishing number of subsets strategy for noise consistency. Fast nonlocal weight calculation is also utilized to reduce the computational cost. Results: Datasets given by the Low Dose CT Grand Challenge were used for the validation, exploiting the training datasets with the regular and quarter dose data. The most important step in this paper was to fine-tune the hyperparameters to provide the best image for diagnosis. Using the regular dose filtered back-projection (FBP) image as ground truth, we could carefully select the hyperparameters by conducting a bias and standard deviation study, and we obtained the best images in a fixed number of iterations. We demonstrated that the proposed method with well selected hyperparameters improved the image quality using quarter dose data. The quarter dose proposed method was compared with the regular dose FBP, quarter dose FBP, and quarter dose l(1)-based 3-D TV method. We confirmed that the quarter dose proposed image was comparable to the regular dose FBP image and was better than images using other quarter dose methods. The reconstructed test images of the accreditation (ACR) CT phantom and 20 patients data were evaluated by radiologists at the Mayo clinic, and this method was awarded first place in the Low Dose CT Grand Challenge. Conclusion: We proposed the iterative CT reconstruction method using a spatially encoded nonlocal penalty and ordered subsets separable quadratic surrogates with the Nesterov's momentum and diminishing number of subsets. The results demonstrated that the proposed method with fine-tuned hyperparameters can significantly improve the image quality and provide accurate diagnostic features at quarter dose. The performance of the proposed method should be further improved for small lesions, and a more thorough evaluation using additional clinical data is required in the future. (C) 2017 American Association of Physicists in Medicine
引用
收藏
页码:e376 / e390
页数:15
相关论文
共 38 条
  • [1] SIMULTANEOUS ALGEBRAIC RECONSTRUCTION TECHNIQUE (SART) - A SUPERIOR IMPLEMENTATION OF THE ART ALGORITHM
    ANDERSEN, AH
    KAK, AC
    [J]. ULTRASONIC IMAGING, 1984, 6 (01) : 81 - 94
  • [2] [Anonymous], INT M FULL 3 DIM IM
  • [3] [Anonymous], IMAGE PROCESSING ALG
  • [4] A non-local algorithm for image denoising
    Buades, A
    Coll, B
    Morel, JM
    [J]. 2005 IEEE COMPUTER SOCIETY CONFERENCE ON COMPUTER VISION AND PATTERN RECOGNITION, VOL 2, PROCEEDINGS, 2005, : 60 - 65
  • [5] Chambolle A, 2004, J MATH IMAGING VIS, V20, P89
  • [6] Adaptive wavelet thresholding for image denoising and compression
    Chang, SG
    Yu, B
    Vetterli, M
    [J]. IEEE TRANSACTIONS ON IMAGE PROCESSING, 2000, 9 (09) : 1532 - 1546
  • [7] Deterministic edge-preserving regularization in computed imaging
    Charbonnier, P
    BlancFeraud, L
    Aubert, G
    Barlaud, M
    [J]. IEEE TRANSACTIONS ON IMAGE PROCESSING, 1997, 6 (02) : 298 - 311
  • [8] Image denoising by sparse 3-D transform-domain collaborative filtering
    Dabov, Kostadin
    Foi, Alessandro
    Katkovnik, Vladimir
    Egiazarian, Karen
    [J]. IEEE TRANSACTIONS ON IMAGE PROCESSING, 2007, 16 (08) : 2080 - 2095
  • [9] Fast nonlocal filtering applied to electron cryomicroscopy
    Darbon, Jerome
    Cunha, Alexandre
    Chan, Tony F.
    Osher, Stanley
    Jensen, Grant J.
    [J]. 2008 IEEE INTERNATIONAL SYMPOSIUM ON BIOMEDICAL IMAGING: FROM NANO TO MACRO, VOLS 1-4, 2008, : 1331 - +
  • [10] First-order methods of smooth convex optimization with inexact oracle
    Devolder, Olivier
    Glineur, Francois
    Nesterov, Yurii
    [J]. MATHEMATICAL PROGRAMMING, 2014, 146 (1-2) : 37 - 75