Constrained TpV Minimization for Enhanced Exploitation of Gradient Sparsity: Application to CT Image Reconstruction.
Exploiting sparsity in the image gradient magnitude has proved to be an effective means for reducing the sampling rate in the projection view angle in computed tomography (CT). Most of the image reconstruction algorithms, developed for this purpose, solve a nonsmooth convex optimization problem involving the image total variation (TV). The TV seminorm is theℓ1norm of the image gradient magnitude, and reducing theℓ1norm is known to encourage sparsity in its argument. Recently, there has been interest in employing nonconvexℓpquasinorms with 0<p<1 for sparsity exploiting image reconstruction, which is potentially more effective thanℓ1because nonconvexℓpis closer toℓ0-a direct measure of sparsity. This paper develops algorithms for constrained minimization of the totalp-variation (TpV),ℓpof the image gradient. Use of the algorithms is illustrated in the context of breast CT-an imaging modality that is still in the research phase and for which constraints on X-ray dose are extremely tight. The TpV-based image reconstruction algorithms are demonstrated on computer simulated data for exploiting gradient magnitude sparsity to reduce the projection view angle sampling. The proposed algorithms are applied to projection data from a realistic breast CT simulation, where the total X-ray dose is equivalent to two-view digital mammography. Following the simulation survey, the algorithms are then demonstrated on a clinical breast CT data set.
Department of Radiology, University of Chicago, Chicago, IL 60637, USA.
Emil Y
EY
Sidky
Theoretical Division T-5, Los Alamos National Laboratory, Los Alamos, NM 87545, USA.
Rick
R
Chartrand
Department of Radiology, University of California Davis Medical Center, Sacramento, CA 95817, USA.
John M
JM
Boone
Department of Radiology, University of Chicago, Chicago, IL 60637, USA.
Xiaochuan
X
Pan
Computed tomography; X-ray tomography; image reconstruction; iterative algorithms; optimization