Constrained [Formula: see text] 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 [Formula: see text] norm of the image gradient magnitude, and reducing the [Formula: see text] norm is known to encourage sparsity in its argument. Recently, there has been interest in employing nonconvex [Formula: see text] quasinorms with [Formula: see text] for sparsity exploiting image reconstruction, which is potentially more effective than [Formula: see text] because nonconvex [Formula: see text] is closer to [Formula: see text] —a direct measure of sparsity. This paper develops algorithms for constrained minimization of the total [Formula: see text]-variation [Formula: see text] , [Formula: see text] of 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 [Formula: see text]-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.