Constrained Total Generalized p-Variation Minimization for Few-View X-Ray Computed Tomography Image Reconstruction

Abstract

Total generalized variation (TGV)-based computed tomography (CT) image reconstruction, which utilizes high-order image derivatives, is superior to total variation-based methods in terms of the preservation of edge information and the suppression of unfavorable staircase effects. However, conventional TGV regularization employs l1-based form, which is not the most direct method for maximizing sparsity prior. In this study, we propose a total generalized p-variation (TGpV) regularization model to improve the sparsity exploitation of TGV and offer efficient solutions to few-view CT image reconstruction problems. To solve the nonconvex optimization problem of the TGpV minimization model, we then present an efficient iterative algorithm based on the alternating minimization of augmented Lagrangian function. All of the resulting subproblems decoupled by variable splitting admit explicit solutions by applying alternating minimization method and generalized p-shrinkage mapping. In addition, approximate solutions that can be easily performed and quickly calculated through fast Fourier transform are derived using the proximal point method to reduce the cost of inner subproblems. The accuracy and efficiency of the simulated and real data are qualitatively and quantitatively evaluated to validate the efficiency and feasibility of the proposed method. Overall, the proposed method exhibits reasonable performance and outperforms the original TGV-based method when applied to few-view problems.

Computed Tomography Sinogram Inpainting With Compound Prior Modelling Both Sinogram and Image Sparsity

“…Other choices, such as total generalized variation [42,43], curvelet transform, or shearlet transform [44], may provide better performance. Moreover, the statistical model or other prior-based model could also be integrated into our proposed model.…”

Image Deblurring under Impulse Noise via Total Generalized Variation and Non-Convex Shrinkage

“…In Equation (16), F −1 2D denotes a two-dimensional inverse Fourier transform, and the division symbol / denotes element-wise matrix division. According to Cramer's rule, || * in Equation (16) refers to the following:…”

“…Zheng et al [15] theoretically compared the advantages of Lp-minimization and L1-minimization, and tested them in a Gaussian noisy setting. Zhang et al [16] applied Lp shrinkage [17] to the reconstruction of tomographic images to obtain an effect superior to the L1-norm constraint. Chen et al [18] applied an Lp-pseudo-norm constraint to the sparse time-frequency representation of the seismic signal spectrum to obtain a time-frequency spectrum with higher time-frequency resolution.…”

Target Tracking Algorithm Based on an Adaptive Feature and Particle Filter

“…In recent years, the L p norm was widely used. Zhang et al [19] proposed an L p norm relaxation to improve the sparsity exploitation of the total generalized variation. Xie et al [20] generalized the nuclear norm minimization to the Schatten L p norm minimization, which achievedd good results both in background subtraction and image denoising.…”