Noise suppression for dual-energy CT via penalized weighted least-square optimization with similarity-based regularization

Reconstruction method suitable for fast CT imaging

Reconstructing computed tomography (CT) images from an extremely limited set of projections is crucial in practical applications. As the available projections significantly decrease, traditional reconstruction and model-based iterative reconstruction methods become constrained. This work aims to seek a reconstruction method applicable to fast CT imaging when available projections are highly sparse. To minimize the time and cost associated with projections acquisition, we propose a deep learning model, X-CTReNet, which parameterizes a nonlinear mapping function from orthogonal projections to CT volumes for 3D reconstruction. The proposed model demonstrates effective capability in inferring CT volumes from two-view projections compared to baseline methods, highlighting the significant potential for drastically reducing projection acquisition in fast CT imaging.

Image Domain Multi-Material Decomposition Noise Suppression Through Basis Transformation and Selective Filtering

Spectral CT can provide material characterization ability to offer more precise material information for diagnosis purposes. However, the material decomposition process generally leads to amplification of noise which significantly limits the utility of the material basis images. To mitigate such problem, an image domain noise suppression method was proposed in this work. The method performs basis transformation of the material basis images based on a singular value decomposition. The noise variances of the original spectral CT images were incorporated in the matrix to be decomposed to ensure that the transformed basis images are statistically uncorrelated. Due to the difference in noise amplitudes in the transformed basis images, a selective filtering method was proposed with the low-noise transformed basis image as guidance. The method was evaluated using both numerical simulation and real clinical dual-energy CT data. Results demonstrated that compared with existing methods, the proposed method performs better in preserving the spatial resolution and the soft tissue contrast while suppressing the image noise. The proposed method is also computationally efficient and can realize real-time noise suppression for clinical spectral CT images.

Projection domain decomposition denoising algorithm based on low rank and similarity-based regularization

BACKGROUND

Projection Domain Decomposition (PDD) is a dual energy reconstruction method which implements the decomposition process before image reconstruction. The advantage of PDD is that it can alleviate beam hardening artifacts and metal artifacts effectively as energy spectra estimation is considered in PDD. However, noise amplification occurs during the decomposition process, which significantly impacts the accuracy of effective atomic number and electron density. Therefore, effective noise reduction techniques are required in PDD.

OBJECTIVE

This study aims to develop a new algorithm capable of minimizing noise while simultaneously preserving edges and fine details.

METHODS

In this study, a denoising algorithm based on low rank and similarity-based regularization (LRSBR) is presented. This algorithm incorporates the low-rank characteristic of tensors into similarity-based regularization (SBR) framework. This method effectively addresses the issue of instability in edge pixels within the SBR algorithm and enhances the structural consistency of dual-energy images.

RESULTS

A series of simulation and practical experiments were conducted on a dual-layer dual-energy CT system. Experiments demonstrate that the proposed method outperforms exiting noise removal methods in Peak Signal-to-noise Ratio (PSNR), Root Mean Square Error (RMSE), and Structural Similarity (SSIM). Meanwhile, there has been a notable enhancement in the visual quality of CT images.

CONCLUSIONS

The proposed algorithm has a significantly improved noise reduction compared to other competing approach in dual-energy CT. Meanwhile, the LRSBR method exhibits outstanding performance in preserving edges and fine structures, making it practical for PDD applications.

Low-Cost Probabilistic 3D Denoising with Applications for Ultra-Low-Radiation Computed Tomography

We propose a pipeline for synthetic generation of personalized Computer Tomography (CT) images, with a radiation exposure evaluation and a lifetime attributable risk (LAR) assessment. We perform a patient-specific performance evaluation for a broad range of denoising algorithms (including the most popular deep learning denoising approaches, wavelets-based methods, methods based on Mumford−Shah denoising, etc.), focusing both on accessing the capability to reduce the patient-specific CT-induced LAR and on computational cost scalability. We introduce a parallel Probabilistic Mumford−Shah denoising model (PMS) and show that it markedly-outperforms the compared common denoising methods in denoising quality and cost scaling. In particular, we show that it allows an approximately 22-fold robust patient-specific LAR reduction for infants and a 10-fold LAR reduction for adults. Using a normal laptop, the proposed algorithm for PMS allows cheap and robust (with a multiscale structural similarity index >90%) denoising of very large 2D videos and 3D images (with over 107 voxels) that are subject to ultra-strong noise (Gaussian and non-Gaussian) for signal-to-noise ratios far below 1.0. The code is provided for open access.

A generalized simultaneous algebraic reconstruction technique (GSART) for dual-energy X-ray computed tomography

BACKGROUND

Dual-energy computed tomography (DECT) is a widely used and actively researched imaging modality that can estimate the physical properties of an object more accurately than single-energy CT (SECT). Recently, iterative reconstruction methods called one-step methods have received attention among various approaches since they can resolve the intermingled limitations of the conventional methods. However, the one-step methods typically have expensive computational costs, and their material decomposition performance is largely affected by the accuracy in the spectral coefficients estimation.

OBJECTIVE

In this study, we aim to develop an efficient one-step algorithm that can effectively decompose into the basis material maps and is less sensitive to the accuracy of the spectral coefficients.

METHODS

By use of a new loss function that employs the non-linear forward model and the weighted squared errors, we propose a one-step reconstruction algorithm named generalized simultaneous algebraic reconstruction technique (GSART). The proposed algorithm was compared with the image-domain material decomposition and other existing one-step reconstruction algorithm.

RESULTS

In both simulation and experimental studies, we demonstrated that the proposed algorithm effectively reduced the beam-hardening artifacts thereby increasing the accuracy in the material decomposition.

CONCLUSIONS

The proposed one-step reconstruction for material decomposition in dual-energy CT outperformed the image-domain approach and the existing one-step algorithm. We believe that the proposed method is a practically very useful addition to the material-selective image reconstruction field.

Generalized deep iterative reconstruction for sparse-view CT imaging

Sparse-view CT is a promising approach in reducing the X-ray radiation dose in clinical CT imaging. However, the CT images reconstructed from the conventional filtered backprojection (FBP) algorithm suffer from severe streaking artifacts. Iterative reconstruction (IR) algorithms have been widely adopted to mitigate these streaking artifacts, but they may prolong the CT imaging time due to the intense data-specific computations. Recently, model-driven deep learning (DL) CT image reconstruction method, which unrolls the iterative optimization procedures into the deep neural network, has shown exciting prospect in improving the image quality and shortening the reconstruction time. In this work, we explore the generalized unrolling scheme for such iterative model to further enhance its performance on sparse-view CT imaging. By using it, the iteration parameters, regularizer term, data-fidelity term and even the mathematical operations are all assumed to be learned and optimized via the network training. Results from the numerical and experimental sparse-view CT imaging demonstrate that the newly proposed network with the maximum generalization provides the best reconstruction performance.

An iterative reconstruction method based on monochromatic images for dual energy CT

PURPOSE

Dual-energy computed tomography (DECT) scans objects using two different X-ray spectra to acquire more information, which is also called dual spectral CT (DSCT) in some articles. Compared to traditional CT, DECT exhibits superior material distinguishability. Therefore, DECT can be widely used in the medical and industrial domains. However, owing to the nonlinearity and ill condition of DECT, studies are underway on DECT reconstruction to obtain high quality images and achieve fast convergence speed. Therefore, in this study, we propose an iterative reconstruction method based on monochromatic images (IRM-MI) to rapidly obtain high-quality images in DECT reconstruction.

METHODS

An IRM-MI is proposed for DECT. The proposed method converts DECT reconstruction problem from the basis material images decomposition to monochromatic images decomposition to significantly improve the convergence speed of DECT reconstruction by changing the coefficient matrix of the original equations to increase the angle of the high- and low-energy projection curves or reduce the condition number of the coefficient matrix. The monochromatic images were then decomposed into basis material images. Furthermore, we conducted numerical experiments to evaluate the performance of the proposed method.

RESULTS

The decomposition results of the simulated data and real data experiments confirmed the effectiveness of the proposed method. Compared to the extended algebraic reconstruction technique (E-ART) method, the proposed method exhibited a significant increase in the convergence speed by increasing the angle of polychromatic projection curves or decreasing the condition number of the coefficient matrix, when choosing the appropriate monochromatic images. Therefore, the proposed method is also advantageous in acquiring high quality and rapidly converged images.

CONCLUSIONS

We developed an iterative reconstruction method based on monochromatic images for the material decomposition for DECT. The numerical experiments using the proposed method validated its capability of decomposing the basis material images. Furthermore, the proposed method achieved faster convergence speed compared to the E-ART method.

Iterative material decomposition for spectral CT using self-supervised Noise2Noise prior

Compared to conventional computed tomography (CT), spectral CT can provide the capability of material decomposition, which can be used in many clinical diagnosis applications. However, the decomposed images can be very noisy due to the dose limit in CT scanning and the noise magnification of the material decomposition process. To alleviate this situation, we proposed an iterative one-step inversion material decomposition algorithm with a Noise2Noise prior. The algorithm estimated material images directly from projection data and used a Noise2Noise prior for denoising. In contrast to supervised deep learning methods, the designed Noise2Noise prior was built based on self-supervised learning and did not need external data for training. In our method, the data consistency term and the Noise2Noise network were alternatively optimized in the iterative framework, respectively, using a separable quadratic surrogate (SQS) and the Adam algorithm. The proposed iterative algorithm was validated and compared to other methods on simulated spectral CT data, preclinical photon-counting CT data and clinical dual-energy CT data. Quantitative analysis showed that our proposed method performs promisingly on noise suppression and structure detail recovery.

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

PURPOSE

Flat-panel detector (FPD) based dual-energy cone-beam computed tomography (DE-CBCT) is a promising imaging technique for dedicated clinical applications. In this paper, we proposed a fully analytical method for fast and effective single-scan DE-CBCT image reconstruction and decomposition.

METHODS

A rotatable Mo filter was inserted between an x-ray source and imaged object to alternately produce low and high-energy x-ray spectra. First, filtered-backprojection (FBP) method was applied on down-sampled projections to reconstruct low and high-energy images. Then, the two images were converted into a vectorized form represented with an amplitude and an argument image. Using amplitude image as a guide, a joint bilateral filter was applied to denoise the argument image. Then, high-quality dual-energy images were recovered from the amplitude image and the denoised argument image. Finally, the recovered dual-energy images were further used for low-noise material decomposition and electron density synthesis. Imaging was conducted on a Catphan^® 600 phantom and an anthropomorphic head phantom. The proposed method was evaluated via comparison with the traditional two-scan method and a commonly used filtering method (HYPR-LR).

RESULTS

On the Catphan^® 600 phantom, the proposed method successfully reduced streaking artifacts and preserved spatial resolution and noise-power-spectrum (NPS) pattern. In the electron density image, the proposed method increased contrast-to-noise ratio (CNR) by more than 2.5 times and achieved <1.2% error for electron density values. On the anthropomorphic head phantom, the proposed method greatly improved the soft-tissue contrast and the fine detail differentiation ability. In the selected ROIs on different human tissues, the differences between the CT number obtained by the proposed method and that by the two-scan method were less than 4 HU. In the material images, the proposed method suppressed noise by over 75.5% compared with two-scan results, and by over 40.4% compared with HYPR-LR results. Implementation of the whole algorithm took 44.5 s for volumetric imaging, including projection preprocessing, FBP reconstruction, joint bilateral filtering, and material decomposition.

CONCLUSIONS

Using down-sampled projections in single-scan DE-CBCT, the proposed method could effectively and efficiently produce high-quality DE-CBCT images and low-noise material decomposition images. This method demonstrated superior performance on spatial resolution enhancement, NPS preservation, noise reduction, and electron density accuracy, indicating better prospect in material differentiation and dose calculation.

Image-domain Material Decomposition for Spectral CT using a Generalized Dictionary Learning

The spectral computed tomography (CT) has huge advantages by providing accurate material information. Unfortunately, due to the instability or overdetermination of material decomposition model, the accuracy of material decomposition can be compromised in practice. Very recently, the dictionary learning based image-domain material decomposition (DLIMD) can obtain high accuracy for material decompositions from reconstructed spectral CT images. This method can explore the correlation of material components to some extent by training a unified dictionary from all material images. In addition, the dictionary learning based prior as a penalty is applied on material components independently, and many parameters would be carefully elaborated in practice. Because the concentration of contrast agent in clinical applications is low, it can result in data inconsistency for dictionary based representation during the iteration process. To avoid the aforementioned limitations and further improve the accuracy of materials, we first construct a generalized dictionary learning based image-domain material decomposition (GDLIMD) model. Then, the material tensor image is unfolded along the mode-1 to enhance the correlation of different materials. Finally, to avoid the data inconsistency of low iodine contrast, a normalization strategy is employed. Both physical phantom and tissue-synthetic phantom experiments demonstrate the proposed GDLIMD method outperforms the DLIMD and direct inversion (DI) methods.

Head and neck multi-organ segmentation on dual-energy CT using dual pyramid convolutional neural networks

Organ delineation is crucial to diagnosis and therapy, while it is also labor-intensive and observer-dependent. Dual energy CT (DECT) provides additional image contrast than conventional single energy CT (SECT), which may facilitate automatic organ segmentation. This work aims to develop an automatic multi-organ segmentation approach using deep learning for head-and-neck region on DECT. We proposed a mask scoring regional convolutional neural network (R-CNN) where comprehensive features are firstly learnt from two independent pyramid networks and are then combined via deep attention strategy to highlight the informative ones extracted from both two channels of low and high energy CT. To perform multi-organ segmentation and avoid misclassification, a mask scoring subnetwork was integrated into the Mask R-CNN framework to build the correlation between the class of potential detected organ's region-of-interest (ROI) and the shape of that organ's segmentation within that ROI. We evaluated our model on DECT images from 127 head-and-neck cancer patients (66 training, 61 testing) with manual contours of 19 organs as training target and ground truth. For large- and mid-sized organs such as brain and parotid, the proposed method successfully achieved average Dice similarity coefficient (DSC) larger than 0.8. For small-sized organs with very low contrast such as chiasm, cochlea, lens and optic nerves, the DSCs ranged between around 0.5 and 0.8. With the proposed method, using DECT images outperforms using SECT in almost all 19 organs with statistical significance in DSC (p<0.05). Meanwhile, by using the DECT, the proposed method is also significantly superior to a recently developed FCN-based method in most of organs in terms of DSC and the 95th percentile Hausdorff distance. Quantitative results demonstrated the feasibility of the proposed method, the superiority of using DECT to SECT, and the advantage of the proposed R-CNN over FCN on the head-and-neck patient study. The proposed method has the potential to facilitate the current head-and-neck cancer radiation therapy workflow in treatment planning.

A material decomposition method for dual-energy CT via dual interactive Wasserstein generative adversarial networks

PURPOSE

Dual-energy computed tomography (DECT) is highly promising for material characterization and identification, whereas reconstructed material-specific images are affected by magnified noise and beam-hardening artifacts. Although various DECT material decomposition methods have been proposed to solve this problem, the quality of the decomposed images is still unsatisfactory, particularly in the image edges. In this study, a data-driven approach using dual interactive Wasserstein generative adversarial networks (DIWGAN) is developed to improve DECT decomposition accuracy and perform edge-preserving images.

METHODS

In proposed DIWGAN, two interactive generators are used to synthesize decomposed images of two basis materials by modeling the spatial and spectral correlations from input DECT reconstructed images, and the corresponding discriminators are employed to distinguish the difference between the generated images and labels. The DECT images reconstructed from high- and low-energy bins are sent to two generators separately, and each generator synthesizes one material-specific image, thereby ensuring the specificity of the network modeling. In addition, the information from different energy bins is exploited through the feature sharing of two generators. During decomposition model training, a hybrid loss function including L₁ loss, edge loss, and adversarial loss is incorporated to preserve the texture and edges in the generated images. Additionally, a selector is employed to define the generator that should be trained in each iteration, which can ensure the modeling ability of two different generators and improve the material decomposition accuracy. The performance of the proposed method is evaluated using digital phantom, XCAT phantom, and real data from a mouse.

RESULTS

On the digital phantom, the regions of bone and soft tissue are strictly and accurately separated using the trained decomposition model. The material densities in different bone and soft-tissue regions are near the ground truth, and the error of material densities is lower than 3 mg/ml. The results from XCAT phantom show that the material-specific images generated by directed matrix inversion and iterative decomposition methods have severe noise and artifacts. Regarding to the learning-based methods, the decomposed images of fully convolutional network (FCN) and butterfly network (Butterfly-Net) still contain varying degrees of artifacts, while proposed DIWGAN can yield high quality images. Compared to Butterfly-Net, the root-mean-square error (RMSE) of soft-tissue images generated by the DIWGAN decreased by 0.01 g/ml, whereas the peak-signal-to-noise ratio (PSNR) and structural similarity (SSIM) of the soft-tissue images reached 31.43 dB and 0.9987, respectively. The mass densities of the decomposed materials are nearest to the ground truth when using the DIWGAN method. The noise standard deviation of the decomposition images reduced by 69%, 60%, 33%, and 21% compared with direct matrix inversion, iterative decomposition, FCN, and Butterfly-Net, respectively. Furthermore, the performance of the mouse data indicates the potential of the proposed material decomposition method in real scanned data.

CONCLUSIONS

A DECT material decomposition method based on deep learning is proposed, and the relationship between reconstructed and material-specific images is mapped by training the DIWGAN model. Results from both the simulation phantoms and real data demonstrate the advantages of this method in suppressing noise and beam-hardening artifacts.

Learning-Based Stopping Power Mapping on Dual-Energy CT for Proton Radiation Therapy

Purpose

Dual-energy computed tomography (DECT) has been used to derive relative stopping power (RSP) maps by obtaining the energy dependence of photon interactions. The DECT-derived RSP maps could potentially be compromised by image noise levels and the severity of artifacts when using physics-based mapping techniques. This work presents a noise-robust learning-based method to predict RSP maps from DECT for proton radiation therapy.

Materials and Methods

The proposed method uses a residual attention cycle-consistent generative adversarial network to bring DECT-to-RSP mapping close to a 1-to-1 mapping by introducing an inverse RSP-to-DECT mapping. To evaluate the proposed method, we retrospectively investigated 20 head-and-neck cancer patients with DECT scans in proton radiation therapy simulation. Ground truth RSP values were assigned by calculation based on chemical compositions and acted as learning targets in the training process for DECT datasets; they were evaluated against results from the proposed method using a leave-one-out cross-validation strategy.

Results

The predicted RSP maps showed an average normalized mean square error of 2.83% across the whole body volume and an average mean error less than 3% in all volumes of interest. With additional simulated noise added in DECT datasets, the proposed method still maintained a comparable performance, while the physics-based stoichiometric method suffered degraded inaccuracy from increased noise level. The average differences from ground truth in dose volume histogram metrics for clinical target volumes were less than 0.2 Gy for D_95% and D_max with no statistical significance. Maximum difference in dose volume histogram metrics of organs at risk was around 1 Gy on average.

Conclusion

These results strongly indicate the high accuracy of RSP maps predicted by our machine-learning–based method and show its potential feasibility for proton treatment planning and dose calculation.

A review on medical imaging synthesis using deep learning and its clinical applications

This paper reviewed the deep learning-based studies for medical imaging synthesis and its clinical application. Specifically, we summarized the recent developments of deep learning-based methods in inter- and intra-modality image synthesis by listing and highlighting the proposed methods, study designs, and reported performances with related clinical applications on representative studies. The challenges among the reviewed studies were then summarized with discussion.

Machine learning in quantitative PET: A review of attenuation correction and low-count image reconstruction methods

The rapid expansion of machine learning is offering a new wave of opportunities for nuclear medicine. This paper reviews applications of machine learning for the study of attenuation correction (AC) and low-count image reconstruction in quantitative positron emission tomography (PET). Specifically, we present the developments of machine learning methodology, ranging from random forest and dictionary learning to the latest convolutional neural network-based architectures. For application in PET attenuation correction, two general strategies are reviewed: 1) generating synthetic CT from MR or non-AC PET for the purposes of PET AC, and 2) direct conversion from non-AC PET to AC PET. For low-count PET reconstruction, recent deep learning-based studies and the potential advantages over conventional machine learning-based methods are presented and discussed. In each application, the proposed methods, study designs and performance of published studies are listed and compared with a brief discussion. Finally, the overall contributions and remaining challenges are summarized.

Reconstruction method for DECT with one half-scan plus a second limited-angle scan using prior knowledge of complementary support set (Pri-CSS)

Dual-energy computed tomography (DECT) has capability to improve material differentiation, but most scanning schemes require two sets of full-scan measurements at different x-ray spectra, limiting its application to imaging system with incomplete scan. In this study, using one half-scan and a second limited-angle scan, we propose a DECT reconstruction method by exploiting the consistent information of gradient images at high- and low-energy spectra, which relaxes the requirement of data acquisition of DECT. Based on the theory of sampling condition analysis, the complementary support set of gradient images plays an important role in image reconstruction because it constitutes the sufficient and necessary condition for accurate CT reconstruction. For DECT, the gradient images of high- and low-energy CT images ideally share the same complementary support set for the same object. Inspired by this idea, we extract the prior knowledge of complementary support set (Pri-CSS) from the gradient image of the first half-scan CT image to promote the second limited-angle CT reconstruction. Pri-CSS will be incorporated into total variation regularization model in the form of constrains. Alternative direction method is applied to iteratively solve the modified optimization model, thereby deriving the proposed algorithm to recover low-energy CT image from limited-angle measurements. The qualitative and quantitative experiments on digital and real data are performed to validate the proposed method. The results show that the proposed method outperforms its counterparts and achieve high reconstruction quality for the designed scanning configuration.

Statistical image-based material decomposition for triple-energy computed tomography using total variation regularization

BACKGROUND

Triple-energy computed tomography (TECT) can obtain x-ray attenuation measurements at three energy spectra, thereby allowing identification of different material compositions with same or very similar attenuation coefficients. This ability is known as material decomposition, which can decompose TECT images into different basis material image. However, the basis material image would be severely degraded when material decomposition is directly performed on the noisy TECT measurements using a matrix inversion method.

OBJECTIVE

To achieve high quality basis material image, we present a statistical image-based material decomposition method for TECT, which uses the penalized weighted least-squares (PWLS) criteria with total variation (TV) regularization (PWLS-TV).

METHODS

The weighted least-squares term involves the noise statistical properties of the material decomposition process, and the TV regularization penalizes differences between local neighboring pixels in a decomposed image, thereby contributing to improving the quality of the basis material image. Subsequently, an alternating optimization method is used to minimize the objective function.

RESULTS

The performance of PWLS-TV is quantitatively evaluated using digital and mouse thorax phantoms. The experimental results show that PWLS-TV material decomposition method can greatly improve the quality of decomposed basis material image compared to the quality of images obtained using the competing methods in terms of suppressing noise and preserving edge and fine structure details.

CONCLUSIONS

The PWLS-TV method can simultaneously perform noise reduction and material decomposition in one iterative step, and it results in a considerable improvement of basis material image quality.

Deep learning-based image quality improvement for low-dose computed tomography simulation in radiation therapy

Low-dose computed tomography (CT) is desirable for treatment planning and simulation in radiation therapy. Multiple rescanning and replanning during the treatment course with a smaller amount of dose than a single conventional full-dose CT simulation is a crucial step in adaptive radiation therapy. We developed a machine learning-based method to improve image quality of low-dose CT for radiation therapy treatment simulation. We used a residual block concept and a self-attention strategy with a cycle-consistent adversarial network framework. A fully convolution neural network with residual blocks and attention gates (AGs) was used in the generator to enable end-to-end transformation. We have collected CT images from 30 patients treated with frameless brain stereotactic radiosurgery (SRS) for this study. These full-dose images were used to generate projection data, which were then added with noise to simulate the low-mAs scanning scenario. Low-dose CT images were reconstructed from this noise-contaminated projection data and were fed into our network along with the original full-dose CT images for training. The performance of our network was evaluated by quantitatively comparing the high-quality CT images generated by our method with the original full-dose images. When mAs is reduced to 0.5% of the original CT scan, the mean square error of the CT images obtained by our method is ∼ 1.6 % , with respect to the original full-dose images. The proposed method successfully improved the noise, contract-to-noise ratio, and nonuniformity level to be close to those of full-dose CT images and outperforms a state-of-the-art iterative reconstruction method. Dosimetric studies show that the average differences of dose-volume histogram metrics are < 0.1    Gy ( p > 0.05 ). These quantitative results strongly indicate that the denoised low-dose CT images using our method maintains image accuracy and quality and are accurate enough for dose calculation in current CT simulation of brain SRS treatment. We also demonstrate the great potential for low-dose CT in the process of simulation and treatment planning.

Low-dose material-specific radiography using monoenergetic photons

Cargo containers constitute the most critical component of global trade: 108 million containers represent the movement of about 95% of the world’s manufactured goods. The steady increase in cargo container shipments has had a profound effect on world security: the threat associated with smuggling of shielded special nuclear material is elevated every year. Containers reaching the borders of the U.S. are currently not radiographically inspected due to time and dose considerations stemming from the use of bremsstrahlung beams for imaging. Bremsstrahlung spectra are low-energy peaked, resulting in low penetration values, especially through dense cargoes. The use of monoenergetic radiography beams could alleviate many of these problems due to higher energy and low background continuum. Using Monte Carlo simulations of a realistic imaging scenario with support from previous experimental measurements, we demonstrate how the use of monoenergetic photon beams in radiography can simultaneously reduce the radiation dose imparted to the cargo and any potential stowaways while increasing image quality. Dual-energy methods are leveraged to calculate material atomic number. Image quality is evaluated by measuring the noise standard deviation, contrast-to-noise ratio, and the pixel error as the dose is decreased.

Development of a deep neural network for generating synthetic dual-energy chest x-ray images with single x-ray exposure

Dual-energy chest radiography (DECR) is a medical imaging technology that can improve diagnostic accuracy. This technique can decompose single-energy chest radiography (SECR) images into separate bone- and soft tissue-only images. This can, however, double the radiation exposure to the patient. To address this limitation, we developed an algorithm for the synthesis of DECR from a SECR through deep learning. To predict high resolution images, we developed a novel deep learning architecture by modifying a conventional U-net to take advantage of the high frequency-dominant information that propagates from the encoding part to the decoding part. In addition, we used the anticorrelated relationship (ACR) of DECR for improving the quality of the predicted images. For training data, 300 pairs of SECR and their corresponding DECR images were used. To test the trained model, 50 DECR images from Yonsei University Severance Hospital and 662 publicly accessible SECRs were used. To evaluate the performance of the proposed method, we compared DECR and predicted images using a structural similarity approach (SSIM). In addition, we quantitatively evaluated image quality calculating the modulation transfer function and coefficient of variation. The proposed model selectively predicted the bone- and soft tissue-only CR images from an SECR image. The strategy for improving the spatial resolution by ACR was effective. Quantitative evaluation showed that the proposed method with ACR showed relatively high SSIM (over 0.85). In addition, predicted images with the proposed ACR model achieved better image quality measures than those of U-net. In conclusion, the proposed method can obtain high-quality bone- and soft tissue-only CR images without the need for additional hardware for double x-ray exposures in clinical practice.

Multi-Energy X-Ray Tomography of an Optical Fiber: The Role of Spatial Averaging

Using a commercial X-ray tomography instrument, we have obtained reconstructions of a graded-index optical fiber with voxels of edge length 1.05 µm at 12 tube voltages. The fiber manufacturer created a graded index in the central region by varying the germanium concentration from a peak value in the center of the core to a very small value at the core-cladding boundary. Operating on 12 tube voltages, we show by a singular value decomposition that there are only two singular vectors with significant weight. Physically, this means scans beyond two tube voltages contain largely redundant information. We concentrate on an analysis of the images associated with these two singular vectors. The first singular vector is dominant and images of the coefficients of the first singular vector at each voxel look are similar to any of the single-energy reconstructions. Images of the coefficients of the second singular vector by itself appear to be noise. However, by averaging the reconstructed voxels in each of several narrow bands of radii, we can obtain values of the second singular vector at each radius. In the core region, where we expect the germanium doping to go from a peak value at the fiber center to zero at the core-cladding boundary, we find that a plot of the two coefficients of the singular vectors forms a line in the two-dimensional space consistent with the dopant decreasing linearly with radial distance from the core center. The coating, made of a polymer rather than silica, is not on this line indicating that the two-dimensional results are sensitive not only to the density but also to the elemental composition.

Optimal virtual monoenergetic image in "TwinBeam" dual-energy CT for organs-at-risk delineation based on contrast-noise-ratio in head-and-neck radiotherapy

PURPOSE

Dual-energy computed tomography (DECT) using TwinBeam CT (TBCT) is a new option for radiation oncology simulators. TBCT scanning provides virtual monoenergetic images which are attractive in treatment planning since lower energies offer better contrast for soft tissues, and higher energies reduce noise. A protocol is needed to achieve optimal performance of this feature. In this study, we investigated the TBCT scan schema with the head-and-neck radiotherapy workflow at our clinic and selected the optimal energy with best contrast-noise-ratio (CNR) in organs-at-risks (OARs) delineation for head-and-neck treatment planning.

METHODS AND MATERIALS

We synthesized monochromatic images from 40 keV to 190 keV at 5 keV increments from data acquired by TBCT. We collected the Hounsfield unit (HU) numbers of OARs (brainstem, mandible, spinal cord, and parotid glands), the HU numbers of marginal regions outside OARs, and the noise levels for each monochromatic image. We then calculated the CNR for the different OARs at each energy level to generate a serial of spectral curves for each OAR. Based on these spectral curves of CNR, the mono-energy corresponding to the max CNR was identified for each OAR of each patient.

RESULTS

Computed tomography scans of ten patients by TBCT were used to test the optimal monoenergetic image for the CNR of OAR. Based on the maximized CNR, the optimal energy values were 78.5 ± 5.3 keV for the brainstem, 78.0 ± 4.2 keV for the mandible, 78.5 ± 5.7 keV for the parotid glands, and 78.5 ± 5.3 keV for the spinal cord. Overall, the optimal energy for the maximum CNR of these OARs in head-and-neck cancer patients was 80 keV.

CONCLUSION

We have proposed a clinically feasible protocol that selects the optimal energy level of the virtual monoenergetic image in TBCT for OAR delineation based on the CNR in head-and-neck OAR. This protocol can be applied in TBCT simulation.

Statistical CT reconstruction using region-aware texture preserving regularization learning from prior normal-dose CT image

In some clinical applications, prior normal-dose CT (NdCT) images are available, and the valuable textures and structure features in them may be used to promote follow-up low-dose CT (LdCT) reconstruction. This study aims to learn texture information from the NdCT images and leverage it for follow-up LdCT image reconstruction to preserve textures and structure features. Specifically, the proposed reconstruction method first learns the texture information from those patches with similar structures in NdCT image, and the similar patches can be clustered by searching context features efficiently from the surroundings of the current patch. Then it utilizes redundant texture information from the similar patches as a priori knowledge to describe specific regions in the LdCT image. The advanced region-aware texture preserving prior is termed as 'RATP'. The main advantage of the PATP prior is that it can properly learn the texture features from available NdCT images and adaptively characterize the region-specific structures in the LdCT image. The experiments using patient data were performed to evaluate the performance of the proposed method. The proposed RATP method demonstrated superior performance in LdCT imaging compared to the filtered back projection (FBP) and statistical iterative reconstruction (SIR) methods using Gaussian regularization, Huber regularization and the original texture preserving regularization.

Single-Scan Dual-Energy CT Using Primary Modulation

Compared with conventional computed tomography (CT), dual-energy CT (DECT) provides better material differentiation but requires projection data acquired with two different effective x-ray spectra, limiting DECT applications to specialized scanners. We propose a hardware-based method, known as PM-DECT, which utilizes primary beam modulation to enable single-scan DECT on a conventional CT scanner. PM-DECT inserts an attenuation sheet with a spatially varying pattern-primary beam modulator-between the x-ray source and imaged object. During a CT scan, the modulator selectively hardens the x-ray beam, thereby increasing the average photon energy at specific detector pixel locations. Thus, PM-DECT simultaneously acquires high and low energy data at each projection angle. From the sparse projection data, high and low energy CT images are jointly reconstructed and simultaneously decomposed into basis materials via an iterative CT reconstruction algorithm with gradient weighting and an improved version of similarity based regularization. Studies on Catphan 600 and anthropomorphic head phantoms demonstrate that PM-DECT retains a high level of spatial resolution compared with conventional CT scans. Electron density values calculated from decomposed images indicate a limited error of 1.12% for PM-DECT. Comparison against a two-scan DECT technique shows that PM-DECT's image reconstruction from sparse data sets contributes only 0.66% error. By granting the opportunity for high-quality single-scan DECT on conventional CT scanners via limited hardware modification, PM-DECT has the potential to liberate DECT from specialized scanners, extending clinical availability, and implementation.

Iterative reconstruction for low dose dual energy CT using information-divergence constrained spectral redundancy information

Dual energy computed tomography (DECT) can improve the capability of differentiating different materials compared with conventional CT. However, due to non-negligible radiation exposure to patients, dose reduction has recently become a critical concern in CT imaging field. In this work, to reduce noise at the same time maintain DECT images quality, we present an iterative reconstruction algorithm for low-dose DECT images where in the objective function of the algorithm consists of a data-fidelity term and a regularization term. The former term is based on alpha-divergence to describe the statistical distribution of the DE sinogram data. And the latter term is based on the redundant information to reflect the prior information of the desired DECT images. For simplicity, the presented algorithm is termed as "AlphaD-aviNLM". To minimize the associative objective function, a modified proximal forward-backward splitting algorithm is proposed. Digital phantom, physical phantom, and patient data were utilized to validate and evaluate the presented AlphaD-aviNLM algorithm. The experimental results characterize the performance of the presented AlphaD-aviNLM algorithm. Speficically, in the digital phantom study, the presented AlphaD-aviNLM algorithm performs better than the PWLS-TV, PWLS-aviNLM, and AlphaD-TV with more than 49%, 34%, and 40% gains for the RMSE metric, 1.3%, 0.4%, and 0.7% gains for the FSIM metric and 13%, 8%, and 11% gains for the PSNR metric. In the physical phantom study, the presented AlphaD-aviNLM algorithm performs better than the PWLS-TV, PWLS-aviNLM, and AlphaD-TV with more than 0.55%, 0.07%, and 0.16% gains for the FSIM metric.

An effective noise reduction method for multi-energy CT images that exploit spatio-spectral features

PURPOSE

To develop and evaluate an image-domain noise reduction method for multi-energy CT (MECT) data.

METHODS

Multi-Energy Non-Local Means (MENLM) is a technique that uses the redundant information in MECT images to achieve noise reduction. In this method, spatio-spectral features are used to determine the similarity between pixels, making the similarity evaluation more robust to image noise. The performance of this MENLM filter was tested on images acquired on a whole-body research photon counting CT system. The impact of filtering on image quality was quantitatively evaluated in phantom studies in terms of image noise level (standard deviation of pixel values), noise power spectrum (NPS), in-plane and cross-plane spatial resolution, CT number accuracy, material decomposition performance, and subjective low-contrast spatial resolution using the American College of Radiology (ACR) CT accreditation phantom. Clinical feasibility was assessed by performing MENLM on contrast-enhanced swine images and unenhanced cadaver head images using clinically relevant doses and dose rates.

RESULTS

The phantom studies demonstrated that the MENLM filter reduced noise substantially and still preserved the shape and peak frequency of the NPS. With 80% noise reduction, MENLM filtering caused no degradation of high-contrast spatial resolution, as illustrated by the modulation transfer function (MTF) and slice sensitivity profile (SSP). CT number accuracy was also maintained for all energy channels, demonstrating that energy resolution was not affected by filtering. Material decomposition performance was improved with MENLM filtering. The subjective evaluation using the ACR phantom demonstrated an improvement in low-contrast performance. MENLM achieved effective noise reduction in both contrast-enhanced swine images and unenhanced cadaver head images, resulting in improved detection of subtle vascular structures and the differentiation of white/gray matter.

CONCLUSION

In MECT, MENLM achieved around 80% noise reduction and greatly improved material decomposition performance and the detection of subtle anatomical/low-contrast features while maintaining spatial and energy resolution. MENLM filtering may improve diagnostic or functional analysis accuracy and facilitate radiation dose and contrast media reduction for MECT.

Statistical image-domain multimaterial decomposition for dual-energy CT

PURPOSE

Dual-energy CT (DECT) enhances tissue characterization because of its basis material decomposition capability. In addition to conventional two-material decomposition from DECT measurements, multimaterial decomposition (MMD) is required in many clinical applications. To solve the ill-posed problem of reconstructing multi-material images from dual-energy measurements, additional constraints are incorporated into the formulation, including volume and mass conservation and the assumptions that there are at most three materials in each pixel and various material types among pixels. The recently proposed flexible image-domain MMD method decomposes pixels sequentially into multiple basis materials using a direct inversion scheme which leads to magnified noise in the material images. In this paper, we propose a statistical image-domain MMD method for DECT to suppress the noise.

METHODS

The proposed method applies penalized weighted least-square (PWLS) reconstruction with a negative log-likelihood term and edge-preserving regularization for each material. The statistical weight is determined by a data-based method accounting for the noise variance of high- and low-energy CT images. We apply the optimization transfer principles to design a serial of pixel-wise separable quadratic surrogates (PWSQS) functions which monotonically decrease the cost function. The separability in each pixel enables the simultaneous update of all pixels.

RESULTS

The proposed method is evaluated on a digital phantom, Catphan©600 phantom and three patients (pelvis, head, and thigh). We also implement the direct inversion and low-pass filtration methods for a comparison purpose. Compared with the direct inversion method, the proposed method reduces noise standard deviation (STD) in soft tissue by 95.35% in the digital phantom study, by 88.01% in the Catphan©600 phantom study, by 92.45% in the pelvis patient study, by 60.21% in the head patient study, and by 81.22% in the thigh patient study, respectively. The overall volume fraction accuracy is improved by around 6.85%. Compared with the low-pass filtration method, the root-mean-square percentage error (RMSE(%)) of electron densities in the Catphan©600 phantom is decreased by 20.89%. As modulation transfer function (MTF) magnitude decreased to 50%, the proposed method increases the spatial resolution by an overall factor of 1.64 on the digital phantom, and 2.16 on the Catphan©600 phantom. The overall volume fraction accuracy is increased by 6.15%.

CONCLUSIONS

We proposed a statistical image-domain MMD method using DECT measurements. The method successfully suppresses the magnified noise while faithfully retaining the quantification accuracy and anatomical structure in the decomposed material images. The proposed method is practical and promising for advanced clinical applications using DECT imaging.

Dual energy CT with one full scan and a second sparse-view scan using structure preserving iterative reconstruction (SPIR)

Conventional dual-energy CT (DECT) reconstruction requires two full-size projection datasets with two different energy spectra. In this study, we propose an iterative algorithm to enable a new data acquisition scheme which requires one full scan and a second sparse-view scan for potential reduction in imaging dose and engineering cost of DECT. A bilateral filter is calculated as a similarity matrix from the first full-scan CT image to quantify the similarity between any two pixels, which is assumed unchanged on a second CT image since DECT scans are performed on the same object. The second CT image from reduced projections is reconstructed by an iterative algorithm which updates the image by minimizing the total variation of the difference between the image and its filtered image by the similarity matrix under data fidelity constraint. As the redundant structural information of the two CT images is contained in the similarity matrix for CT reconstruction, we refer to the algorithm as structure preserving iterative reconstruction (SPIR). The proposed method is evaluated on both digital and physical phantoms, and is compared with the filtered-backprojection (FBP) method, the conventional total-variation-regularization-based algorithm (TVR) and prior-image-constrained-compressed-sensing (PICCS). SPIR with a second 10-view scan reduces the image noise STD by a factor of one order of magnitude with same spatial resolution as full-view FBP image. SPIR substantially improves over TVR on the reconstruction accuracy of a 10-view scan by decreasing the reconstruction error from 6.18% to 1.33%, and outperforms TVR at 50 and 20-view scans on spatial resolution with a higher frequency at the modulation transfer function value of 10% by an average factor of 4. Compared with the 20-view scan PICCS result, the SPIR image has 7 times lower noise STD with similar spatial resolution. The electron density map obtained from the SPIR-based DECT images with a second 10-view scan has an average error of less than 1%.