One-dimensional scatter-subtraction method for brachytherapy dose calculation near bounded heterogeneities

Jeffrey F. Williamson, Zuofena Li, John Wong

Research output: Contribution to journalArticle

Abstract

Neglect of tissue and applicator heterogeneities in brachytherapy treatment planning is due in part to lack of accurate, general, and fast three-dimensional (3D) dose-computation algorithms. A novel dose-calculation algorithm that accounts for the lateral dimensions and location of the heterogeneity as well as its thickness has been developed. For simple 2D, water-equivalent density heterogeneities, the algorithm is shown to be applicable to a wide range of photon energies and is 500–1000 times more efficient than Monte Carlo photon-transport calculations. The model is based upon reducing the bounded 2D heterogeneity problem to two 1D problems by dividing the scattering volume into two regions: A cone-shaped region that subtends the heterogeneity with its apex at the source and the complementary cone that contributes scatter dose by diffusion around the heterogeneity. The input data consist of precalculated scatter-to-primary ratios (SPRs) for collimated isotropic point sources. The central-axis “mini-beam” problem for a slab heterogeneity is solved by a simple 1D scatter integration model that accounts for both the thickness of the heterogeneity and its location relative to the point of interest. The scatter contribution arising outside the mini-beam is modeled as the difference in SPRs corrected for transmission through the heterogeneity. The algorithm agrees, on average, with sample Monte Carlo calculations within 1% to 7% for 125I, 192Ir, and 100 keV point sources along the axes of water-equivalent cylindrical heterogeneities (p=0-12.6 g/cm3, 3.6, and 24 mm diameters, and dose-perturbation factors of 0.44-1.33 relative to the homogeneous case). The potential of generalizing the scatter-subtraction approach to encompass 3D heterogeneities, those consisting of high-atomic number media, and those of irregular shape, is discussed.

Original languageEnglish (US)
Pages (from-to)233-244
Number of pages12
JournalMedical Physics
Volume20
Issue number1
DOIs
StatePublished - 1993
Externally publishedYes

Fingerprint

Brachytherapy
Photons
Water
Therapeutics

Keywords

  • brachytherapy dosimetry
  • heterogeneity corrections
  • Monte Carlo simulation

ASJC Scopus subject areas

  • Biophysics
  • Radiology Nuclear Medicine and imaging

Cite this

One-dimensional scatter-subtraction method for brachytherapy dose calculation near bounded heterogeneities. / Williamson, Jeffrey F.; Li, Zuofena; Wong, John.

In: Medical Physics, Vol. 20, No. 1, 1993, p. 233-244.

Research output: Contribution to journalArticle

@article{4e9efd4ebb194d769965d9233a697926,
title = "One-dimensional scatter-subtraction method for brachytherapy dose calculation near bounded heterogeneities",
abstract = "Neglect of tissue and applicator heterogeneities in brachytherapy treatment planning is due in part to lack of accurate, general, and fast three-dimensional (3D) dose-computation algorithms. A novel dose-calculation algorithm that accounts for the lateral dimensions and location of the heterogeneity as well as its thickness has been developed. For simple 2D, water-equivalent density heterogeneities, the algorithm is shown to be applicable to a wide range of photon energies and is 500–1000 times more efficient than Monte Carlo photon-transport calculations. The model is based upon reducing the bounded 2D heterogeneity problem to two 1D problems by dividing the scattering volume into two regions: A cone-shaped region that subtends the heterogeneity with its apex at the source and the complementary cone that contributes scatter dose by diffusion around the heterogeneity. The input data consist of precalculated scatter-to-primary ratios (SPRs) for collimated isotropic point sources. The central-axis “mini-beam” problem for a slab heterogeneity is solved by a simple 1D scatter integration model that accounts for both the thickness of the heterogeneity and its location relative to the point of interest. The scatter contribution arising outside the mini-beam is modeled as the difference in SPRs corrected for transmission through the heterogeneity. The algorithm agrees, on average, with sample Monte Carlo calculations within 1{\%} to 7{\%} for 125I, 192Ir, and 100 keV point sources along the axes of water-equivalent cylindrical heterogeneities (p=0-12.6 g/cm3, 3.6, and 24 mm diameters, and dose-perturbation factors of 0.44-1.33 relative to the homogeneous case). The potential of generalizing the scatter-subtraction approach to encompass 3D heterogeneities, those consisting of high-atomic number media, and those of irregular shape, is discussed.",
keywords = "brachytherapy dosimetry, heterogeneity corrections, Monte Carlo simulation",
author = "Williamson, {Jeffrey F.} and Zuofena Li and John Wong",
year = "1993",
doi = "10.1118/1.597090",
language = "English (US)",
volume = "20",
pages = "233--244",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "1",

}

TY - JOUR

T1 - One-dimensional scatter-subtraction method for brachytherapy dose calculation near bounded heterogeneities

AU - Williamson, Jeffrey F.

AU - Li, Zuofena

AU - Wong, John

PY - 1993

Y1 - 1993

N2 - Neglect of tissue and applicator heterogeneities in brachytherapy treatment planning is due in part to lack of accurate, general, and fast three-dimensional (3D) dose-computation algorithms. A novel dose-calculation algorithm that accounts for the lateral dimensions and location of the heterogeneity as well as its thickness has been developed. For simple 2D, water-equivalent density heterogeneities, the algorithm is shown to be applicable to a wide range of photon energies and is 500–1000 times more efficient than Monte Carlo photon-transport calculations. The model is based upon reducing the bounded 2D heterogeneity problem to two 1D problems by dividing the scattering volume into two regions: A cone-shaped region that subtends the heterogeneity with its apex at the source and the complementary cone that contributes scatter dose by diffusion around the heterogeneity. The input data consist of precalculated scatter-to-primary ratios (SPRs) for collimated isotropic point sources. The central-axis “mini-beam” problem for a slab heterogeneity is solved by a simple 1D scatter integration model that accounts for both the thickness of the heterogeneity and its location relative to the point of interest. The scatter contribution arising outside the mini-beam is modeled as the difference in SPRs corrected for transmission through the heterogeneity. The algorithm agrees, on average, with sample Monte Carlo calculations within 1% to 7% for 125I, 192Ir, and 100 keV point sources along the axes of water-equivalent cylindrical heterogeneities (p=0-12.6 g/cm3, 3.6, and 24 mm diameters, and dose-perturbation factors of 0.44-1.33 relative to the homogeneous case). The potential of generalizing the scatter-subtraction approach to encompass 3D heterogeneities, those consisting of high-atomic number media, and those of irregular shape, is discussed.

AB - Neglect of tissue and applicator heterogeneities in brachytherapy treatment planning is due in part to lack of accurate, general, and fast three-dimensional (3D) dose-computation algorithms. A novel dose-calculation algorithm that accounts for the lateral dimensions and location of the heterogeneity as well as its thickness has been developed. For simple 2D, water-equivalent density heterogeneities, the algorithm is shown to be applicable to a wide range of photon energies and is 500–1000 times more efficient than Monte Carlo photon-transport calculations. The model is based upon reducing the bounded 2D heterogeneity problem to two 1D problems by dividing the scattering volume into two regions: A cone-shaped region that subtends the heterogeneity with its apex at the source and the complementary cone that contributes scatter dose by diffusion around the heterogeneity. The input data consist of precalculated scatter-to-primary ratios (SPRs) for collimated isotropic point sources. The central-axis “mini-beam” problem for a slab heterogeneity is solved by a simple 1D scatter integration model that accounts for both the thickness of the heterogeneity and its location relative to the point of interest. The scatter contribution arising outside the mini-beam is modeled as the difference in SPRs corrected for transmission through the heterogeneity. The algorithm agrees, on average, with sample Monte Carlo calculations within 1% to 7% for 125I, 192Ir, and 100 keV point sources along the axes of water-equivalent cylindrical heterogeneities (p=0-12.6 g/cm3, 3.6, and 24 mm diameters, and dose-perturbation factors of 0.44-1.33 relative to the homogeneous case). The potential of generalizing the scatter-subtraction approach to encompass 3D heterogeneities, those consisting of high-atomic number media, and those of irregular shape, is discussed.

KW - brachytherapy dosimetry

KW - heterogeneity corrections

KW - Monte Carlo simulation

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

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

U2 - 10.1118/1.597090

DO - 10.1118/1.597090

M3 - Article

C2 - 8455505

AN - SCOPUS:0027462867

VL - 20

SP - 233

EP - 244

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 1

ER -