Dose optimization with first-order total-variation minimization for dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT).

PURPOSE A new treatment scheme coined as dense angularly sampled and sparse intensity modulated radiation therapy (DASSIM-RT) has recently been proposed to bridge the gap between IMRT and VMAT. By increasing the angular sampling of radiation beams while eliminating dispensable segments of the incident fields, DASSIM-RT is capable of providing improved conformity in dose distributions while maintaining high delivery efficiency. The fact that DASSIM-RT utilizes a large number of incident beams represents a major computational challenge for the clinical applications of this powerful treatment scheme. The purpose of this work is to provide a practical solution to the DASSIM-RT inverse planning problem. METHODS The inverse planning problem is formulated as a fluence-map optimization problem with total-variation (TV) minimization. A newly released L1-solver, template for first-order conic solver (TFOCS), was adopted in this work. TFOCS achieves faster convergence with less memory usage as compared with conventional quadratic programming (QP) for the TV form through the effective use of conic forms, dual-variable updates, and optimal first-order approaches. As such, it is tailored to specifically address the computational challenges of large-scale optimization in DASSIM-RT inverse planning. Two clinical cases (a prostate and a head and neck case) are used to evaluate the effectiveness and efficiency of the proposed planning technique. DASSIM-RT plans with 15 and 30 beams are compared with conventional IMRT plans with 7 beams in terms of plan quality and delivery efficiency, which are quantified by conformation number (CN), the total number of segments and modulation index, respectively. For optimization efficiency, the QP-based approach was compared with the proposed algorithm for the DASSIM-RT plans with 15 beams for both cases. RESULTS Plan quality improves with an increasing number of incident beams, while the total number of segments is maintained to be about the same in both cases. For the prostate patient, the conformation number to the target was 0.7509, 0.7565, and 0.7611 with 80 segments for IMRT with 7 beams, and DASSIM-RT with 15 and 30 beams, respectively. For the head and neck (HN) patient with a complicated target shape, conformation numbers of the three treatment plans were 0.7554, 0.7758, and 0.7819 with 75 segments for all beam configurations. With respect to the dose sparing to the critical structures, the organs such as the femoral heads in the prostate case and the brainstem and spinal cord in the HN case were better protected with DASSIM-RT. For both cases, the delivery efficiency has been greatly improved as the beam angular sampling increases with the similar or better conformal dose distribution. Compared with conventional quadratic programming approaches, first-order TFOCS-based optimization achieves far faster convergence and smaller memory requirements in DASSIM-RT. CONCLUSIONS The new optimization algorithm TFOCS provides a practical and timely solution to the DASSIM-RT or other inverse planning problem requiring large memory space. The new treatment scheme is shown to outperform conventional IMRT in terms of dose conformity to both the targetand the critical structures, while maintaining high delivery efficiency.

[1]  Lei Xing,et al.  Bridging the gap between IMRT and VMAT: dense angularly sampled and sparse intensity modulated radiation therapy. , 2011, Medical physics.

[2]  Tae-Suk Suh,et al.  Inverse planning for IMRT with nonuniform beam profiles using total-variation regularization (TVR). , 2011, Medical physics.

[3]  Emmanuel J. Candès,et al.  Templates for convex cone problems with applications to sparse signal recovery , 2010, Math. Program. Comput..

[4]  Emmanuel J. Candès,et al.  NESTA: A Fast and Accurate First-Order Method for Sparse Recovery , 2009, SIAM J. Imaging Sci..

[5]  Lei Zhu,et al.  Compressed sensing based cone-beam computed tomography reconstruction with a first-order methoda). , 2010, Medical physics.

[6]  Andrew Jackson,et al.  Volumetric modulated arc therapy: planning and evaluation for prostate cancer cases. , 2010, International journal of radiation oncology, biology, physics.

[7]  Steve B. Jiang,et al.  GPU-based fast cone beam CT reconstruction from undersampled and noisy projection data via total variation. , 2010, Medical physics.

[8]  Tae-Suk Suh,et al.  Beam's-eye-view Dosimetrics-guided inverse planning for aperture-modulated arc therapy. , 2009, International journal of radiation oncology, biology, physics.

[9]  Lei Zhu,et al.  Search for IMRT inverse plans with piecewise constant fluence maps using compressed sensing techniques. , 2009, Medical physics.

[10]  Steve Webb,et al.  Single-Arc IMRT? , 2009, Physics in medicine and biology.

[11]  Lei Xing,et al.  Using total-variation regularization for intensity modulated radiation therapy inverse planning with field-specific numbers of segments , 2008, Physics in medicine and biology.

[12]  Yin Zhang,et al.  Fixed-Point Continuation for l1-Minimization: Methodology and Convergence , 2008, SIAM J. Optim..

[13]  E. Sidky,et al.  Image reconstruction in circular cone-beam computed tomography by constrained, total-variation minimization , 2008, Physics in medicine and biology.

[14]  Jie Tang,et al.  Prior image constrained compressed sensing (PICCS): a method to accurately reconstruct dynamic CT images from highly undersampled projection data sets. , 2008, Medical physics.

[15]  Karl Otto,et al.  Volumetric modulated arc therapy: IMRT in a single gantry arc. , 2007, Medical physics.

[16]  Lei Dong,et al.  Long-term results of the M. D. Anderson randomized dose-escalation trial for prostate cancer. , 2008, International journal of radiation oncology, biology, physics.

[17]  Mário A. T. Figueiredo,et al.  Gradient Projection for Sparse Reconstruction: Application to Compressed Sensing and Other Inverse Problems , 2007, IEEE Journal of Selected Topics in Signal Processing.

[18]  Y. Nesterov Gradient methods for minimizing composite objective function , 2007 .

[19]  Karl Otto,et al.  Direct aperture optimization for IMRT using Monte Carlo generated beamlets. , 2006, Medical physics.

[20]  Joos V Lebesque,et al.  Intensity-modulated radiotherapy of breast cancer using direct aperture optimization. , 2006, Radiotherapy and oncology : journal of the European Society for Therapeutic Radiology and Oncology.

[21]  Georges Noël,et al.  Conformity index: a review. , 2006, International journal of radiation oncology, biology, physics.

[22]  Emmanuel J. Candès,et al.  Near-Optimal Signal Recovery From Random Projections: Universal Encoding Strategies? , 2004, IEEE Transactions on Information Theory.

[23]  E. Candès,et al.  Robust uncertainty principles: exact signal reconstruction from highly incomplete frequency information , 2004, IEEE Transactions on Information Theory.

[24]  Stephen P. Boyd,et al.  Convex Optimization , 2004, Algorithms and Theory of Computation Handbook.

[25]  Arvind Kumar,et al.  A Column Generation Approach to Radiation Therapy Treatment Planning Using Aperture Modulation , 2005, SIAM J. Optim..

[26]  S. Sahni,et al.  Optimal leaf sequencing with elimination of tongue-and-groove underdosage. , 2004, Physics in medicine and biology.

[27]  C. Cotrutz,et al.  Segment-based dose optimization using a genetic algorithm. , 2003, Physics in medicine and biology.

[28]  Cedric X. Yu,et al.  Guidance document on delivery, treatment planning, and clinical implementation of IMRT: report of the IMRT Subcommittee of the AAPM Radiation Therapy Committee. , 2003, Medical physics.

[29]  S Webb,et al.  Use of a quantitative index of beam modulation to characterize dose conformality: illustration by a comparison of full beamlet IMRT, few-segment IMRT (fsIMRT) and conformal unmodulated radiotherapy. , 2003, Physics in medicine and biology.

[30]  Stefan G Scheib,et al.  Quantifying the degree of conformity in radiosurgery treatment planning. , 2003, International journal of radiation oncology, biology, physics.

[31]  Arthur L. Boyer,et al.  The physics of intensity-modulated radiation therapy , 2002 .

[32]  D M Shepard,et al.  Direct aperture optimization: a turnkey solution for step-and-shoot IMRT. , 2002, Medical physics.

[33]  Uwe Oelfke,et al.  Acceleration of intensity-modulated radiotherapy dose calculation by importance sampling of the calculation matrices. , 2002, Medical physics.

[34]  P. Xia,et al.  Intensity-modulated radiotherapy in the treatment of nasopharyngeal carcinoma: an update of the UCSF experience. , 2001, International journal of radiation oncology, biology, physics.

[35]  K M Ayyangar,et al.  Leaf sequencing techniques for MLC-based IMRT. , 2001, Medical dosimetry : official journal of the American Association of Medical Dosimetrists.

[36]  P S Cho,et al.  Leaf sequencing with secondary beam blocking under leaf positioning constraints for continuously modulated radiotherapy beams. , 2001, Medical physics.

[37]  S L Hancock,et al.  Role of beam orientation optimization in intensity-modulated radiation therapy. , 2001, International journal of radiation oncology, biology, physics.

[38]  R. Garcia,et al.  Evaluation dosimétrique d'une radiothérapie conformationnelle: le facteur de conformation , 2000 .

[39]  B. Chauvet,et al.  [Dosimetric evaluation of conformal radiotherapy: conformity factor]. , 2000, Cancer radiotherapie : journal de la Societe francaise de radiotherapie oncologique.

[40]  A L Boyer,et al.  Synchronizing dynamic multileaf collimators for producing two-dimensional intensity-modulated fields with minimum beam delivery time. , 1999, International journal of radiation oncology, biology, physics.

[41]  T Bortfeld,et al.  Optimized planning using physical objectives and constraints. , 1999, Seminars in radiation oncology.

[42]  T Knöös,et al.  Volumetric and dosimetric evaluation of radiation treatment plans: radiation conformity index. , 1998, International journal of radiation oncology, biology, physics.

[43]  P. Xia,et al.  Multileaf collimator leaf sequencing algorithm for intensity modulated beams with multiple static segments. , 1998, Medical physics.

[44]  I. Kawrakow Improved modeling of multiple scattering in the Voxel Monte Carlo model. , 1997, Medical physics.

[45]  M. Moerland,et al.  A conformation number to quantify the degree of conformality in brachytherapy and external beam irradiation: application to the prostate. , 1997, International journal of radiation oncology, biology, physics.

[46]  I. Kawrakow,et al.  3D electron dose calculation using a Voxel based Monte Carlo algorithm (VMC). , 1996, Medical physics.