GPU-based Low Dose CT Reconstruction via Edge-preserving Total Variation Regularization

Abstract

High radiation dose in CT scans increases a lifetime risk of cancer and has become a major clinical concern. Recently, iterative reconstruction algorithms with Total Variation (TV) regularization have been developed to reconstruct CT images from highly undersampled data acquired at low mAs levels in order to reduce the imaging dose. Nonetheless, TV regularization may lead to over-smoothed images and lost edge information. To solve this problem, in this work we develop an iterative CT reconstruction algorithm with edge-preserving TV regularization to reconstruct CT images from highly undersampled data obtained at low mAs levels. The CT image is reconstructed by minimizing an energy consisting of an edge-preserving TV norm and a data fidelity term posed by the x-ray projections. The edge-preserving TV term is proposed to preferentially perform smoothing only on non-edge part of the image in order to avoid over-smoothing, which is realized by introducing a penalty weight to the original total variation norm. Our iterative algorithm is implemented on GPU to improve its speed. We test our reconstruction algorithm on a digital NCAT phantom, a physical chest phantom, and a Catphan phantom. Reconstruction results from a conventional FBP algorithm and a TV regularization method without edge preserving penalty are also presented for comparison purpose. The experimental results illustrate that both TV-based algorithm and our edge-preserving TV algorithm outperform the conventional FBP algorithm in suppressing the streaking artifacts and image noise under the low dose context. Our edge-preserving algorithm is superior to the TV-based algorithm in that it can preserve more information of fine structures and therefore maintain acceptable spatial resolution.

21 pages, 6 figures, 2 tables

References in corpus (6)

• Accurate image reconstruction from few-views and limited-angle data in divergent-beam CT
• Implementation and evaluation of various demons deformable image registration algorithms on GPU
• Development of a GPU-based Monte Carlo dose calculation code for coupled electron-photon transport
• GPU-based ultra fast IMRT plan optimization
• GPU-based ultra fast dose calculation using a finite pencil beam model
• GPU-based ultra-fast direct aperture optimization for online adaptive radiation therapy

Cited by in corpus (26)

• Low Dose CT Image Denoising Using a Generative Adversarial Network with Wasserstein Distance and Perceptual Loss
• Sharpness-aware Low dose CT denoising using conditional generative adversarial network
• Structure-sensitive Multi-scale Deep Neural Network for Low-Dose CT Denoising
• CTformer: Convolution-free Token2Token Dilated Vision Transformer for Low-dose CT Denoising
• A comprehensive study on the relationship between image quality and imaging dose in low-dose cone beam CT
• Learning Sub-Sampling and Signal Recovery with Applications in Ultrasound Imaging
• A Variational Reconstruction Method for Undersampled Dynamic X-ray Tomography based on Physical Motion Models
• Neural Networks-based Regularization for Large-Scale Medical Image Reconstruction
• Systematic Review on Learning-based Spectral CT
• Noise2Context: Context-assisted Learning 3D Thin-layer Low Dose CT Without Clean Data
• A function space framework for structural total variation regularization with applications in inverse problems
• A very fast iterative algorithm for TV-regularized image reconstruction with applications to low-dose and few-view CT
• Sparse dynamic tomography. A shearlet-based approach for iodine perfusion in plant stems
• Medical Imaging Synthesis using Deep Learning and its Clinical Applications: A Review
• Unsupervised/Semi-supervised Deep Learning for Low-dose CT Enhancement
• TransCT: Dual-path Transformer for Low Dose Computed Tomography
• ROI-Wise Material Decomposition in Spectral Photon-Counting CT
• Projected Newton Method for noise constrained Tikhonov regularization
• TED-net: Convolution-free T2T Vision Transformer-based Encoder-decoder Dilation network for Low-dose CT Denoising
• A Re-weighted Joint Spatial-Radon Domain CT Image Reconstruction Model for Metal Artifact Reduction
• Regularization with optimal space-time priors
• Learning Priors in High-frequency Domain for Inverse Imaging Reconstruction
• Missing Cone Artifacts Removal in ODT using Unsupervised Deep Learning in Projection Domain
• Learning to solve TV regularized problems with unrolled algorithms
• Comparison of Sinogram-based Iterative Reconstruction with Compressed Sensing Techniques in X-ray CT
• An Algorithmic Theory of Dependent Regularizers, Part 1: Submodular Structure