A new Monte Carlo-based treatment plan optimization approach for intensity modulated radiation therapy

Yongbao Li, Zhen Tian, Feng Shi, Ting Song, Zhaoxia Wu, Yaqiang Liu, Steve Jiang, Xun Jia

Research output: Contribution to journalArticle

7 Citations (Scopus)

Abstract

Intensity-modulated radiation treatment (IMRT) plan optimization needs beamlet dose distributions. Pencil-beam or superposition/convolution type algorithms are typically used because of their high computational speed. However, inaccurate beamlet dose distributions may mislead the optimization process and hinder the resulting plan quality. To solve this problem, the Monte Carlo (MC) simulation method has been used to compute all beamlet doses prior to the optimization step. The conventional approach samples the same number of particles from each beamlet. Yet this is not the optimal use of MC in this problem. In fact, there are beamlets that have very small intensities after solving the plan optimization problem. For those beamlets, it may be possible to use fewer particles in dose calculations to increase efficiency. Based on this idea, we have developed a new MC-based IMRT plan optimization framework that iteratively performs MC dose calculation and plan optimization. At each dose calculation step, the particle numbers for beamlets were adjusted based on the beamlet intensities obtained through solving the plan optimization problem in the last iteration step. We modified a GPU-based MC dose engine to allow simultaneous computations of a large number of beamlet doses. To test the accuracy of our modified dose engine, we compared the dose from a broad beam and the summed beamlet doses in this beam in an inhomogeneous phantom. Agreement within 1% for the maximum difference and 0.55% for the average difference was observed. We then validated the proposed MC-based optimization schemes in one lung IMRT case. It was found that the conventional scheme required 106 particles from each beamlet to achieve an optimization result that was 3% difference in fluence map and 1% difference in dose from the ground truth. In contrast, the proposed scheme achieved the same level of accuracy with on average 1.2×105 particles per beamlet. Correspondingly, the computation time including both MC dose calculations and plan optimizations was reduced by a factor of 4.4, from 494 to 113 s, using only one GPU card.

Original languageEnglish (US)
Pages (from-to)2903-2919
Number of pages17
JournalPhysics in Medicine and Biology
Volume60
Issue number7
DOIs
StatePublished - Apr 7 2015

Fingerprint

Radiotherapy
Radiation
Monte Carlo Method
Therapeutics
Lung

Keywords

  • intensity modulated radiation therapy
  • Monte Carlo
  • optimization

ASJC Scopus subject areas

  • Radiology Nuclear Medicine and imaging
  • Radiological and Ultrasound Technology

Cite this

A new Monte Carlo-based treatment plan optimization approach for intensity modulated radiation therapy. / Li, Yongbao; Tian, Zhen; Shi, Feng; Song, Ting; Wu, Zhaoxia; Liu, Yaqiang; Jiang, Steve; Jia, Xun.

In: Physics in Medicine and Biology, Vol. 60, No. 7, 07.04.2015, p. 2903-2919.

Research output: Contribution to journalArticle

Li, Yongbao ; Tian, Zhen ; Shi, Feng ; Song, Ting ; Wu, Zhaoxia ; Liu, Yaqiang ; Jiang, Steve ; Jia, Xun. / A new Monte Carlo-based treatment plan optimization approach for intensity modulated radiation therapy. In: Physics in Medicine and Biology. 2015 ; Vol. 60, No. 7. pp. 2903-2919.
@article{1772ac44522645e68c004f6228836e00,
title = "A new Monte Carlo-based treatment plan optimization approach for intensity modulated radiation therapy",
abstract = "Intensity-modulated radiation treatment (IMRT) plan optimization needs beamlet dose distributions. Pencil-beam or superposition/convolution type algorithms are typically used because of their high computational speed. However, inaccurate beamlet dose distributions may mislead the optimization process and hinder the resulting plan quality. To solve this problem, the Monte Carlo (MC) simulation method has been used to compute all beamlet doses prior to the optimization step. The conventional approach samples the same number of particles from each beamlet. Yet this is not the optimal use of MC in this problem. In fact, there are beamlets that have very small intensities after solving the plan optimization problem. For those beamlets, it may be possible to use fewer particles in dose calculations to increase efficiency. Based on this idea, we have developed a new MC-based IMRT plan optimization framework that iteratively performs MC dose calculation and plan optimization. At each dose calculation step, the particle numbers for beamlets were adjusted based on the beamlet intensities obtained through solving the plan optimization problem in the last iteration step. We modified a GPU-based MC dose engine to allow simultaneous computations of a large number of beamlet doses. To test the accuracy of our modified dose engine, we compared the dose from a broad beam and the summed beamlet doses in this beam in an inhomogeneous phantom. Agreement within 1{\%} for the maximum difference and 0.55{\%} for the average difference was observed. We then validated the proposed MC-based optimization schemes in one lung IMRT case. It was found that the conventional scheme required 106 particles from each beamlet to achieve an optimization result that was 3{\%} difference in fluence map and 1{\%} difference in dose from the ground truth. In contrast, the proposed scheme achieved the same level of accuracy with on average 1.2×105 particles per beamlet. Correspondingly, the computation time including both MC dose calculations and plan optimizations was reduced by a factor of 4.4, from 494 to 113 s, using only one GPU card.",
keywords = "intensity modulated radiation therapy, Monte Carlo, optimization",
author = "Yongbao Li and Zhen Tian and Feng Shi and Ting Song and Zhaoxia Wu and Yaqiang Liu and Steve Jiang and Xun Jia",
year = "2015",
month = "4",
day = "7",
doi = "10.1088/0031-9155/60/7/2903",
language = "English (US)",
volume = "60",
pages = "2903--2919",
journal = "Physics in Medicine and Biology",
issn = "0031-9155",
publisher = "IOP Publishing Ltd.",
number = "7",

}

TY - JOUR

T1 - A new Monte Carlo-based treatment plan optimization approach for intensity modulated radiation therapy

AU - Li, Yongbao

AU - Tian, Zhen

AU - Shi, Feng

AU - Song, Ting

AU - Wu, Zhaoxia

AU - Liu, Yaqiang

AU - Jiang, Steve

AU - Jia, Xun

PY - 2015/4/7

Y1 - 2015/4/7

N2 - Intensity-modulated radiation treatment (IMRT) plan optimization needs beamlet dose distributions. Pencil-beam or superposition/convolution type algorithms are typically used because of their high computational speed. However, inaccurate beamlet dose distributions may mislead the optimization process and hinder the resulting plan quality. To solve this problem, the Monte Carlo (MC) simulation method has been used to compute all beamlet doses prior to the optimization step. The conventional approach samples the same number of particles from each beamlet. Yet this is not the optimal use of MC in this problem. In fact, there are beamlets that have very small intensities after solving the plan optimization problem. For those beamlets, it may be possible to use fewer particles in dose calculations to increase efficiency. Based on this idea, we have developed a new MC-based IMRT plan optimization framework that iteratively performs MC dose calculation and plan optimization. At each dose calculation step, the particle numbers for beamlets were adjusted based on the beamlet intensities obtained through solving the plan optimization problem in the last iteration step. We modified a GPU-based MC dose engine to allow simultaneous computations of a large number of beamlet doses. To test the accuracy of our modified dose engine, we compared the dose from a broad beam and the summed beamlet doses in this beam in an inhomogeneous phantom. Agreement within 1% for the maximum difference and 0.55% for the average difference was observed. We then validated the proposed MC-based optimization schemes in one lung IMRT case. It was found that the conventional scheme required 106 particles from each beamlet to achieve an optimization result that was 3% difference in fluence map and 1% difference in dose from the ground truth. In contrast, the proposed scheme achieved the same level of accuracy with on average 1.2×105 particles per beamlet. Correspondingly, the computation time including both MC dose calculations and plan optimizations was reduced by a factor of 4.4, from 494 to 113 s, using only one GPU card.

AB - Intensity-modulated radiation treatment (IMRT) plan optimization needs beamlet dose distributions. Pencil-beam or superposition/convolution type algorithms are typically used because of their high computational speed. However, inaccurate beamlet dose distributions may mislead the optimization process and hinder the resulting plan quality. To solve this problem, the Monte Carlo (MC) simulation method has been used to compute all beamlet doses prior to the optimization step. The conventional approach samples the same number of particles from each beamlet. Yet this is not the optimal use of MC in this problem. In fact, there are beamlets that have very small intensities after solving the plan optimization problem. For those beamlets, it may be possible to use fewer particles in dose calculations to increase efficiency. Based on this idea, we have developed a new MC-based IMRT plan optimization framework that iteratively performs MC dose calculation and plan optimization. At each dose calculation step, the particle numbers for beamlets were adjusted based on the beamlet intensities obtained through solving the plan optimization problem in the last iteration step. We modified a GPU-based MC dose engine to allow simultaneous computations of a large number of beamlet doses. To test the accuracy of our modified dose engine, we compared the dose from a broad beam and the summed beamlet doses in this beam in an inhomogeneous phantom. Agreement within 1% for the maximum difference and 0.55% for the average difference was observed. We then validated the proposed MC-based optimization schemes in one lung IMRT case. It was found that the conventional scheme required 106 particles from each beamlet to achieve an optimization result that was 3% difference in fluence map and 1% difference in dose from the ground truth. In contrast, the proposed scheme achieved the same level of accuracy with on average 1.2×105 particles per beamlet. Correspondingly, the computation time including both MC dose calculations and plan optimizations was reduced by a factor of 4.4, from 494 to 113 s, using only one GPU card.

KW - intensity modulated radiation therapy

KW - Monte Carlo

KW - optimization

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

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

U2 - 10.1088/0031-9155/60/7/2903

DO - 10.1088/0031-9155/60/7/2903

M3 - Article

VL - 60

SP - 2903

EP - 2919

JO - Physics in Medicine and Biology

JF - Physics in Medicine and Biology

SN - 0031-9155

IS - 7

ER -