Fast iterative algorithms for three-dimensional inverse treatment planning

L. Xing, Russell J Hamilton, D. Spelbring, C. A. Pelizzari, G. T Y Chen, A. L. Boyer

Research output: Contribution to journalArticle

59 Citations (Scopus)

Abstract

Three types of iterative algorithms, algebraic inverse treatment planning (AITP), simultaneous iterative inverse treatment planning (SIITP), and iterative least-square inverse treatment planning (ILSITP), differentiated according to their updating sequences, were generalized to three dimension with true beam geometry and dose model. A rapid ray-tracing approach was developed to optimize the primary beam components. Instead of recalculating the dose matrix at each iteration, the dose distribution was generated by scaling up or down the dose matrix elements of the previous iteration. This significantly increased the calculation speed. The iterative algorithms started with an initial intensity profile for each beam, specified by a two-dimensional pixel beam map of M elements. The calculation volume was divided into N voxels, and the calculation was done by repeatedly comparing the calculated and desired doses and adjusting the values of the beam map elements to minimize an objective function. In AITP, the iteration is performed voxel by voxel. For each voxel, the dose discrepancy was evaluated and the contributing pencil beams were updated. In ILSITP and SIITP, the iteration proceeded pencil beam by pencil beam instead of voxel by voxel. In all cases, the iteration procedure was repeated until the best possible dose distribution was achieved. The algorithms were applied to two examples and the results showed that the iterative techniques were able to produce superior isodose distributions.

Original languageEnglish (US)
Pages (from-to)1845-1849
Number of pages5
JournalMedical Physics
Volume25
Issue number10
DOIs
StatePublished - Oct 1998
Externally publishedYes

Fingerprint

Least-Squares Analysis

Keywords

  • Intensity modulation
  • Inverse problem
  • Iterative method
  • Treatment planning

ASJC Scopus subject areas

  • Biophysics

Cite this

Xing, L., Hamilton, R. J., Spelbring, D., Pelizzari, C. A., Chen, G. T. Y., & Boyer, A. L. (1998). Fast iterative algorithms for three-dimensional inverse treatment planning. Medical Physics, 25(10), 1845-1849. https://doi.org/10.1118/1.598374

Fast iterative algorithms for three-dimensional inverse treatment planning. / Xing, L.; Hamilton, Russell J; Spelbring, D.; Pelizzari, C. A.; Chen, G. T Y; Boyer, A. L.

In: Medical Physics, Vol. 25, No. 10, 10.1998, p. 1845-1849.

Research output: Contribution to journalArticle

Xing, L, Hamilton, RJ, Spelbring, D, Pelizzari, CA, Chen, GTY & Boyer, AL 1998, 'Fast iterative algorithms for three-dimensional inverse treatment planning', Medical Physics, vol. 25, no. 10, pp. 1845-1849. https://doi.org/10.1118/1.598374
Xing, L. ; Hamilton, Russell J ; Spelbring, D. ; Pelizzari, C. A. ; Chen, G. T Y ; Boyer, A. L. / Fast iterative algorithms for three-dimensional inverse treatment planning. In: Medical Physics. 1998 ; Vol. 25, No. 10. pp. 1845-1849.
@article{2230741e4b5340a6be9595016704f440,
title = "Fast iterative algorithms for three-dimensional inverse treatment planning",
abstract = "Three types of iterative algorithms, algebraic inverse treatment planning (AITP), simultaneous iterative inverse treatment planning (SIITP), and iterative least-square inverse treatment planning (ILSITP), differentiated according to their updating sequences, were generalized to three dimension with true beam geometry and dose model. A rapid ray-tracing approach was developed to optimize the primary beam components. Instead of recalculating the dose matrix at each iteration, the dose distribution was generated by scaling up or down the dose matrix elements of the previous iteration. This significantly increased the calculation speed. The iterative algorithms started with an initial intensity profile for each beam, specified by a two-dimensional pixel beam map of M elements. The calculation volume was divided into N voxels, and the calculation was done by repeatedly comparing the calculated and desired doses and adjusting the values of the beam map elements to minimize an objective function. In AITP, the iteration is performed voxel by voxel. For each voxel, the dose discrepancy was evaluated and the contributing pencil beams were updated. In ILSITP and SIITP, the iteration proceeded pencil beam by pencil beam instead of voxel by voxel. In all cases, the iteration procedure was repeated until the best possible dose distribution was achieved. The algorithms were applied to two examples and the results showed that the iterative techniques were able to produce superior isodose distributions.",
keywords = "Intensity modulation, Inverse problem, Iterative method, Treatment planning",
author = "L. Xing and Hamilton, {Russell J} and D. Spelbring and Pelizzari, {C. A.} and Chen, {G. T Y} and Boyer, {A. L.}",
year = "1998",
month = "10",
doi = "10.1118/1.598374",
language = "English (US)",
volume = "25",
pages = "1845--1849",
journal = "Medical Physics",
issn = "0094-2405",
publisher = "AAPM - American Association of Physicists in Medicine",
number = "10",

}

TY - JOUR

T1 - Fast iterative algorithms for three-dimensional inverse treatment planning

AU - Xing, L.

AU - Hamilton, Russell J

AU - Spelbring, D.

AU - Pelizzari, C. A.

AU - Chen, G. T Y

AU - Boyer, A. L.

PY - 1998/10

Y1 - 1998/10

N2 - Three types of iterative algorithms, algebraic inverse treatment planning (AITP), simultaneous iterative inverse treatment planning (SIITP), and iterative least-square inverse treatment planning (ILSITP), differentiated according to their updating sequences, were generalized to three dimension with true beam geometry and dose model. A rapid ray-tracing approach was developed to optimize the primary beam components. Instead of recalculating the dose matrix at each iteration, the dose distribution was generated by scaling up or down the dose matrix elements of the previous iteration. This significantly increased the calculation speed. The iterative algorithms started with an initial intensity profile for each beam, specified by a two-dimensional pixel beam map of M elements. The calculation volume was divided into N voxels, and the calculation was done by repeatedly comparing the calculated and desired doses and adjusting the values of the beam map elements to minimize an objective function. In AITP, the iteration is performed voxel by voxel. For each voxel, the dose discrepancy was evaluated and the contributing pencil beams were updated. In ILSITP and SIITP, the iteration proceeded pencil beam by pencil beam instead of voxel by voxel. In all cases, the iteration procedure was repeated until the best possible dose distribution was achieved. The algorithms were applied to two examples and the results showed that the iterative techniques were able to produce superior isodose distributions.

AB - Three types of iterative algorithms, algebraic inverse treatment planning (AITP), simultaneous iterative inverse treatment planning (SIITP), and iterative least-square inverse treatment planning (ILSITP), differentiated according to their updating sequences, were generalized to three dimension with true beam geometry and dose model. A rapid ray-tracing approach was developed to optimize the primary beam components. Instead of recalculating the dose matrix at each iteration, the dose distribution was generated by scaling up or down the dose matrix elements of the previous iteration. This significantly increased the calculation speed. The iterative algorithms started with an initial intensity profile for each beam, specified by a two-dimensional pixel beam map of M elements. The calculation volume was divided into N voxels, and the calculation was done by repeatedly comparing the calculated and desired doses and adjusting the values of the beam map elements to minimize an objective function. In AITP, the iteration is performed voxel by voxel. For each voxel, the dose discrepancy was evaluated and the contributing pencil beams were updated. In ILSITP and SIITP, the iteration proceeded pencil beam by pencil beam instead of voxel by voxel. In all cases, the iteration procedure was repeated until the best possible dose distribution was achieved. The algorithms were applied to two examples and the results showed that the iterative techniques were able to produce superior isodose distributions.

KW - Intensity modulation

KW - Inverse problem

KW - Iterative method

KW - Treatment planning

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

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

U2 - 10.1118/1.598374

DO - 10.1118/1.598374

M3 - Article

C2 - 9800690

AN - SCOPUS:0031732963

VL - 25

SP - 1845

EP - 1849

JO - Medical Physics

JF - Medical Physics

SN - 0094-2405

IS - 10

ER -