Hybrid scatter correction for CT imaging

The purpose of this study was to develop and evaluate the hybrid scatter correction algorithm (HSC) for CT imaging. Therefore, two established ways to perform scatter correction, i.e. physical scatter correction based on Monte Carlo simulations and a convolution-based scatter correction algorithm, were combined in order to perform an object-dependent, fast and accurate scatter correction. Based on a reconstructed CT volume, patient-specific scatter intensity is estimated by a coarse Monte Carlo simulation that uses a reduced amount of simulated photons in order to reduce the simulation time. To further speed up the Monte Carlo scatter estimation, scatter intensities are simulated only for a fraction of all projections. In a second step, the high noise estimate of the scatter intensity is used to calibrate the open parameters in a convolution-based algorithm which is then used to correct measured intensities for scatter. Furthermore, the scatter-corrected intensities are used in order to reconstruct a scatter-corrected CT volume data set. To evaluate the scatter reduction potential of HSC, we conducted simulations in a clinical CT geometry and measurements with a flat detector CT system. In the simulation study, HSC-corrected images were compared to scatter-free reference images. For the measurements, no scatter-free reference image was available. Therefore, we used an image corrected with a low-noise Monte Carlo simulation as a reference. The results show that the HSC can significantly reduce scatter artifacts. Compared to the reference images, the error due to scatter artifacts decreased from 100% for uncorrected images to a value below 20% for HSC-corrected images for both the clinical (simulated data) and the flat detector CT geometry (measurement). Compared to a low-noise Monte Carlo simulation, with the HSC the number of photon histories can be reduced by about a factor of 100 per projection without losing correction accuracy. Furthermore, it was sufficient to calibrate the parameters in the convolution model at an angular increment of about 20°. The reduction of the simulated photon histories together with the reduced amount of simulated Monte Carlo scatter projections decreased the total runtime of the scatter correction by about two orders of magnitude for the cases investigated here when using the HSC instead of a low-noise Monte Carlo simulation for scatter correction.