A CT-based Monte Carlo Dose Calculations for Proton Therapy Using a New Interface Program
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.
[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.
[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.
@article{"International Journal of Engineering, Mathematical and Physical Sciences:56826", author = "A. Esmaili Torshabi and A. Terakawa and K. Ishii and H. Yamazaki and S. Matsuyama and Y. Kikuchi and M. Nakhostin and H. Sabet and A. Ishizaki and W. Yamashita and T. Togashi and J. Arikawa and H. Akiyama and K. Koyata", title = "A CT-based Monte Carlo Dose Calculations for Proton Therapy Using a New Interface Program", 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", volume = "3", number = "5", pages = "341-6", }