Use of the FLUKA Monte Carlo code for 3D patient-specific dosimetry on PET-CT and SPECT-CT images

F. Botta, A. Mairani, Robert Hobbs, A. Vergara Gil, M. Pacilio, K. Parodi, M. Cremonesi, M. A Coca Pérez, A. Di Dia, M. Ferrari, F. Guerriero, G. Battistoni, G. Pedroli, G. Paganelli, L. A Torres Aroche, George Sgouros

Research output: Contribution to journalArticle

Abstract

Patient-specific absorbed dose calculation for nuclear medicine therapy is a topic of increasing interest. 3D dosimetry at the voxel level is one of the major improvements for the development of more accurate calculation techniques, as compared to the standard dosimetry at the organ level. This study aims to use the FLUKA Monte Carlo code to perform patient-specific 3D dosimetry through direct Monte Carlo simulation on PET-CT and SPECT-CT images. To this aim, dedicated routines were developed in the FLUKA environment. Two sets of simulations were performed on model and phantom images. Firstly, the correct handling of PET and SPECT images was tested under the assumption of homogeneous water medium by comparing FLUKA results with those obtained with the voxel kernel convolution method and with other Monte Carlo-based tools developed to the same purpose (the EGS-based 3D-RD software and the MCNP5-based MCID). Afterwards, the correct integration of the PET/SPECT and CT information was tested, performing direct simulations on PET/CT images for both homogeneous (water) and non-homogeneous (water with air, lung and bone inserts) phantoms. Comparison was performed with the other Monte Carlo tools performing direct simulation as well. The absorbed dose maps were compared at the voxel level. In the case of homogeneous water, by simulating 108 primary particles a 2% average difference with respect to the kernel convolution method was achieved; such difference was lower than the statistical uncertainty affecting the FLUKA results. The agreement with the other tools was within 3-4%, partially ascribable to the differences among the simulation algorithms. Including the CT-based density map, the average difference was always within 4% irrespective of the medium (water, air, bone), except for a maximum 6% value when comparing FLUKA and 3D-RD in air. The results confirmed that the routines were properly developed, opening the way for the use of FLUKA for patient-specific, image-based dosimetry in nuclear medicine.

Original languageEnglish (US)
Pages (from-to)8099-8120
Number of pages22
JournalPhysics in Medicine and Biology
Volume58
Issue number22
DOIs
StatePublished - Nov 21 2013

Fingerprint

Water
Nuclear Medicine
Air
Bone and Bones
Single-Photon Emission-Computed Tomography
Uncertainty
Software
Single Photon Emission Computed Tomography Computed Tomography
Positron Emission Tomography Computed Tomography
Lung
Therapeutics

ASJC Scopus subject areas

  • Radiology Nuclear Medicine and imaging
  • Radiological and Ultrasound Technology

Cite this

Use of the FLUKA Monte Carlo code for 3D patient-specific dosimetry on PET-CT and SPECT-CT images. / Botta, F.; Mairani, A.; Hobbs, Robert; Vergara Gil, A.; Pacilio, M.; Parodi, K.; Cremonesi, M.; Pérez, M. A Coca; Di Dia, A.; Ferrari, M.; Guerriero, F.; Battistoni, G.; Pedroli, G.; Paganelli, G.; Aroche, L. A Torres; Sgouros, George.

In: Physics in Medicine and Biology, Vol. 58, No. 22, 21.11.2013, p. 8099-8120.

Research output: Contribution to journalArticle

Botta, F, Mairani, A, Hobbs, R, Vergara Gil, A, Pacilio, M, Parodi, K, Cremonesi, M, Pérez, MAC, Di Dia, A, Ferrari, M, Guerriero, F, Battistoni, G, Pedroli, G, Paganelli, G, Aroche, LAT & Sgouros, G 2013, 'Use of the FLUKA Monte Carlo code for 3D patient-specific dosimetry on PET-CT and SPECT-CT images', Physics in Medicine and Biology, vol. 58, no. 22, pp. 8099-8120. https://doi.org/10.1088/0031-9155/58/22/8099
Botta, F. ; Mairani, A. ; Hobbs, Robert ; Vergara Gil, A. ; Pacilio, M. ; Parodi, K. ; Cremonesi, M. ; Pérez, M. A Coca ; Di Dia, A. ; Ferrari, M. ; Guerriero, F. ; Battistoni, G. ; Pedroli, G. ; Paganelli, G. ; Aroche, L. A Torres ; Sgouros, George. / Use of the FLUKA Monte Carlo code for 3D patient-specific dosimetry on PET-CT and SPECT-CT images. In: Physics in Medicine and Biology. 2013 ; Vol. 58, No. 22. pp. 8099-8120.
@article{fc9cecddf80e44108abc14617c9e4c9b,
title = "Use of the FLUKA Monte Carlo code for 3D patient-specific dosimetry on PET-CT and SPECT-CT images",
abstract = "Patient-specific absorbed dose calculation for nuclear medicine therapy is a topic of increasing interest. 3D dosimetry at the voxel level is one of the major improvements for the development of more accurate calculation techniques, as compared to the standard dosimetry at the organ level. This study aims to use the FLUKA Monte Carlo code to perform patient-specific 3D dosimetry through direct Monte Carlo simulation on PET-CT and SPECT-CT images. To this aim, dedicated routines were developed in the FLUKA environment. Two sets of simulations were performed on model and phantom images. Firstly, the correct handling of PET and SPECT images was tested under the assumption of homogeneous water medium by comparing FLUKA results with those obtained with the voxel kernel convolution method and with other Monte Carlo-based tools developed to the same purpose (the EGS-based 3D-RD software and the MCNP5-based MCID). Afterwards, the correct integration of the PET/SPECT and CT information was tested, performing direct simulations on PET/CT images for both homogeneous (water) and non-homogeneous (water with air, lung and bone inserts) phantoms. Comparison was performed with the other Monte Carlo tools performing direct simulation as well. The absorbed dose maps were compared at the voxel level. In the case of homogeneous water, by simulating 108 primary particles a 2{\%} average difference with respect to the kernel convolution method was achieved; such difference was lower than the statistical uncertainty affecting the FLUKA results. The agreement with the other tools was within 3-4{\%}, partially ascribable to the differences among the simulation algorithms. Including the CT-based density map, the average difference was always within 4{\%} irrespective of the medium (water, air, bone), except for a maximum 6{\%} value when comparing FLUKA and 3D-RD in air. The results confirmed that the routines were properly developed, opening the way for the use of FLUKA for patient-specific, image-based dosimetry in nuclear medicine.",
author = "F. Botta and A. Mairani and Robert Hobbs and {Vergara Gil}, A. and M. Pacilio and K. Parodi and M. Cremonesi and P{\'e}rez, {M. A Coca} and {Di Dia}, A. and M. Ferrari and F. Guerriero and G. Battistoni and G. Pedroli and G. Paganelli and Aroche, {L. A Torres} and George Sgouros",
year = "2013",
month = "11",
day = "21",
doi = "10.1088/0031-9155/58/22/8099",
language = "English (US)",
volume = "58",
pages = "8099--8120",
journal = "Physics in Medicine and Biology",
issn = "0031-9155",
publisher = "IOP Publishing Ltd.",
number = "22",

}

TY - JOUR

T1 - Use of the FLUKA Monte Carlo code for 3D patient-specific dosimetry on PET-CT and SPECT-CT images

AU - Botta, F.

AU - Mairani, A.

AU - Hobbs, Robert

AU - Vergara Gil, A.

AU - Pacilio, M.

AU - Parodi, K.

AU - Cremonesi, M.

AU - Pérez, M. A Coca

AU - Di Dia, A.

AU - Ferrari, M.

AU - Guerriero, F.

AU - Battistoni, G.

AU - Pedroli, G.

AU - Paganelli, G.

AU - Aroche, L. A Torres

AU - Sgouros, George

PY - 2013/11/21

Y1 - 2013/11/21

N2 - Patient-specific absorbed dose calculation for nuclear medicine therapy is a topic of increasing interest. 3D dosimetry at the voxel level is one of the major improvements for the development of more accurate calculation techniques, as compared to the standard dosimetry at the organ level. This study aims to use the FLUKA Monte Carlo code to perform patient-specific 3D dosimetry through direct Monte Carlo simulation on PET-CT and SPECT-CT images. To this aim, dedicated routines were developed in the FLUKA environment. Two sets of simulations were performed on model and phantom images. Firstly, the correct handling of PET and SPECT images was tested under the assumption of homogeneous water medium by comparing FLUKA results with those obtained with the voxel kernel convolution method and with other Monte Carlo-based tools developed to the same purpose (the EGS-based 3D-RD software and the MCNP5-based MCID). Afterwards, the correct integration of the PET/SPECT and CT information was tested, performing direct simulations on PET/CT images for both homogeneous (water) and non-homogeneous (water with air, lung and bone inserts) phantoms. Comparison was performed with the other Monte Carlo tools performing direct simulation as well. The absorbed dose maps were compared at the voxel level. In the case of homogeneous water, by simulating 108 primary particles a 2% average difference with respect to the kernel convolution method was achieved; such difference was lower than the statistical uncertainty affecting the FLUKA results. The agreement with the other tools was within 3-4%, partially ascribable to the differences among the simulation algorithms. Including the CT-based density map, the average difference was always within 4% irrespective of the medium (water, air, bone), except for a maximum 6% value when comparing FLUKA and 3D-RD in air. The results confirmed that the routines were properly developed, opening the way for the use of FLUKA for patient-specific, image-based dosimetry in nuclear medicine.

AB - Patient-specific absorbed dose calculation for nuclear medicine therapy is a topic of increasing interest. 3D dosimetry at the voxel level is one of the major improvements for the development of more accurate calculation techniques, as compared to the standard dosimetry at the organ level. This study aims to use the FLUKA Monte Carlo code to perform patient-specific 3D dosimetry through direct Monte Carlo simulation on PET-CT and SPECT-CT images. To this aim, dedicated routines were developed in the FLUKA environment. Two sets of simulations were performed on model and phantom images. Firstly, the correct handling of PET and SPECT images was tested under the assumption of homogeneous water medium by comparing FLUKA results with those obtained with the voxel kernel convolution method and with other Monte Carlo-based tools developed to the same purpose (the EGS-based 3D-RD software and the MCNP5-based MCID). Afterwards, the correct integration of the PET/SPECT and CT information was tested, performing direct simulations on PET/CT images for both homogeneous (water) and non-homogeneous (water with air, lung and bone inserts) phantoms. Comparison was performed with the other Monte Carlo tools performing direct simulation as well. The absorbed dose maps were compared at the voxel level. In the case of homogeneous water, by simulating 108 primary particles a 2% average difference with respect to the kernel convolution method was achieved; such difference was lower than the statistical uncertainty affecting the FLUKA results. The agreement with the other tools was within 3-4%, partially ascribable to the differences among the simulation algorithms. Including the CT-based density map, the average difference was always within 4% irrespective of the medium (water, air, bone), except for a maximum 6% value when comparing FLUKA and 3D-RD in air. The results confirmed that the routines were properly developed, opening the way for the use of FLUKA for patient-specific, image-based dosimetry in nuclear medicine.

UR - http://www.scopus.com/inward/record.url?scp=84896503871&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=84896503871&partnerID=8YFLogxK

U2 - 10.1088/0031-9155/58/22/8099

DO - 10.1088/0031-9155/58/22/8099

M3 - Article

C2 - 24200697

AN - SCOPUS:84896503871

VL - 58

SP - 8099

EP - 8120

JO - Physics in Medicine and Biology

JF - Physics in Medicine and Biology

SN - 0031-9155

IS - 22

ER -