Abstract:PURPOSE:Dual-Energy Computed Tomography (DECT) is of great interest in medical imaging, security inspection, and nondestructive testing. Most DECT reconstruction methods focus on producing two material images with different linear attenuation coefficients. However, the ability to reconstruct three or more basis materials is clinically and industrially important. Under the assumption that there are at most three materials in each pixel, there are a few methods that estimate multiple material images from DECT measurements by enforcing sum-to-one and a box constraint ([0 1]) derived from both the volume and mass conservation assumption. The recently proposed image-domain multimaterial decomposition (MMD) method introduces edge-preserving regularization for each material image. It enforces the assumption that there are at most three materials in each pixel using a time-consuming loop over all possible material triplets. However, this method neglects relations among material images. We propose a new image-domain MMD model for DECT that considers the prior information that different material images have common or complementary edges and encourages sparsity of material composition in each pixel using regularization.METHOD:The proposed PWLS-TNV-ℓ0 method uses penalized weighted least-square (PWLS) reconstruction with three regularization terms. The first term is total nuclear variation (TNV) that accounts for the image property that basis material images share common or complementary boundaries and each material image is piecewise constant. The second term is an ℓ0 norm that encourages each pixel containing a small subset of material types out of several possible materials. The third term is a characteristic function based on sum-to-one and a box constraint derived from the volume and mass conservation assumption. We apply the Alternating Direction Method of Multipliers (ADMM) to optimize the cost function of the PWLS-TNV-ℓ0 method.RESULT:We evaluated the proposed method on a simulated digital phantom, Catphan©600 phantom and patient's pelvis data. We implemented two existing image-domain MMD methods for DECT, the Direct Inversion and the PWLS-EP-LOOP method. We initialized the PWLS-TNV-ℓ0 method and the PWLS-EP-LOOP method with the results of the Direct Inversion method and compared performance of the proposed method with that of the PWLS-EP-LOOP method. The proposed method lowers the bias of decomposed material fractions by 84.47% in the digital phantom study, by 99.50% in the Catphan©600 phantom study, and by 99.64% in the pelvis patient study, respectively, compared to the PWLS-EP-LOOP method. The proposed method reduces noise standard deviation (STD) by 52.21% in the Catphan©600 phantom study, and by 16.74% in the patient's pelvis study, compared to the PWLS-EP-LOOP method. The proposed method increases volume fraction accuracy by 6.04%,20.55%, and 13.46% for the digital phantom, the Catphan©600 phantom, and the patient's pelvis study, respectively, compared to the PWLS-EP-LOOP method. Compared with the PWLS-EP-LOOP method, the root mean square percentage error [RMSE(%)] of electron densities in the Catphan©600 phantom is decreased by about 7.39%.CONCLUSIONS:We proposed an image-domain MMD method, PWLS-TNV-ℓ0 , for DECT. The PWLS-TNV-ℓ0 method takes low rank property of material image gradients, sparsity of material composition and mass and volume conservation into consideration. The proposed method suppresses noise, reduces cross contamination, and improves accuracy in the decomposed material images, compared to the PWLS-EP-LOOP method.

High-Resolution Model-based Material Decomposition for Multi-Layer Flat-Panel Detectors.

High-resolution model-based material decomposition in dual-layer flat-panel CBCT

Model-based Material Decomposition with System Blur Modeling

Material decomposition from photon-counting CT using a convolutional neural network and energy-integrating CT training labels

Model-based material decomposition with a penalized nonlinear least-squares CT reconstruction algorithm

Accurate Multi-Material Decomposition in Dual-Energy CT: A Phantom Study

Approaches for Three Material Decomposition Using a Triple-Layer Flat-Panel Detector.

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

Dual-energy head cone-beam CT using a dual-layer flat-panel detector: physics-based material decomposition

Image Domain Multi-Material Decomposition Using Single Energy CT.

MMD‐Net: Image domain multi‐material decomposition network for dual‐energy CT imaging

Image‐domain Multimaterial Decomposition for Dual‐energy CT Based on Prior Information of Material Images

Material decomposition with a prototype photon-counting detector CT system: expanding a stoichiometric dual-energy CT method via energy bin optimization and K-edge imaging

Iterative Clustering Material Decomposition Aided by Empirical Spectral Correction for High-Resolution Photon-Counting Detectors in Micro-CT

Deep learning-based material decomposition of iodine and calcium in mobile photon counting detector CT

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

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

Optimized Spatial-Spectral CT for Multi-Material Decomposition

Dynamic Material Decomposition Method for MeV Dual-Energy X-ray CT

Robust Multimaterial Decomposition of Spectral CT Using Convolutional Neural Networks

Material decomposition using dual-energy CT with unsupervised learning