Monte Carlo calculation based on hydrogen composition of the tissue for MV photon radiotherapy

被引:8
作者
Demol, Benjamin [1 ,2 ,3 ]
Viard, Romain [2 ]
Reynaert, Nick [1 ]
机构
[1] Ctr Oscar Lambret, Dept Radiotherapy, F-59000 Lille, France
[2] AQUILAB SAS, Loos Les Lille, France
[3] IEMN, UMR CNRS 8520, Villeneuve Dascq, France
来源
JOURNAL OF APPLIED CLINICAL MEDICAL PHYSICS | 2015年 / 16卷 / 05期
关键词
Monte Carlo; hydrogen content; stoichiometric calibration; megavoltage photon radiotherapy; magnetic resonance imaging; MRI-ONLY RADIOTHERAPY; ELECTRON-DENSITY; MAGNETIC-RESONANCE; RADIATION-THERAPY; DOSE CALCULATIONS; PSEUDO-CT; ECHO-TIME; SEQUENCES; BRAIN; UNITS;
D O I
10.1120/jacmp.v16i5.5586
中图分类号
R8 [特种医学]; R445 [影像诊断学];
学科分类号
1002 ; 100207 ; 1009 ;
摘要
The purpose of this study was to demonstrate that Monte Carlo treatment planning systems require tissue characterization (density and composition) as a function of CT number. A discrete set of tissue classes with a specific composition is introduced. In the current work we demonstrate that, for megavoltage photon radiotherapy, only the hydrogen content of the different tissues is of interest. This conclusion might have an impact on MRI-based dose calculations and on MVCT calibration using tissue substitutes. A stoichiometric calibration was performed, grouping tissues with similar atomic composition into 15 dosimetrically equivalent subsets. To demonstrate the importance of hydrogen, a new scheme was derived, with correct hydrogen content, complemented by oxygen (all elements differing from hydrogen are replaced by oxygen). Mass attenuation coefficients and mass stopping powers for this scheme were calculated and compared to the original scheme. Twenty-five CyberKnife treatment plans were recalculated by an in-house developed Monte Carlo system using tissue density and hydrogen content derived from the CT images. The results were compared to Monte Carlo simulations using the original stoichiometric calibration. Between 300 keV and 3 MeV, the relative difference of mass attenuation coefficients is under 1% within all subsets. Between 10 keV and 20 MeV, the relative difference of mass stopping powers goes up to 5% in hard bone and remains below 2% for all other tissue subsets. Dose-volume histograms (DVHs) of the treatment plans present no visual difference between the two schemes. Relative differences of dose indexes D-98, D-95, D-50, D-05, D-02, and D-mean were analyzed and a distribution centered around zero and of standard deviation below 2% (3 sigma) was established. On the other hand, once the hydrogen content is slightly modified, important dose differences are obtained. Monte Carlo dose planning in the field of megavoltage photon radiotherapy is fully achievable using only hydrogen content of tissues, a conclusion that might impact MRI dose calculation, but can also help selecting the optimal tissue substitutes when calibrating MVCT devices.
引用
收藏
页码:117 / 130
页数:14
相关论文
共 27 条
[1]   GEANT4-a simulation toolkit [J].
Agostinelli, S ;
Allison, J ;
Amako, K ;
Apostolakis, J ;
Araujo, H ;
Arce, P ;
Asai, M ;
Axen, D ;
Banerjee, S ;
Barrand, G ;
Behner, F ;
Bellagamba, L ;
Boudreau, J ;
Broglia, L ;
Brunengo, A ;
Burkhardt, H ;
Chauvie, S ;
Chuma, J ;
Chytracek, R ;
Cooperman, G ;
Cosmo, G ;
Degtyarenko, P ;
Dell'Acqua, A ;
Depaola, G ;
Dietrich, D ;
Enami, R ;
Feliciello, A ;
Ferguson, C ;
Fesefeldt, H ;
Folger, G ;
Foppiano, F ;
Forti, A ;
Garelli, S ;
Giani, S ;
Giannitrapani, R ;
Gibin, D ;
Cadenas, JJG ;
González, I ;
Abril, GG ;
Greeniaus, G ;
Greiner, W ;
Grichine, V ;
Grossheim, A ;
Guatelli, S ;
Gumplinger, P ;
Hamatsu, R ;
Hashimoto, K ;
Hasui, H ;
Heikkinen, A ;
Howard, A .
NUCLEAR INSTRUMENTS & METHODS IN PHYSICS RESEARCH SECTION A-ACCELERATORS SPECTROMETERS DETECTORS AND ASSOCIATED EQUIPMENT, 2003, 506 (03) :250-303
[2]   Patch-based generation of a pseudo CT from conventional MRI sequences for MRI-only radiotherapy of the brain [J].
Andreasen, Daniel ;
Van Leemput, Koen ;
Hansen, Rasmus H. ;
Andersen, Jon A. L. ;
Edmund, Jens M. .
MEDICAL PHYSICS, 2015, 42 (04) :1596-1605
[3]  
[Anonymous], 2011, PIRS794REVB NRCC
[4]  
Berger M., NIST, PML, DOI [DOI 10.18434/T48G6X, 10.18434/T48G6X]
[5]  
Berger M. J., 2005, ESTAR, PSTAR, and ASTAR: Computer Programs for Calculating Stopping-Power and Range Tables for Electrons, Protons, and Helium Ions
[6]   On the alpha particles of radium, and their loss of range in passing through various atoms and molecules [J].
Bragg, W. H. ;
Elder, M. A. .
PHILOSOPHICAL MAGAZINE, 1905, 10 (55-60) :318-340
[7]   An Atlas-Based Electron Density Mapping Method for Magnetic Resonance Imaging (MRI)-Alone Treatment Planning and Adaptive MRI-Based Prostate Radiation Therapy [J].
Dowling, Jason A. ;
Lambert, Jonathan ;
Parker, Joel ;
Salvado, Olivier ;
Fripp, Jurgen ;
Capp, Anne ;
Wratten, Chris ;
Denham, James W. ;
Greer, Peter B. .
INTERNATIONAL JOURNAL OF RADIATION ONCOLOGY BIOLOGY PHYSICS, 2012, 83 (01) :E5-E11
[8]   A voxel-based investigation for MRI-only radiotherapy of the brain using ultra short echo times [J].
Edmund, Jens M. ;
Kjer, Hans M. ;
Van Leemput, Koen ;
Hansen, Rasmus H. ;
Andersen, Jon A. L. ;
Andreasen, Daniel .
PHYSICS IN MEDICINE AND BIOLOGY, 2014, 59 (23) :7501-7519
[9]   A unifying probabilistic Bayesian approach to derive electron density from MRI for radiation therapy treatment planning [J].
Gudur, Madhu Sudhan Reddy ;
Hara, Wendy ;
Le, Quynh-Thu ;
Wang, Lei ;
Xing, Lei ;
Li, Ruijiang .
PHYSICS IN MEDICINE AND BIOLOGY, 2014, 59 (21) :6595-6606
[10]   X-RAY ATTENUATION COEFFICIENTS OF ELEMENTS AND MIXTURES [J].
JACKSON, DF ;
HAWKES, DJ .
PHYSICS REPORTS-REVIEW SECTION OF PHYSICS LETTERS, 1981, 70 (03) :169-233