Monte Carlo calculation of radiation dose in CT examinations using phantom and patient tomographic models

By using a voxel-based Monte Carlo simulation technique, we developed and validated a method to calculate radiation-absorbed dose in the computed tomography (CT) examinations from the images of phantoms and patients. The ionising radiation transport was simulated using the EGS4 code system. The geometry of the X-ray beam (focus-to-axis distance, field of view, collimation, and primary and beam-shaper filtration) and the X-ray spectral distribution (HiSpeed LX/i) were included in the simulation. Each axial CT image was reduced to a 256 x 256 matrix and stacked in a volume. The patient images were segmented before the simulation of radiation transport by using four categories of materials, such as air, lung, muscle and bone. To test the voxel-based method, the values of the radiation dose derived from a simulated CT exposure were calculated and compared with those obtained from the measurements performed within the dosimetry phantoms. To complete the scope of the work, series of CT scans of the trunk of an anthropomorphic phantom and patients were simulated to calculate the average dose in each 1-cm-wide transverse slice (ADS). The comparison between the simulated and measured dose data for the CT indices showed a difference of <5% in all the cases. The estimated mean values of ADS from the chest, abdomen and pelvis of the anthropomorphic phantom were approximately 1.7-2 times the weighted CT dose index (CTDI(w)) value, whereas the mean ADS values for these anatomical areas were 1.3-2 times the CTDI(w) of patients. The voxel-based simulation method provided a technique for estimating the individual patient doses in the CT examinations.