Dual-energy CT imaging with limited-angular-range data

Image-domain material decomposition for dual-energy CT using unsupervised learning with data-fidelity loss

BACKGROUND

Dual-energy computed tomography (DECT) and material decomposition play vital roles in quantitative medical imaging. However, the decomposition process may suffer from significant noise amplification, leading to severely degraded image signal-to-noise ratios (SNRs). While existing iterative algorithms perform noise suppression using different image priors, these heuristic image priors cannot accurately represent the features of the target image manifold. Although deep learning-based decomposition methods have been reported, these methods are in the supervised-learning framework requiring paired data for training, which is not readily available in clinical settings.

PURPOSE

This work aims to develop an unsupervised-learning framework with data-measurement consistency for image-domain material decomposition in DECT.

METHODS

The proposed framework combines iterative decomposition and deep learning-based image prior in a generative adversarial network (GAN) architecture. In the generator module, a data-fidelity loss is introduced to enforce the measurement consistency in material decomposition. In the discriminator module, the discriminator is trained to differentiate the low-noise material-specific images from the high-noise images. In this scheme, paired images of DECT and ground-truth material-specific images are not required for the model training. Once trained, the generator can perform image-domain material decomposition with noise suppression in a single step.

RESULTS

In the simulation studies of head and lung digital phantoms, the proposed method reduced the standard deviation (SD) in decomposed images by 97% and 91% from the values in direct inversion results. It also generated decomposed images with structural similarity index measures (SSIMs) greater than 0.95 against the ground truth. In the clinical head and lung patient studies, the proposed method suppressed the SD by 95% and 93% compared to the decomposed images of matrix inversion.

CONCLUSIONS

Since the invention of DECT, noise amplification during material decomposition has been one of the biggest challenges, impeding its quantitative use in clinical practice. The proposed method performs accurate material decomposition with efficient noise suppression. Furthermore, the proposed method is within an unsupervised-learning framework, which does not require paired data for model training and resolves the issue of lack of ground-truth data in clinical scenarios.

Accurate Reconstruction of Multiple Basis Images Directly From Dual Energy CT Data

OBJECTIVE

We develop optimization-based algorithms to accurately reconstruct multiple ( 2) basis images directly from dual-energy (DE) data in CT.

METHODS

In medical and industrial CT imaging, some basis materials such as bone, metals, and contrast agents of interest are confined often spatially within regions in the image. Exploiting this observation, we develop an optimization-based algorithm to reconstruct, directly from DE data, basis-region images from which multiple ( 2) basis images and virtual monochromatic images (VMIs) can be obtained over the entire image array.

RESULTS

We conduct experimental studies using simulated and real DE data in CT, and evaluate basis images and VMIs obtained in terms of visual inspection and quantitative metrics. The study results reveal that the algorithm developed can accurately and robustly reconstruct multiple ( 2) basis images directly from DE data.

CONCLUSIONS

The developed algorithm can yield accurate multiple ( 2) basis images, VMIs, and physical quantities of interest from DE data in CT.

SIGNIFICANCE

The work may provide insights into the development of practical procedures for reconstructing multiple basis images, VMIs, and physical quantities from DE data in applications. The work can be extended to reconstruct multiple basis images in multi-spectral or photon-counting CT.

Directional TV algorithm for image reconstruction from sparse-view projections in EPR imaging

Objective.Electron paramagnetic resonance (EPR) imaging is an advanced in vivo oxygen imaging modality. The main drawback of EPR imaging is the long scanning time. Sparse-view projections collection is an effective fast scanning pattern. However, the commonly-used filtered back projection (FBP) algorithm is not competent to accurately reconstruct images from sparse-view projections because of the severe streak artifacts. The aim of this work is to develop an advanced algorithm for sparse reconstruction of 3D EPR imaging.Methods.The optimization based algorithms including the total variation (TV) algorithm have proven to be effective in sparse reconstruction in EPR imaging. To further improve the reconstruction accuracy, we propose the directional TV (DTV) model and derive its Chambolle-Pock solving algorithm.Results.After the algorithm correctness validation on simulation data, we explore the sparse reconstruction capability of the DTV algorithm via a simulated six-sphere phantom and two real bottle phantoms filled with OX063 trityl solution and scanned by an EPR imager with a magnetic field strength of 250 G.Conclusion.Both the simulated and real data experiments show that the DTV algorithm is superior to the existing FBP and TV-type algorithms and a deep learning based method according to visual inspection and quantitative evaluations in sparse reconstruction of EPR imaging.Significance.These insights gained in this work may be used in the development of fast EPR imaging workflow of practical significance.

Directional TV algorithm for fast EPR imaging

Precise radiation guided by oxygen images has demonstrated superiority over the traditional radiation methods. Electron paramagnetic resonance (EPR) imaging has proven to be the most advanced oxygen imaging modality. However, the main drawback of EPR imaging is the long scan time. For each projection, we usually need to collect the projection many times and then average them to achieve high signal-to-noise ratio (SNR). One approach to fast scan is to reduce the repeating time for each projection. While the projections would be noisy and thus the traditional commonly-use filtered backprojection (FBP) algorithm would not be capable of accurately reconstructing images. Optimization-based iterative algorithms may accurately reconstruct images from noisy projections for they may incorporate prior information into optimization models. Based on the total variation (TV) algorithms for EPR imaging, in this work, we propose a directional TV (DTV) algorithm to further improve the reconstruction accuracy. We construct the DTV constrained, data divergence minimization (DTVcDM) model, derive its Chambolle-Pock (CP) solving algorithm, validate the correctness of the whole algorithm, and perform evaluations via simulated and real data. The experimental results show that the DTV algorithm outperforms the existing TV and FBP algorithms in fast EPR imaging. Compared to the standard FBP algorithm, the proposed algorithm may achieve 10 times of acceleration.

A model-based direct inversion network (MDIN) for dual spectral computed tomography

Objective. Dual spectral computed tomography (DSCT) is a very challenging problem in the field of imaging. Due to the nonlinearity of its mathematical model, the images reconstructed by the conventional CT usually suffer from the beam hardening artifacts. Additionally, several existing DSCT methods rely heavily on the information of the spectra, which is often not readily available in applications. To address this problem, in this study, we aim to develop a novel approach to improve the DSCT reconstruction performance.Approach. A model-based direct inversion network (MDIN) is proposed for DSCT, which can directly predict the basis material images from the collected polychromatic projections. The all operations are performed in the network, requiring neither the conventional algorithms nor the information of the spectra. It can be viewed as an approximation to the inverse procedure of DSCT imaging model. The MDIN is composed of projection pre-decomposition module (PD-module), domain transformation layer (DT-layer), and image post-decomposition module (ID-module). The PD-module first performs the pre-decomposition on the polychromatic projections that consists of a series of stacked one-dimensional convolution layers. The DT-layer is designed to obtain the preliminary decomposed results, which has the characteristics of sparsely connected and learnable parameters. And the ID-module uses a deep neural network to further decompose the reconstructed results of the DT-layer so as to achieve higher-quality basis material images.Main results. Numerical experiments demonstrate that the proposed MDIN has significant advantages in substance decomposition, artifact reduction and noise suppression compared to other methods in the DSCT reconstruction.Significance. The proposed method has a flexible applicability, which can be extended to other CT problems, such as multi-spectral CT and low dose CT.

Prototyping optimization-based image reconstructions from limited-angular-range data in dual-energy CT

Image reconstruction from data collected over full-angular range (FAR) in dual-energy CT (DECT) is well-studied. There exists interest in DECT with advanced scan configurations in which data are collected only over limited-angular ranges (LARs) for meeting unique workflow needs in certain practical imaging applications, and thus in the algorithm development for image reconstruction from such LAR data. The objective of the work is to investigate and prototype image reconstructions in DECT with LAR scans. We investigate and prototype optimization programs with various designs of constraints on the directional-total-variations (DTVs) of virtual monochromatic images and/or basis images, and derive the DTV algorithms to numerically solve the optimization programs for achieving accurate image reconstruction from data collected in a slew of different LAR scans. Using simulated and real data acquired with low- and high-kV spectra over LARs, we conduct quantitative studies to demonstrate and evaluate the optimization programs and their DTV algorithms developed. As the results of the numerical studies reveal, while the DTV algorithms yield images of visual quality and quantitative accuracy comparable to that of the existing algorithms from FAR data, the former reconstruct images with improved visualization, reduced artifacts, and also enhanced quantitative accuracy when applied to LAR data in DECT. Optimization-based, one-step algorithms, including the DTV algorithms demonstrated, can be developed for quantitative image reconstruction from spectral data collected over LARs of extents that are considerably smaller than the FAR in DECT. The theoretical and numerical results obtained can be exploited for prototyping designs of optimization-based reconstructions and LAR scans in DECT, and they may also yield insights into the development of reconstruction procedures in practical DECT applications. The approach and algorithms developed can naturally be applied to investigating image reconstruction from LAR data in multi-spectral and photon-counting CT.

A dense and U-shaped transformer with dual-domain multi-loss function for sparse-view CT reconstruction

OBJECTIVE

CT image reconstruction from sparse-view projections is an important imaging configuration for low-dose CT, as it can reduce radiation dose. However, the CT images reconstructed from sparse-view projections by traditional analytic algorithms suffer from severe sparse artifacts. Therefore, it is of great value to develop advanced methods to suppress these artifacts. In this work, we aim to use a deep learning (DL)-based method to suppress sparse artifacts.

METHODS

Inspired by the good performance of DenseNet and Transformer architecture in computer vision tasks, we propose a Dense U-shaped Transformer (D-U-Transformer) to suppress sparse artifacts. This architecture exploits the advantages of densely connected convolutions in capturing local context and Transformer in modelling long-range dependencies, and applies channel attention to fusion features. Moreover, we design a dual-domain multi-loss function with learned weights for the optimization of the model to further improve image quality.

RESULTS

Experimental results of our proposed D-U-Transformer yield performance improvements on the well-known Mayo Clinic LDCT dataset over several representative DL-based models in terms of artifact suppression and image feature preservation. Extensive internal ablation experiments demonstrate the effectiveness of the components in the proposed model for sparse-view computed tomography (SVCT) reconstruction.

SIGNIFICANCE

The proposed method can effectively suppress sparse artifacts and achieve high-precision SVCT reconstruction, thus promoting clinical CT scanning towards low-dose radiation and high-quality imaging. The findings of this work can be applied to denoising and artifact removal tasks in CT and other medical images.

A dual-energy CT reconstruction method based on anchor network from dual quarter scans

Compared with conventional single-energy computed tomography (CT), dual-energy CT (DECT) provides better material differentiation but most DECT imaging systems require dual full-angle projection data at different X-ray spectra. Relaxing the requirement of data acquisition is an attractive research to promote the applications of DECT in wide range areas and reduce the radiation dose as low as reasonably achievable. In this work, we design a novel DECT imaging scheme with dual quarter scans and propose an efficient method to reconstruct the desired DECT images from the dual limited-angle projection data. We first study the characteristics of limited-angle artifacts under dual quarter scans scheme, and find that the negative and positive artifacts of DECT images are complementarily distributed in image domain because the corresponding X-rays of high- and low-energy scans are symmetric. Inspired by this finding, a fusion CT image is generated by integrating the limited-angle DECT images of dual quarter scans. This strategy enhances the true image information and suppresses the limited-angle artifacts, thereby restoring the image edges and inner structures. Utilizing the capability of neural network in the modeling of nonlinear problem, a novel Anchor network with single-entry double-out architecture is designed in this work to yield the desired DECT images from the generated fusion CT image. Experimental results on the simulated and real data verify the effectiveness of the proposed method. This work enables DECT on imaging configurations with half-scan and largely reduces scanning angles and radiation doses.

Dual-Energy CT in Cardiothoracic Imaging: Current Developments

This article describes the technical principles and clinical applications of dual-energy computed tomography (DECT) in the context of cardiothoracic imaging with a focus on current developments and techniques. Since the introduction of DECT, different vendors developed distinct hard and software approaches for generating multi-energy datasets and multiple DECT applications that were developed and clinically investigated for different fields of interest. Benefits for various clinical settings, such as oncology, trauma and emergency radiology, as well as musculoskeletal and cardiovascular imaging, were recently reported in the literature. State-of-the-art applications, such as virtual monoenergetic imaging (VMI), material decomposition, perfused blood volume imaging, virtual non-contrast imaging (VNC), plaque removal, and virtual non-calcium (VNCa) imaging, can significantly improve cardiothoracic CT image workflows and have a high potential for improvement of diagnostic accuracy and patient safety.

4D-image reconstruction directly from limited-angular-range data in continuous-wave electron paramagnetic resonance imaging

OBJECTIVE

We investigate and develop optimization-based algorithms for accurate reconstruction of four-dimensional (4D)-spectral-spatial (SS) images directly from data collected over limited angular ranges (LARs) in continuous-wave (CW) electron paramagnetic resonance imaging (EPRI).

METHODS

Basing on a discrete-to-discrete data model devised in CW EPRI employing the Zeeman-modulation (ZM) scheme for data acquisition, we first formulate the image reconstruction problem as a convex, constrained optimization program that includes a data fidelity term and also constraints on the individual directional total variations (DTVs) of the 4D-SS image. Subsequently, we develop a primal-dual-based DTV algorithm, simply referred to as the DTV algorithm, to solve the constrained optimization program for achieving image reconstruction from data collected in LAR scans in CW-ZM EPRI.

RESULTS

We evaluate the DTV algorithm in simulated- and real-data studies for a variety of LAR scans of interest in CW-ZM EPRI, and visual and quantitative results of the studies reveal that 4D-SS images can be reconstructed directly from LAR data, which are visually and quantitatively comparable to those obtained from data acquired in the standard, full-angular-range (FAR) scan in CW-ZM EPRI.

CONCLUSION

An optimization-based DTV algorithm is developed for accurately reconstructing 4D-SS images directly from LAR data in CW-ZM EPRI. Future work includes the development and application of the optimization-based DTV algorithm for reconstructions of 4D-SS images from FAR and LAR data acquired in CW EPRI employing schemes other than the ZM scheme.

SIGNIFICANCE

The DTV algorithm developed may be exploited potentially for enabling and optimizing CW EPRI with minimized imaging time and artifacts by acquiring data in LAR scans.

TIME-Net: Transformer-Integrated Multi-Encoder Network for limited-angle artifact removal in dual-energy CBCT

Dual-energy cone-beam computed tomography (DE-CBCT) is a promising imaging technique with foreseeable clinical applications. DE-CBCT images acquired with two different spectra can provide material-specific information. Meanwhile, the anatomical consistency and energy-domain correlation result in significant information redundancy, which could be exploited to improve image quality. In this context, this paper develops the Transformer-Integrated Multi-Encoder Network (TIME-Net) for DE-CBCT to remove the limited-angle artifacts. TIME-Net comprises three encoders (image encoder, prior encoder, and transformer encoder), two decoders (low- and high-energy decoders), and one feature fusion module. Three encoders extract various features for image restoration. The feature fusion module compresses these features into more compact shared features and feeds them to the decoders. Two decoders perform differential learning for DE-CBCT images. By design, TIME-Net could obtain high-quality DE-CBCT images using two complementary quarter-scans, holding great potential to reduce radiation dose and shorten the acquisition time. Qualitative and quantitative analyses based on simulated data and real rat data have demonstrated the promising performance of TIME-Net in artifact removal, subtle structure restoration, and reconstruction accuracy preservation. Two clinical applications, virtual non-contrast (VNC) imaging and iodine quantification, have proved the potential utility of the DE-CBCT images provided by TIME-Net.

Dual spectral limited-angle CT imaging regularized by edge-preserving diffusion and smoothing

Limited-angle CT scan is an effective way for nondestructive inspection of planar objects, and various methods have been proposed accordingly. When the scanned object contains high-absorption material, such as metal, existing methods may fail due to the beam hardening of X-rays. In order to overcome this problem, we adopt a dual spectral limited-angle CT scan and propose a corresponding image reconstruction algorithm, which takes the polychromatic property of the X-ray into consideration, makes basis material images free of beam hardening artifacts and metal artifacts, and then helps depress the limited-angle artifacts. Experimental results on both simulated PCB data and real data demonstrate the effectiveness of the proposed algorithm.

Accurate Image Reconstruction in Dual-Energy CT with Limited-Angular-Range Data Using a Two-Step Method

Dual-energy CT (DECT) with scans over limited-angular ranges (LARs) may allow reductions in scan time and radiation dose and avoidance of possible collision between the moving parts of a scanner and the imaged object. The beam-hardening (BH) and LAR effects are two sources of image artifacts in DECT with LAR data. In this work, we investigate a two-step method to correct for both BH and LAR artifacts in order to yield accurate image reconstruction in DECT with LAR data. From low- and high-kVp LAR data in DECT, we first use a data-domain decomposition (DDD) algorithm to obtain LAR basis data with the non-linear BH effect corrected for. We then develop and tailor a directional-total-variation (DTV) algorithm to reconstruct from the LAR basis data obtained basis images with the LAR effect compensated for. Finally, using the basis images reconstructed, we create virtual monochromatic images (VMIs), and estimate physical quantities such as iodine concentrations and effective atomic numbers within the object imaged. We conduct numerical studies using two digital phantoms of different complexity levels and types of structures. LAR data of low- and high-kVp are generated from the phantoms over both single-arc (SA) and two-orthogonal-arc (TOA) LARs ranging from 14∘ to 180∘. Visual inspection and quantitative assessment of VMIs obtained reveal that the two-step method proposed can yield VMIs in which both BH and LAR artifacts are reduced, and estimation accuracy of physical quantities is improved. In addition, concerning SA and TOA scans with the same total LAR, the latter is shown to yield more accurate images and physical quantity estimations than the former. We investigate a two-step method that combines the DDD and DTV algorithms to correct for both BH and LAR artifacts in image reconstruction, yielding accurate VMIs and estimations of physical quantities, from low- and high-kVp LAR data in DECT. The results and knowledge acquired in the work on accurate image reconstruction in LAR DECT may give rise to further understanding and insights into the practical design of LAR scan configurations and reconstruction procedures for DECT applications.

Image reconstruction from data over two orthogonal arcs of limited-angular ranges

PURPOSE

Computed tomography (CT) scanning over limited-angular ranges (LARs) is of practical interest in possible reduction of imaging dose and time and in design of non-standard scans. This work aims to investigate image reconstruction for two non-overlapping arcs of LARs, and to demonstrate that they may allow more accurate image reconstruction than may a single arc of LAR.

METHODS

We consider a configuration with two non-overlapping arcs of LARs α₁ and α₂ , whose centers are separated by 90°, and refer to it as a two-orthogonal-arc configuration. Data are generated from a chest phantom with two-orthogonal-arc configurations over total angular coverage α_τ = α₁ + α₂ ranging from 18° to 180°, and images are reconstructed subsequently by use of the directional-total-variation (DTV) algorithm. For comparison, we also consider image reconstruction for a single-arc configuration of angular range α_τ . Quantitative metrics such as the normalized root-mean-square-error (nRMSE) are used for evaluation of image reconstruction accuracy.

RESULTS

Visual inspection and quantitative analysis of images reconstructed reveal that a two-orthogonal-arc configuration generally yields more accurate image reconstruction than does its single-arc counterpart. As total angular range α_τ increases, the DTV algorithm yields image reconstruction with enhanced accuracy, as expected. Also, if α_τ remains constant, the two-orthogonal-arc configuration with α₁ = α₂ generally leads to image reconstruction more accurate than those of two-orthogonal-arc configurations with α₁ ≠ α₂ , as the nRMSE of the former can be lower than that of the latter for up to more than one order of magnitude.

CONCLUSIONS

Appropriately designed two-orthogonal-arc configurations may be exploited for improving image-reconstruction accuracy in CT imaging with reduced angular coverage. This study may yield insights into the design of innovative CT scans for lowering scan time and radiation dose, and/or for avoiding scan collision in, e.g., C-arm CT.

Dual-energy CT imaging over non-overlapping, orthogonal arcs of limited-angular ranges

BACKGROUND

Interest exists in dual-energy computed tomography (DECT) imaging with scanning arcs of limited-angular ranges (LARs) for reducing scan time and radiation dose, and for enabling scan configurations of C-arm CT that can avoid possible collision between the rotating X-ray tube/detector and the imaged subject.

OBJECTIVE

In this work, we investigate image reconstruction for a type of configurations of practical DECT interest, referred to as the two-orthogonal-arc configuration, in which low- and high-kVp data are collected over two non-overlapping arcs of equal LAR α, ranging from 30° to 90°, separated by 90°. The configuration can readily be implemented, e.g., on CT with dual sources separated by 90° or with the slow-kVp-switching technique.

METHODS

The directional-total-variation (DTV) algorithm developed previously for image reconstruction in conventional, single-energy CT is tailored to enable image reconstruction in DECT with two-orthogonal-arc configurations.

RESULTS

Performing visual inspection and quantitative analysis of monochromatic images obtained and effective atomic numbers estimated, we observe that the monochromatic images of the DTV algorithm from LAR data are with substantially reduced LAR artifacts, which are observed otherwise in those of existing algorithms, and thus visually correlate reasonably well, in terms of metrics PCC and nMI, with their reference images obtained from full-angular-range data. In addition, effective atomic numbers estimated from LAR data of DECT with two-orthogonal-arc configurations are in reasonable agreement, with relative errors up to ∼  10%, with those estimated from full-angular-range data in DECT.

CONCLUSIONS

The results acquired in the work may yield insights into the design of LAR configurations of practical dual-energy application relevance in diagnostic CT or C-arm CT imaging.