×

zbMATH — the first resource for mathematics

Using eigenstructure of the Hessian to reduce the dimension of the intensity modulated radiation therapy optimization problem. (English) Zbl 1104.92030
Summary: Optimization is of vital importance when performing intensity modulated radiation therapy to treat cancer tumors. The optimization problem is typically large-scale with a nonlinear objective function and bounds on the variables, and we solve it using a quasi-Newton sequential quadratic programming method. This study investigates the effect on the optimal solution, and hence treatment outcome, when solving an approximate optimization problem of lower dimension. Through a spectral decomposition, eigenvectors and eigenvalues of an approximation to the Hessian are computed. An approximate optimization problem of reduced dimension is formulated by introducing eigenvector weights as optimization parameters, where only eigenvectors corresponding to large eigenvalues are included.
The approach is evaluated on a clinical prostate case. Compared to bixel weight optimization, eigenvector weight optimization with few parameters results in faster initial decline in the objective function, but with inferior final solution. Another approach, which combines eigenvector weights and bixel weights as variables, gives lower final objective values than what bixel weight optimization does. However, this advantage comes at the expense of the pre-computational time for the spectral decomposition.

MSC:
92C50 Medical applications (general)
90C90 Applications of mathematical programming
90C20 Quadratic programming
Keywords:
IMRT
Software:
NPSOL
PDF BibTeX XML Cite
Full Text: DOI
References:
[1] Ahnesjö, A. and M.M. Aspradakis. (1999). ”Dose Calculations for External Photon Beams in Radiotherapy.” Physics in Medicine and Biology, 44(11), R99–R155. · doi:10.1088/0031-9155/44/11/201
[2] Alber, M., G. Meedt, F. Nüsslin, and R. Reemtsen. (2002). ”On the Degeneracy of the IMRT Optimization Problem.” Medical Physics, 29(11), 2584–2589. · doi:10.1118/1.1500402
[3] Bednarz, G., D. Michalski, C. Houser, M.S. Huq, Y. Xiao, P.R. Anne, and J.M. Galvin. (2002). ”The use of Mixed-Integer Programming for Inverse Treatment Planning with Pre-Defined Field Segments.” Physics in Medicine and Biology, 47(13), 2235–2245. · doi:10.1088/0031-9155/47/13/304
[4] Carlsson, F. and A. Forsgren. (2006). ”Iterative Regularization in Intensity-Modulated Radiation Therapy Optimization.” Medical Physics, 33(1), 225–234. · doi:10.1118/1.2148918
[5] Deasy, J.O. (1997). ”Multiple Local Minima in Radiotherapy Optimization Problems with Dose-Volume Constraints.” Medical Physics, 24(7), 1157–1161. · doi:10.1118/1.598017
[6] Gill, P.E., W. Murray, M.A. Saunders, and M.H. Wright. (1986). ”User’s Guide for NPSOL (Version 4.0): a Fortran Package for Nonlinear Programming.” Stanford University SOL 86-2.
[7] Golub, G.H. and C.F. Van Loan. (1996). Matrix Computations. 3rd edn. The Johns Hopkins University Press, Baltimore, Maryland, 1996. ISBN 0-8018-5414-8. · Zbl 0865.65009
[8] Küfer, K-.H., A. Scherrer, M. Monz, F. Alonso, H. Trinkaus, T. Bortfeld, and C. Thieke. (2003). ”Intensity-Modulated Radiotherapy–A Large Scale Multi-Criteria Programming Problem.” OR Spectrum, 25, 223–249. · Zbl 1018.92015 · doi:10.1007/s00291-003-0125-7
[9] Löf, J. (2000). Development of a General Framework for Optimization of Radiation Therapy. PhD thesis, Stockholm University.
[10] Markman, J., D. Low, A. Beavis, and J. Deasy. (2002) ”Beyond Bixels: Generalizing the Optimization Parameters for Intensity Modulated Radiation Therapy.” Medical Physics, 29(10), 2298–2304. · doi:10.1118/1.1508799
[11] Rowbottom, C.G. and S. Webb. (2002). ”Configuration Space Analysis of Common Cost Functions in Radiotherapy Beam-Weight Optimization Algorithms.” Physics in Medicine and Biology, 47(1), 65–77. · doi:10.1088/0031-9155/47/1/305
[12] Wu, Q. and R. Mohan. (2002) ”Multiple Local Minima in IMRT Optimization Based on Dose-Volume Criteria.” Medical Physics, 29(7), 1514–1527. · doi:10.1118/1.1485059
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. It attempts to reflect the references listed in the original paper as accurately as possible without claiming the completeness or perfect precision of the matching.