Using total-variation regularization for intensity modulated radiation therapy inverse planning with field-specific numbers of segments

Currently, there are two types of treatment planning algorithms for intensity modulated radiation therapy (IMRT). The beamlet-based algorithm generates beamlet intensity maps with high complexity, resulting in large numbers of segments in the delivery after a leaf-sequencing algorithm is applied. The segment-based direct aperture optimization (DAO) algorithm includes the physical constraints of the deliverable apertures in the calculation, and achieves a conformal dose distribution using a small number of segments. However, the number of segments is pre-fixed in most of the DAO approaches, and the typical random search scheme in the optimization is computationally intensive. A regularization-based algorithm is proposed to overcome the drawbacks of the DAO method. Instead of smoothing the beamlet intensity maps as in many existing methods, we include a total-variation term in the optimization objective function to reduce the number of signal levels of the beam intensity maps. An aperture rectification algorithm is then applied to generate a significantly reduced number of deliverable apertures. As compared to the DAO algorithm, our method has an efficient form of quadratic optimization, with an additional advantage of optimizing field-specific numbers of segments based on the modulation complexity. The proposed approach is evaluated using two clinical cases. Under the condition that the clinical acceptance criteria of the treatment plan are satisfied, for the prostate patient, the total number of segments for five fields is reduced from 61 using the Eclipse planning system to 35 using the proposed algorithm; for the head and neck patient, the total number of segments for seven fields is reduced from 107 to 28. The head and neck result is also compared to that using an equal number of four segments for each field. The comparison shows that using field-specific numbers of segments achieves a much improved dose distribution.

[1]  L. Rudin,et al.  Nonlinear total variation based noise removal algorithms , 1992 .

[2]  T. Bortfeld,et al.  X-ray field compensation with multileaf collimators. , 1994, International journal of radiation oncology, biology, physics.

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

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

[5]  S. Webb,et al.  Inverse planning with constraints to generate smoothed intensity-modulated beams. , 1998, Physics in medicine and biology.

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

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

[8]  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.

[9]  A L Boyer,et al.  Optimization of importance factors in inverse planning. , 1999, Physics in medicine and biology.

[10]  U Oelfke,et al.  Inverse planning for x-ray rotation therapy: a general solution of the inverse problem. , 1999, Physics in medicine and biology.

[11]  F Nüsslin,et al.  Intensity modulated photon beams subject to a minimal surface smoothing constraint. , 2000, Physics in medicine and biology.

[12]  L Xing,et al.  Computer verification of fluence map for intensity modulated radiation therapy. , 2000, Medical physics.

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

[14]  M. Phillips,et al.  Reduction of computational dimensionality in inverse radiotherapy planning using sparse matrix operations. , 2001, Physics in medicine and biology.

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

[16]  J Yang,et al.  Smoothing intensity-modulated beam profiles to improve the efficiency of delivery. , 2001, Medical physics.

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

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

[19]  Lijun Ma,et al.  Smoothing intensity-modulated treatment delivery under hardware constraints. , 2002, Medical physics.

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

[21]  Lei Xing,et al.  Incorporating leaf transmission and head scatter corrections into step-and-shoot leaf sequences for IMRT. , 2003, International journal of radiation oncology, biology, physics.

[22]  H. Romeijn,et al.  A novel linear programming approach to fluence map optimization for intensity modulated radiation therapy treatment planning. , 2003, Physics in medicine and biology.

[23]  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.

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

[25]  Ping Xia,et al.  A new smoothing procedure to reduce delivery segments for static MLC-based IMRT planning. , 2004, Medical physics.

[26]  Yair Censor,et al.  Inherent smoothness of intensity patterns for intensity modulated radiation therapy generated by simultaneous projection algorithms. , 2004, Physics in medicine and biology.

[27]  Yair Censor,et al.  The dose–volume constraint satisfaction problem for inverse treatment planning with field segments , 2004, Physics in medicine and biology.

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

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

[30]  Jari P. Kaipio,et al.  Parallelized Bayesian inversion for three-dimensional dental X-ray imaging , 2006, IEEE Transactions on Medical Imaging.

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

[32]  Sebastiaan Breedveld,et al.  Fast, multiple optimizations of quadratic dose objective functions in IMRT , 2006, Physics in medicine and biology.

[33]  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.

[34]  K. T. Block,et al.  Undersampled radial MRI with multiple coils. Iterative image reconstruction using a total variation constraint , 2007, Magnetic resonance in medicine.

[35]  S. Webb,et al.  Direct-aperture optimization applied to selection of beam orientations in intensity-modulated radiation therapy , 2007, Physics in medicine and biology.

[36]  Alan Nichol,et al.  Direct aperture optimization for online adaptive radiation therapy. , 2007, Medical physics.

[37]  Benedick A Fraass,et al.  Reduction of IMRT beam complexity through the use of beam modulation penalties in the objective function. , 2007, Medical physics.