Absence of multiple local minima effects in intensity modulated optimization with dose-volume constraints

Jorge Llacer, Jospeh O. Deasy, Thomas R. Bortfeld, Timothy D. Solberg, Claus Promberger

Research output: Contribution to journalArticle

51 Citations (Scopus)

Abstract

This paper reports on the analysis of intensity modulated radiation treatment optimization problems in the presence of non-convex feasible parameter spaces caused by the specification of dose-volume constraints for the organs-at-risk (OARs). The main aim was to determine whether the presence of those non-convex spaces affects the optimization of clinical cases in any significant way. This was done in two phases: (1) Using a carefully designed two-dimensional mathematical phantom that exhibits two ontrollable minima and with randomly initialized beamlet weights, we developed a methodology for exploring the nature of the convergence characteristics of quadratic cost function optimizations (deterministic or stochastic). The methodology is based on observing the statistical behaviour of the residual cost at the end of optimizations in which the stopping criterion is progressively more demanding and carrying out those optimizations to very small error changes per iteration. (2) Seven clinical cases were then analysed with dose-volume constraints that are stronger than originally used in the clinic. The clinical cases are two prostate cases differently posed, a meningioma case, two head-and-neck cases, a spleen case and a spine case. Of the 14 different sets of optimizations (with and without the specification of maximum doses allowed for the OARs), 12 fail to show any effect due to the existence of non-convex feasible spaces. The remaining two sets of optimizations show evidence of multiple minima in the solutions, but those minima are very close to each other in cost and the resulting treatment plans are practically identical, as measured by the quality of the dose-volume histograms (DVHs). We discuss the differences between fluence maps resulting from those similar treatment plans. We provide a possible reason for the observed results and conclude that, although the study is necessarily limited, the annealing characteristics of a simulated annealing method may not be justified in clinical optimization in the presence of dose-volume constraints. The results of optimizations by the Newton gradient (NG) method with a quadratic cost function are reported in detail. An adaptive simulated annealing method, optimizing the same function, and the dynamically penalized likelihood method, optimizing a log likelihood function, have also been used in the study. The results of the latter two methods have only been discussed briefly, as they yielded the same conclusions as the NG method.

Original languageEnglish (US)
Pages (from-to)183-210
Number of pages28
JournalPhysics in Medicine and Biology
Volume48
Issue number2
DOIs
StatePublished - Jan 21 2003

Fingerprint

dosage
optimization
Organs at Risk
Costs and Cost Analysis
costs
Gradient methods
simulated annealing
Simulated annealing
Likelihood Functions
organs
Cost functions
newton
specifications
Meningioma
Health Care Costs
methodology
Specifications
Prostate
gradients
Spine

ASJC Scopus subject areas

  • Biomedical Engineering
  • Physics and Astronomy (miscellaneous)
  • Radiology Nuclear Medicine and imaging
  • Radiological and Ultrasound Technology

Cite this

Absence of multiple local minima effects in intensity modulated optimization with dose-volume constraints. / Llacer, Jorge; Deasy, Jospeh O.; Bortfeld, Thomas R.; Solberg, Timothy D.; Promberger, Claus.

In: Physics in Medicine and Biology, Vol. 48, No. 2, 21.01.2003, p. 183-210.

Research output: Contribution to journalArticle

Llacer, Jorge ; Deasy, Jospeh O. ; Bortfeld, Thomas R. ; Solberg, Timothy D. ; Promberger, Claus. / Absence of multiple local minima effects in intensity modulated optimization with dose-volume constraints. In: Physics in Medicine and Biology. 2003 ; Vol. 48, No. 2. pp. 183-210.
@article{ed3abf529748483091e24ee72003c31e,
title = "Absence of multiple local minima effects in intensity modulated optimization with dose-volume constraints",
abstract = "This paper reports on the analysis of intensity modulated radiation treatment optimization problems in the presence of non-convex feasible parameter spaces caused by the specification of dose-volume constraints for the organs-at-risk (OARs). The main aim was to determine whether the presence of those non-convex spaces affects the optimization of clinical cases in any significant way. This was done in two phases: (1) Using a carefully designed two-dimensional mathematical phantom that exhibits two ontrollable minima and with randomly initialized beamlet weights, we developed a methodology for exploring the nature of the convergence characteristics of quadratic cost function optimizations (deterministic or stochastic). The methodology is based on observing the statistical behaviour of the residual cost at the end of optimizations in which the stopping criterion is progressively more demanding and carrying out those optimizations to very small error changes per iteration. (2) Seven clinical cases were then analysed with dose-volume constraints that are stronger than originally used in the clinic. The clinical cases are two prostate cases differently posed, a meningioma case, two head-and-neck cases, a spleen case and a spine case. Of the 14 different sets of optimizations (with and without the specification of maximum doses allowed for the OARs), 12 fail to show any effect due to the existence of non-convex feasible spaces. The remaining two sets of optimizations show evidence of multiple minima in the solutions, but those minima are very close to each other in cost and the resulting treatment plans are practically identical, as measured by the quality of the dose-volume histograms (DVHs). We discuss the differences between fluence maps resulting from those similar treatment plans. We provide a possible reason for the observed results and conclude that, although the study is necessarily limited, the annealing characteristics of a simulated annealing method may not be justified in clinical optimization in the presence of dose-volume constraints. The results of optimizations by the Newton gradient (NG) method with a quadratic cost function are reported in detail. An adaptive simulated annealing method, optimizing the same function, and the dynamically penalized likelihood method, optimizing a log likelihood function, have also been used in the study. The results of the latter two methods have only been discussed briefly, as they yielded the same conclusions as the NG method.",
author = "Jorge Llacer and Deasy, {Jospeh O.} and Bortfeld, {Thomas R.} and Solberg, {Timothy D.} and Claus Promberger",
year = "2003",
month = "1",
day = "21",
doi = "10.1088/0031-9155/48/2/304",
language = "English (US)",
volume = "48",
pages = "183--210",
journal = "Physics in Medicine and Biology",
issn = "0031-9155",
publisher = "IOP Publishing Ltd.",
number = "2",

}

TY - JOUR

T1 - Absence of multiple local minima effects in intensity modulated optimization with dose-volume constraints

AU - Llacer, Jorge

AU - Deasy, Jospeh O.

AU - Bortfeld, Thomas R.

AU - Solberg, Timothy D.

AU - Promberger, Claus

PY - 2003/1/21

Y1 - 2003/1/21

N2 - This paper reports on the analysis of intensity modulated radiation treatment optimization problems in the presence of non-convex feasible parameter spaces caused by the specification of dose-volume constraints for the organs-at-risk (OARs). The main aim was to determine whether the presence of those non-convex spaces affects the optimization of clinical cases in any significant way. This was done in two phases: (1) Using a carefully designed two-dimensional mathematical phantom that exhibits two ontrollable minima and with randomly initialized beamlet weights, we developed a methodology for exploring the nature of the convergence characteristics of quadratic cost function optimizations (deterministic or stochastic). The methodology is based on observing the statistical behaviour of the residual cost at the end of optimizations in which the stopping criterion is progressively more demanding and carrying out those optimizations to very small error changes per iteration. (2) Seven clinical cases were then analysed with dose-volume constraints that are stronger than originally used in the clinic. The clinical cases are two prostate cases differently posed, a meningioma case, two head-and-neck cases, a spleen case and a spine case. Of the 14 different sets of optimizations (with and without the specification of maximum doses allowed for the OARs), 12 fail to show any effect due to the existence of non-convex feasible spaces. The remaining two sets of optimizations show evidence of multiple minima in the solutions, but those minima are very close to each other in cost and the resulting treatment plans are practically identical, as measured by the quality of the dose-volume histograms (DVHs). We discuss the differences between fluence maps resulting from those similar treatment plans. We provide a possible reason for the observed results and conclude that, although the study is necessarily limited, the annealing characteristics of a simulated annealing method may not be justified in clinical optimization in the presence of dose-volume constraints. The results of optimizations by the Newton gradient (NG) method with a quadratic cost function are reported in detail. An adaptive simulated annealing method, optimizing the same function, and the dynamically penalized likelihood method, optimizing a log likelihood function, have also been used in the study. The results of the latter two methods have only been discussed briefly, as they yielded the same conclusions as the NG method.

AB - This paper reports on the analysis of intensity modulated radiation treatment optimization problems in the presence of non-convex feasible parameter spaces caused by the specification of dose-volume constraints for the organs-at-risk (OARs). The main aim was to determine whether the presence of those non-convex spaces affects the optimization of clinical cases in any significant way. This was done in two phases: (1) Using a carefully designed two-dimensional mathematical phantom that exhibits two ontrollable minima and with randomly initialized beamlet weights, we developed a methodology for exploring the nature of the convergence characteristics of quadratic cost function optimizations (deterministic or stochastic). The methodology is based on observing the statistical behaviour of the residual cost at the end of optimizations in which the stopping criterion is progressively more demanding and carrying out those optimizations to very small error changes per iteration. (2) Seven clinical cases were then analysed with dose-volume constraints that are stronger than originally used in the clinic. The clinical cases are two prostate cases differently posed, a meningioma case, two head-and-neck cases, a spleen case and a spine case. Of the 14 different sets of optimizations (with and without the specification of maximum doses allowed for the OARs), 12 fail to show any effect due to the existence of non-convex feasible spaces. The remaining two sets of optimizations show evidence of multiple minima in the solutions, but those minima are very close to each other in cost and the resulting treatment plans are practically identical, as measured by the quality of the dose-volume histograms (DVHs). We discuss the differences between fluence maps resulting from those similar treatment plans. We provide a possible reason for the observed results and conclude that, although the study is necessarily limited, the annealing characteristics of a simulated annealing method may not be justified in clinical optimization in the presence of dose-volume constraints. The results of optimizations by the Newton gradient (NG) method with a quadratic cost function are reported in detail. An adaptive simulated annealing method, optimizing the same function, and the dynamically penalized likelihood method, optimizing a log likelihood function, have also been used in the study. The results of the latter two methods have only been discussed briefly, as they yielded the same conclusions as the NG method.

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

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

U2 - 10.1088/0031-9155/48/2/304

DO - 10.1088/0031-9155/48/2/304

M3 - Article

C2 - 12587904

AN - SCOPUS:0037458222

VL - 48

SP - 183

EP - 210

JO - Physics in Medicine and Biology

JF - Physics in Medicine and Biology

SN - 0031-9155

IS - 2

ER -