Bi‐objective optimization of catheter positions for high‐dose‐rate prostate brachytherapy

• Published in: Medical physics (Lancaster) Vol. 47; no. 12; pp. 6077 - 6086
• Main Authors: Meer, Marjolein C., Bosman, Peter A.N., Niatsetski, Yury, Alderliesten, Tanja, Wieringen, Niek, Pieters, Bradley R., Bel, Arjan
• Format: Journal Article
• Language: English
• Published: United States John Wiley and Sons Inc 01-12-2020
• Subjects: bi‐objective optimization; Brachytherapy; catheter positions; Catheters; HDR brachytherapy; Humans; Male; prostate neoplasms; Prostatic Neoplasms - radiotherapy; Radiotherapy Dosage; Radiotherapy Planning, Computer-Assisted; Retrospective Studies; THERAPEUTIC INTERVENTIONS; treatment planning; bi-objective optimization; treatment planning; prostate neoplasms; catheter positions; HDR brachytherapy
• ISSN: 0094-2405; 2473-4209
• DOI: 10.1002/mp.14505

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 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. The P‐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.