UC Davis

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 ℓ₁ norm of the image gradient magnitude, and reducing the ℓ₁ norm is known to encourage sparsity in its argument. Recently, there has been interest in employing nonconvex ℓ_p quasinorms with 0<p<1 for sparsity exploiting image reconstruction, which is potentially more effective than ℓ₁ because nonconvex ℓ_p is closer to ℓ₀-a direct measure of sparsity. This paper develops algorithms for constrained minimization of the total p-variation (TpV), ℓ_p 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 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.