Purpose: To quantify the radiation dose by organ received by the patient during the acquisition of CT images using a Monte Carlo method and GPU parallelization with the purpose of including dose data in the patient medical history.
Methods: The Rando phantom was prepared locating inside MOSFET dosemeters to obtain 3D absolute dose points, which are then compared with the values obtained by the MC-GPU simulation.The authors have developed a MATLAB interface, which performs the required simulations automatically from the CT images acquired, including the positions where the experimental values have been obtained. The MATLAB code developed, called CT-dose, automatically creates a MC-GPU voxelized model from the CT images taking into account the Hounsfield numbers. After that, considering the DICOM data the program creates also automatically the rest of the input needed in the MC-GPU simulation and it calls MC-GPU to simulate. Finally, CT-dose delivers a table with the dose data by organ and a 3D visualization of the dose distribution, which can be represented in the program interface or using the PARAVIEW software.
Results: The simulation results in the points where the MOSFET dosimeters were located inside the phantom shows good agreement with respect to experimental values obtained. Due to the punctual characterization of the dose, one could extrapolate that the dose distribution calculation should be accurate. However, more experimental validation with depth dose curves would be necessary for the final validation of the calculations.
Conclusion: This work presents an interface that could be useful to calculate the organ dose data of the patients after CT scans with an accurate Monte Carlo methodology. The parallelization of the MC-GPU code in graphical processors allows the possibility of using this methodology in real cases, which is interesting considering for example the EURATOM 2013/59 Directive.