iCycle: Integrated, multicriterial beam angle, and profile optimization for generation of coplanar and noncoplanar IMRT plans

Application and progress of artificial intelligence in radiation therapy dose prediction

Radiation therapy (RT) nowadays is a main treatment modality of cancer. To ensure the therapeutic efficacy of patients, accurate dose distribution is often required, which is a time-consuming and labor-intensive process. In addition, due to the differences in knowledge and experience among participants and diverse institutions, the predicted dose are often inconsistent. In last several decades, artificial intelligence (AI) has been applied in various aspects of RT, several products have been implemented in clinical practice and confirmed superiority. In this paper, we will review the research of AI in dose prediction, focusing on the progress in deep learning (DL).

Lexicographic optimization-based planning for stereotactic radiosurgery of brain metastases

AIM

To validate a fully-automated lexicographic optimization-planning system (mCycle, Elekta) for single-(SL) and multiple-(ML, up to 4 metastases) lesions in intracranial stereotactic radiosurgery (SRS, 21 Gy, single fraction).

METHODS

A pre-determined priority list, Wish-List (WL), represents a dialogue between planner and clinician, establishing strict constraints and pursuing objectives. In order to satisfy the clinical protocol without manual intervention, four patients were required to tweak and fine-tune each WL (SLp, MLp) for coplanar arcs. Thirty-five testing plans (20 SLp, 15 MLp) were automatically re-planned (mCP). Automatic and manual plans were compared including dose constraints, conformality, modulation complexity score (MCS), delivery time, and local gamma analysis (2%/2 mm). To ensure plan clinical acceptability, two radiation oncologists conducted an independent blind plan choice.

RESULTS

Each WL-tuning took 3 days. Estimated median manual plans and mCP calculation time were 8 and 3 h, respectively. Significant increases in SLp and MLp target coverage and conformity were registered. mCP showed a not significant and clinically acceptable higher median brain V12Gy. SLp registered a -5.8% MU decrease with comparable median delivery time (MP 2.0 min, mCP 1.9 min) while MLp showed a +9.8% MU increase and longer delivery time (MP 3.5 min, mCP 4.4 min). mCP MCS resulted significantly higher without affecting gamma passing rates. At blind choice, mCP were preferred in the majority of cases.

CONCLUSIONS

Lexicographic optimization produced acceptable SRS plans with coplanar arcs significantly reducing the overall planning time in cases with up to 4 brain metastases. These planning improvements suggest further investigations by setting high-quality non-coplanar arc plans as a reference.

Validation of Fully Automated Robust Multicriterial Treatment Planning for Head and Neck Cancer IMPT

PURPOSE

Our purpose was to compare robust intensity modulated proton therapy (IMPT) plans, automatically generated with wish-list-based multicriterial optimization as implemented in Erasmus-iCycle, with manually created robust clinical IMPT plans for patients with head and neck cancer.

METHODS AND MATERIALS

Thirty-three patients with head and neck cancer were retrospectively included. All patients were previously treated with a manually created IMPT plan with 7000 cGy dose prescription to the primary tumor (clinical target volume [CTV]7000) and 5425 cGy dose prescription to the bilateral elective volumes (CTV5425). Plans had a 4-beam field configuration and were generated with scenario-based robust optimization (21 scenarios, 3-mm setup error, and ±3% density uncertainty for the CTVs). Three clinical plans were used to configure the Erasmus-iCycle wish-list for automated generation of robust IMPT plans for the other 30 included patients, in line with clinical planning requirements. Automatically and manually generated IMPT plans were compared for (robust) target coverage, organ-at-risk (OAR) doses, and normal tissue complication probabilities (NTCP). No manual fine-tuning of automatically generated plans was performed.

RESULTS

For all automatically generated plans, voxel-wise minimum D98% values for the CTVs were within clinical constraints and similar to manual plans. All investigated OAR parameters were favorable in the automatically generated plans (all P < .001). Median reductions in mean dose to OARs went up to 667 cGy for the inferior pharyngeal constrictor muscle, and median reductions in D0.03cm3 in serial OARs ranged up to 1795 cGy for the spinal cord surface. The observed lower mean dose in parallel OARs resulted in statistically significant lower NTCP for xerostomia (grade ≥2: 34.4% vs 38.0%; grade ≥3: 9.0% vs 10.2%) and dysphagia (grade ≥2: 11.8% vs 15.0%; grade ≥3: 1.8% vs 2.8%).

CONCLUSIONS

Erasmus-iCycle was able to produce IMPT dose distributions fully automatically with similar (robust) target coverage and improved OAR doses and NTCPs compared with clinical manual planning, with negligible hands-on planning workload.

Reducing the lateral dose penumbra in IMPT by incorporating transmission pencil beams

OBJECTIVE

In intensity-modulated proton therapy (IMPT), Bragg peaks result in steep distal dose fall-offs, while the lateral IMPT dose fall-off is often less steep than in photon therapy. High-energy pristine transmission ('shoot through') pencil beams have no Bragg peak in the patient, but show a sharp lateral penumbra at the target level. We investigated whether combining Bragg peaks with Transmission pencil beams ('IMPT&TPB') could improve head-and-neck plans by exploiting the steep lateral dose fall-off of transmission pencil beams.

APPROACH

Our system for automated multi-criteria IMPT plan optimisation was extended for combined optimisation of BPs and TPBs. The system generates for each patient a Pareto-optimal plan using a generic 'wish-list' with prioritised planning objectives and hard constraints. For eight nasopharynx cancer patients (NPC) and eight oropharynx cancer (OPC) patients, the IMPT&TPB plan was compared to the competing conventional IMPT plan with only Bragg peaks, which was generated with the same optimiser, but without transmission pencil beams.

MAIN RESULTS

Clinical OAR and target constraints were met in all plans. By allowing transmission pencil beams in the optimisation, on average 14 of the 25 investigated OAR plan parameters significantly improved for NPC, and 9 of the 17 for OPC, while only one OPC parameter showed small but significant deterioration. Non-significant differences were found in the remaining parameters. In NPC, cochlea Dmean reduced by up to 17.5 Gy and optic nerve D2% by up to 11.1 Gy.

CONCLUSION

Compared to IMPT, IMPT&TPB resulted in comparable target coverage with overall superior OAR sparing, the latter originating from steeper dose fall-offs close to OARs.

Dosimetric impact of bone marrow sparing for robustly optimized IMPT for locally advanced cervical cancer

BACKGROUND AND PURPOSE

To investigate the trade-off between bone marrow sparing (BMS) and dose to organs at risk (OARs) for intensity modulated proton therapy (IMPT) for women with locally advanced cervical cancer (LACC).

MATERIALS AND METHODS

Twenty LACC patients were retrospectively included. IMPT plans were created for each patient using automated treatment planning. These plans progressively reduced bone marrow mean doses by steps of 1 GyRBE, while constraining target coverage and conformality. The relation between bone marrow dose and bladder, small bowel, rectum, and sigmoid doses was evaluated.

RESULTS

A total of 140 IMPT plans were created. Plans without BMS had an average [range] bone marrow mean dose of 17.3 [14.7-21.6] GyRBE , which reduced to 12.0 [10.0-14.0] GyRBE with maximum BMS. The mean OAR dose [range] increased modestly for 1 GyRBE BMS: 0.2 [0.0 - 0.6] GyRBE for bladder, 0.3 [-0.2 - 0.7] GyRBE for rectum, 0.4 [0.1 - 0.8] GyRBE for small bowel, and 0.2 [-0.2 - 0.4] GyRBE for sigmoid. Moreover, for maximum BMS, mean OAR doses [range] escalated by 3.3 [0.1 - 6.7] GyRBE for bladder, 5.8 [1.8 - 12.4] GyRBE for rectum, 3.9 [1.6 - 5.9] GyRBE for small bowel, and 2.7 [0.6 - 5.9] GyRBE for sigmoid.

CONCLUSION

Achieving 1 GyRBE BMS for IMPT is feasible for LACC patients with limited dosimetric impact on other OARs. While further bone marrow dose reduction is possible for some patients, it may increase OAR doses substantially for others. Hence, we recommend a personalized approach when introducing BMS into clinical IMPT treatment planning to carefully assess individual patient benefits and risks.

Dosimetric evaluation of different planning strategies for hypofractionated whole-breast irradiation technique

Purpose.The dose hotspot areas in hypofractionated whole-breast irradiation (WBI) greatly increase the risk of acute skin toxicity because of the anatomical peculiarities of the breast. In this study, we presented several novel planning strategies that integrate multiple sub-planning target volumes (sub-PTVs), field secondary placement, and RapidPlan models for right-sided hypofractionated WBI.Methods.A total of 35 cases of WBI with a dose of 42.5 Gy for PTVs using tangential intensity-modulated radiotherapy (IMRT) were selected. Both PTVs were planned for simultaneous treatment using the original manual multiple sub-PTV plan (OMMP) and the original manual single-PTV plan (OMSP). The manual field secondary placement multiple sub-PTV plan (m-FSMP) with multiple objects on the original PTV and the manual field secondary placement single-objective plan (m-FSSP) were initially planned, which were distribution-based of V105 (volume receiving 105% of the prescription dose). In addition, two RapidPlan-based plans were developed, including the RapidPlan-based multiple sub-PTVs plan (r-FSMP) and the RapidPlan-based single-PTV plan (r-FSSP). Dosimetric parameters of the plans were compared, and V105 was evaluated using multivariate analysis to determine how it was related to the volume of PTV and the interval of lateral beam angles (ILBA).Results.The lowest mean V105 (5.64 ± 6.5%) of PTV was observed in m-FSMP compared to other manual plans. Upon validation, r-FSSP demonstrated superior dosimetric quality for OAR compared to the two other manual planning methods, except for V5(the volume of ipsilateral lung receiving 5 Gy) of the ipsilateral lung. While r-FSMP showed no significant difference (p = 0.06) compared to r-FSSP, it achieved the lowest V105 value (4.3 ± 4.5%), albeit with a slight increase in the dose to some OARs. Multivariate GEE linear regression showed that V105 is significantly correlated with target volume and ILBA.Conclusions.m-FSMP and r-FSMP can substantially enhance the homogeneity index (HI) and reduce V105, thereby minimizing the risk of acute skin toxicities, even though there may be a slight dose compromise for certain OARs.

Penalty weight tuning in high dose rate brachytherapy using multi-objective Bayesian optimization

Objective.Treatment plan optimization in high dose rate brachytherapy often requires manual fine-tuning of penalty weights for each objective, which can be time-consuming and dependent on the planner's experience. To automate this process, this study used a multi-criteria approach called multi-objective Bayesian optimization with q-noisy expected hypervolume improvement as its acquisition function (MOBO-qNEHVI).Approach.The treatment plans of 13 prostate cancer patients were retrospectively imported to a research treatment planning system, RapidBrachyMTPS, where fast mixed integer optimization (FMIO) performs dwell time optimization given a set of penalty weights to deliver 15 Gy to the target volume. MOBO-qNEHVI was used to find patient-specific Pareto optimal penalty weight vectors that yield clinically acceptable dose volume histogram metrics. The relationship between the number of MOBO-qNEHVI iterations and the number of clinically acceptable plans per patient (acceptance rate) was investigated. The performance time was obtained for various parameter configurations.Main results.MOBO-qNEHVI found clinically acceptable treatment plans for all patients. With increasing the number of MOBO-qNEHVI iterations, the acceptance rate grew logarithmically while the performance time grew exponentially. Fixing the penalty weight of the tumour volume to maximum value, adding the target dose as a parameter, initiating MOBO-qNEHVI with 25 parallel sampling of FMIO, and running 6 MOBO-qNEHVI iterations found solutions that delivered 15 Gy to the hottest 95% of the clinical target volume while respecting the dose constraints to the organs at risk. The average acceptance rate for each patient was 89.74% ± 8.11%, and performance time was 66.6 ± 12.6 s. The initiation took 22.47 ± 7.57 s, and each iteration took 7.35 ± 2.45 s to find one Pareto solution.Significance.MOBO-qNEHVI combined with FMIO can automatically explore the trade-offs between treatment plan objectives in a patient specific manner within a minute. This approach can reduce the dependency of plan quality on planner's experience and reduce dose to the organs at risk.

Multi-institutional evaluation of a Pareto navigation guided automated radiotherapy planning solution for prostate cancer

BACKGROUND

Current automated planning solutions are calibrated using trial and error or machine learning on historical datasets. Neither method allows for the intuitive exploration of differing trade-off options during calibration, which may aid in ensuring automated solutions align with clinical preference. Pareto navigation provides this functionality and offers a potential calibration alternative. The purpose of this study was to validate an automated radiotherapy planning solution with a novel multi-dimensional Pareto navigation calibration interface across two external institutions for prostate cancer.

METHODS

The implemented 'Pareto Guided Automated Planning' (PGAP) methodology was developed in RayStation using scripting and consisted of a Pareto navigation calibration interface built upon a 'Protocol Based Automatic Iterative Optimisation' planning framework. 30 previous patients were randomly selected by each institution (IA and IB), 10 for calibration and 20 for validation. Utilising the Pareto navigation interface automated protocols were calibrated to the institutions' clinical preferences. A single automated plan (VMATAuto) was generated for each validation patient with plan quality compared against the previously treated clinical plan (VMATClinical) both quantitatively, using a range of DVH metrics, and qualitatively through blind review at the external institution.

RESULTS

PGAP led to marked improvements across the majority of rectal dose metrics, with Dmean reduced by 3.7 Gy and 1.8 Gy for IA and IB respectively (p < 0.001). For bladder, results were mixed with low and intermediate dose metrics reduced for IB but increased for IA. Differences, whilst statistically significant (p < 0.05) were small and not considered clinically relevant. The reduction in rectum dose was not at the expense of PTV coverage (D98% was generally improved with VMATAuto), but was somewhat detrimental to PTV conformality. The prioritisation of rectum over conformality was however aligned with preferences expressed during calibration and was a key driver in both institutions demonstrating a clear preference towards VMATAuto, with 31/40 considered superior to VMATClinical upon blind review.

CONCLUSIONS

PGAP enabled intuitive adaptation of automated protocols to an institution's planning aims and yielded plans more congruent with the institution's clinical preference than the locally produced manual clinical plans.

A fast and robust constraint-based online re-optimization approach for automated online adaptive intensity modulated proton therapy in head and neck cancer

Objective. In head-and-neck cancer intensity modulated proton therapy, adaptive radiotherapy is currently restricted to offline re-planning, mitigating the effect of slow changes in patient anatomies. Daily online adaptations can potentially improve dosimetry. Here, a new, fully automated online re-optimization strategy is presented. In a retrospective study, this online re-optimization approach was compared to our trigger-based offline re-planning (offlineTBre-planning) schedule, including extensive robustness analyses.Approach. The online re-optimization method employs automated multi-criterial re-optimization, using robust optimization with 1 mm setup-robustness settings (in contrast to 3 mm for offlineTBre-planning). Hard planning constraints and spot addition are used to enforce adequate target coverage, avoid prohibitively large maximum doses and minimize organ-at-risk doses. For 67 repeat-CTs from 15 patients, fraction doses of the two strategies were compared for the CTVs and organs-at-risk. Per repeat-CT, 10.000 fractions with different setup and range robustness settings were simulated using polynomial chaos expansion for fast and accurate dose calculations.Main results. For 14/67 repeat-CTs, offlineTBre-planning resulted in <50% probability ofD98%≥ 95% of the prescribed dose (Dpres) in one or both CTVs, which never happened with online re-optimization. With offlineTBre-planning, eight repeat-CTs had zero probability of obtainingD98%≥ 95%Dpresfor CTV7000, while the minimum probability with online re-optimization was 81%. Risks of xerostomia and dysphagia grade ≥ II were reduced by 3.5 ± 1.7 and 3.9 ± 2.8 percentage point [mean ± SD] (p< 10-5for both). In online re-optimization, adjustment of spot configuration followed by spot-intensity re-optimization took 3.4 min on average.Significance. The fast online re-optimization strategy always prevented substantial losses of target coverage caused by day-to-day anatomical variations, as opposed to the clinical trigger-based offline re-planning schedule. On top of this, online re-optimization could be performed with smaller setup robustness settings, contributing to improved organs-at-risk sparing.

SISS-MCO: large scale sparsity-induced spot selection for fast and fully-automated robust multi-criteria optimisation of proton plans

Objective.Intensity modulated proton therapy (IMPT) is an emerging treatment modality for cancer. However, treatment planning for IMPT is labour-intensive and time-consuming. We have developed a novel approach for multi-criteria optimisation (MCO) of robust IMPT plans (SISS-MCO) that is fully automated and fast, and we compare it for head and neck, cervix, and prostate tumours to a previously published method for automated robust MCO (IPBR-MCO, van de Water 2013).Approach.In both auto-planning approaches, the applied automated MCO of spot weights was performed with wish-list driven prioritised optimisation (Breedveld 2012). In SISS-MCO, spot weight MCO was applied once for every patient after sparsity-induced spot selection (SISS) for pre-selection of the most relevant spots from a large input set of candidate spots. IPBR-MCO had several iterations of spot re-sampling, each followed by MCO of the weights of the current spots.Main results.Compared to the published IPBR-MCO, the novel SISS-MCO resulted in similar or slightly superior plan quality. Optimisation times were reduced by a factor of 6 i.e. from 287 to 47 min. Numbers of spots and energy layers in the final plans were similar.Significance.The novel SISS-MCO automatically generated high-quality robust IMPT plans. Compared to a published algorithm for automated robust IMPT planning, optimisation times were reduced on average by a factor of 6. Moreover, SISS-MCO is a large scale approach; this enables optimisation of more complex wish-lists, and novel research opportunities in proton therapy.

Dosimetrically motivated beam-angle optimization for non-coplanar arc radiotherapy with and without dynamic collimator rotation

BACKGROUND

Non-coplanar techniques have shown to improve the achievable dose distribution compared to standard coplanar techniques for multiple treatment sites but finding optimal beam directions is challenging. Dynamic collimator trajectory radiotherapy (colli-DTRT) is a new intensity modulated radiotherapy technique that uses non-coplanar partial arcs and dynamic collimator rotation.

PURPOSE

To solve the beam angle optimization (BAO) problem for colli-DTRT and non-coplanar VMAT (NC-VMAT) by determining the table-angle and the gantry-angle ranges of the partial arcs through iterative 4π fluence map optimization (FMO) and beam direction elimination.

METHODS

BAO considers all available beam directions sampled on a gantry-table map with the collimator angle aligned to the superior-inferior axis (colli-DTRT) or static (NC-VMAT). First, FMO is performed, and beam directions are scored based on their contributions to the objective function. The map is thresholded to remove the least contributing beam directions, and arc candidates are formed by adjacent beam directions with the same table angle. Next, FMO and arc candidate trimming, based on objective function penalty score, is performed iteratively until a desired total gantry angle range is reached. Direct aperture optimization on the final set of colli-DTRT or NC-VMAT arcs generates deliverable plans. colli-DTRT and NC-VMAT plans were created for seven clinically-motivated cases with targets in the head and neck (two cases), brain, esophagus, lung, breast, and prostate. colli-DTRT and NC-VMAT were compared to coplanar VMAT plans as well as to class-solution non-coplanar VMAT plans for the brain and head and neck cases. Dosimetric validation was performed for one colli-DTRT (head and neck) and one NC-VMAT (breast) plan using film measurements.

RESULTS

Target coverage and conformity was similar for all techniques. colli-DTRT and NC-VMAT plans had improved dosimetric performance compared to coplanar VMAT for all treatment sites except prostate where all techniques were equivalent. For the head and neck and brain cases, mean dose reduction-in percentage of the prescription dose-to parallel organs was on average 0.7% (colli-DTRT), 0.8% (NC-VMAT) and 0.4% (class-solution) compared to VMAT. The reduction in D2% for the serial organs was on average 1.7% (colli-DTRT), 2.0% (NC-VMAT) and 0.9% (class-solution). For the esophagus, lung, and breast cases, mean dose reduction to parallel organs was on average 0.2% (colli-DTRT) and 0.3% (NC-VMAT) compared to VMAT. The reduction in D2% for the serial organs was on average 1.3% (colli-DTRT) and 0.9% (NC-VMAT). Estimated delivery times for colli-DTRT and NC-VMAT were below 4 min for a full gantry angle range of 720°, including transitions between arcs, except for the brain case where multiple arcs covered the whole table angle range. These times are in the same order as the class-solution for the head and neck and brain cases. Total optimization times were 25%-107% longer for colli-DTRT, including BAO, compared to VMAT.

CONCLUSIONS

We successfully developed dosimetrically motivated BAO for colli-DTRT and NC-VMAT treatment planning. colli-DTRT and NC-VMAT are applicable to multiple treatment sites, including body sites, with beneficial or equivalent dosimetric performances compared to coplanar VMAT and reasonable delivery times.

Script-based automatic radiotherapy planning for cervical cancer

BACKGROUND

This study aimed to develop fully automated script-based radiotherapy treatment plans for cervical cancer patients, and evaluate them against clinically accepted plans, as validation before clinical implementation.

MATERIAL AND METHODS

In this retrospective planning study, treatment plans for 25 locally advanced cervical cancer (LACC) patients with up to three dose levels were included. Fully automated plans were created using an in-house developed Python script in RayStation, and compared to clinically accepted manually made plans. Quantitatively, relevant dose statistics were compared, and average dose volume histograms (DVHs) were analyzed. Qualitatively, a blinded plan comparison was conducted between the clinical and automatic plans. The accuracy of treatment plan delivery was verified with the Delta4 Phantom+.

RESULTS

The quantitative evaluation showed that target coverage was acceptable for all the automatic and clinical plans. The automatic plans were significantly more conformal than the clinical plans; median of 1.03 vs. 1.12. Mean doses to almost all organs at risk (OARs) were reduced in the automatic plans, with a median reduction of between 0.6 Gy and 1.9 Gy. In the blinded plan comparison, the automatic plans were the preferred plans or of equal quality as the clinical plans in 99% of the cases. In addition, plan delivery was excellent, with a mean gamma passing rate of 99.8%. Complete script-based plans were generated in 30-45 min; about four to ten times faster than manually made plans.

CONCLUSION

The automatic plans had acceptable target coverage, lower doses to almost all OARs, more conformal dose distributions, and were predominantly preferred by the clinicians. Based on these results, our institution has implemented the script for clinical use.

Cochlear-optimized treatment planning in photon and proton radiosurgery for vestibular schwannoma patients

Objective

To investigate the potential to reduce the cochlear dose with robotic photon radiosurgery or intensity-modulated proton therapy planning for vestibular schwannomas.

Materials and Methods

Clinically delivered photon radiosurgery treatment plans were compared to five cochlear-optimized plans: one photon and four proton plans (total of 120). A 1x12 Gy dose was prescribed. Photon plans were generated with Precision (Cyberknife, Accuray) with no PTV margin for set-up errors. Proton plans were generated using an in-house automated multi-criterial planning system with three or nine-beam arrangements, and applying 0 or 3 mm robustness for set-up errors during plan optimization and evaluation (and 3 % range robustness). The sample size was calculated based on a reduction of cochlear Dmean > 1.5 Gy(RBE) from the clinical plans, and resulted in 24 patients.

Results

Compared to the clinical photon plans, a reduction of cochlear Dmean > 1.5 Gy(RBE) could be achieved in 11/24 cochlear-optimized photon plans, 4/24 and 6/24 cochlear-optimized proton plans without set-up robustness for three and nine-beam arrangement, respectively, and in 0/24 proton plans with set-up robustness. The cochlea could best be spared in cases with a distance between tumor and cochlea. Using nine proton beams resulted in a reduced dose to most organs at risk.

Conclusion

Cochlear dose reduction is possible in vestibular schwannoma radiosurgery while maintaining tumor coverage, especially when the tumor is not adjacent to the cochlea. With current set-up robustness, proton therapy is capable of providing lower dose to organs at risk located distant to the tumor, but not for organs adjacent to it. Consequently, photon plans provided better cochlear sparing than proton plans.

Improving knowledge-based treatment planning for lung cancer radiotherapy with automatic multi-criteria optimized training plans

BACKGROUND

Knowledge-based planning (KBP) is a method for automated radiotherapy treatment planning where appropriate optimization objectives for new patients are predicted based on a library of training plans. KBP can save time and improve organ at-risk sparing and inter-patient consistency compared to manual planning, but its performance depends on the quality of the training plans. We used another system for automated planning, which generates multi-criteria optimized (MCO) plans based on a wish list, to create training plans for the KBP model, to allow seamless integration of knowledge from a new system into clinical routine. Model performance was compared for KBP models trained with manually created and automatic MCO treatment plans.

MATERIAL AND METHODS

Two RapidPlan models with the same 30 locally advanced non-small cell lung cancer patients included were created, one containing manually created clinical plans (RP_CLIN) and one containing fully automatic multi-criteria optimized plans (RP_MCO). For 15 validation patients, model performance was compared in terms of dose-volume parameters and normal tissue complication probabilities, and an oncologist performed a blind comparison of the clinical (CLIN), RP_CLIN, and RP_MCO plans.

RESULTS

The heart and esophagus doses were lower for RP_MCO compared to RP_CLIN, resulting in an average reduction in the risk of 2-year mortality by 0.9 percentage points and the risk of acute esophageal toxicity by 1.6 percentage points with RP_MCO. The oncologist preferred the RP_MCO plan for 8 patients and the CLIN plan for 7 patients, while the RP_CLIN plan was not preferred for any patients.

CONCLUSION

RP_MCO improved OAR sparing compared to RP_CLIN and was selected for implementation in the clinic. Training a KBP model with clinical plans may lead to suboptimal output plans, and making an extra effort to optimize the library plans in the KBP model creation phase can improve the plan quality for many future patients.

Influence of optimisation parameters on directly deliverable Pareto fronts explored for prostate cancer

PURPOSE

In inverse radiotherapy treatment planning, the Pareto front is the set of optimal solutions to the multi-criteria problem of adequately irradiating the planning target volume (PTV) while reducing dose to organs at risk (OAR). The Pareto front depends on the chosen optimisation parameters whose influence (clinically relevant versus not clinically relevant) is investigated in this paper.

METHODS

Thirty-one prostate cancer patients treated at our clinic were randomly selected. We developed an in-house Python script that controlled the commercial treatment planning system RayStation to calculate directly deliverable Pareto fronts. We calculated reference Pareto fronts for a given set of objective functions, varying the PTV coverage and the mean dose of the primary OAR (rectum) and fixing the mean doses of the secondary OARs (bladder and femoral heads). We calculated the fronts for different sets of objective functions and different mean doses to secondary OARs. We compared all fronts using a specific metric (clinical distance measure).

RESULTS

The in-house script was validated for directly deliverable Pareto front calculations in two and three dimensions. The Pareto fronts depended on the choice of objective functions and fixed mean doses to secondary OARs, whereas the parameters most influencing the front and leading to clinically relevant differences were the dose gradient around the PTV, the weight of the PTV objective function, and the bladder mean dose.

CONCLUSIONS

Our study suggests that for multi-criteria optimisation of prostate treatments using external therapy, dose gradient around the PTV and bladder mean dose are the most influencial parameters.

Updating approach for lexicographic optimization-based planning to improve cervical cancer plan quality

BACKGROUND

To investigate the capability of a not-yet commercially available fully automated lexicographic optimization (LO) planning algorithm, called mCycle (Elekta AB, Stockholm, Sweden), to further improve the plan quality of an already-validated Wish List (WL) pushing on the organs-at-risk (OAR) sparing without compromising target coverage and plan delivery accuracy.

MATERIAL AND METHODS

Twenty-four mono-institutional consecutive cervical cancer Volumetric-Modulated Arc Therapy (VMAT) plans delivered between November 2019 and April 2022 (50 Gy/25 fractions) have been retrospectively selected. In mCycle the LO planning algorithm was combined with the a-priori multi-criterial optimization (MCO). Two versions of WL have been defined to reproduce manual plans (WL01), and to improve the OAR sparing without affecting minimum target coverage and plan delivery accuracy (WL02). Robust WLs have been tuned using a subset of 4 randomly selected patients. The remaining plans have been automatically re-planned by using the designed WLs. Manual plans (MP) and mCycle plans (mCP01 and mCP02) were compared in terms of dose distributions, complexity, delivery accuracy, and clinical acceptability. Two senior physicians independently performed a blind clinical evaluation, ranking the three competing plans. Furthermore, a previous defined global quality index has been used to gather into a single score the plan quality evaluation.

RESULTS

The WL tweaking requests 5 and 3 working days for the WL01 and the WL02, respectively. The re-planning took in both cases 3 working days. mCP01 best performed in terms of target coverage (PTV V95% (%): MP 98.0 [95.6-99.3], mCP01 99.2 [89.7-99.9], mCP02 96.9 [89.4-99.5]), while mCP02 showed a large OAR sparing improvement, especially in the rectum parameters (e.g., Rectum D50% (Gy): MP 41.7 [30.2-47.0], mCP01 40.3 [31.4-45.8], mCP02 32.6 [26.9-42.6]). An increase in plan complexity has been registered in mCPs without affecting plan delivery accuracy. In the blind comparisons, all automated plans were considered clinically acceptable, and mCPs were preferred over MP in 90% of cases. Globally, automated plans registered a plan quality score at least comparable to MP.

CONCLUSIONS

This study showed the flexibility of the Lexicographic approach in creating more demanding Wish Lists able to potentially minimize toxicities in RT plans.

A unified path seeking algorithm for IMRT and IMPT beam orientation optimization

Objective. Fully automated beam orientation optimization (BOO) for intensity-modulated radiotherapy and intensity modulated proton therapy (IMPT) is gaining interest, since achieving optimal plan quality for an unknown number of fixed beam arrangements is tedious. Fast group sparsity-based optimization methods have been proposed to find the optimal orientation, but manual tuning is required to eliminate the exact number of beams from a large candidate set. Here, we introduce a fast, automated gradient descent-based path-seeking algorithm (PathGD), which performs fluence map optimization for sequentially added beams, to visualize the dosimetric benefit of one added field at a time.Approach. Several configurations of 2-4 proton and 5-15 photon beams were selected for three head-and-neck patients using PathGD, which was compared to group sparsity-regularized BOO solved with the fast iterative shrinkage-thresholding algorithm (GS-FISTA), and manually selected IMPT beams or one coplanar photon VMAT arc (MAN). Once beams were chosen, all plans were compared on computational efficiency, dosimetry, and for proton plans, robustness.Main results. With each added proton beam, Clinical Target Volume (CTV) and organs at risk (OAR) dosimetric cost improved on average across plans by [1.1%, 13.6%], and for photons, [0.6%, 2.0%]. Comparing algorithms, beam selection for PathGD was faster than GS-FISTA on average by 35%, and PathGD matched the CTV coverage of GS-FISTA plans while reducing OAR mean and maximum dose in all structures by an average of 13.6%. PathGD was able to improve CTV [Dmax, D95%] by [2.6%, 5.2%] and reduced worst-case [max, mean] dose in OARs by [11.1%, 13.1%].Significance. The benefit of a path-seeking algorithm is the beam-by-beam analysis of dosimetric cost. PathGD was shown to be most efficient and dosimetrically desirable amongst group sparsity and manual BOO methods, and highlights the sensitivity of beam addition for IMPT in particular.

Clinical benefit of range uncertainty reduction in proton treatment planning based on dual-energy CT for neuro-oncological patients

OBJECTIVE

Several studies have shown that dual-energy CT (DECT) can lead to improved accuracy for proton range estimation. This study investigated the clinical benefit of reduced range uncertainty, enabled by DECT, in robust optimisation for neuro-oncological patients.

METHODS

DECT scans for 27 neuro-oncological patients were included. Commercial software was applied to create stopping-power ratio (SPR) maps based on the DECT scan. Two plans were robustly optimised on the SPR map, keeping the beam and plan settings identical to the clinical plan. One plan was robustly optimised and evaluated with a range uncertainty of 3% (as used clinically; denoted 3%-plan); the second plan applied a range uncertainty of 2% (2%-plan). Both plans were clinical acceptable and optimal. The dose-volume histogram parameters were compared between the two plans. Two experienced neuro-radiation oncologists determined the relevant dose difference for each organ-at-risk (OAR). Moreover, the OAR toxicity levels were assessed.

RESULTS

For 24 patients, a dose reduction >0.5/1 Gy (relevant dose difference depending on the OAR) was seen in one or more OARs for the 2%-plan; e.g. for brainstem D0.03cc in 10 patients, and hippocampus D40% in 6 patients. Furthermore, 12 patients had a reduction in toxicity level for one or two OARs, showing a clear benefit for the patient.

CONCLUSION

Robust optimisation with reduced range uncertainty allows for reduction of OAR toxicity, providing a rationale for clinical implementation. Based on these results, we have clinically introduced DECT-based proton treatment planning for neuro-oncological patients, accompanied with a reduced range uncertainty of 2%.

ADVANCES IN KNOWLEDGE

This study shows the clinical benefit of range uncertainty reduction from 3% to 2% in robustly optimised proton plans. A dose reduction to one or more OARs was seen for 89% of the patients, and 44% of the patients had an expected toxicity level decrease.

Novel in-house knowledge-based automated planning system for lung cancer treated with intensity-modulated radiotherapy

PURPOSE

The goal of this study was to propose a knowledge-based planning system which could automatically design plans for lung cancer patients treated with intensity-modulated radiotherapy (IMRT).

METHODS AND MATERIALS

From May 2018 to June 2020, 612 IMRT treatment plans of lung cancer patients were retrospectively selected to construct a planning database. Knowledge-based planning (KBP) architecture named αDiar was proposed in this study. It consisted of two parts separated by a firewall. One was the in-hospital workstation, and the other was the search engine in the cloud. Based on our previous study, A‑Net in the in-hospital workstation was used to generate predicted virtual dose images. A search engine including a three-dimensional convolutional neural network (3D CNN) was constructed to derive the feature vectors of dose images. By comparing the similarity of the features between virtual dose images and the clinical dose images in the database, the most similar feature was found. The optimization parameters (OPs) of the treatment plan corresponding to the most similar feature were assigned to the new plan, and the design of a new treatment plan was automatically completed. After αDiar was developed, we performed two studies. The first retrospective study was conducted to validate whether this architecture was qualified for clinical practice and involved 96 patients. The second comparative study was performed to investigate whether αDiar could assist dosimetrists in improving the quality of planning for the patients. Two dosimetrists were involved and designed plans for only one trial with and without αDiar; 26 patients were involved in this study.

RESULTS

The first study showed that about 54% (52/96) of the automatically generated plans would achieve the dosimetric constraints of the Radiation Therapy Oncology Group (RTOG) and about 93% (89/96) of the automatically generated plans would achieve the dosimetric constraints of the National Comprehensive Cancer Network (NCCN). The second study showed that the quality of treatment planning designed by junior dosimetrists was improved with the help of αDiar.

CONCLUSIONS

Our results showed that αDiar was an effective tool to improve planning quality. Over half of the patients' plans could be designed automatically. For the remaining patients, although the automatically designed plans did not fully meet the clinical requirements, their quality was also better than that of manual plans.

Per-fraction planning to enhance optimization degrees of freedom compared to the conventional single-plan approach

Objective. In conventional radiotherapy, a single treatment plan is generated pre-treatment, and delivered in daily fractions. In this study, we propose to generate different treatment plans for all fractions ('Per-fraction' planning) to reduce cumulative organs at risk (OAR) doses. Per-fraction planning was compared to the 'Conventional' single-plan approach for non-coplanar 4 × 9.5 Gy prostate stereotactic body radiation therapy (SBRT).Approach. An in-house application for fully automated, non-coplanar multi-criterial treatment planning with integrated beam angle and fluence optimization was used for plan generations. For the Conventional approach, a single 12-beam non-coplanar IMRT plan with individualized beam angles was generated for each of the 20 included patients. In Per-fraction planning, four fraction plans were generated for each patient. For each fraction, a different set of patient-specific 12-beam configurations could be automatically selected. Per-fraction plans were sequentially generated by adding dose to already generated fraction plan(s). For each fraction, the cumulative- and fraction dose were simultaneously optimized, allowing some minor constraint violations in fraction doses, but not in cumulative.Main results. In the Per-fraction approach, on average 32.9 ± 3.1 [29;39] unique beams per patient were used. PTV doses in the separate Per-fraction plans were acceptable and highly similar to those in Conventional plans, while also fulfilling all OAR hard constraints. When comparing total cumulative doses, Per-fraction planning showed improved bladder sparing for all patients with reductions in Dmean of 22.6% (p= 0.0001) and in D1cc of 2.0% (p= 0.0001), reductions in patient volumes receiving 30% and 50% of the prescribed dose of 54.7% and 6.3%, respectively, and a 3.1% lower rectum Dmean (p= 0.007). Rectum D1cc was 4.1% higher (p= 0.0001) and Urethra dose was similar.Significance. In this proof-of-concept paper, Per-fraction planning resulted in several dose improvements in healthy tissues compared to the Conventional single-plan approach, for similar PTV dose. By keeping the number of beams per fraction the same as in Conventional planning, reported dosimetric improvements could be obtained without increase in fraction durations. Further research is needed to explore the full potential of the Per-fraction planning approach.

The impact of bone marrow sparing on organs at risk dose for cervical cancer: a Pareto front analysis

Background and purpose

To quantify the increase in bladder and rectum dose of a bone marrow sparing (BMS) VMAT strategy for primary treatment of locally advanced cervical cancer (LACC).

Materials and methods

Twenty patients with stage IB-IVA cervical cancer were selected for this study. The whole Pelvic Bones (PB) was taken as substitute for bone marrow. For every patient, Pareto-optimal plans were generated to explore the trade-off between rectum, bladder, and PB mean dose. The PB mean dose was decreased in steps of 1 Gy. For each step, the increase in rectum and bladder mean dose was quantified. The increase in mean dose of other OAR compared to no BMS was constrained to 1 Gy.

Results

In total, 931 plans of 19 evaluable patients were analyzed. The average [range] mean dose of PB without BMS was 22.8 [20.7-26.2] Gy. When maximum BMS was applied, the average reduction in mean PB dose was 5.4 [3.0-6.8] Gy resulting in an average mean PB dose of 17.5 [15.8-19.8] Gy. For <1 Gy increase in both the bladder and the rectum mean dose, the PB mean dose could be decreased by >2 Gy, >3 Gy, >4 Gy, and >5 Gy for 19/19, 13/19, 5/19, and 1/19 patients, respectively.

Conclusion

Based on the comprehensive three-dimensional Pareto front analysis, we conclude that 2-5 Gy BMS can be implemented without a clinically relevant increase in mean dose to other OAR. If BMS is too dominant, it results in a large increase in mean dose to other OAR. Therefore, we recommend implementing moderate BMS for the treatment of LACC patients with VMAT.

Online adaptive planning methods for intensity-modulated radiotherapy

Online adaptive radiation therapy aims at adapting a patient's treatment plan to their current anatomy to account for inter-fraction variations before daily treatment delivery. As this process needs to be accomplished while the patient is immobilized on the treatment couch, it requires time-efficient adaptive planning methods to generate a quality daily treatment plan rapidly. The conventional planning methods do not meet the time requirement of online adaptive radiation therapy because they often involve excessive human intervention, significantly prolonging the planning phase. This article reviews the planning strategies employed by current commercial online adaptive radiation therapy systems, research on online adaptive planning, and artificial intelligence's potential application to online adaptive planning.

Multi-isocenter VMAT craniospinal irradiation using feasibility dose-volume histogram-guided auto-planning technique

This study aims to propose a novel treatment planning methodology for multi-isocenter volumetric modulated arc therapy (VMAT) craniospinal irradiation (CSI) using the special feasibility dose-volume histogram (FDVH)-guided auto-planning (AP) technique. Three different multi-isocenter VMAT -CSI plans were created, including manually based plans (MUPs), conventional AP plans (CAPs) and FDVH-guided AP plans (FAPs). The CAPs and FAPs were specially designed by combining multi-isocenter VMAT and AP techniques in the Pinnacle treatment planning system. Specially, the personalized optimization parameters for FAPs were generated using the FDVH function implemented in PlanIQ software, which provides the ideal organs at risk (OARs) sparing for the specific anatomical geometry based on the valuable assumption of the dose fall-off. Compared to MUPs, CAPs and FAPs significantly reduced the dose for most of the OARs. FAPs achieved the best homogeneity index (0.092 ± 0.013) and conformity index (0.980 ± 0.011), while CAPs were slightly inferior to the FAPs but superior to the MUPs. As opposed to MUPs, FAPs delivered a lower dose to OARs, whereas the difference between FAPs and CAPs was not statistically significant except for the optic chiasm and inner ear_L. The two AP approaches had similar MUs, which were significantly lower than the MUPs. The planning time of FAPs (145.00 ± 10.25 min) was slightly lower than that of CAPs (149.83 ± 14.37 min) and was substantially lower than that of MUPs (157.92 ± 16.11 min) with P < 0.0167. Overall, introducing the multi-isocenter AP technique into VMAT-CSI yielded positive outcomes and may play an important role in clinical CSI planning in the future.

Learning image representations for content-based image retrieval of radiotherapy treatment plans

Objective.In this work, we propose a content-based image retrieval (CBIR) method for retrieving dose distributions of previously planned patients based on anatomical similarity. Retrieved dose distributions from this method can be incorporated into automated treatment planning workflows in order to streamline the iterative planning process. As CBIR has not yet been applied to treatment planning, our work seeks to understand which current machine learning models are most viable in this context.Approach.Our proposed CBIR method trains a representation model that produces latent space embeddings of a patient's anatomical information. The latent space embeddings of new patients are then compared against those of previous patients in a database for image retrieval of dose distributions. All source code for this project is available on github.Main results.The retrieval performance of various CBIR methods is evaluated on a dataset consisting of both publicly available image sets and clinical image sets from our institution. This study compares various encoding methods, ranging from simple autoencoders to more recent Siamese networks like SimSiam, and the best performance was observed for the multitask Siamese network.Significance.Our current results demonstrate that excellent image retrieval performance can be obtained through slight changes to previously developed Siamese networks. We hope to integrate CBIR into automated planning workflow in future works.

Built-in wavelet-induced smoothness to reduce plan complexity in intensity modulated radiation therapy (IMRT)

Objective.Reducing plan complexity in intensity modulated radiation therapy (IMRT) to ensure dosimetric accuracy and delivery efficiency of the radiation treatment plans. We propose a novel approach by representing the beamlet intensities using an incomplete wavelet basis that explicitly excludes fluctuating intensity maps from the decision space (explicit hard constraint). This technique provides a built-in wavelet-induced smoothness that improves both dosimetric plan quality and delivery efficiency.Approach.The beamlet intensity maps need to be especially smooth in the leaf travel direction (referred to as theX-direction). We treat the intensity map of each beam as a 2D image and represent it using the wavelets corresponding to low-frequency changes in theX-direction (i.e. approximation and horizontal). The absence of wavelets corresponding to high-frequency changes (i.e. vertical and diagonal) induces built-in smoothness. We still utilize a regularization term in the objective function to promote smoothness in theY-direction (perpendicular to theX-direction) and further possible smoothness in theX-direction. This technique has been tested on three patient cases of different disease sites (paraspinal, lung, prostate) and all final evaluations and comparisons have been performed on an FDA-approved commercial treatment planning system (Varian Eclipse^TM).Main results.Wavelet-induced smoothness reduced monitor units by about 10%, 45%, and 14% for paraspinal, lung, and prostate cases, respectively. It also improved organ at risk sparing, especially on the complex paraspinal case where it resulted in about 7%, 13%, and 14% less mean dose to esophagus, lung, and cord, respectively. Moreover, built-in wavelet-induced smoothness desensitizes the results to changing the weight associated to the regularization term, and thereby mitigates the weight fine-tuning difficulty.Significance.Fluence smoothness is often achieved by smoothing the beamlet intensity maps using a proper regularization term in the objective function aiming at disincentivizing fluctuation in the beamlet intensities (implicit soft constraint). This work reports a novel application of wavelets in imposing an explicit smoothness hard constraint in the search space using an incomplete wavelet basis. This idea has been successfully applied to exclude complex and clinically irrelevant radiation plans from the search space. The code and pertained models along with a sample dataset are released on our LowDimRT GitHub (https://github.com/PortPy-Project/LowDimRT).

Pencil-beam Delivery Pattern Optimization Increases Dose Rate for Stereotactic FLASH Proton Therapy

PURPOSE

FLASH dose rates >40 Gy/s are readily available in proton therapy (PT) with cyclotron-accelerated beams and pencil-beam scanning (PBS). The PBS delivery pattern will affect the local dose rate, as quantified by the PBS dose rate (PBS-DR), and therefore needs to be accounted for in FLASH-PT with PBS, but it is not yet clear how. Our aim was to optimize patient-specific scan patterns for stereotactic FLASH-PT of early-stage lung cancer and lung metastases, maximizing the volume irradiated with PBS-DR >40 Gy/s of the organs at risk voxels irradiated to >8 Gy (FLASH coverage).

METHODS AND MATERIALS

Plans to 54 Gy/3 fractions with 3 equiangular coplanar 244 MeV proton shoot-through transmission beams for 20 patients were optimized with in-house developed software. Planning target volume-based planning with a 5 mm margin was used. Planning target volume ranged from 4.4 to 84 cc. Scan-pattern optimization was performed with a Genetic Algorithm, run in parallel for 20 independent populations (islands). Mapped crossover, inversion, swap, and shift operators were applied to achieve (local) optimality on each island, with migration between them for global optimality. The cost function was chosen to maximize the FLASH coverage per beam at >8 Gy, >40 Gy/s, and 40 nA beam current. The optimized patterns were evaluated on FLASH coverage, PBS-DR distribution, and population PBS-DR-volume histograms, compared with standard line-by-line scanning. Robustness against beam current variation was investigated.

RESULTS

The optimized patterns have a snowflake-like structure, combined with outward swirling for larger targets. A population median FLASH coverage of 29.0% was obtained for optimized patterns compared with 6.9% for standard patterns, illustrating a significant increase in FLASH coverage for optimized patterns. For beam current variations of 5 nA, FLASH coverage varied between -6.1%-point and 2.2%-point for optimized patterns.

CONCLUSIONS

Significant improvements on the PBS-DR and, hence, on FLASH coverage and potential healthy-tissue sparing are obtained by sequential scan-pattern optimization. The optimizer is flexible and may be further fine-tuned, based on the exact conditions for FLASH.

Patient anatomy-specific trade-offs between sub-clinical disease coverage and normal tissue dose reduction in head-and-neck cancer

PURPOSE

Risk of subclinical disease decreases with increasing distance from the GTV in head- and-neck squamous cell carcinoma (HNSCC). Depending on individual patient anatomy, OAR sparing could be improved by reducing target coverage in regions with low risk of subclinical spread. Using automated multi-criteria optimization, we investigate patient-specific optimal trade-offs between target periphery coverage and OAR sparing.

METHODS

VMAT plans for 39 HNSCC patients were retrospectively created following our clinical three-target-level protocol: high-risk (PTV1), intermediate-risk (PTV2, 5 mm expansion from PTV1), and elective (PTV3). A baseline plan fulfilling clinical constraints (D 99 % ≥95 % for all PTVs) was compared to three plans with reduced PTV2 coverage (goals: PTV2 D 99 % ≥90 % or 85 %, or no PTV2) at the outer edge of PTV2. Plans were compared on PTV D 99 %, OAR D mean, and NTCP (xerostomia/dysphagia).

RESULTS

Trade-offs between PTV2 coverage and OAR doses varied considerably between patients. For plans with PTV2 D 99 % -goal 90 %, median PTV2 D 99 % was 91.5 % resulting in xerostomia (≥grade 4) and dysphagia (≥grade 2) NTCP decrease of median [maximum] 1.9 % [5.3 %] and 1.1 % [4.1 %], respectively, compared to nominal PTV2 D 99 % -goal 95 %. For PTV2 D 99 % -goal 85 % median PTV D 99 % was 87 % with NTCP improvements of 4.6 % [9.9 %] and 1.5 % [5.4 %]. For no-margin plans, PTV2 D 99 % decreased to 83.3 % with NTCP reductions of 5.1 % [10.2 %] and 1.4 % [6.1 %].

CONCLUSION

Clinically relevant, patient-specific reductions in OARs and NTCP were observed at limited cost in target under-coverage at the outermost PTV edge. Given the observed inter-patient variations, individual evaluation is warranted to determine whether trade- offs would benefit a specific patient.

Feasibility and dosimetric evaluation of single- and multi-isocentre stereotactic body radiation therapy for multiple liver metastases

Objectives

Single-isocentre volumetric-modulated arc therapy (VMAT) stereotactic body radiation therapy (SBRT) improves treatment efficiency and patient compliance for patients with multiple liver metastases (MLM). However, the potential increase in dose spillage to normal liver tissue using a single-isocentre technique has not yet been studied. We comprehensively evaluated the quality of single- and multi-isocentre VMAT-SBRT for MLM and propose a RapidPlan-based automatic planning (AP) approach for MLM SBRT.

Methods

A total of 30 patients with MLM (two or three lesions) were selected for this retrospective study. We manually replanned all patients treated with MLM SBRT by using the single-isocentre (MUS) and multi-isocentre (MUM) techniques. Then, we randomly selected 20 MUS and MUM plans for training to generate the single-isocentre RapidPlan model (RPS) and the multi-isocentre RapidPlan model (RPM). Finally, we used data from the remaining 10 patients to validate RPS and RPM.

Results

Compared with MUS, MUM reduced the mean dose delivered to the right kidney by 0.3 Gy. The mean liver dose (MLD) was 2.3 Gy higher for MUS compared with MUM. However, the monitor units, delivery time, and V20Gy of normal liver (liver-gross tumour volume) for MUM were significantly higher than for MUS. Based on validation, RPS and RPM slightly improved the MLD, V20Gy, normal tissue complications, and dose sparing to the right and left kidneys and spinal cord compared with manual plans (MUS vs RPS and MUM vs RPM), but RPS and RPM significantly increased monitor units and delivery time.

Conclusions

The single-isocentre VMAT-SBRT approach could be used for MLM to reduce treatment time and patient comfort at the cost of a small increase in the MLD. Compared with the manual plans, RapidPlan-based plans, especially RPS, have slightly improved quality.

Automation of pencil beam scanning proton treatment planning for intracranial tumours

PURPOSE

To evaluate the feasibility of comprehensive automation of an intra-cranial proton treatment planning.

MATERIALS AND METHODS

Class solution (CS) beam configuration selection allows the user to identify predefined beam configuration based on target localization; automatic CS (aCS) will then explore all the possible CS beam geometries. Ten patients, already used for the evaluation of the automatic selection of the beam configuration, have been also employed to training an algorithm based on the computation of a benchmark dose exploit automatic general planning solution (GPS) optimization with a wish list approach for the planning optimization. An independent cohort of ten patients has been then used for the evaluation step between the clinical and the GPS plan in terms of dosimetric quality of plans and the time needed to generate a plan.

RESULTS

The definition of a beam configuration requires on average 22 min (range 9-29 min). The average time for GPS plan generation is 18 min (range 7-26 min). Median dose differences (GPS-Manual) for each OAR constraints are: brainstem -1.60 Gy, left cochlea -1.22 Gy, right cochlea -1.42 Gy, left eye 0.55 Gy, right eye -2.33 Gy, optic chiasm -1.87 Gy, left optic nerve -4.45 Gy, right optic nerve -2.48 Gy and optic tract -0.31 Gy. Dosimetric CS and aCS plan evaluation shows a slightly worsening of the OARs values except for the optic tract and optic chiasm for both CS and aCS, where better results have been observed.

CONCLUSION

This study has shown the feasibility and implementation of the automatic planning system for intracranial tumors. The method developed in this work is ready to be implemented in a clinical workflow.

Comprehensive dosimetric and clinical evaluation of lexicographic optimization-based planning for cervical cancer

AIM

In this study, a not yet commercially available fully-automated lexicographic optimization (LO) planning algorithm, called mCycle (Elekta AB, Stockholm, Sweden), was validated for cervical cancer.

MATERIAL AND METHODS

Twenty-four mono-institutional consecutive treatment plans (50 Gy/25 fx) delivered between November 2019 and April 2022 were retrospectively selected. The automatic re-planning was performed by mCycle, implemented in the Monaco TPS research version (v5.59.13), in which the LO and Multicriterial Optimization (MCO) are coupled with Monte Carlo calculation. mCycle optimization follows an a priori assigned priority list, the so-called Wish List (WL), representing a dialogue between the radiation oncologist and the planner, setting hard constraints and following objectives. The WL was tuned on a patient subset according to the institution's clinical protocol to obtain an optimal plan in a single optimization. This robust WL was then used to automatically re-plan the remaining patients. Manual plans (MP) and mCycle plans (mCP) were compared in terms of dose distributions, complexity (modulation complexity score, MCS), and delivery accuracy (perpendicular diode matrices, gamma analysis-passing ratio, PR). Their clinical acceptability was assessed through the blind choice of two radiation oncologists. Finally, a global quality score index (SI) was defined to gather into a single number the plan evaluation process.

RESULTS

The WL tuning requested four patients. The 20 automated re-planning tasks took three working days. The median optimization and calculation time can be estimated at 4 h and just over 1 h per MP and mCP, respectively. The dose comparison showed a comparable organ-at-risk spare. The planning target volume coverage increased (V95%: MP 98.0% [95.6-99.3]; mCP 99.2%[89.7-99.9], p >0.05). A significant increase has been registered in MCS (MP 0.29 [0.24-0.34]; mCP 0.26 [0.23-0.30], p <0.05) without affecting delivery accuracy (PR (3%/3mm): MP 97.0% [92.7-99.2]; mCP 97.1% [95.0-98.6], p >0.05). In the blind choice, all mCP results were clinically acceptable and chosen over MP in more than 75% of cases. The median SI score was 0.69 [0.41-0.84] and 0.73 [0.51-0.82] for MP and mCP, respectively (p >0.05).

CONCLUSIONS

mCycle plans were comparable to clinical manual plans, more complex but accurately deliverable and registering a similar SI. Automated plans outperformed manual plans in blinded clinical choice.

An online adaptive plan library approach for intensity modulated proton therapy for head and neck cancer

BACKGROUND AND PURPOSE

In intensity modulated proton therapy (IMPT), the impact of setup errors and anatomical changes is commonly mitigated by robust optimization with population-based setup robustness (SR) settings and offline replanning. In this study we propose and evaluate an alternative approach based on daily plan selection from patient-specific pre-treatment established plan libraries (PLs). Clinical implementation of the PL strategy would be rather straightforward compared to daily online re-planning.

MATERIALS AND METHODS

For 15 head-and-neck cancer patients, the planning CT was used to generate a PL with 5 plans, robustly optimized for increasing SR: 0, 1, 2, 3, 5 mm, and 3% range robustness. Repeat CTs (rCTs) and realistic setup and range uncertainty distributions were used for simulation of treatment courses for the PL approach, treatments with fixed SR (fSR₃) and a trigger-based offline adaptive schedule for 3 mm SR (fSR₃OfA). Daily plan selection in the PL approach was based only on recomputed dose to the CTV on the rCT.

RESULTS

Compared to using fSR₃ and fSR₃OfA, the risk of xerostomia grade ≥ II & III and dysphagia ≥ grade III were significantly reduced with the PL. For 6/15 patients the risk of xerostomia and/or dysphagia ≥ grade II could be reduced by > 2% by using PL. For the other patients, adherence to target coverage constraints was often improved. fSR₃OfA resulted in significantly improved coverage compared to PL for selected patients.

CONCLUSION

The proposed PL approach resulted in overall reduced NTCPs compared to fSR₃ and fSR₃OfA at limited cost in target coverage.

Neural network-assisted automated image registration for MRI-guided adaptive brachytherapy in cervical cancer

PURPOSE

In image-guided adaptive brachytherapy (IGABT) a quantitative evaluation of the dosimetric changes between fractions due to anatomical variations, can be implemented via rigid registration of images from subsequent fractions based on the applicator as a reference structure. With available treatment planning systems (TPS), this is a manual and time-consuming process. The aim of this retrospective study was to automate this process. A neural network (NN) was trained to predict the applicator structure from MR images. The resulting segmentation was used to automatically register MR-volumes.

MATERIAL AND METHODS

DICOM images and plans of 56 patients treated for cervical cancer with high dose-rate (HDR) brachytherapy were used in the study. A 2D and a 3D NN were trained to segment applicator structures on clinical T2-weighted MRI datasets. Different rigid registration algorithms were investigated and compared. To evaluate a fully automatic registration workflow, the NN-predicted applicator segmentations (AS) were used for rigid image registration with the best performing algorithm. The DICE coefficient and mean distance error between dwell positions (MDE) were used to evaluate segmentation and registration performance.

RESULTS

The mean DICE coefficient for the predicted AS was 0.70 ± 0.07 and 0.58 ± 0.04 for the 3D NN and 2D NN, respectively. Registration algorithms achieved MDE errors from 8.1 ± 3.7 mm (worst) to 0.7 ± 0.5 mm (best), using ground-truth AS. Using the predicted AS from the 3D NN together with the best registration algorithm, an MDE of 2.7 ± 1.4 mm was achieved.

CONCLUSION

Using a combination of deep learning models and state of the art image registration techniques has been demonstrated to be a promising solution for automatic image registration in IGABT. In combination with auto-contouring of organs at risk, the auto-registration workflow from this study could become part of an online-dosimetric interfraction evaluation workflow in the future.

Clinical implementation and validation of an automated adaptive workflow for proton therapy

Background and purpose

Treatment quality of proton therapy can be monitored by repeat-computed tomography scans (reCTs). However, manual re-delineation of target contours can be time-consuming. To improve the workflow, we implemented an automated reCT evaluation, and assessed if automatic target contour propagation would lead to the same clinical decision for plan adaptation as the manual workflow.

Materials and methods

This study included 79 consecutive patients with a total of 250 reCTs which had been manually evaluated. To assess the feasibility of automated reCT evaluation, we propagated the clinical target volumes (CTVs) deformably from the planning-CT to the reCTs in a commercial treatment planning system. The dose-volume-histogram parameters were extracted for manually re-delineated (CTV_manual) and deformably mapped target contours (CTV_auto). It was compared if CTV_manual and CTV_auto both satisfied/failed the clinical constraints. Duration of the reCT workflows was also recorded.

Results

In 92% (N = 229) of the reCTs correct flagging was obtained. Only 4% (N = 9) of the reCTs presented with false negatives (i.e., at least one clinical constraint failed for CTV_manual, but all constraints were satisfied for CTV_auto), while 5% (N = 12) of the reCTs led to a false positive. Only for one false negative reCT a plan adaption was made in clinical practice, i.e., only one adaptation would have been missed, suggesting that automated reCT evaluation was possible. Clinical introduction hereof led to a time reduction of 49 h (from 65 to 16 h).

Conclusion

Deformable target contour propagation was clinically acceptable. A script-based automatic reCT evaluation workflow has been introduced in routine clinical practice.

Personalized mid-course FDG-PET based adaptive treatment planning for non-small cell lung cancer using machine learning and optimization

Objective. Traditional radiotherapy (RT) treatment planning of non-small cell lung cancer (NSCLC) relies on population-wide estimates of organ tolerance to minimize excess toxicity. The goal of this study is to develop a personalized treatment planning based on patient-specific lung radiosensitivity, by combining machine learning and optimization.Approach. Sixty-nine non-small cell lung cancer patients with baseline and mid-treatment [18]F-fluorodeoxyglucose (FDG)-PET images were retrospectively analyzed. A probabilistic Bayesian networks (BN) model was developed to predict the risk of radiation pneumonitis (RP) at three months post-RT using pre- and mid-treatment FDG information. A patient-specific dose modifying factor (DMF), as a surrogate for lung radiosensitivity, was estimated to personalize the normal tissue toxicity probability (NTCP) model. This personalized NTCP was then integrated into a NTCP-based optimization model for RT adaptation, ensuring tumor coverage and respecting patient-specific lung radiosensitivity. The methodology was employed to adapt the treatment planning of fifteen NSCLC patients.Main results. The magnitude of the BN predicted risks corresponded with the RP severity. Average predicted risk for grade 1-4 RP were 0.18, 0.42, 0.63, and 0.76, respectively (p< 0.001). The proposed model yielded an average area under the receiver-operating characteristic curve (AUROC) of 0.84, outperforming the AUROCs of LKB-NTCP (0.77), and pre-treatment BN (0.79). Average DMF for the radio-tolerant (RP grade = 1) and radiosensitive (RP grade ≥ 2) groups were 0.8 and 1.63,p< 0.01. RT personalization resulted in five dose escalation strategies (average mean tumor dose increase = 6.47 Gy, range = [2.67-17.5]), and ten dose de-escalation (average mean lung dose reduction = 2.98 Gy [0.8-5.4]), corresponding to average NTCP reduction of 15% [4-27].Significance. Personalized FDG-PET-based mid-treatment adaptation of NSCLC RT could significantly lower the RP risk without compromising tumor control. The proposed methodology could help the design of personalized clinical trials for NSCLC patients.

Patient selection, inter-fraction plan robustness and reduction of toxicity risk with deep inspiration breath hold in intensity-modulated radiotherapy of locally advanced non-small cell lung cancer

Background

State-of-the-art radiotherapy of locally advanced non-small cell lung cancer (LA-NSCLC) is performed with intensity-modulation during free breathing (FB). Previous studies have found encouraging geometric reproducibility and patient compliance of deep inspiration breath hold (DIBH) radiotherapy for LA-NSCLC patients. However, dosimetric comparisons of DIBH with FB are sparse, and DIBH is not routinely used for this patient group. The objective of this simulation study was therefore to compare DIBH and FB in a prospective cohort of LA-NSCLC patients treated with intensity-modulated radiotherapy (IMRT).

Methods

For 38 LA-NSCLC patients, 4DCTs and DIBH CTs were acquired for treatment planning and during the first and third week of radiotherapy treatment. Using automated planning, one FB and one DIBH IMRT plan were generated for each patient. FB and DIBH was compared in terms of dosimetric parameters and NTCP. The treatment plans were recalculated on the repeat CTs to evaluate robustness. Correlations between ΔNTCPs and patient characteristics that could potentially predict the benefit of DIBH were explored.

Results

DIBH reduced the median D_mean to the lungs and heart by 1.4 Gy and 1.1 Gy, respectively. This translated into reductions in NTCP for radiation pneumonitis grade ≥2 from 20.3% to 18.3%, and for 2-year mortality from 51.4% to 50.3%. The organ at risk sparing with DIBH remained significant in week 1 and week 3 of treatment, and the robustness of the target coverage was similar for FB and DIBH. While the risk of radiation pneumonitis was consistently reduced with DIBH regardless of patient characteristics, the ability to reduce the risk of 2-year mortality was evident among patients with upper and left lower lobe tumors but not right lower lobe tumors.

Conclusion

Compared to FB, DIBH allowed for smaller target volumes and similar target coverage. DIBH reduced the lung and heart dose, as well as the risk of radiation pneumonitis and 2-year mortality, for 92% and 74% of LA-NSCLC patients, respectively. However, the advantages varied considerably between patients, and the ability to reduce the risk of 2-year mortality was dependent on tumor location. Evaluation of repeat CTs showed similar robustness of the dose distributions with each technique.

Fast VMAT planning for prostate radiotherapy: dosimetric validation of a deep learning-based initial segment generation method

Objective. To develop and evaluate a deep learning based fast volumetric modulated arc therapy (VMAT) plan generation method for prostate radiotherapy. Approach. A customized 3D U-Net was trained and validated to predict initial segments at 90 evenly distributed control points of an arc, linked to our research treatment planning system (TPS) for segment shape optimization (SSO) and segment weight optimization (SWO). For 27 test patients, the VMAT plans generated based on the deep learning prediction (VMATDL) were compared with VMAT plans generated with a previously validated automated treatment planning method (VMATref). For all test cases, the deep learning prediction accuracy, plan dosimetric quality, and the planning efficiency were quantified and analyzed. Main results. For all 27 test cases, the resulting plans were clinically acceptable. The V 95% for the PTV2 was greater than 99%, and the V 107% was below 0.2%. Statistically significant difference in target coverage was not observed between the VMATref and VMATDL plans (P = 0.3243 > 0.05). The dose sparing effect to the OARs between the two groups of plans was similar. Small differences were only observed for the Dmean of rectum and anus. Compared to the VMATref, the VMATDL reduced 29.3% of the optimization time on average. Significance. A fully automated VMAT plan generation method may result in significant improvement in prostate treatment planning efficiency. Due to the clinically acceptable dosimetric quality and high efficiency, it could potentially be used for clinical planning application and real-time adaptive therapy application after further validation.

Automated multi-criterial planning with beam angle optimization to establish non-coplanar VMAT class solutions for nasopharyngeal carcinoma

PURPOSE

Complexity in selecting optimal non-coplanar beam setups and prolonged delivery times may hamper the use of non-coplanar treatments for nasopharyngeal carcinoma (NPC). Automated multi-criterial planning with integrated beam angle optimization was used to define non-coplanar VMAT class solutions (CSs), each consisting of a coplanar arc and additional 1 or 2 fixed, non-coplanar partial arcs.

METHODS

Automated planning was used to generate a coplanar VMAT plan with 5 complementary computer-optimized non-coplanar IMRT beams (VMAT+5) for each of the 20 included patients. Subsequently, the frequency distribution of the 100 patient-specific non-coplanar IMRT beam directions was used to select non-coplanar arcs for supplementing coplanar VMAT. A second investigated CS with only one non-coplanar arc consisted of coplanar VMAT plus a partial arc at 90° couch angle (VMATCS90). Plans generated with the two VMATCSs were compared to coplanar VMAT.

RESULTS

VMAT+5 analysis resulted in VMATCS60: two partial non-coplanar arcs at couch angles 60° and -60° to complement coplanar VMAT. Compared to coplanar VMAT, the non-coplanar VMATCS60 and VMATCS90 yielded substantial average dose reductions in OARs associated with xerostomia and dysphagia, i.e., parotids, submandibular glands, oral cavity and swallowing muscles (p < 0.05) for the same PTV coverage and without violating hard constraints. Impact of non-coplanar treatment and superiority of either VMACS60 and VMATCS90 was highly patient dependent.

CONCLUSIONS

Compared to coplanar VMAT, dose to OARs was substantially reduced with a CS with one or two non-coplanar arcs. Preferences for coplanar or one or two additional arcs are highly patient-specific, balancing plan quality and treatment time.

Reflections on beam configuration optimization for intensity-modulated proton therapy

Presumably, intensity-modulated proton radiotherapy (IMPT) is the most powerful form of proton radiotherapy. In the current state of the art, IMPT beam configurations (i.e. the number of beams and their directions) are, in general, chosen subjectively based on prior experience and practicality. Beam configuration optimization (BCO) for IMPT could, in theory, significantly enhance IMPT’s therapeutic potential. However, BCO is complex and highly computer resource-intensive. Some algorithms for BCO have been developed for intensity-modulated photon therapy (IMRT). They are rarely used clinically mainly because the large number of beams typically employed in IMRT renders BCO essentially unnecessary. Moreover, in the newer form of IMRT, volumetric modulated arc therapy, there are no individual static beams. BCO is of greater importance for IMPT because it typically employs a very small number of beams (2-4) and, when the number of beams is small, BCO is critical for improving plan quality. However, the unique properties and requirements of protons, particularly in IMPT, make BCO challenging. Protons are more sensitive than photons to anatomic changes, exhibit variable relative biological effectiveness along their paths, and, as recently discovered, may spare the immune system. Such factors must be considered in IMPT BCO, though doing so would make BCO more resource intensive and make it more challenging to extend BCO algorithms developed for IMRT to IMPT. A limited amount of research in IMPT BCO has been conducted; however, considerable additional work is needed for its further development to make it truly effective and computationally practical. This article aims to provide a review of existing BCO algorithms, most of which were developed for IMRT, and addresses important requirements specific to BCO for IMPT optimization that necessitate the modification of existing approaches or the development of new effective and efficient ones.

Ultra-fast pencil beam scanning proton therapy for locally advanced non-small-cell lung cancers: field delivery within a single breath-hold

PURPOSE

The use of motion mitigation techniques such as breath-hold can reduce the dosimetric uncertainty of lung cancer proton therapy. We studied the feasibility of pencil beam scanning (PBS) proton therapy field delivery within a single breath-hold at PSI's Gantry 2.

METHODS

In PBS proton therapy, the delivery time for a field is determined by the beam-on time and the dead time between proton spots (the time required to change the energy and/or lateral position). We studied ways to reduce beam-on and lateral scanning time, without sacrificing dosimetric plan quality, aiming at a single field delivery time of 15 seconds at maximum. We tested this approach on 10 lung cases with varying target volumes. To reduce the beam-on time, we increased the beam current at the isocenter by developing new beam optics for PSI's PROSCAN beamline and Gantry 2. To reduce the dead time between the spots, we used spot-reduced plan optimization.

RESULTS

We found that it is possible to achieve conventional fractionated (2 Gy(RBE)/fraction) and hypofractionated (6 Gy(RBE)/fraction) field delivery times within a single breath-hold (<15 sec) for a variety non-small-cell lung cancer cases.

CONCLUSION

In summary, the combination of spot reduction and improved beam line transmission is a promising approach for the treatment of mobile tumours within clinically achievable breath-hold durations.

Influence of Beam Angle on Normal Tissue Complication Probability of Knowledge-Based Head and Neck Cancer Proton Planning

Knowledge-based planning solutions have brought significant improvements in treatment planning. However, the performance of a proton-specific knowledge-based planning model in creating knowledge-based plans (KBPs) with beam angles differing from those used to train the model remains unexplored. We used a previously validated RapidPlanPT model and scripting to create nine KBPs, one with default and eight with altered beam angles, for 10 recent oropharynx cancer patients. The altered-angle plans were compared against the default-angle ones in terms of grade 2 dysphagia and xerostomia normal tissue complication probability (NTCP), mean doses of several organs at risk, and dose homogeneity index (HI). As KBP could be suboptimal, a proof of principle automatic iterative optimizer (AIO) was added with the aim of reducing the plan NTCP. There were no statistically significant differences in NTCP or HI between default- and altered-angle KBPs, and the altered-angle plans showed a <1% reduction in NTCP. AIO was able to reduce the sum of grade 2 NTCPs in 66/90 cases with mean a reduction of 3.5 ± 1.8%. While the altered-angle plans saw greater benefit from AIO, both default- and altered-angle plans could be improved, indicating that the KBP model alone was not completely optimal to achieve the lowest NTCP. Overall, the data showed that the model was robust to the various beam arrangements within the range described in this analysis.

Meta-optimization for fully automated radiation therapy treatment planning

Objective. Radiation therapy treatment planning is a time-consuming process involving iterative adjustments of hyperparameters. To automate the treatment planning process, we propose a meta-optimization framework, called MetaPlanner (MP).Approach. Our MP algorithm automates planning by performing meta-optimization of treatment planning hyperparameters. The algorithm uses a derivative-free method (i.e. parallel Nelder-Mead simplex search) to search for weight configurations that minimize a meta-scoring function. Meta-scoring is performed by constructing a tier list of the relevant considerations (e.g. dose homogeneity, conformity, spillage, and OAR sparing) to mimic the clinical decision-making process. Additionally, we have made our source code publicly available via github.Main results. The proposed MP method is evaluated on two datasets (21 prostate cases and 6 head and neck cases) collected as part of clinical workflow. MP is applied to both IMRT and VMAT planning and compared to a baseline of manual VMAT plans. MP in both IMRT and VMAT scenarios has comparable or better performance than manual VMAT planning for all evaluated metrics.Significance. Our proposed MP provides a general framework for fully automated treatment planning that produces high quality treatment plans. Our MP method promises to substantially reduce the workload of treatment planners while maintaining or improving plan quality.

Artificial Intelligence in Radiation Therapy

Artificial intelligence (AI) has great potential to transform the clinical workflow of radiotherapy. Since the introduction of deep neural networks, many AI-based methods have been proposed to address challenges in different aspects of radiotherapy. Commercial vendors have started to release AI-based tools that can be readily integrated to the established clinical workflow. To show the recent progress in AI-aided radiotherapy, we have reviewed AI-based studies in five major aspects of radiotherapy including image reconstruction, image registration, image segmentation, image synthesis, and automatic treatment planning. In each section, we summarized and categorized the recently published methods, followed by a discussion of the challenges, concerns, and future development. Given the rapid development of AI-aided radiotherapy, the efficiency and effectiveness of radiotherapy in the future could be substantially improved through intelligent automation of various aspects of radiotherapy.

TBS-BAO: fully automated beam angle optimization for IMRT guided by a total-beam-space reference plan

Properly selected beam angles contribute to the quality of radiotherapy treatment plans. However, the beam angle optimization (BAO) problem is difficult to solve to optimality due to its non-convex discrete nature with many local minima. In this study, we propose TBS-BAO, a novel approach for solving the BAO problem, and test it for non-coplanar robotic CyberKnife radiotherapy for prostate cancer. First, an ideal Pareto-optimal reference dose distribution is automatically generated usinga priorimulti-criterial fluence map optimization (FMO) to generate a plan that includes all candidate beams (total-beam-space, TBS). Then, this ideal dose distribution is reproduced as closely as possible in a subsequent segmentation/beam angle optimization step (SEG/BAO), while limiting the number of allowed beams to a user-selectable preset value. SEG/BAO aims at a close reproduction of the ideal dose distribution. For each of 33 prostate SBRT patients, 18 treatment plans with different pre-set numbers of allowed beams were automatically generated with the proposed TBS-BAO. For each patient, the TBS-BAO plans were then compared to a plan that was automatically generated with an alternative BAO method (Erasmus-iCycle) and to a high-quality manually generated plan. TBS-BAO was able to automatically generate plans with clinically feasible numbers of beams (∼25), with a quality highly similar to corresponding 91-beam ideal reference plans. Compared to the alternative Erasmus-iCycle BAO approach, similar plan quality was obtained for 25-beam segmented plans, while computation times were reduced from 10.7 hours to 4.8/1.5 hours, depending on the applied pencil-beam resolution in TBS-BAO. 25-beam TBS-BAO plans had similar quality as manually generated plans with on average 48 beams, while delivery times reduced from 22.3 to 18.4/18.1 min. TBS reference plans could effectively steer the discrete non-convex BAO.

Comparing Multi-Objective Local Search Algorithms for the Beam Angle Selection Problem

In intensity-modulated radiation therapy, treatment planners aim to irradiate the tumour according to a medical prescription while sparing surrounding organs at risk as much as possible. Although this problem is inherently a multi-objective optimisation (MO) problem, most of the models in the literature are single-objective ones. For this reason, a large number of single-objective algorithms have been proposed in the literature to solve such single-objective models rather than multi-objective ones. Further, a difficulty that one has to face when solving the MO version of the problem is that the algorithms take too long before converging to a set of (approximately) non-dominated points. In this paper, we propose and compare three different strategies, namely random PLS (rPLS), judgement-function-guided PLS (jPLS) and neighbour-first PLS (nPLS), to accelerate a previously proposed Pareto local search (PLS) algorithm to solve the beam angle selection problem in IMRT. A distinctive feature of these strategies when compared to the PLS algorithms in the literature is that they do not evaluate their entire neighbourhood before performing the dominance analysis. The rPLS algorithm randomly chooses the next non-dominated solution in the archive and it is used as a baseline for the other implemented algorithms. The jPLS algorithm first chooses the non-dominated solution in the archive that has the best objective function value. Finally, the nPLS algorithm first chooses the solutions that are within the neighbourhood of the current solution. All these strategies prevent us from evaluating a large set of BACs, without any major impairment in the obtained solutions’ quality. We apply our algorithms to a prostate case and compare the obtained results to those obtained by the PLS from the literature. The results show that algorithms proposed in this paper reach a similar performance than PLS and require fewer function evaluations.

Implementation of automatic plan optimization in Italy: Status and perspectives

PURPOSE

To investigate and report on the diffusion and clinical use of automated radiotherapy planning systems in Italy and to assess the perspectives of the community of Italian medical physicists involved in radiotherapy on the use of these tools.

MATERIALS AND METHODS

A survey of medical physicists (one per Institute) of 175 radiotherapy centers in Italy was conducted between February 21st and April 1st, 2021. The information collected included the institute's characteristics, plan activity, availability/use of automatic tools and related issues regarding satisfaction, criticisms, expectations, and perceived professional modifications. Responses were analysed, including the impact of a few variables such as the institute type and experience.

RESULTS

125 of the centers (71%) answered the survey, with regional variability (range: 47%-100%); among these, 49% have a TPS with some automatic option. Clinical use of automatic planning is present in 33% of the centers, with 13% applying it in >50% of their plans. Among the 125 responding centres the most used systems are Pinnacle (16%), Raystation (9%) and Eclipse (4%). The majority of participants consider the use of automated techniques to be beneficial, while only 1% do not see any advantage; 83% of respondents see the possibility of enriching their professional role as a potential benefit, while 3% see potential threats.

CONCLUSIONS

Our survey shows that 49% of the responding centres have an automatic planning solution although clinically used in only 33% of the cases. Most physicists consider the use of automated techniques to be beneficial and show a prevalently positive attitude.

A data-driven approach to optimal beam/arc angle selection for liver stereotactic body radiation therapy treatment planning

BACKGROUND

Stereotactic body radiation therapy (SBRT) for liver cancer has shown promising therapeutic effects. Effective treatment relies not only on the precise delivery provided by image-guided radiation therapy (IGRT) but also high dose gradient formed around the treatment volume to spare functional liver tissue, which is highly dependent on the beam/arc angle selection. In this study, we aim to develop a decision support model to learn human planner's beam navigation approach for beam angle/arc angle selection for liver SBRT.

METHODS

A total of 27 liver SBRT/HIGRT patients (10 IMRT, 17 VMAT/DCA) were included in this study. A dosimetric budget index was defined for each beam angle/control point considering dose penetration through the patient body and liver tissue. Optimal beam angle setting (beam angles for IMRT and start/terminal angle for VMAT/DCA) was determined by minimizing the loss function defined as the sum of total dosimetric budget index and beam span penalty function. Leave-one-out validation was exercised on all 27 cases while weighting coefficients in the loss function was tuned in nested cross validation. To compare the efficacy of the model, a model plan was generated using automatically generated beam setting while retaining the original optimization constraints in the clinical plan. Model plan was normalized to the same planning target volume (PTV) V100% as the clinical plans. Dosimetric endpoints including PTV D98%, D2%, liver V20Gy and total MU were compared between two plan groups. Wilcoxon Signed-Rank test was performed with the null hypothesis being that no difference exists between two plan groups.

RESULTS

Beam setting prediction was instantaneous. Mean PTV D98% was 91.3% and 91.3% (P=0.566), while mean PTV D2% was 107.9% and 108.1% (P=0.164) for clinical plan and model plan respectively. Liver V20Gy showed no significant difference (P=0.590) with 23.3% for clinical plan and 23.4% for the model plan. Total MU is comparable (P=0.256) between the clinical plan (avg. 2,389.6) and model plan (avg. 2,319.6).

CONCLUSIONS

The evidence driven beam setting model yielded similar plan quality as hand-crafted clinical plan. It is capable of capturing human's knowledge in beam selection decision making. This model could facilitate decision making for beam angle selection while eliminating lengthy trial-and-error process of adjusting beam setting during liver SBRT treatment planning.

Conic formulation of fluence map optimization problems

The convexity of objectives and constraints in fluence map optimization (FMO) for radiation therapy has been extensively studied. Next to convexity, there is another important characteristic of optimization functions and problems, which has thus far not been considered in FMO literature: conic representation. Optimization problems that are conically representable using quadratic, exponential and power cones are solvable with advanced primal-dual interior-point algorithms. These algorithms guarantee an optimal solution in polynomial time and have good performance in practice. In this paper, we construct conic representations for most FMO objectives and constraints. This paper is the first that shows that FMO problems containing multiple biological evaluation criteria can be solved in polynomial time. For fractionation-corrected functions for which no exact conic reformulation is found, we provide an accurate approximation that is conically representable. We present numerical results on the TROTS data set, which demonstrate very stable numerical performance for solving FMO problems in conic form. With ongoing research in the optimization community, improvements in speed can be expected, which makes conic optimization a promising alternative for solving FMO problems.

Evaluation of alternative parameter settings for dose restoration and full plan adaptation in IMPT for prostate cancer

BACKGROUND/PURPOSE

Intensity-modulated proton therapy is highly sensitive to anatomical variations. A dose restoration method and a full plan adaptation method have been developed earlier, both requiring several parameter settings. This study evaluates the validity of the previously selected settings by systematically comparing them to alternatives.

MATERIALS/METHODS

The dose restoration method takes a prior plan and uses an energy-adaptation followed by a spot-intensity re-optimization to restore the plan to its initial state. The full adaptation method uses an energy-adaptation followed by the addition of new spots and a spot-intensity optimization to fit the new anatomy. We varied: 1) The margins and robustness settings of the prior plan, 2) the spot-addition sample size, i.e. the number of added spots, 3) the spot-addition stopping criterion, and 4) the spot-intensity optimization approach. The last three were evaluated only for the full plan adaptation. Evaluations were done on 88 CT scans of 11 prostate cancer patients. Dose was prescribed as 55 Gy(RBE) to the lymph nodes and seminal vesicles with a boost to 74 Gy(RBE) to the prostate.

RESULTS

For the dose restoration method, changing the applied CTV-to-PTV margins and plan robustness in the prior plans yielded insufficient target coverage or increased OAR doses. For the full plan adaptation, more spot-addition iterations and using a different optimization approach resulted in lower OAR doses compared to the default settings while maintaining target coverage. However, the calculation times increased by up to 20 times, making these variations infeasible for online-adaptation.

CONCLUSION

We recommend maintaining the default setting for the dose restoration approach. For the full plan adaptation we recommend to focus on fine-tuning the optimization-parameters, and apart from this using the default settings.