A hybrid approach for fast simulation of x-ray computed tomography

With the advent of powerful computers the use of Monte Carlo (MC) techniques for radiation transport simulation has become one of the most popular methods for modeling x-ray computed tomography (CT) scanners. One significant problem in the use of MC calculations is the presence of statistical uncertainties in the estimates. A simple but not practical way to decrease statistical uncertainties is to run MC simulations for sufficiently long time (large number of histories) and use efficient variance reduction techniques. In this paper, we propose a hybrid approach combining MC and analytical simulations to speed up x-ray CT modeling. In addition, an analytic method for modeling of beam hardening effect (BHE) in CT is proposed for fast generation of raw data for use in validation of BHE correction algorithms. In this method, the contribution of the primary component to the projections is calculated through analytical simulations using ray-tracing methods whereas the contribution of the scatter component is calculated using pure MC simulations. The final projections are determined by appropriate combination of both simulation results. The hybrid simulator and analytic BHE modelling were validated through comparison with experimental measurements and pure MC simulation of different phantoms performed on clinical CT scanners. In general, there is good agreement between the simulated and measured projections and reconstructed images. This work has demonstrated the feasibility of using the hybrid approach for fast and accurate modeling of CT scanners thus allowing to evaluate the effect of physical, geometrical and other design parameters on CT scanner performance and image quality. The results also validate the proposed method for analytical modeling of BHE which is a useful tool for fast generation of raw CT datasets.