Purpose. Recently, we introduced a bi-objective optimization approach based on dose-volume indices to automatically create clinically good HDR prostate brachytherapy plans. To calculate dose-volume... Show morePurpose. Recently, we introduced a bi-objective optimization approach based on dose-volume indices to automatically create clinically good HDR prostate brachytherapy plans. To calculate dose-volume indices, a reconstruction algorithm is used to determine the 3D organ shape from 2D contours, inevitably containing settings that influence the result. We augment the optimization approach to quickly find plans that are robust to differences in 3D reconstruction. Methods. Studied reconstruction settings were: interpolation between delineated organ contours, overlap between contours, and organ shape at the top and bottom contour. Two options for each setting yields 8 possible 3D organ reconstructions per patient, over which the robust model defines minimax optimization. For the original model, settings were based on our treatment planning system. Both models were tested on data of 26 patients and compared by re-evaluating selected optimized plans both in the original model (1 organ reconstruction, the difference determines the cost), and in the robust model (8 organ reconstructions, the difference determines the benefit). Results. Robust optimization increased the run time from 3 to 6 min. The median cost for robust optimization as observed in the original model was -0.25% in the dose-volume indices with a range of [-0.01%, -1.03%]. The median benefit of robust optimization as observed in the robust model was 0.93% with a range of [0.19%, 4.16%]. For 4 patients, selected plans that appeared good when optimized in the original model, violated the clinical protocol with more than 1% when considering different settings. This was not the case for robustly optimized plans. Conclusions. Plans of high quality, irrespective of 3D organ reconstruction settings, can be obtained using our robust optimization approach. With its limited effect on total runtime, our approach therefore offers a way to account for dosimetry uncertainties that result from choices in organ reconstruction settings that is viable in clinical practice. Show less
Meer, M.C. van der; Bosman, P.A.N.; Niatsetski, Y.; Alderliesten, T.; Wieringen, N. van; Pieters, B.R.; Bel, A. 2020
Purpose Bi-objective simultaneous optimization of catheter positions and dwell times for high-dose-rate (HDR) prostate brachytherapy, based directly on dose-volume indices, has shown promising... Show morePurpose Bi-objective simultaneous optimization of catheter positions and dwell times for high-dose-rate (HDR) prostate brachytherapy, based directly on dose-volume indices, has shown promising results. However, optimization with the state-of-the-art evolutionary algorithm MO-RV-GOMEA so far required several hours of runtime, and resulting catheter positions were not always clinically feasible. The aim of this study is to extend the optimization model and apply GPU parallelization to achieve clinically acceptable computation times. The resulting optimization procedure is compared with a previously introduced method based solely on geometric criteria, the adapted Centroidal Voronoi Tessellations (CVT) algorithm. Methods Bi-objective simultaneous optimization was performed with a GPU-parallelized version of MO-RV-GOMEA. This optimization of catheter positions and dwell times was retrospectively applied to the data of 26 patients previously treated with HDR prostate brachytherapy for 8-16 catheters (steps of 2). Optimization of catheter positions using CVT was performed in seconds, after which optimization of only the dwell times using MO-RV-GOMEA was performed in 1 min. Results Simultaneous optimization of catheter positions and dwell times using MO-RV-GOMEA was performed in 5 min. For 16 down to 8 catheters (steps of 2), MO-RV-GOMEA found plans satisfying the planning-aims for 20, 20, 18, 14, and 11 out of the 26 patients, respectively. CVT achieved this for 19, 17, 13, 9, and 2 patients, respectively. TheP-value for the difference between MO-RV-GOMEA and CVT was 0.023 for 16 catheters, 0.005 for 14 catheters, and <0.001 for 12, 10, and 8 catheters. Conclusions With bi-objective simultaneous optimization on a GPU, high-quality catheter positions can now be obtained within 5 min, which is clinically acceptable, but slower than CVT. For 16 catheters, the difference between MO-RV-GOMEA and CVT is clinically irrelevant. For 14 catheters and less, MO-RV-GOMEA outperforms CVT in finding plans satisfying all planning-aims. Show less
Maree, S.C.; Bosman, P.A.N.; Wieringen, N. van; Niatsetski, Y.; Pieters, B.R.; Bel, A.; Alderliesten, T. 2020
We present an automatic bi-objective parameter-tuning approach for inverse planning methods for high-dose-rate prostate brachytherapy, which aims to overcome the difficult and time-consuming manual... Show moreWe present an automatic bi-objective parameter-tuning approach for inverse planning methods for high-dose-rate prostate brachytherapy, which aims to overcome the difficult and time-consuming manual parameter tuning that is currently required to obtain patient-specific high-quality treatment plans. We modelled treatment planning as a bi-objective optimization problem, in which dose-volume-based planning criteria related to target coverage are explicitly separated from organ-sparing criteria. When this model is optimized, a large set of high-quality plans with different trade-offs can be obtained. This set can be visualized as an insightful patient-specific trade-off curve. In our parameter-tuning approach, the parameters of inverse planning methods are automatically tuned, aimed to maximize the two objectives of the bi-objective planning model. By generating trade-off curves for different inverse planning methods, their maximally achievable plan quality can be insightfully compared. Automatic parameter tuning furthermore allows to construct standard parameter sets (class solutions) representing different trade-offs in a principled way, which can be directly used in current clinical practice. In this work, we considered the inverse planning methods IPSA and HIPO. Thirty-nine previously treated prostate cancer patients were included. We compared automatic parameter tuning, random parameter sampling, and the maximally achievable plan quality obtained by directly optimizing the bi-objective planning model with the state-of-the-art optimization software GOMEA. We showed that for each patient, a set of plans with a wide range of trade-offs could be obtained using automatic parameter tuning for both IPSA and HIPO. By tuning HIPO, better trade-offs were obtained than by tuning IPSA. For most patients, automatic tuning of HIPO resulted in plans close to the maximally achievable plan quality obtained by optimizing the bi-objective planning model directly. Automatic parameter tuning was shown to improve plan quality significantly compared to random parameter sampling. Finally, from the automatically-tuned plans, three class solutions were successfully constructed representing different trade-offs. Show less