Morteza Salehjahromi1, Qian Wang1, Yanbo Zhang2, Lars A Gjesteby3, Dan Harrison3, Ge Wang3, Peter M Edic4, Hengyong Yu1. 1. Department of Electrical and Computer Engineering, University of Massachusetts Lowell, Lowell, MA, 01854, USA. 2. Ping An Technology, U.S. Research Lab, Palo Alto, CA, 94306, USA. 3. Department of Biomedical Engineering, Rensselaer Polytechnic Institute, Troy, NY, 12180, USA. 4. GE Research, Niskayuna, NY, 12309, USA.
Abstract
PURPOSE: In computed tomography (CT), miscalibrated or imperfect detector elements produce stripe artifacts in the sinogram. The stripe artifacts in Radon space are responsible for concentric ring artifacts in the reconstructed images. In this work, a novel optimization model is proposed to remove the ring artifacts in an iterative reconstruction procedure. METHOD: In the proposed optimization model, a novel ring total variation (RTV) regularization is developed to penalize the ring artifacts in the image domain. Moreover, to correct the sinogram, a new correcting vector is proposed to compensate for malfunctioning of detectors in the projection domain. The optimization problem is solved by using the alternating minimization scheme (AMS). In each iteration, the fidelity term along with the RTV regularization is solved using the alternating direction method of multipliers (ADMM) to find the image, and then the correcting coefficient vector is updated for certain detectors according to the obtained image. Because the sinogram and the image are simultaneously updated, the proposed method basically performs in both image and sinogram domains. RESULTS: The proposed method is evaluated using both simulated and physical phantom datasets containing different ring artifact patterns. In the simulated datasets, the Shepp-Logan phantom, a real chest scan image and a noisy low-contrast phantom are considered for the performance evaluation of our method. We compare the quantitative root mean square error (RMSE) and structural similarity (SSIM) results of our algorithm with wavelet-Fourier sinogram filtering method by Munch et al., the ring artifact reduction method by Brun et al., and the TV-based ring correction method by Paleo and Mirone. Our proposed method is also evaluated using a physical phantom dataset where strong ring artifacts are manifest due to the miscalibration of a large number of detectors. Our proposed method outperforms the competing methods in terms of both qualitative and quantitative evaluation results. CONCLUSION: The experimental results in both simulated and physical phantom datasets show that the proposed method achieves the state-of-the-art ring artifact reduction performance in terms of RMSE, SSIM, and subjective visual quality.
PURPOSE: In computed tomography (CT), miscalibrated or imperfect detector elements produce stripe artifacts in the sinogram. The stripe artifacts in Radon space are responsible for concentric ring artifacts in the reconstructed images. In this work, a novel optimization model is proposed to remove the ring artifacts in an iterative reconstruction procedure. METHOD: In the proposed optimization model, a novel ring total variation (RTV) regularization is developed to penalize the ring artifacts in the image domain. Moreover, to correct the sinogram, a new correcting vector is proposed to compensate for malfunctioning of detectors in the projection domain. The optimization problem is solved by using the alternating minimization scheme (AMS). In each iteration, the fidelity term along with the RTV regularization is solved using the alternating direction method of multipliers (ADMM) to find the image, and then the correcting coefficient vector is updated for certain detectors according to the obtained image. Because the sinogram and the image are simultaneously updated, the proposed method basically performs in both image and sinogram domains. RESULTS: The proposed method is evaluated using both simulated and physical phantom datasets containing different ring artifact patterns. In the simulated datasets, the Shepp-Logan phantom, a real chest scan image and a noisy low-contrast phantom are considered for the performance evaluation of our method. We compare the quantitative root mean square error (RMSE) and structural similarity (SSIM) results of our algorithm with wavelet-Fourier sinogram filtering method by Munch et al., the ring artifact reduction method by Brun et al., and the TV-based ring correction method by Paleo and Mirone. Our proposed method is also evaluated using a physical phantom dataset where strong ring artifacts are manifest due to the miscalibration of a large number of detectors. Our proposed method outperforms the competing methods in terms of both qualitative and quantitative evaluation results. CONCLUSION: The experimental results in both simulated and physical phantom datasets show that the proposed method achieves the state-of-the-art ring artifact reduction performance in terms of RMSE, SSIM, and subjective visual quality.
Keywords:
alternating direction method of multipliers; alternating minimization scheme; computed tomography; iterative reconstruction; ring artifact; ring total variation