Dataset Information

Superiorization of projection algorithms for linearly constrained inverse radiotherapy treatment planning.

ABSTRACT:

Objective

We apply the superiorization methodology to the constrained intensity-modulated radiation therapy (IMRT) treatment planning problem. Superiorization combines a feasibility-seeking projection algorithm with objective function reduction: The underlying projection algorithm is perturbed with gradient descent steps to steer the algorithm towards a solution with a lower objective function value compared to one obtained solely through feasibility-seeking.

Approach

Within the open-source inverse planning toolkit matRad, we implement a prototypical algorithmic framework for superiorization using the well-established Agmon, Motzkin, and Schoenberg (AMS) feasibility-seeking projection algorithm and common nonlinear dose optimization objective functions. Based on this prototype, we apply superiorization to intensity-modulated radiation therapy treatment planning and compare it with (i) bare feasibility-seeking (i.e., without any objective function) and (ii) nonlinear constrained optimization using first-order derivatives. For these comparisons, we use the TG119 water phantom, the head-and-neck and the prostate patient of the CORT dataset.

Main results

Bare feasibility-seeking with AMS confirms previous studies, showing it can find solutions that are nearly equivalent to those found by the established piece-wise least-squares optimization approach. The superiorization prototype solved the linearly constrained planning problem with similar dosimetric performance to that of a general-purpose nonlinear constrained optimizer while showing smooth convergence in both constraint proximity and objective function reduction.

Significance

Superiorization is a useful alternative to constrained optimization in radiotherapy inverse treatment planning. Future extensions with other approaches to feasibility-seeking, e.g., with dose-volume constraints and more sophisticated perturbations, may unlock its full potential for high performant inverse treatment planning.

SUBMITTER: Barkmann F

PROVIDER: S-EPMC10685292 | biostudies-literature | 2023

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Superiorization of projection algorithms for linearly constrained inverse radiotherapy treatment planning.

Barkmann Florian F Censor Yair Y Wahl Niklas N

Frontiers in oncology 20231026

<h4>Objective</h4>We apply the superiorization methodology to the constrained intensity-modulated radiation therapy (IMRT) treatment planning problem. Superiorization combines a feasibility-seeking projection algorithm with objective function reduction: The underlying projection algorithm is perturbed with gradient descent steps to steer the algorithm towards a solution with a lower objective function value compared to one obtained solely through feasibility-seeking.<h4>Approach</h4>Within the o ...[more]

PMID: 38033492

Similar Datasets

Project description:Four-dimensional computed tomography (4D-CT) is commonly used to account for respiratory motion of target volumes in radiotherapy to the thorax. From the 4D-CT acquisition, a maximum-intensity projection (MIP) image set can be created and used to help define the tumor motion envelope or the internal gross tumor volume (iGTV). The purpose of this study was to quantify the differences in automatically contoured target volumes for usage in the delivery of stereotactic body radiation therapy using MIP data sets generated from one of the four methods: (1) 4D-CT phase-binned (PB) based on retrospective phase calculations, (2) 4D-CT phase-corrected phase-binned (PC-PB) based on motion extrema, (3) 4D-CT amplitude-binned (AB), and (4) cine CT built from all available images.MIP image data sets using each of the four methods were generated for a cohort of 28 patients who had prior thoracic 4D-CT scans that exhibited lung tumor motion of at least 1 cm. Each MIP image set was automatically contoured on commercial radiation treatment planning system. Margins were added to the iGTV to observe differences in the final simulated planning target volumes (PTVs).For all patients, the iGTV measured on the MIP generated from the entire cine CT data set (iGTVcine) was the largest. Expressed as a percentage of iGTVcine, 4D-CT iGTV (all sorting methods) ranged from 83.8% to 99.1%, representing differences in the absolute volume ranging from 0.02 to 4.20 cm3; the largest average and range of 4D-CT iGTV measurements was from the PC-PB data set. Expressed as a percentage of PTVcine (expansions applied to iGTVeine), the 4D-CT PTV ranged from 87.6% to 99.6%, representing differences in the absolute volume ranging from 0.08 to 7.42 cm3. Regions of the measured respiratory waveform corresponding to a rapid change of phase or amplitude showed an increased susceptibility to the selection of identical images for adjacent bins. Duplicate image selection was most common in the AB implementation, followed by the PC-PB method. The authors also found that the image associated with the minimum amplitude measurement did not always correlate with the image that showed maximum tumor motion extent.The authors identified cases in which the MIP generated from a 4D-CT sorting process under-represented the iGTV by more than 10% or up to 4.2 cm3 when compared to the iGTVcine. They suggest utilization of a MIP generated from the full cine CT data set to ensure maximum inclusive tumor extent.

Project description:PURPOSE:To develop and implement a fully automated approach to intensity modulated radiation therapy (IMRT) treatment planning. METHOD:The optimization algorithm is developed based on a hierarchical constrained optimization technique and is referred internally at our institution as expedited constrained hierarchical optimization (ECHO). Beamlet contributions to regions-of-interest are precomputed and captured in the influence matrix. Planning goals are of two classes: hard constraints that are strictly enforced from the first step (e.g., maximum dose to spinal cord), and desirable goals that are sequentially introduced in three constrained optimization problems (better planning target volume (PTV) coverage, lower organ at risk (OAR) doses, and smoother fluence map). After solving the optimization problems using external commercial optimization engines, the optimal fluence map is imported into an FDA-approved treatment planning system (TPS) for leaf sequencing and accurate full dose calculation. The dose-discrepancy between the optimization and TPS dose calculation is then calculated and incorporated into optimization by a novel dose correction loop technique using Lagrange multipliers. The correction loop incorporates the leaf sequencing and scattering effects into optimization to improve the plan quality and reduce the calculation time. The resultant optimal fluence map is again imported into TPS for leaf sequencing and final dose calculation for plan evaluation and delivery. The workflow is automated using application program interface (API) scripting, requiring user interaction solely to prepare the contours and beam arrangement prior to launching the ECHO plug-in from the TPS. For each site, parameters and objective functions are chosen to represent clinical priorities. The first site chosen for clinical implementation was metastatic paraspinal lesions treated with stereotactic body radiotherapy (SBRT). As a first step, 75 ECHO paraspinal plans were generated retrospectively and compared with clinically treated plans generated by planners using VMAT (volumetric modulated arc therapy) with 4 to 6 partial arcs. Subsequently, clinical deployment began in April, 2017. RESULTS:In retrospective study, ECHO plans were found to be dosimetrically superior with respect to tumor coverage, plan conformity, and OAR sparing. For example, the average PTV D95%, cord and esophagus max doses, and Paddick Conformity Index were improved, respectively, by 1%, 6%, 14%, and 15%, at a negligible 3% cost of the average skin D10cc dose. CONCLUSION:Hierarchical constrained optimization is a powerful and flexible tool for automated IMRT treatment planning. The dosimetric correction step accurately accounts for detailed dosimetric multileaf collimator and scattering effects. The system produces high-quality, Pareto optimal plans and avoids the time-consuming trial-and-error planning process.

Project description:PurposeTo propose and validate a fully automated multicriterial treatment planning solution for a CyberKnife® equipped with an InCiseTM 2 multileaf collimator.MethodsThe AUTO BAO plans are generated using fully automated prioritized multicriterial optimization (AUTO MCO) of pencil-beam fluence maps with integrated noncoplanar beam angle optimization (BAO), followed by MLC segment generation. Both the AUTO MCO and segmentation algorithms have been developed in-house. AUTO MCO generates for each patient a single, high-quality Pareto-optimal IMRT plan. The segmentation algorithm then accurately mimics the AUTO MCO 3D dose distribution, while considering all candidate beams simultaneously, rather than replicating the fluence maps. Pencil-beams, segment dose depositions, and final dose calculations are performed with a stand-alone version of the clinical dose calculation engine. For validation, AUTO BAO plans were generated for 33 prostate SBRT patients and compared to reference plans (REF) that were manually generated with the commercial treatment planning system (TPS), in absence of time pressure. REF plans were also compared to AUTO RB plans, for which fluence map optimization was performed for the beam angle configuration used in the REF plan, and the segmentation could use all these beams or only a subset, depending on the dosimetry.ResultsAUTO BAO plans were clinically acceptable and dosimetrically similar to REF plans, but had on average reduced numbers of beams ((beams in AUTO BAO)/(beams in REF) (relative improvement): 24.7/48.3 (-49%)), segments (59.5/98.9 (-40%)), and delivery times (17.1/22.3 min. (-23%)). Dosimetry of AUTO RB and REF were also similar, but AUTO RB used on average fewer beams (38.0/48.3 (-21%)) and had on average shorter delivery times (18.6/22.3 min. (-17%)). Delivered Monitor Units (MU) were similar for all three planning approaches.ConclusionsA new, vendor-independent optimization workflow for fully automated generation of deliverable high-quality CyberKnife® plans was proposed, including BAO. Compared to manual planning with the commercial TPS, fraction delivery times were reduced by 5.3 min. (-23%) due to large reductions in beam and segment numbers.

Dataset Information

Superiorization of projection algorithms for linearly constrained inverse radiotherapy treatment planning.

Objective

Approach

Main results

Significance

Publications

Superiorization of projection algorithms for linearly constrained inverse radiotherapy treatment planning.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets