X-Ray Scatter Estimation Using Deep Splines

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 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.

A practical and robust method for beam blocker-based cone beam CT scatter correction

Objective. In the traditional beam-blocker based cone beam CT (CBCT) scatter correction, the scatter measured in the region shaded by lead strips was multiplied by a correction factor to directly represent the scatter in the unblocked region. The correction factor optimization is a tedious process and lacks an objective stop criterion. To skip the optimization process, an indirect scatter estimation method was developed and validated in phantom imaging.Approach.A beam-blocker made of lead strips was mounted between the x-ray source and object for scatter estimation. The primary signal between lead strips in the blocked region was first calculated by subtracting the measured scatter, and then used to calculate the scatter signal in the unblocked region corresponding to the same attenuation path. The calculated scatter signal was smoothed via local filtration and used to correct the measured projection in the unblocked region. Finally, the CBCT was reconstructed via Feldkamp-Davis-Kress algorithm. A Catphan and a head phantom were used to verify the performance of the proposed method in both full- and half-blocker scenarios, and with and without a bow-tie filter.Main Results. For scans without the bow-tie filter, the CT number error was reduced to 3.97±2.27 and 5.51±3.90 HU in the full- and half-blocker scenarios, respectively, for the Catphan, and to 4.01±2.18 and 7.97 ± 4.05 HU for the head phantom. When the bow-tie filter was applied, the CT number error was reduced to 2.29±1.42 and 6.72±0.77 HU in the full- and half-blocker scenarios, respectively, for the Catphan, and 2.35±1.25 and 4.96 ± 1.89 HU for the head phantom.Significance. The proposed method effectively avoids the influence of the inserted beam blocker itself on the scatter intensity estimation, and proves a more practical and robust way for the beam-blocker based scatter correction in CBCT scanning.

Deep learning for x-ray scatter correction in dedicated breast CT

BACKGROUND

Accurate correction of x-ray scatter in dedicated breast computed tomography (bCT) imaging may result in improved visual interpretation and is crucial to achieve quantitative accuracy during image reconstruction and analysis.

PURPOSE

To develop a deep learning (DL) model to correct for x-ray scatter in bCT projection images.

METHODS

A total of 115 patient scans acquired with a bCT clinical system were segmented into the major breast tissue types (skin, adipose, and fibroglandular tissue). The resulting breast phantoms were divided into training (n = 110) and internal validation cohort (n = 5). Training phantoms were augmented by a factor of four by random translation of the breast in the image field of view. Using a previously validated Monte Carlo (MC) simulation algorithm, 12 primary and scatter bCT projection images with a 30-degree step were generated from each phantom. For each projection, the thickness map and breast location in the field of view were also calculated. A U-Net based DL model was developed to estimate the scatter signal based on the total input simulated image and trained single-projection-wise, with the thickness map and breast location provided as additional inputs. The model was internally validated using MC-simulated projections and tested using an external data set of 10 phantoms derived from images acquired with a different bCT system. For this purpose, the mean relative difference (MRD) and mean absolute error (MAE) were calculated. To test for accuracy in reconstructed images, a full bCT acquisition was mimicked with MC-simulations and then assessed by calculating the MAE and the structural similarity (SSIM). Subsequently, scatter was estimated and subtracted from the bCT scans of three patients to obtain the scatter-corrected image. The scatter-corrected projections were reconstructed and compared with the uncorrected reconstructions by evaluating the correction of the cupping artifact, increase in image contrast, and contrast-to-noise ratio (CNR).

RESULTS

The mean MRD and MAE across all cases (min, max) for the internal validation set were 0.04% (-1.1%, 1.3%) and 2.94% (2.7%, 3.2%), while for the external test set they were -0.64% (-1.6%, 0.2%) and 2.84% (2.3%, 3.5%), respectively. For MC-simulated reconstruction slices, the computed SSIM was 0.99 and the MAE was 0.11% (range: 0%, 0.35%) with a single outlier slice of 2.06%. For the three patient bCT reconstructed images, the correction increased the contrast by a mean of 25% (range: 20%, 30%), and reduced the cupping artifact. The mean CNR increased by 0.32 after scatter correction, which was not found to be significant (95% confidence interval: [-0.01, 0.65], p = 0.059). The time required to correct the scatter in a single bCT projection was 0.2 s on an NVIDIA GeForce GTX 1080 GPU.

CONCLUSION

The developed DL model could accurately estimate scatter in bCT projection images and could enhance contrast and correct for cupping artifact in reconstructed patient images without significantly affecting the CNR. The time required for correction would allow its use in daily clinical practice, and the reported accuracy will potentially allow quantitative reconstructions.

Learning-based occupational x-ray scatter estimation

Objective. During x-ray-guided interventional procedures, the medical staff is exposed to scattered ionizing radiation caused by the patient. To increase the staff’s awareness of the invisible radiation and monitor dose online, computational scatter estimation methods are convenient. However, such methods are usually based on Monte Carlo (MC) simulations, which are inherently computationally expensive. Yet, in the interventional environment, immediate feedback to the personnel is desirable. Approach. In this work, we propose deep neural networks to mitigate the computational effort of MC simulations. Our learning-based models consider detailed models of the (outer) patient shape and (inner) anatomy, additional objects in the room, and the x-ray tube spectrum to cover imaging settings encountered in real interventional settings. We investigate two cases of scatter prediction. First, we employ network architectures to estimate the full three-dimensional (3D) scatter distribution. Second, we investigate the prediction of two-dimensional (2D) intensity projections that facilitate the intra-procedural visualization. Main results. Depending on the dimensionality of the estimated scatter distribution and the network architecture, the mean relative error of each network is in the range of 12% and 14% compared to MC simulations. However, 3D scatter distributions can be estimated within 60 ms and 2D distributions within 15 ms. Significance. Overall, our method is suitable to support the online assessment of scattered ionizing radiation in the interventional environment and can help to lower the occupational radiation risk.