Calculation of scatter in cone beam CT : Steps towards a virtual tomograph

Abstract

Scattered photons—shortly scatter—are generated by interaction processes when photon beams interact with matter. In diagnostic radiology, they deteriorate image quality since they add an undesirable signal that lowers the contrast in projection radiography and causes cupping and streak artefacts in computed tomography (CT). Scatter is one of the most detrimental factors in cone beam CT owing to irradiation geometries using wide beams. It cannot be fully eliminated, nevertheless its amount can be lowered via scatter reduction techniques (air gaps, antiscatter grids, collimators) and its effect on medical images can be suppressed via scatter correction algorithms. Aim: Develop a tool—a virtual tomograph—that simulates projections and performs image reconstructions similarly to a real CT scanner. Use this tool to evaluate the effect of scatter on projections and reconstructed images in cone beam CT. Propose improvements in CT scanner design and image reconstruction algorithms. Methods: A software toolkit (CTmod) based on the application development framework ROOT was written to simulate primary and scatter projections using analytic and Monte Carlo methods, respectively. It was used to calculate the amount of scatter in cone beam CT for anthropomorphic voxel phantoms and water cylinders. Configurations with and without bowtie filters, antiscatter grids, and beam hardening corrections were investigated. Filtered back-projection was used to reconstruct images. Automatic threshold segmentation of volumetric CT data of anthropomorphic phantoms with known tissue compositions was tested to evaluate its usability in an iterative image reconstruction algorithm capable of performing scatter correction. Results: It was found that computer speed was the limiting factor for the deployment of this method in clinical CT scanners. It took several hours to calculate a single projection depending on the complexity of the geometry, number of simulated detector elements, and statistical precision. Data calculated using the CTmod code confirmed the already known facts that the amount of scatter is almost linearly proportional to the beam width, the scatter-to-primary ratio (SPR) can be larger than 1 for body-size objects, and bowtie filters can decrease the SPR in certain regions of projections. Ideal antiscatter grids significantly lowered the amount of scatter. The beneficial effect of classical antiscatter grids in cone beam CT with flat panel imagers was not confirmed by other researchers nevertheless new grid designs are still being tested. A simple formula estimating the effect of scatter on the quality of reconstructed images was suggested and tested. Conclusions: It was shown that computer simulations could calculate the amount of scatter in diagnostic radiology. The Monte Carlo method was too slow for a routine use in contemporary clinical practice nevertheless it could be used to optimize CT scanner design and, with some enhancements, it could become a part of an image reconstruction algorithm that performs scatter correction