PND-Net: Physics-Inspired Non-Local Dual-Domain Network for Metal Artifact Reduction.

Metal artifacts caused by the presence of metallic implants tremendously degrade the quality of reconstructed computed tomography (CT) images and therefore affect the clinical diagnosis or reduce the accuracy of organ delineation and dose calculation in radiotherapy. Although various deep learning methods have been proposed for metal artifact reduction (MAR), most of them aim to restore the corrupted sinogram within the metal trace, which removes beam hardening artifacts but ignores other components of metal artifacts. In this paper, based on the physical property of metal artifacts which is verified via Monte Carlo (MC) simulation, we propose a novel physics-inspired non-local dual-domain network (PND-Net) for MAR in CT imaging. Specifically, we design a novel non-local sinogram decomposition network (NSD-Net) to acquire the weighted artifact component and develop an image restoration network (IR-Net) to reduce the residual and secondary artifacts in the image domain. To facilitate the generalization and robustness of our method on clinical CT images, we employ a trainable fusion network (F-Net) in the artifact synthesis path to achieve unpaired learning. Furthermore, we design an internal consistency loss to ensure the data fidelity of anatomical structures in the image domain and introduce the linear interpolation sinogram as prior knowledge to guide sinogram decomposition. NSD-Net, IR-Net, and F-Net are jointly trained so that they can benefit from one another. Extensive experiments on simulation and clinical data demonstrate that our method outperforms state-of-the-art MAR methods.

Adaptive scatter kernel deconvolution modeling for cone-beam CT scatter correction via deep reinforcement learning.

BACKGROUND: Scattering photons can seriously contaminate cone-beam CT (CBCT) image quality with severe artifacts and substantial degradation of CT value accuracy, which is a major concern limiting the widespread application of CBCT in the medical field. The scatter kernel deconvolution (SKD) method commonly used in clinic requires a Monte Carlo (MC) simulation to determine numerous quality-related kernel parameters, and it cannot realize intelligent scatter kernel parameter optimization, causing limited accuracy of scatter estimation. PURPOSE: Aiming at improving the scatter estimation accuracy of the SKD algorithm, an intelligent scatter correction framework integrating the SKD with deep reinforcement learning (DRL) scheme is proposed. METHODS: Our method firstly builds a scatter kernel model to iteratively convolve with raw projections, and then the deep Q-network of the DRL scheme is introduced to intelligently interact with the scatter kernel to achieve a projection adaptive parameter optimization. The potential of the proposed framework is demonstrated on CBCT head and pelvis simulation data and experimental CBCT measurement data. Furthermore, we have implemented the U-net based scatter estimation approach for comparison. RESULTS: The simulation study demonstrates that the mean absolute percentage error (MAPE) of the proposed method is less than 9.72% and the peak signal-to-noise ratio (PSNR) is higher than 23.90 dB, while for the conventional SKD algorithm, the minimum MAPE is 17.92% and the maximum PSNR is 19.32 dB. In the measurement study, we adopt a hardware-based beam stop array algorithm to obtain the scatter-free projections as a comparison baseline, and our method can achieve superior performance with MAPE < 17.79% and PSNR > 16.34 dB. CONCLUSIONS: In this paper, we propose an intelligent scatter correction framework that integrates the physical scatter kernel model with DRL algorithm, which has the potential to improve the accuracy of the clinical scatter correction method to obtain better CBCT imaging quality.

A transmission measurement-based spectrum estimation method incorporating X-ray tube voltage fluctuation.

Background: The energy spectrum is the property of the X-ray tube that describes the energy fluence per unit interval of photon energy. The existing indirect methods for estimating the spectrum ignore the influence caused by the voltage fluctuation of the X-ray tube. Methods: In this work, we propose a method for estimating the X-ray energy spectrum more accurately by including the voltage fluctuation of the X-ray tube. It expresses the spectrum as the weighted summation of a set of model spectra within a certain voltage fluctuation range. The difference between the raw projection and the estimated projection is considered as the objective function for obtaining the corresponding weight of each model spectrum. The equilibrium optimizer (EO) algorithm is used to find the weight combination that minimizes the objective function. Finally, the estimated spectrum is obtained. We refer to the proposed method as the poly-voltage method. The method is mainly aimed at the cone-beam computed tomography (CBCT) system. Results: The model spectra mixture evaluation and projection evaluation showed that the reference spectrum can be combined by multiple model spectra. They also showed that it is appropriate to choose about 10% of the preset voltage as the voltage range of the model spectra, which can match the reference spectrum and projection quite well. The phantom evaluation showed that the beam-hardening artifact can be corrected using the estimated spectrum via the poly-voltage method, and the poly-voltage method provides not only the accurate reprojection but also an accurate spectrum. The normalized root mean square error (NRMSE) index between the spectrum generated via the poly-voltage method and the reference spectrum could be kept within 3% according to above evaluations. There existed a 1.77% percentage error between the estimated scatter of polymethyl methacrylate (PMMA) phantom using the two spectra generated via the poly-voltage method and the single-voltage method, and it could be considered for scatter simulation. Conclusions: Our proposed poly-voltage method could estimate the spectrum more accurately for both ideal and more realistic voltage spectra, and it is robust to the different modes of voltage pulse.

Deep learning based de-overlapping correction of projections from a flat-panel micro array X-ray source: Simulation study.

PURPOSE: Flat-panel X-ray source is an experimental X-ray emitter with target application of static computer tomography (CT), which can save imaging space and time. However, the X-ray cone beams emitted by the densely arranged micro-ray sources are overlapped, causing serious structural overlapping and visual blur in the projection results. Traditional deoverlapping methods can hardly solve this problem well. METHOD: We converted the overlapping cone beam projections to parallel beam projections through a U-like neural network and selected structural similarity (SSIM) loss as the loss function. In this study, we converted three kinds of overlapping cone beam projections of the Shepp-Logan, line-pairs, and abdominal data with two overlapping levels to corresponding parallel beam projections. Training completed, we tested the model using the test set data that was not used at the training phase, and evaluated the difference between the test set conversion results and their corresponding parallel beams through three indicators: mean squared error (MSE), peak signal-to-noise ratio (PSNR) and SSIM. In addition, projections from head phantoms were applied for generalization test. RESULT: In the Shepp-Logan low-overlapping task, we obtained a MSE of 1.624×10-5, a PSNR of 47.892 dB, and a SSIM of 0.998 which are the best results of the six experiments. For the most challenging abdominal task, the MSE, PSNR, and SSIM are 1.563×10-3, 28.0586 dB, and 0.983, respectively. In more generalized data, the model also achieved good results. CONCLUSION: This study proves the feasibility of utilizing the end-to-end U-net for deblurring and deoverlapping in the flat-panel X-ray source domain.

A correlated sampling-based Monte Carlo simulation for fast CBCT iterative scatter correction.

BACKGROUND: In recent years, cone-beam computed tomography (CBCT) has played an important role in medical imaging. However, the applications of CBCT are limited due to the severe scatter contamination. Conventional Monte Carlo (MC) simulation can provide accurate scatter estimation for scatter correction, but the expensive computational cost has always been the bottleneck of MC method in clinical application. PURPOSE: In this work, an MC simulation method combined with a variance reduction technique called correlated sampling is proposed for fast iterative scatter correction. METHODS: Correlated sampling exploits correlation between similar simulation systems to reduce the variance of interest quantities. Specifically, conventional MC simulation is first performed on the scatter-contaminated CBCT to generate the initial scatter signal. In the subsequent correction iterations, scatter estimation is then updated by applying correlated MC sampling to the latest corrected CBCT images by reusing the random number sequences of the task-related photons in conventional MC. Afterward, the corrected projections obtained by subtracting the scatter estimation from raw projections are utilized for FDK reconstruction. These steps are repeated until an adequate scatter correction is obtained. The performance of the proposed framework is evaluated by the accuracy of the scatter estimation, the quality of corrected CBCT images and efficiency. RESULTS: Overall, the difference in mean absolute percentage error between scatter estimation with and without correlated sampling is 0.25% for full-fan case and 0.34% for half-fan case, respectively. In simulation studies, scatter artifacts are substantially eliminated, where the mean absolute error value is reduced from 15 to 2 HU in full-fan case and from 53 to 13 HU in half-fan case. Scatter-to-primary ratio is reduced to 0.02 for full-fan and 0.04 for half-fan, respectively. In phantom study, the contrast-to-noise ratio (CNR) is increased by a factor of 1.63, and the contrast is increased by a factor of 1.77. As for clinical studies, the CNR is improved by 11% and 14% for half-fan and full-fan, respectively. The contrast after correction is increased by 19% for half-fan and 44% for full-fan. Furthermore, root mean square error is also effectively reduced, especially from 78 to 4 HU for full-fan. Experimental results demonstrate that the figure of merit is improved between 23 and 43 folds when using correlated sampling. The proposed method takes less than 25 s for the whole iterative scatter correction process. CONCLUSIONS: The proposed correlated sampling-based MC simulation method can achieve fast and accurate scatter correction for CBCT, making it suitable for real-time clinical use.

Multiresolution residual deep neural network for improving pelvic CBCT image quality.

PURPOSE: Cone-beam computed tomography (CBCT) is frequently used for accurate image-guided radiation therapy. However, the poor CBCT image quality prevents its further clinical use. Thus, it is important to improve the HU accuracy and structure preservation of CBCT images. METHODS: In this study, we proposed a novel method to generate synthetic CT (sCT) images from CBCT images. A multiresolution residual deep neural network (RDNN) was adopted for image regression from CBCT images to planning CT (pCT) images. At the coarse level, RDNN was first trained with a large amount of lower resolution images, which can make the network focus on coarse information and prevent overfitting problems. More fine information was obtained gradually by fine-tuning the coarse model using fewer number of higher resolution images. Our model was optimized by using aligned pCT and CBCT image pairs of a particular body region of 153 prostate cancer patients treated in our hospital (120 for training and 33 for testing). Five-fold cross-validation was used to tune the hyperparameters and the testing data were used to evaluate the performance of the final models. RESULTS: The mean absolute error (MAE) between CBCT and pCT on the testing data was 352.56 HU, while the MAE between the sCT and pCT images was 52.18 HU for our proposed multiresolution RDNN model, which reduced the MAE by 85.20% (p < 0.01). In addition, the average structural similarity index measure between the sCT and CBCT was 19.64% (p = 0.01) higher than that of pCT and CBCT. CONCLUSIONS: The sCT images generated using our proposed multiresolution RDNN have higher HU accuracy and structural fidelity, which may promote the further applications of CBCT images in the clinic for structure segmentation, dose calculation, and adaptive radiotherapy planning.