A CT-based Monte Carlo Dose Calculations for Proton Therapy Using a New Interface Program

Authors: A. Esmaili Torshabi, A. Terakawa, K. Ishii, H. Yamazaki, S. Matsuyama, Y. Kikuchi, M. Nakhostin, H. Sabet, A. Ishizaki, W. Yamashita, T. Togashi, J. Arikawa, H. Akiyama, K. Koyata

Abstract:

The purpose of this study is to introduce a new interface program to calculate a dose distribution with Monte Carlo method in complex heterogeneous systems such as organs or tissues in proton therapy. This interface program was developed under MATLAB software and includes a friendly graphical user interface with several tools such as image properties adjustment or results display. Quadtree decomposition technique was used as an image segmentation algorithm to create optimum geometries from Computed Tomography (CT) images for dose calculations of proton beam. The result of the mentioned technique is a number of nonoverlapped squares with different sizes in every image. By this way the resolution of image segmentation is high enough in and near heterogeneous areas to preserve the precision of dose calculations and is low enough in homogeneous areas to reduce the number of cells directly. Furthermore a cell reduction algorithm can be used to combine neighboring cells with the same material. The validation of this method has been done in two ways; first, in comparison with experimental data obtained with 80 MeV proton beam in Cyclotron and Radioisotope Center (CYRIC) in Tohoku University and second, in comparison with data based on polybinary tissue calibration method, performed in CYRIC. These results are presented in this paper. This program can read the output file of Monte Carlo code while region of interest is selected manually, and give a plot of dose distribution of proton beam superimposed onto the CT images.

Keywords: Monte Carlo, CT images, Quadtree decomposition, Interface program, Proton beam

Digital Object Identifier (DOI): doi.org/10.5281/zenodo.1334351

References:

[1] S. Chiavassa, M. Bardies, F. Guiraud-Vitaux, D. Bruel, J. R. Jourdian, D. Franck, and I. Aubineau-Laniece, "OEDIPE: a personalized dosimetric tool associating voxel-based models with MCNPX," Cancer Biother., vol. 20, no. 3, pp. 325-332, 2005.
[2] K. A. Van Riper, "A CT AND MRI SCAN TO MCNP INPUT CONVERSION PROGRAM," Radiat Prot. Dosim., vol. 115, no. 1-4, pp. 513-516, Dec. 2005.
[3] C. Oliveira, H. Yoriyaz, F. C. Oliveira, and L. M. Ferreira, "Monte Carlo simulation for dose distribution calculations in a CT-based phantom at the Portuguese gamma irradiation facility," Nucl. Inst. and Meth. in Phys. Res. B, vol. 213, pp. 662-665, Jan. 2004.
[4] W. Ai-dong, W. Yi-can, T. Sheng-xiang, and Z. Jiang-hui, "Effect of CT Image-based Voxel Size On Monte Carlo Dose Calculation," in Proc. 27th Annu. Conf. Engineering in Medicine and Biology, Shanghai, 2005,pp. 6449-6451.
[5] N. Kanematsu, N. Matsufuji, R. Kohno, S. Minohara, and T. Kanai, "A CT calibration method based on the Polybinary tissue model for radiotherapy treatment planning," Phys. Med. Biol., vol. 48, pp. 1053- 1064, 2003.
[6] International Commission on Radiation Units and Measurements, "Tissue Substitutes in Radiation Dosimetry and Measurement," ICRU rep. no. 44 Maryland, Bethesda, 1989.
[7] R. Kohno, A. Nohtomi, Y. Takada, T. Terunuma, T. Sakae, and K. Matsumoto, "A compensation method of an imaging plate response to clinical proton beams," Nucl. Inst. and Meth. in Phys. Res. A, vol. 481, pp. 669-674, 2004.