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>

IVDST: A Fast Algorithm for Atomic Norm Minimization in Line Spectral Estimation

A Coordinate Descent Approach to Atomic Norm Denoising

Fast Two-Dimensional Atomic Norm Minimization in Spectrum Estimation and Denoising.

Fast X-ray CT image reconstruction using the linearized augmented Lagrangian method with ordered subsets

Accelerated singular value thresholding for matrix completion.

When atomic norm meets the G-filter: A general framework for line spectral estimation

New Reweighted Atomic Norm Minimization Approach for Line Spectral Estimation

Super-Resolution Compressed Sensing for Line Spectral Estimation: an Iterative Reweighted Approach

A Novel Gradient Descent Least-Squares (GDLSs) Algorithm for Efficient Gridless Line Spectrum Estimation With Applications in Tomographic SAR Imaging

Enhancing Sparsity and Resolution Via Reweighted Atomic Norm Minimization

The variant of the iterative shrinkage-thresholding algorithm for minimization of the l 1 over l ∞ norms

Separation-Free Super-Resolution from Compressed Measurements is Possible: an Orthonormal Atomic Norm Minimization Approach

Line Spectral Analysis Using the G-Filter: An Atomic Norm Minimization Approach

A Fast Algorithm for Sparse Reconstruction Based on Shrinkage, Subspace Optimization, and Continuation

Sparse Clutter Estimation for STAP Based on Decouple Atomic Norm Minimization

Gridless sparse ISAR imaging using modified atomic norm minimization

Few Views Image Reconstruction Using Alternating Direction Method Via ‐norm Minimization

Fast Linearized Bregman Iteration for Compressive Sensing and Sparse Denoising

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

A Primal Douglas–Rachford Splitting Method for the Constrained Minimization Problem in Compressive Sensing

Robust Depth Linear Error Decomposition with Double Total Variation and Nuclear Norm for Dynamic MRI Reconstruction