Analysis of the bias induced by voxel and unstructured mesh Monte Carlo models for the MCNP6 code in orthovoltage applications

Lorenzo Isolan, Marco De Pietri, Mauro Iori, Andrea Botti, Elisabetta Cagni, Marco Sumini

Research output: Contribution to journalArticle

Abstract

The use of Monte Carlo methods for the set up of Treatment Planning Systems (TPS) in radiotherapy applications is a current standard. The most advanced modeling techniques aim at directly link the output of CT scans to a patient-specific model build up instead of standard phantoms and look-up tables. This link represents a critical step since even the most accurate segmentation and organ volume definition must be translated into a suitable input for the Monte Carlo code. During that step, the segmented volume is usually mapped on a regular geometrical lattice (voxels). A more sophisticated option appears to be a volume description based on an Unstructured Mesh (UM) geometry typical of current finite element codes. In this paper, we compared the two approaches analyzing the bias induced by the different choices. Starting from anonymous patient DICOM files coming from a CT scan, suitable segmentation and volume definition have been carried out and voxel and UM-based equivalent models for the Monte Carlo code MCNP6 have been built. Some computational phantoms, covering some significant portion of the human body (head and lower limb), have been used as a benchmark of the dose distribution obtained from X-ray sources commonly used in orthovoltage applications. The specific clinical application has been chosen because, despite being used on a growing number of patients and with multi-Gy doses, treatment planning is still based on poorly characterized dose mapping systems such as look-up tables or low-resolution computational phantoms for MC codes. Also, experimental measurements on phantom slabs irradiated by an X-ray source were carried out preliminarily to validate the simulated radiation source. As shown in the results, the UM computational phantoms (built through the Simpleware SCAN-IP™ tool) can reduce the bias induced by the regularity of the classical cubic voxel geometry, give a more accurate description of volumes and complex surfaces and, thanks to an optimized discretization of the volumes, are also able to reduce the computational work. For the same UM models, various comparisons with different voxel sizes have been produced to investigate the dose distributions and evaluate the voxelization effects. The comparison shows a convergence pattern of the voxel model to the UM one. It has been possible to verify that the most significant bias occurs where the dose gradient is higher, along the beam borders and at tissue interfaces. The simulations showed that the UM can be really effective and reliable to compute the dose distributions within computational anthropomorphic phantoms obtained from the patient’s CT scans.

Original languageEnglish
JournalRadiation Effects and Defects in Solids
DOIs
Publication statusPublished - Jan 1 2019

Fingerprint

Computerized tomography
mesh
dosage
Digital Imaging and Communications in Medicine (DICOM)
Planning
X rays
planning
Geometry
Radiotherapy
Dosimetry
Monte Carlo methods
human body
Tissue
geometry
limbs
radiation sources
files
borders
Radiation
regularity

Keywords

  • MCNP
  • Monte Carlo
  • orthovoltage
  • radiotherapy
  • unstructured mesh

ASJC Scopus subject areas

  • Radiation
  • Nuclear and High Energy Physics
  • Materials Science(all)
  • Condensed Matter Physics

Cite this

@article{b53133de4b8d44409841c0ab5284c351,
title = "Analysis of the bias induced by voxel and unstructured mesh Monte Carlo models for the MCNP6 code in orthovoltage applications",
abstract = "The use of Monte Carlo methods for the set up of Treatment Planning Systems (TPS) in radiotherapy applications is a current standard. The most advanced modeling techniques aim at directly link the output of CT scans to a patient-specific model build up instead of standard phantoms and look-up tables. This link represents a critical step since even the most accurate segmentation and organ volume definition must be translated into a suitable input for the Monte Carlo code. During that step, the segmented volume is usually mapped on a regular geometrical lattice (voxels). A more sophisticated option appears to be a volume description based on an Unstructured Mesh (UM) geometry typical of current finite element codes. In this paper, we compared the two approaches analyzing the bias induced by the different choices. Starting from anonymous patient DICOM files coming from a CT scan, suitable segmentation and volume definition have been carried out and voxel and UM-based equivalent models for the Monte Carlo code MCNP6 have been built. Some computational phantoms, covering some significant portion of the human body (head and lower limb), have been used as a benchmark of the dose distribution obtained from X-ray sources commonly used in orthovoltage applications. The specific clinical application has been chosen because, despite being used on a growing number of patients and with multi-Gy doses, treatment planning is still based on poorly characterized dose mapping systems such as look-up tables or low-resolution computational phantoms for MC codes. Also, experimental measurements on phantom slabs irradiated by an X-ray source were carried out preliminarily to validate the simulated radiation source. As shown in the results, the UM computational phantoms (built through the Simpleware SCAN-IP™ tool) can reduce the bias induced by the regularity of the classical cubic voxel geometry, give a more accurate description of volumes and complex surfaces and, thanks to an optimized discretization of the volumes, are also able to reduce the computational work. For the same UM models, various comparisons with different voxel sizes have been produced to investigate the dose distributions and evaluate the voxelization effects. The comparison shows a convergence pattern of the voxel model to the UM one. It has been possible to verify that the most significant bias occurs where the dose gradient is higher, along the beam borders and at tissue interfaces. The simulations showed that the UM can be really effective and reliable to compute the dose distributions within computational anthropomorphic phantoms obtained from the patient’s CT scans.",
keywords = "MCNP, Monte Carlo, orthovoltage, radiotherapy, unstructured mesh",
author = "Lorenzo Isolan and {De Pietri}, Marco and Mauro Iori and Andrea Botti and Elisabetta Cagni and Marco Sumini",
year = "2019",
month = "1",
day = "1",
doi = "10.1080/10420150.2019.1596104",
language = "English",
journal = "Radiation Effects and Defects in Solids",
issn = "1042-0150",
publisher = "Taylor and Francis Ltd.",

}

TY - JOUR

T1 - Analysis of the bias induced by voxel and unstructured mesh Monte Carlo models for the MCNP6 code in orthovoltage applications

AU - Isolan, Lorenzo

AU - De Pietri, Marco

AU - Iori, Mauro

AU - Botti, Andrea

AU - Cagni, Elisabetta

AU - Sumini, Marco

PY - 2019/1/1

Y1 - 2019/1/1

N2 - The use of Monte Carlo methods for the set up of Treatment Planning Systems (TPS) in radiotherapy applications is a current standard. The most advanced modeling techniques aim at directly link the output of CT scans to a patient-specific model build up instead of standard phantoms and look-up tables. This link represents a critical step since even the most accurate segmentation and organ volume definition must be translated into a suitable input for the Monte Carlo code. During that step, the segmented volume is usually mapped on a regular geometrical lattice (voxels). A more sophisticated option appears to be a volume description based on an Unstructured Mesh (UM) geometry typical of current finite element codes. In this paper, we compared the two approaches analyzing the bias induced by the different choices. Starting from anonymous patient DICOM files coming from a CT scan, suitable segmentation and volume definition have been carried out and voxel and UM-based equivalent models for the Monte Carlo code MCNP6 have been built. Some computational phantoms, covering some significant portion of the human body (head and lower limb), have been used as a benchmark of the dose distribution obtained from X-ray sources commonly used in orthovoltage applications. The specific clinical application has been chosen because, despite being used on a growing number of patients and with multi-Gy doses, treatment planning is still based on poorly characterized dose mapping systems such as look-up tables or low-resolution computational phantoms for MC codes. Also, experimental measurements on phantom slabs irradiated by an X-ray source were carried out preliminarily to validate the simulated radiation source. As shown in the results, the UM computational phantoms (built through the Simpleware SCAN-IP™ tool) can reduce the bias induced by the regularity of the classical cubic voxel geometry, give a more accurate description of volumes and complex surfaces and, thanks to an optimized discretization of the volumes, are also able to reduce the computational work. For the same UM models, various comparisons with different voxel sizes have been produced to investigate the dose distributions and evaluate the voxelization effects. The comparison shows a convergence pattern of the voxel model to the UM one. It has been possible to verify that the most significant bias occurs where the dose gradient is higher, along the beam borders and at tissue interfaces. The simulations showed that the UM can be really effective and reliable to compute the dose distributions within computational anthropomorphic phantoms obtained from the patient’s CT scans.

AB - The use of Monte Carlo methods for the set up of Treatment Planning Systems (TPS) in radiotherapy applications is a current standard. The most advanced modeling techniques aim at directly link the output of CT scans to a patient-specific model build up instead of standard phantoms and look-up tables. This link represents a critical step since even the most accurate segmentation and organ volume definition must be translated into a suitable input for the Monte Carlo code. During that step, the segmented volume is usually mapped on a regular geometrical lattice (voxels). A more sophisticated option appears to be a volume description based on an Unstructured Mesh (UM) geometry typical of current finite element codes. In this paper, we compared the two approaches analyzing the bias induced by the different choices. Starting from anonymous patient DICOM files coming from a CT scan, suitable segmentation and volume definition have been carried out and voxel and UM-based equivalent models for the Monte Carlo code MCNP6 have been built. Some computational phantoms, covering some significant portion of the human body (head and lower limb), have been used as a benchmark of the dose distribution obtained from X-ray sources commonly used in orthovoltage applications. The specific clinical application has been chosen because, despite being used on a growing number of patients and with multi-Gy doses, treatment planning is still based on poorly characterized dose mapping systems such as look-up tables or low-resolution computational phantoms for MC codes. Also, experimental measurements on phantom slabs irradiated by an X-ray source were carried out preliminarily to validate the simulated radiation source. As shown in the results, the UM computational phantoms (built through the Simpleware SCAN-IP™ tool) can reduce the bias induced by the regularity of the classical cubic voxel geometry, give a more accurate description of volumes and complex surfaces and, thanks to an optimized discretization of the volumes, are also able to reduce the computational work. For the same UM models, various comparisons with different voxel sizes have been produced to investigate the dose distributions and evaluate the voxelization effects. The comparison shows a convergence pattern of the voxel model to the UM one. It has been possible to verify that the most significant bias occurs where the dose gradient is higher, along the beam borders and at tissue interfaces. The simulations showed that the UM can be really effective and reliable to compute the dose distributions within computational anthropomorphic phantoms obtained from the patient’s CT scans.

KW - MCNP

KW - Monte Carlo

KW - orthovoltage

KW - radiotherapy

KW - unstructured mesh

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

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

U2 - 10.1080/10420150.2019.1596104

DO - 10.1080/10420150.2019.1596104

M3 - Article

JO - Radiation Effects and Defects in Solids

JF - Radiation Effects and Defects in Solids

SN - 1042-0150

ER -