Abstract:Background: Dual-energy computed tomography (DECT) is a promising technique, which can provide unique capability for material quantification. The iterative reconstruction of material maps requires spectral information and its accuracy is affected by spectral mismatch. Simultaneously estimating the spectra and reconstructing material maps avoids extra workload on spectrum estimation and the negative impact of spectral mismatch. However, existing methods are not satisfactory in image detail preservation, edge retention, and convergence rate. The purpose of this paper was to mine the similarity between the reconstructed images and the material images to improve the imaging quality, and to design an effective iteration strategy to improve the convergence efficiency. Methods: The material-image subspace decomposition-based iterative reconstruction (MISD-IR) with spectrum estimation was proposed for DECT. MISD-IR is an optimized model combining spectral estimation and material reconstruction with fast convergence speed and promising noise suppression capability. We proposed to reconstruct the material maps based on extended simultaneous algebraic reconstruction techniques and estimation of the spectrum with model spectral. To stabilize the iteration and alleviate the influence of errors, we introduced a weighted proximal operator based on the block coordinate descending algorithm (WP-BCD). Furthermore, the reconstructed computed tomography (CT) images were introduced to suppress the noise based on subspace decomposition, which relies on non-local regularization to prevent noise accumulation. Results: In numerical experiments, the results of MISD-IR were closer to the ground truth compared with other methods. In real scanning data experiments, the results of MISD-IR showed sharper edges and details. Compared with other one-step iterative methods in the experiment, the running time of MISD-IR was reduced by 75%. Conclusions: The proposed MISD-IR can achieve accurate material decomposition (MD) without known energy spectrum in advance, and has good retention of image edges and details. Compared with other one-step iterative methods, it has high convergence efficiency.

Material Discrimination By Multi-Spectral Ct Based On Image Total Variation And Tensor Dictionary

Contrast-Medium Anisotropy-Aware Tensor Total Variation Model for Robust Cerebral Perfusion CT Reconstruction with Low-Dose Scans

Dictionary learning based image-domain material decomposition for spectral CT

DLIMD: Dictionary Learning based Image-domain Material Decomposition for spectral CT

Dynamic-dual-energy spectral CT for improving multi-material decomposition in image-domain.

Reconstruction and Multi-Material Decomposition of Spectral Ct in Dynamic-Threshold-Based Counting and Integrating Modes

Image Domain Multi-Material Decomposition Using Single Energy CT.

A Post-Processed Multi-Material Decomposition Method Based on Dynamic Dual-energy Detectors

Sparse-View Spectral CT Reconstruction Based on Tensor Decomposition and Total Generalized Variation

Spectral CT Reconstruction with Weighted Non-Local Total-Variation Minimization

Multi-Material Decomposition for Single Energy CT Using Material Sparsity Constraint

A TV-constrained Decomposition Method for Spectral CT

Robust Multimaterial Decomposition of Spectral CT Using Convolutional Neural Networks

Optimized Spatial-Spectral CT for Multi-Material Decomposition

MISD-IR: material-image subspace decomposition-based iterative reconstruction with spectrum estimation for dual-energy computed tomography

Reweighted Anisotropic Total Variation Minimization for Limited-Angle CT Reconstruction

Dual-energy CT Reconstruction Based on Dictionary Learning and Total Variation Constraint

Image-domain multi-material decomposition for dual-energy CT based on correlation and sparsity of material images

Image-domain Material Decomposition Based on Spectra Recovery in Spectral CT with Photon-counting Detectors

Application of Deep Leaming in Multi-Material Decomposition of Spectral CT

Low-dose spectral CT reconstruction using L0 image gradient and tensor dictionary