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

Monte Carlo tree search -based non-coplanar trajectory design for station parameter optimized radiation therapy (SPORT)

An important yet challenging problem in LINAC-based rotational arc radiation therapy is the design of beam trajectory, which requires simultaneous consideration of delivery efficiency and final dose distribution. In this work, we propose a novel trajectory selection strategy by developing a Monte Carlo tree search (MCTS) algorithm during the beam trajectory selection process. To search through the vast number of possible trajectories, the MCTS algorithm was implemented. In this approach, a candidate trajectory is explored by starting from a leaf node and sequentially examining the next level of linked nodes with consideration of geometric and physical constraints. The maximum Upper Confidence Bounds for Trees, which is a function of average objective function value and the number of times the node under testing has been visited, was employed to intelligently select the trajectory. For each candidate trajectory, we run an inverse fluence map optimization with an infinity norm regularization. The ranking of the plan as measured by the corresponding objective function value was then fed back to update the statistics of the nodes on the trajectory. The method was evaluated with a chest wall and a brain case, and the results were compared with the coplanar and noncoplanar 4pi beam configurations. For both clinical cases, the MCTS method found effective and easy-to-deliver trajectories within an hour. As compared with the coplanar plans, it offers much better sparing of the OARs while maintaining the PTV coverage. The quality of the MCTS-generated plan is found to be comparable to the 4pi plans. Artificial intelligence based on MCTS is valuable to facilitate the design of beam trajectory and paves the way for future clinical use of non-coplanar treatment delivery.

Fraction-variant beam orientation optimization for non-coplanar IMRT

Conventional beam orientation optimization (BOO) algorithms for IMRT assume that the same set of beam angles is used for all treatment fractions. In this paper we present a BOO formulation based on group sparsity that simultaneously optimizes non-coplanar beam angles for all fractions, yielding a fraction-variant (FV) treatment plan. Beam angles are selected by solving a multi-fraction fluence map optimization problem involving 500-700 candidate beams per fraction, with an additional group sparsity term that encourages most candidate beams to be inactive. The optimization problem is solved using the fast iterative shrinkage-thresholding algorithm. Our FV BOO algorithm is used to create five-fraction treatment plans for digital phantom, prostate, and lung cases as well as a 30-fraction plan for a head and neck case. A homogeneous PTV dose coverage is maintained in all fractions. The treatment plans are compared with fraction-invariant plans that use a fixed set of beam angles for all fractions. The FV plans reduced OAR mean dose and D 2 values on average by 3.3% and 3.8% of the prescription dose, respectively. Notably, mean OAR dose was reduced by 14.3% of prescription dose (rectum), 11.6% (penile bulb), 10.7% (seminal vesicle), 5.5% (right femur), 3.5% (bladder), 4.0% (normal left lung), 15.5% (cochleas), and 5.2% (chiasm). D 2 was reduced by 14.9% of prescription dose (right femur), 8.2% (penile bulb), 12.7% (proximal bronchus), 4.1% (normal left lung), 15.2% (cochleas), 10.1% (orbits), 9.1% (chiasm), 8.7% (brainstem), and 7.1% (parotids). Meanwhile, PTV homogeneity defined as D 95/D 5 improved from .92 to .95 (digital phantom), from .95 to .98 (prostate case), and from .94 to .97 (lung case), and remained constant for the head and neck case. Moreover, the FV plans are dosimetrically similar to conventional plans that use twice as many beams per fraction. Thus, FV BOO offers the potential to reduce delivery time for non-coplanar IMRT.

A comprehensive formulation for volumetric modulated arc therapy planning

PURPOSE

Volumetric modulated arc therapy (VMAT) is a widely employed radiation therapy technique, showing comparable dosimetry to static beam intensity modulated radiation therapy (IMRT) with reduced monitor units and treatment time. However, the current VMAT optimization has various greedy heuristics employed for an empirical solution, which jeopardizes plan consistency and quality. The authors introduce a novel direct aperture optimization method for VMAT to overcome these limitations.

METHODS

The comprehensive VMAT (comVMAT) planning was formulated as an optimization problem with an L2-norm fidelity term to penalize the difference between the optimized dose and the prescribed dose, as well as an anisotropic total variation term to promote piecewise continuity in the fluence maps, preparing it for direct aperture optimization. A level set function was used to describe the aperture shapes and the difference between aperture shapes at adjacent angles was penalized to control MLC motion range. A proximal-class optimization solver was adopted to solve the large scale optimization problem, and an alternating optimization strategy was implemented to solve the fluence intensity and aperture shapes simultaneously. Single arc comVMAT plans, utilizing 180 beams with 2° angular resolution, were generated for a glioblastoma multiforme case, a lung (LNG) case, and two head and neck cases-one with three PTVs (H&N3PTV) and one with foue PTVs (H&N4PTV)-to test the efficacy. The plans were optimized using an alternating optimization strategy. The plans were compared against the clinical VMAT (clnVMAT) plans utilizing two overlapping coplanar arcs for treatment.

RESULTS

The optimization of the comVMAT plans had converged within 600 iterations of the block minimization algorithm. comVMAT plans were able to consistently reduce the dose to all organs-at-risk (OARs) as compared to the clnVMAT plans. On average, comVMAT plans reduced the max and mean OAR dose by 6.59% and 7.45%, respectively, of the prescription dose. Reductions in max dose and mean dose were as high as 14.5 Gy in the LNG case and 15.3 Gy in the H&N3PTV case. PTV coverages measured by D95, D98, and D99 were within 0.25% of the prescription dose. By comprehensively optimizing all beams, the comVMAT optimizer gained the freedom to allow some selected beams to deliver higher intensities, yielding a dose distribution that resembles a static beam IMRT plan with beam orientation optimization.

CONCLUSIONS

The novel nongreedy VMAT approach simultaneously optimizes all beams in an arc and then directly generates deliverable apertures. The single arc VMAT approach thus fully utilizes the digital Linac's capability in dose rate and gantry rotation speed modulation. In practice, the new single VMAT algorithm generates plans superior to existing VMAT algorithms utilizing two arcs.

Dose domain regularization of MLC leaf patterns for highly complex IMRT plans

PURPOSE

The advent of automated beam orientation and fluence optimization enables more complex intensity modulated radiation therapy (IMRT) planning using an increasing number of fields to exploit the expanded solution space. This has created a challenge in converting complex fluences to robust multileaf collimator (MLC) segments for delivery. A novel method to regularize the fluence map and simplify MLC segments is introduced to maximize delivery efficiency, accuracy, and plan quality.

METHODS

In this work, we implemented a novel approach to regularize optimized fluences in the dose domain. The treatment planning problem was formulated in an optimization framework to minimize the segmentation-induced dose distribution degradation subject to a total variation regularization to encourage piecewise smoothness in fluence maps. The optimization problem was solved using a first-order primal-dual algorithm known as the Chambolle-Pock algorithm. Plans for 2 GBM, 2 head and neck, and 2 lung patients were created using 20 automatically selected and optimized noncoplanar beams. The fluence was first regularized using Chambolle-Pock and then stratified into equal steps, and the MLC segments were calculated using a previously described level reducing method. Isolated apertures with sizes smaller than preset thresholds of 1-3 bixels, which are square units of an IMRT fluence map from MLC discretization, were removed from the MLC segments. Performance of the dose domain regularized (DDR) fluences was compared to direct stratification and direct MLC segmentation (DMS) of the fluences using level reduction without dose domain fluence regularization.

RESULTS

For all six cases, the DDR method increased the average planning target volume dose homogeneity (D95/D5) from 0.814 to 0.878 while maintaining equivalent dose to organs at risk (OARs). Regularized fluences were more robust to MLC sequencing, particularly to the stratification and small aperture removal. The maximum and mean aperture sizes using the DDR were consistently larger than those from DMS for all tested number of segments.

CONCLUSIONS

The fluence map to MLC segmentation conversion problem was formulated as a secondary optimization problem in the dose domain to minimize the smoothness-regularized dose discrepancy. The large scale optimization problem was solved using a primal-dual algorithm that transformed complicated fluences into maps that were more robust to the MLC segmentation and sequencing, affording fewer and larger segments with minimal degradation to dose distribution.

Feasibility of a unified approach to intensity-modulated radiation therapy and volume-modulated arc therapy optimization and delivery

PURPOSE

To study the feasibility of unified intensity-modulated arc therapy (UIMAT) which combines intensity-modulated radiotherapy (IMRT) and volumetric-modulated arc therapy (VMAT) optimization and delivery to produce superior radiation treatment plans, both in terms of dose distribution and efficiency of beam delivery when compared with either VMAT or IMRT alone.

METHODS

An inverse planning algorithm for UIMAT was prototyped within the pinnacle treatment planning system (Philips Healthcare). The IMRT and VMAT deliveries are unified within the same arc, with IMRT being delivered at specific gantry angles within the arc. Optimized gantry angles for the IMRT and VMAT phases are assigned automatically by the inverse optimization algorithm. Optimization of the IMRT and VMAT phases is done simultaneously using a direct aperture optimization algorithm. Five treatment plans each for prostate, head and neck, and lung were generated using a unified optimization technique and compared with clinical IMRT or VMAT plans. Delivery verification was performed with an ArcCheck phantom (Sun Nuclear) on a Varian TrueBeam linear accelerator (Varian Medical Systems).

RESULTS

In this prototype implementation, the UIMAT plans offered the same target dose coverage while reducing mean doses to organs at risk by 8.4% for head-and-neck cases, 5.7% for lung cases, and 3.5% for prostate cases, compared with the VMAT or IMRT plans. In addition, UIMAT can be delivered with similar efficiency as VMAT.

CONCLUSIONS

In this proof-of-concept work, a novel radiation therapy optimization and delivery technique that interlaces VMAT or IMRT delivery within the same arc has been demonstrated. Initial results show that unified VMAT/IMRT has the potential to be superior to either standard IMRT or VMAT.

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.

Feasibility study of volumetric modulated arc therapy with constant dose rate for endometrial cancer

To investigate the feasibility, efficiency, and delivery accuracy of volumetric modulated arc therapy with constant dose rate (VMAT-CDR) for whole-pelvic radiotherapy (WPRT) of endometrial cancer. The nine-field intensity-modulated radiotherapy (IMRT), VMAT with variable dose-rate (VMAT-VDR), and VMAT-CDR plans were created for 9 patients with endometrial cancer undergoing WPRT. The dose distribution of planning target volume (PTV), organs at risk (OARs), and normal tissue (NT) were compared. The monitor units (MUs) and treatment delivery time were also evaluated. For each VMAT-CDR plan, a dry run was performed to assess the dosimetric accuracy with MatriXX from IBA. Compared with IMRT, the VMAT-CDR plans delivered a slightly greater V20 of the bowel, bladder, pelvis bone, and NT, but significantly decreased the dose to the high-dose region of the rectum and pelvis bone. The MUs decreased from 1105 with IMRT to 628 with VMAT-CDR. The delivery time also decreased from 9.5 to 3.2 minutes. The average gamma pass rate was 95.6% at the 3%/3mm criteria with MatriXX pretreatment verification for 9 patients. VMAT-CDR can achieve comparable plan quality with significant shorter delivery time and smaller number of MUs compared with IMRT for patients with endometrial cancer undergoing WPRT. It can be accurately delivered and be an alternative to IMRT on the linear accelerator without VDR capability.