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.

Feasibility study of block-matching and 3D filtering denoising algorithm in multi-material decomposition technique for dual-energy computed tomography

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

Noise Suppression in Image-Domain Multi-Material Decomposition for Dual-Energy CT

Image Domain Multi-Material Decomposition Using Single Energy CT.

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

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

Fast and effective single‐scan dual‐energy cone‐beam CT reconstruction and decomposition denoising based on dual‐energy vectorization

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

Quantitative material decomposition using spectral computed tomography with an energy-resolved photon-counting detector

Feasibility study of three-material decomposition in dual-energy cone-beam CT imaging with deep learning

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

Iterative dual energy material decomposition from spatial mismatched raw data sets

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

Material decomposition using dual-energy CT with unsupervised learning

Image-Domain Multimaterial Decomposition for Dual-Energy CT Based on Dictionary Learning and Relative Total Variation

Material decomposition for simulated dual-energy breast computed tomography via hybrid optimization method

O ct 2 01 7 Image-domain multi-material decomposition for dual-energy CT based on correlation and sparsity of material images

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

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

Dual-Energy Ct Reconstruction Using Deep Mutual-Domain Knowledge For Basis Decomposition And Denoising

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