A General-Thresholding Solution for [Formula Omitted] Regularized CT Reconstruction

It is well known that [Formula Omitted] minimization can be used to recover sufficiently sparse unknown signals in the compressive sensing field. The [Formula Omitted] regularization method, a generalized version between the well-known [Formula Omitted] regularization and the [Formula Omitted] regularization, has been proposed for a sparser solution. In this paper, we derive several quasi-analytic thresholding representations for the [Formula Omitted] regularization. The derived representations are exact matches for the well-known soft-threshold filtering for the [Formula Omitted] regularization and the hard-threshold filtering for the [Formula Omitted] regularization. The error bounds of the approximate general formulas are analyzed. The general-threshold representation formulas are incorporated into an iterative thresholding framework for a fast solution of an [Formula Omitted] regularized computed tomography (CT) reconstruction. A series of simulated and realistic data experiments are conducted to evaluate the performance of the proposed general-threshold filtering algorithm for CT reconstruction, and it is also compared with the well-known reweighted approach. Compared with the reweighted algorithm, the proposed general-threshold filtering algorithm can substantially reduce the necessary view number for an accurate reconstruction of the Shepp-Logan phantom. In addition, the proposed general-threshold filtering algorithm performs well in terms of image quality, reconstruction accuracy, convergence speed, and sensitivity to parameters.