Abstract:<section class="article-section__content"><h3 class="article-section__sub-title section1"> Purpose</h3><p>Spectral computed tomography (CT) is proposed by extending the conventional CT along the energy dimension. One newly implementation is to employ an energy‐discriminating photon counting detector (PCD), which can distinguish photon energy and divide a whole X‐ray spectrum into several energy bins with appropriate post‐processing steps. The state‐of‐the‐art PCD‐based spectral CT has superior energy‐resolution and material distinguishability, and it further has a great potential in both medical and industrial applications. To improve the reconstruction quality and decomposition accuracy, in this work, we propose an optimization‐based spectral CT reconstruction method with an innovational sparsity constraint.</p></section><section class="article-section__content"><h3 class="article-section__sub-title section1"> Methods</h3><p>We first employ a locally linear transform to the reconstructed channel images, and the structural similarity along the spectral dimension is effectively converted to a one‐dimensional (1D) gradient sparsity. Then, combining the prior knowledge of piecewise constant in the spatial domain (e.g. a two‐dimensional (2D) gradient sparsity feature), we unify both spectral and spatial dimensions and establish a joint three‐dimensional (3D) gradient sparsity. In addition, we use the <i>L</i><sub>0</sub>‐norm to measure the proposed sparsity and incorporate it as a smoothness constraint to concretize a general optimization framework. Furthermore, we develop the corresponding iterative algorithm to solve the optimization problem. </p></section><section class="article-section__content"><h3 class="article-section__sub-title section1"> Results</h3><p>Both visual results and quantitative indexes of numerical simulations and phantom experiments demonstrate the proposed method outperform the conventional filtered backprojection (FBP), total variation (TV), 2D <i>L</i><sub>0</sub>‐norm (<i>L </i><sub>0</sub>), and TV with low rank (TVLR) based methods. From the image and ROI comparisons, we find the proposed method performs well in noise suppression, detail maintenance, and decomposition accuracy. However, the FBP suffers severe noise, the TV and <i>L </i><sub>0</sub> are difficult to work consistently among different energy bins, and the TVLR fails to avoid gray value shift. The image quality assessments, such as peak signal‐to‐noise ratio (PSNR), normal mean absolute deviation (NMAD) and structural similarity (SSIM), also consistently indicate the proposed method can effectively removing noise and keeping fine structures in both channel‐wise reconstructions and material decompositions. </p></section><section class="article-section__content"><h3 class="article-section__sub-title section1"> Conclusions</h3><p>By employing a locally linear transform, the structural similarity among spectral channel images is converted to a 1D gradient sparsity, and the gray value shift is effectively avoided when the difference measurement is minimized. The 3D <i>L</i><sub>0</sub>‐norm jointly and uniformly measures the gradient sparsity in both spectral and spatial dimensions. The cooperation of locally linear transform and 3D <i>L</i><sub>0</sub>‐norm well reinforces the global sparse features and keeps the correlation along spectral dimension without bringing gray‐value distortions. The corresponding constraint optimization model is fast and stably solved by using an alternative direction technique. Both numerical simulations and phantom experiments confirm the superior performance of the proposed method in noise suppression, structure maintenance, and accurate decomposition. </p></section>

Accelerated Gradient Methods for Total-Variation-based CT Image Reconstruction

Reweighted Anisotropic Total Variation Minimization for Limited-Angle CT Reconstruction

A nonlinear total variation based computed tomography (CT) image reconstruction method using gradient reinforcement

Reconstruction of Sparse-View X-ray Computed Tomography Based on Adaptive Total Variation Minimization

Constrained Tpv Minimization for Enhanced Exploitation of Gradient Sparsity: Application to Ct Image Reconstruction

A low-complexity 2-point step size gradient projection method with selective function evaluations for smoothed total variation based CBCT reconstructions

Sparse-view cone-beam computed tomography iterative reconstruction based on new multi-gradient direction total variation

GPU-based Low Dose CT Reconstruction via Edge-preserving Total Variation Regularization

Template-based CT Reconstruction with Optimal Transport and Total Generalized Variation

A fast iterative reconstruction method based on the selective total variation for sparse angular CT

A Fast Regularized Iterative Algorithm for Fan-Beam CT Reconstruction.

A very fast iterative algorithm for TV-regularized image reconstruction with applications to low-dose and few-view CT

Constrained <formula formulatype="inline"><tex Notation="TeX">${\rm T}p{\rm V}$</tex> </formula> Minimization for Enhanced Exploitation of Gradient Sparsity: Application to CT Image Reconstruction

Low-dose CT reconstruction via edge-preserving total variation regularization

3D Anisotropic Total Variation Method for Limited-Angle Ct Reconstruction

A Reweighted Total Variation Minimization Method for Few View CT Reconstruction in the Instant CT

Sparse-view x-ray CT reconstruction via total generalized variation regularization

Locally linear transform based three‐dimensional gradient ‐norm minimization for spectral CT reconstruction

Image Reconstruction in Cone Beam Computed Tomography Using Controlled Gradient Sparsity

A Generalized Accelerated Proximal Gradient Approach for Total-Variation-Based Image Restoration

Superiorized algorithm for reconstruction of CT images from sparse-view and limited-angle polyenergetic data