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>

Multi-Energy Ct Reconstruction Based on Low Rank and Sparsity with the Split-Bregman Method (Mlrss)

Sparsity-constrained SENSE reconstruction: An efficient implementation using a fast composite splitting algorithm

Patch-wise non-local low-rank for few-view multi-energy CT reconstruction

Spectral CT reconstruction via low-rank representation and structure preserving regularization

A tensor PRISM algorithm for multi-energy CT reconstruction and comparative studies.

Learned Multi-layer Residual Sparsifying Transform Model for Low-dose CT Reconstruction

Limited-angle Multi-Energy CT Using Joint Clustering Prior and Sparsity Regularization

Refined Locally Linear Transform-Based Spectral-Domain Gradient Sparsity and Its Applications in Spectral CT Reconstruction

MVMS-RCN: A Dual-Domain Unfolding CT Reconstruction with Multi-sparse-view and Multi-scale Refinement-correction

Low-Rank Plus Sparse Reconstruction Using Dictionary Learning for 3D-Mri

Sparse/Low Rank Constrained Reconstruction for Dynamic PET Imaging.

A prior image constraint robust principal component analysis reconstruction method for sparse segmental multi-energy computed tomography

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

Low-Rank Tensor Subspace Decomposition with Weighted Group Sparsity for the Acceleration of Non-Cartesian Dynamic MRI.

Nonconvex ${L_ {1/2}} $-Regularized Nonlocal Self-similarity Denoiser for Compressive Sensing based CT Reconstruction

Weighted Tensor Low-Rankness and Learnable Analysis Sparse Representation Model for Texture Preserving Low-Dose CT Reconstruction

Hybrid plug-and-play CT image restoration using nonconvex low-rank group sparsity and deep denoiser priors

Sparse-View CT Reconstruction Based on a Hybrid Domain Model with Multi-Level Wavelet Transform

Photon-counting spectral CT reconstruction with sparse and double low-rank components fusion

A dual-domain deep learning-based reconstruction method for fully 3D sparse data helical CT

Dual-Domain Reconstruction Network Incorporating Multi-Level Wavelet Transform and Recurrent Convolution for Sparse View Computed Tomography Imaging