Beam's-Eye-View Dosimetrics–Guided Inverse Planning for Aperture-Modulated Arc Therapy

A Prospective 4π Radiation Therapy Clinical Study in Recurrent High-Grade Glioma Patients

PURPOSE

To evaluate the feasibility, safety, dosimetric benefits, delivery efficiency, and patient comfort in the clinical implementation of 4π radiation therapy.

METHODS AND MATERIALS

Eleven patients with recurrent high-grade glioma were recruited for the trial. 4π plans integrating beam orientation and fluence-map optimization were created using an in-house column-generation algorithm. The collision-free beam solution space throughout the 4π steradian was determined using a computer-aided-design model of the Varian TrueBeam system and a human subject. Twenty beams were optimized for each case and imported into Eclipse for intensity modulated radiation therapy planning. Beam orientations with neighboring couch kicks were merged for increased delivery efficiency, generating plans with an average of 16 beam orientations. Volumetric modulated arc therapy (VMAT) plans with 3-4 arcs were also generated for each case, and the plan achieving superior dosimetric quality was selected for treatment. Patient comfort was surveyed after every fraction. Multiple 2-dimensional X-ray images were obtained to measure intrafractional motion.

RESULTS

Of 11 patients, 9 were treated with 4π. Mean and maximum organ at risk doses were equal or significantly lower (P < .05) with 4π than with VMAT. Particularly substantial dose reduction of 2.92 Gy in the average accumulated brainstem maximum dose enabled treatments that would otherwise not satisfy safe dose constraints with VMAT. One patient was not treated because neither plan met the dosimetric criteria. The other was treated with VMAT owing to comparable dosimetry resulting from a planning target volume located in a separate co-plane superior to organs at risk. Treatments were well tolerated, with an average patient comfort score of 8.6/10. Intrafractional motion was <1.5 mm for all delivered fractions, and the average delivery time was 34.1 minutes.

CONCLUSIONS

The feasibility, safety, dosimetric benefits, delivery efficiency, and patient comfort of 4π radiation therapy have been clinically demonstrated with a prospective clinical trial. The results elucidate the potential and challenges of wider clinical implementations.

Trajectory optimization in radiotherapy using sectioning (TORUS)

PURPOSE

A challenging problem in trajectory optimization for radiotherapy is properly handling the synchronization of the medical accelerators dynamic delivery. The initial coarse sampling of control points implemented in a Progressive Resolution Optimization type approach (VMAT) routinely results in MLC aperture forming contention issues as the sampling resolution increases. This work presents an approach to optimize continuous, beam-on radiation trajectories through exploration of the anatomical topology present in the patient and formation of a novel dual-metric graph optimization problem.

METHODS

This work presents a new perspective on trajectory optimization in radiotherapy using the concept of sectioning (TORUS). TORUS avoids degradation of 3D dose optimization quality by mapping the connectedness of target regions from the BEV perspective throughout the space of deliverable coordinates. This connectedness information is then incorporated into a graph optimization problem to define ideal trajectories. The unique usage of two distance functions in this graph optimization permits the TORUS algorithm to generate efficient dynamic trajectories for delivery while maximizing the angular flux through all PTV voxels. 3D dose optimization is performed for trajectories using a commercial TPS progressive resolution optimizer.

RESULTS

The TORUS algorithm is applied to three example treatments: chest wall, scalp, and the TG-119 C-shape phantom. When static collimator coplanar trajectories are generated for the chest wall and scalp cases, the TORUS trajectories are found to outperform both 7 field IMRT and 2 arc VMAT plans in delivery time, organ at risk sparing, conformality, and homogeneity. For the TG-119 phantom, when static couch and collimator non-coplanar trajectories are optimized, TORUS trajectories have superior sparing of the central core avoidance with shorter delivery times, with similar dose conformality and homogeneity.

CONCLUSIONS

The TORUS algorithm is able to automatically generate trajectories having improved plan quality and delivery time over standard IMRT and VMAT treatments. TORUS offers an exciting and promising avenue forward toward increasing dynamic capabilities in radiation delivery.

Beam’s-eye-view dosimetrics (BEVD) guided rotational station parameter optimized radiation therapy (SPORT) planning based on reweighted total-variation minimization

Conventional VMAT optimizes aperture shapes and weights at uniformly sampled stations, which is a generalization of the concept of a control point. Recently, rotational station parameter optimized radiation therapy (SPORT) has been proposed to improve the plan quality by inserting beams to the regions that demand additional intensity modulations, thus formulating nonuniform beam sampling. This work presents a new rotational SPORT planning strategy based on reweighted total-variation (TV) minimization (min.), using beam’s-eye-view dosimetrics (BEVD) guided beam selection. The convex programming based reweighted TV min. assures the simplified fluence-map, which facilitates single-aperture selection at each station for single-arc delivery. For the rotational arc treatment planning and non-uniform beam angle setting, the mathematical model needs to be modified by additional penalty term describing the fluence-map similarity and by determination of appropriate angular weighting factors. The proposed algorithm with additional penalty term is capable of achieving more efficient and deliverable plans adaptive to the conventional VMAT and SPORT planning schemes by reducing the dose delivery time about 5 to 10 s in three clinical cases (one prostate and two head-and-neck (HN) cases with a single and multiple targets). The BEVD guided beam selection provides effective and yet easy calculating methodology to select angles for denser, non-uniform angular sampling in SPORT planning. Our BEVD guided SPORT treatment schemes improve the dose sparing to femoral heads in the prostate and brainstem, parotid glands and oral cavity in the two HN cases, where the mean dose reduction of those organs ranges from 0.5 to 2.5 Gy. Also, it increases the conformation number assessing the dose conformity to the target from 0.84, 0.75 and 0.74 to 0.86, 0.79 and 0.80 in the prostate and two HN cases, while preserving the delivery efficiency, relative to conventional single-arc VMAT plans.

Total-variation regularization based inverse planning for intensity modulated arc therapy

Intensity modulated arc therapy (IMAT) delivers conformal dose distributions through continuous gantry rotation with constant or variable speed while modulating the field aperture shape and weight. The enlarged angular space and machine delivery constraints make inverse planning of IMAT more intractable as compared to its counterpart of fixed gantry IMRT. Currently, IMAT inverse planning is being done using two extreme methods: the first one computes in beamlet domain with a subsequent arc leaf sequencing, and the second proceeds in machine parameter domain with entire emphasis placed on a pre-determined delivery method without exploring potentially better alternative delivery schemes. Towards truly optimizing the IMAT treatment on a patient specific basis, in this work we propose a total-variation based inverse planning framework for IMAT, which takes advantage of the useful features of the above two existing approaches while avoiding their shortcomings. A quadratic optimization algorithm has been implemented to demonstrate the performance and advantage of the proposed approach. Applications of the technique to a prostate case and a head and neck case indicate that the algorithm is capable of generating IMAT plans with patient specific numbers of arcs efficiently. Superior dose distributions and delivery time are achieved with a maximum number of apertures of three for each field. As compared to conventional beamlet-based algorithms, our method regularizes the field modulation complexity during optimization, and permits us to obtain the best possible plan with a pre-set modulation complexity of fluences. As illustrated in both prostate and head-and-neck case studies, the proposed method produces more favorable dose distributions than the segment-based algorithms, by optimally accommodating the clinical need of intensity modulation levels for each individual field. On a more fundamental level, our formulation preserves the convexity of optimization and makes the search of the global optimal solution possible with a deterministic method.

Nonisocentric treatment strategy for breast radiation therapy: a proof of concept study

PURPOSE

To propose a nonisocentric treatment strategy as a special form of station parameter optimized radiation therapy, to improve sparing of critical structures while preserving target coverage in breast radiation therapy.

METHODS AND MATERIALS

To minimize the volume of exposed lung and heart in breast irradiation, we propose a novel nonisocentric treatment scheme by strategically placing nonconverging beams with multiple isocenters. As its name suggests, the central axes of these beams do not intersect at a single isocenter as in conventional breast treatment planning. Rather, the isocenter locations and beam directions are carefully selected, in that each beam is only responsible for a certain subvolume of the target, so as to minimize the volume of irradiated normal tissue. When put together, the beams will provide an adequate coverage of the target and expose only a minimal amount of normal tissue to radiation. We apply the nonisocentric planning technique to 2 previously treated clinical cases (breast and chest wall).

RESULTS

The proposed nonisocentric technique substantially improved sparing of the ipsilateral lung. Compared with conventional isocentric plans using 2 tangential beams, the mean lung dose was reduced by 38% and 50% using the proposed technique, and the volume of the ipsilateral lung receiving ≥ 20 Gy was reduced by a factor of approximately 2 and 3 for the breast and chest wall cases, respectively. The improvement in lung sparing is even greater compared with volumetric modulated arc therapy.

CONCLUSIONS

A nonisocentric implementation of station parameter optimized radiation therapy has been proposed for breast radiation therapy. The new treatment scheme overcomes the limitations of existing approaches and affords a useful tool for conformal breast radiation therapy, especially in cases with extreme chest wall curvature.

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.

Bridging the gap between IMRT and VMAT: dense angularly sampled and sparse intensity modulated radiation therapy

PURPOSE

To propose an alternative radiation therapy (RT) planning and delivery scheme with optimal angular beam sampling and intrabeam modulation for improved dose distribution while maintaining high delivery efficiency.

METHODS

In the proposed approach, coined as dense angularly sampled and sparse intensity modulated RT (DASSIM-RT), a large number of beam angles are used to increase the angular sampling, leading to potentially more conformal dose distributions as compared to conventional IMRT. At the same time, intensity modulation of the incident beams is simplified to eliminate the dispensable segments, compensating the increase in delivery time caused by the increased number of beams and facilitating the plan delivery. In a sense, the proposed approach shifts and transforms, in an optimal fashion, some of the beam segments in conventional IMRT to the added beams. For newly available digital accelerators, the DASSIM-RT delivery can be made very efficient by concatenating the beams so that they can be delivered sequentially without operator's intervention. Different from VMAT, the level of intensity modulation in DASSIS-RT is field specific and optimized to meet the need of each beam direction. Three clinical cases (a head and neck (HN) case, a pancreas case, and a lung case) are used to evaluate the proposed RT scheme. DASSIM-RT, VMAT, and conventional IMRT plans are compared quantitatively in terms of the conformality index (CI) and delivery efficiency.

RESULTS

Plan quality improves generally with the number and intensity modulation of the incident beams. For a fixed number of beams or fixed level of intensity modulation, the improvement saturates after the intensity modulation or number of beams reaches to a certain level. An interplay between the two variables is observed and the saturation point depends on the values of both variables. For all the cases studied here, the CI of DASSIM-RT with 15 beams and 5 intensity levels (0.90, 0.79, and 0.84 for the HN, pancreas, and lung cases, respectively) is similar with that of conventional IMRT with seven beams and ten intensity levels (0.88, 0.79, and 0.83) and is higher than that of single-arc VMAT (0.75, 0.75, and 0.82). It is also found that the DASSIM-RT plans generally have better sparing of organs-at-risk than IMRT plans. It is estimated that the dose delivery time of DASSIM-RT with 15 beams and 5 intensity levels is about 4.5, 4.4, and 4.2 min for the HN, pancreas, and lung case, respectively, similar to that of IMRT plans with 7 beams and 10 intensity levels.

CONCLUSION

DASSIS-RT bridges the gap between IMRT and VMAT and allows optimal sampling of angular space and intrabeam modulation, thus it provides improved conformity in dose distributions while maintaining high delivery efficiency.

Choreographing couch and collimator in volumetric modulated arc therapy

PURPOSE

To design and optimize trajectory-based, noncoplanar subarcs for volumetric modulated arc therapy (VMAT) deliverable on both Varian TrueBEAM system and traditional accelerators; and to investigate their potential advantages for treating central nervous system (CNS) tumors.

METHODS AND MATERIALS

To guide the computerized selection of beam trajectories consisting of simultaneous couch, gantry, and collimator motion, a score function was implemented to estimate the geometric overlap between targets and organs at risk for each couch/gantry angle combination. An initial set of beam orientations is obtained as a function of couch and gantry angle, according to a minimum search of the score function excluding zones of collision. This set is grouped into multiple continuous and extended subarcs subject to mechanical limitations using a hierarchical clustering algorithm. After determination of couch/gantry trajectories, a principal component analysis finds the collimator angle at each beam orientation that minimizes residual target-organ at risk overlaps. An in-house VMAT optimization algorithm determines the optimal multileaf collimator position and monitor units for control points within each subarc. A retrospective study of 10 CNS patients compares the proposed method of VMAT trajectory with dynamic gantry, leaves, couch, and collimator motion (Tra-VMAT); a standard noncoplanar VMAT with no couch/collimator motion within subarcs (Std-VMAT); and noncoplanar intensity-modulated radiotherapy (IMRT) plans that were clinically used.

RESULTS

Tra-VMAT provided improved target dose conformality and lowered maximum dose to brainstem, optic nerves, and chiasm by 7.7%, 1.1%, 2.3%, and 1.7%, respectively, compared with Std-VMAT. Tra-VMAT provided higher planning target volume minimum dose and reduced maximum dose to chiasm, optic nerves, and cochlea by 6.2%, 1.3%, 6.3%, and 8.4%, respectively, and reduced cochlea mean dose by 8.7%, compared with IMRT. Tra-VMAT averaged beam-on time was comparable to Std-VMAT but significantly (45%) less than IMRT.

CONCLUSION

Optimized couch, gantry, and collimator trajectories may be integrated into VMAT with improved mechanical flexibility and may provide better dosimetric properties and improved efficiency in the treatment of CNS tumors.

Inverse planning for four-dimensional (4D) volumetric modulated arc therapy

PURPOSE

To develop a 4D volumetric modulated arc therapy (VMAT) inverse planning framework.

METHODS

4D VMAT inverse planning aims to derive an aperture and weight modulated arc therapy treatment plan that optimizes the accumulated dose distribution from all gantry angles and breathing phases. Under an assumption that the gantry rotation and patient breathing are synchronized (i.e., there is a functional relationship between the phase of the patient breathing cycle and the beam angle), the authors compute the contribution from different respiration phases through the registration of the phased CT images. The accumulative dose distribution is optimized by iteratively adjusting the aperture shape and weight of each beam through the minimization of the planning objective function. For comparison, traditional 3D VMAT plans are also performed for the two cases and the performance of the proposed technique is demonstrated.

RESULTS

A framework for 4D VMAT inverse planning has been proposed. With the consideration of the extra dimension of time in VMAT, a tighter target margin can be achieved with a full duty cycle, which is otherwise not achievable simultaneously by either 3D VMAT optimization or gated VMAT.

CONCLUSIONS

The 4D VMAT planning formulism proposed here provides useful insight on how the "time" dimension can be exploited in rotational arc therapy to maximally compensate for the intrafraction organ motion.

A method for optimizing LINAC treatment geometry for volumetric modulated arc therapy of multiple brain metastases

PURPOSE

Volumetric modulated arc therapy (VMAT) is a rotational delivery technique in which MLC shapes, dose rate, and gantry rotation speed are optimized to produce conformal dose distributions. The aim of this work is to develop a beam projection method for deriving the optimal table and collimator angles for multilesion treatment planning.

METHODS

The method consists of four steps. The first step is to define the vector of beam-eye-view (BEV)-Y-axis in the treatment planning CT coordinates. The second step is to project each target onto the BEV-Y-axis vector. In the third step, the best table and collimator angle are found with a brute-force optimization technique that minimizes MLC leaf sharing between lesions. The fourth step is to generate an optimized VMAT plan with appropriate table/collimator angles and evaluate the plan quality.

RESULTS

The authors tested the method on three example cases with targets of various locations in the brain and sizes ranging from 1.18 to 17.86 cm(3). Applying the optimized geometric parameter to generate VMAT plan, a reduction of the 12 Gy volume was more than 6.1% for all cases; the plan homogeneity (D2%-D95%) was improved from 5.88 +/- 1.21 to 5.21 +/- 0.93 Gy vs a VMAT plan with the manufacturer recommended table and collimator angles.

CONCLUSIONS

The authors conclude that the use of the projection method minimizes the sharing of MLC leaves between lesions and improves the plan quality for multilesion VMAT delivery.

A unified framework for 3D radiation therapy and IMRT planning: plan optimization in the beamlet domain by constraining or regularizing the fluence map variations

The purpose of this work is to demonstrate that physical constraints on fluence gradients in 3D radiation therapy (RT) planning can be incorporated into beamlet optimization explicitly by direct constraint on the spatial variation of the fluence maps or implicitly by using total-variation regularization (TVR). The former method forces the fluence to vary in accordance with the known form of a wedged field and latter encourages the fluence to take the known form of the wedged field by requiring the derivatives of the fluence maps to be piece-wise constant. The performances of the proposed methods are evaluated by using a brain cancer case and a head and neck case. It is found that both approaches are capable of providing clinically sensible 3D RT solutions with monotonically varying fluence maps. For currently available 3D RT delivery schemes based on the use of customized physical or dynamic wedges, constrained optimization seems to be more useful because the optimized fields are directly deliverable. Working in the beamlet domain provides a natural way to model the spatial variation of the beam fluence. The proposed methods take advantage of the fact that 3D RT is a special form of intensity-modulated radiation therapy (IMRT) and finds the optimal plan by searching for fields with a certain type of spatial variation. The approach provides a unified framework for 3D CRT and IMRT plan optimization.

Arc-modulated radiation therapy based on linear models

This paper reports an inverse arc-modulated radiation therapy planning technique based on linear models. It is implemented with a two-step procedure. First, fluence maps for 36 fixed-gantry beams are generated using a linear model-based intensity-modulated radiation therapy (IMRT) optimization algorithm. The 2D fluence maps are decomposed into 1D fluence profiles according to each leaf pair position. Second, a mixed integer linear model is used to construct the leaf motions of an arc delivery that reproduce the 1D fluence profile previously derived from the static gantry IMRT optimization. The multi-leaf collimator (MLC) sequence takes into account the starting and ending leaf positions in between the neighbouring apertures, such that the MLC segments of the entire treatment plan are deliverable in a continuous arc. Since both steps in the algorithm use linear models, implementation is simple and straightforward. Details of the algorithm are presented, and its conceptual correctness is verified with clinical cases representing prostate and head-and-neck treatments.