Iterative regularization in intensity-modulated radiation therapy optimization

被引:16
作者
Carlsson, F
Forsgren, A
机构
[1] RaySearch Labs, SE-11134 Stockholm, Sweden
[2] Royal Inst Technol KTH, Dept Math, SE-10044 Stockholm, Sweden
关键词
intensity-modulated radiation therapy; quasi-Newton method; conjugate gradient method; regularization; iterative regularization;
D O I
10.1118/1.2148918
中图分类号
R8 [特种医学]; R445 [影像诊断学];
学科分类号
1002 ; 100207 ; 1009 ;
摘要
A common way to solve intensity-modulated radiation therapy (IMRT) optimization problems is to use a beamlet-based approach. The approach is usually employed in a three-step manner: first a beamlet-weight optimization problem is solved, then the fluence profiles are converted into stepand-shoot segments, and finally postoptimization of the segment weights is performed. A drawback of beamlet-based approaches is that beamlet-weight optimization problems are ill-conditioned and have to be regularized in order to produce smooth fluence profiles that are suitable for conversion. The purpose of this paper is twofold: first, to explain the suitability of solving beamlet-based IMRT problems by a BFGS quasi-Newton sequential quadratic programming method with diagonal initial Hessian estimate, and second, to empirically show that beamlet-weight optimization problems should be solved in relatively few iterations when using this optimization method. The explanation of the suitability is based on viewing the optimization method as an iterative regularization method. In iterative regularization, the optimization problem is solved approximately by iterating long enough to obtain a solution close to the optimal one, but terminating before too much noise occurs. Iterative regularization requires an optimization method that initially proceeds in smooth directions and makes rapid initial progress. Solving ten beamlet-based IMRT problems with dose-volume objectives and bounds on the beamlet-weights, we find that the considered optimization method fulfills the requirements for performing iterative regularization. After segment-weight optimization, the treatments obtained using 35 beamlet-weight iterations outperform the treatments obtained using 100 beamlet-weight iterations, both in terms of objective value and of target uniformity. We conclude that iterating too long may in fact deteriorate the quality of the deliverable plan. (c) 2006 American Association of Physicists in Medicine.
引用
收藏
页码:225 / 234
页数:10
相关论文
共 34 条
[1]   Dose calculations for external photon beams in radiotherapy [J].
Ahnesjö, A ;
Aspradakis, MM .
PHYSICS IN MEDICINE AND BIOLOGY, 1999, 44 (11) :R99-R155
[2]   Intensity modulated photon beams subject to a minimal surface smoothing constraint [J].
Alber, M ;
Nüsslin, F .
PHYSICS IN MEDICINE AND BIOLOGY, 2000, 45 (05) :N49-N52
[3]   On the degeneracy of the IMRT optimization problem [J].
Alber, M ;
Meedt, G ;
Nüsslin, F ;
Reemtsen, R .
MEDICAL PHYSICS, 2002, 29 (11) :2584-2589
[4]   Optimization of intensity modulated radiotherapy under constraints for static and dynamic MLC delivery [J].
Alber, M ;
Nüsslin, F .
PHYSICS IN MEDICINE AND BIOLOGY, 2001, 46 (12) :3229-3239
[5]  
ALBER M, IN PRESS OPTIM METHO
[6]   A limited-memory, quasi-Newton preconditioner for nonnegatively constrained image reconstruction [J].
Bardsley, JM .
JOURNAL OF THE OPTICAL SOCIETY OF AMERICA A-OPTICS IMAGE SCIENCE AND VISION, 2004, 21 (05) :724-731
[7]   A nonnegatively constrained convex programming method for image reconstruction [J].
Bardsley, JM ;
Vogel, CR .
SIAM JOURNAL ON SCIENTIFIC COMPUTING, 2003, 25 (04) :1326-1343
[8]   METHODS OF IMAGE-RECONSTRUCTION FROM PROJECTIONS APPLIED TO CONFORMATION RADIOTHERAPY [J].
BORTFELD, T ;
BURKELBACH, J ;
BOESECKE, R ;
SCHLEGEL, W .
PHYSICS IN MEDICINE AND BIOLOGY, 1990, 35 (10) :1423-1434
[9]   OPTIMIZATION OF STATIONARY AND MOVING BEAM RADIATION-THERAPY TECHNIQUES [J].
BRAHME, A .
RADIOTHERAPY AND ONCOLOGY, 1988, 12 (02) :129-140
[10]  
CARLSSON F, 2004, TRITAMAT2004OSI