Particle distribution estimation is an important issue in medical diagnosis. In particular, photon scattering in some medical devices extremely degrades image quality and causes measurement inaccuracy. The Monte Carlo (MC) algorithm is regarded as the most accurate particle estimation approach but is still time-consuming, even with graphic processing unit (GPU) acceleration. The goal of this work is to develop an automatic scatter estimation framework for high-efficiency photon distribution estimation. Specifically, a GPU-based MC simulation initially yields a raw scatter signal with a low photon number to hasten scatter generation. In the proposed method, assume that the scatter signal follows Poisson distribution, where an optimization objective function fused with sparse feature penalty is modeled. Then, an over-relaxation algorithm is deduced mathematically to solve this objective function. For optimizing the parameters in the over-relaxation algorithm, the deep Q -network in the deep reinforcement learning scheme is built to intelligently interact with the over-relaxation algorithm to accurately and rapidly estimate a scatter signal with the large range of photon numbers. Experimental results demonstrated that our proposed framework can achieve superior performance with structural similarity > 0.94 , peak signal-to-noise ratio > 26.55 dB , and relative absolute error < 5.62 % , and the lowest computation time for one scatter image generation can be within 2 s.
BackgroundScattering 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.PurposeAiming 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.MethodsOur 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.ResultsThe 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.ConclusionsIn 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.
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. Scatterto-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 fullfan. 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.
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.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.