Parameterless Gene-pool Optimal Mixing Evolutionary Algorithms.

When it comes to solving optimization problems with evolutionary algorithms (EAs) in a reliable and scalable manner, detecting and exploiting linkage information, i.e., dependencies between variables, can be key. In this article, we present the latest version of, and propose substantial enhancements to, the Gene-pool Optimal Mixing Evoutionary Algorithm (GOMEA): an EA explicitly designed to estimate and exploit linkage information. We begin by performing a largescale search over several GOMEA design choices to understand what matters most and obtain a generally best-performing version of the algorithm. Next, we introduce a novel version of GOMEA, called CGOMEA, where linkage-based variation is further improved by filtering solution mating based on conditional dependencies. We compare our latest version of GOMEA, the newly introduced CGOMEA, and another contending linkage-aware EA, DSMGA-II, in an extensive experimental evaluation, involving a benchmark set of 9 black-box problems that can only be solved efficiently if their inherent dependency structure is unveiled and exploited. Finally, in an attempt to make EAs more usable and resilient to parameter choices, we investigate the performance of different automatic population management schemes for GOMEA and CGOMEA, de facto making the EAs parameterless. Our results show that GOMEA and CGOMEA significantly outperform the original GOMEA and DSMGA-II on most problems, setting a new state of the art for the field.

Achieving Highly Scalable Evolutionary Real-Valued Optimization by Exploiting Partial Evaluations.

It is known that to achieve efficient scalability of an Evolutionary Algorithm (EA), dependencies (also known as linkage) must be properly taken into account during variation. In a Gray-Box Optimization (GBO) setting, exploiting prior knowledge regarding these dependencies can greatly benefit optimization. We specifically consider the setting where partial evaluations are possible, meaning that the partial modification of a solution can be efficiently evaluated. Such problems are potentially very difficult, for example, non-separable, multimodal, and multiobjective. The Gene-pool Optimal Mixing Evolutionary Algorithm (GOMEA) can effectively exploit partial evaluations, leading to a substantial improvement in performance and scalability. GOMEA was recently shown to be extendable to real-valued optimization through a combination with the real-valued estimation of distribution algorithm AMaLGaM. In this article, we definitively introduce the Real-Valued GOMEA (RV-GOMEA), and introduce a new variant, constructed by combining GOMEA with what is arguably the best-known real-valued EA, the Covariance Matrix Adaptation Evolution Strategies (CMA-ES). Both variants of GOMEA are compared to L-BFGS and the Limited Memory CMA-ES (LM-CMA-ES). We show that both variants of RV-GOMEA achieve excellent performance and scalability in a GBO setting, which can be orders of magnitude better than that of EAs unable to efficiently exploit the GBO setting.

Sensitivity of dose-volume indices to computation settings in high-dose-rate prostate brachytherapy treatment plan evaluation.

PURPOSE: To investigate the variation in computed dose-volume (DV) indices for high-dose-rate (HDR) prostate brachytherapy that can result from typical differences in computation settings in treatment planning systems (TPSs). METHODS: Five factors were taken into account: number of dose-calculation points, radioactive source description, interpolation between delineated contours, intersections between delineated organ contours, and organ shape at the top and bottom contour using either full or partial slice thickness. Using in-house developed software, the DV indices of the treatment plans of 26 patients were calculated with different settings, and compared to a baseline setting that closely followed the default settings of the TPS used in our medical center. Studied organs were prostate and seminal vesicles, denoted as targets, and bladder, rectum, and urethra, denoted as organs at risk (OARs), which were delineated on MRI scans with a 3.3 mm slice thickness. RESULTS: When sampling a fixed number of points in each organ, in order to achieve a width of the 95% confidence interval over all patients of the DV indices of 1% or less, only 32,000 points had to be sampled per target, but 256,000 points had to be sampled per OAR. For the remaining factors, DV indices changed up to 0.4% for rectum, 1.3% for urethra, and 2.6% for prostate. DV indices of the bladder changed especially if the high-dose-region was (partly) located at the most caudal contour, up to 8.5%, and DV indices of the vesicles changed especially if there were few delineated contours, up to 9.8%, both due to the use of full slice thickness for the top and bottom contour. CONCLUSIONS: The values of DV indices used in prostate HDR brachytherapy treatment planning are influenced by the computation settings in a TPS, especially at the most caudal part of the bladder, as well as in the seminal vesicles.

How do patient characteristics and anatomical features correlate to accuracy of organ dose reconstruction for Wilms' tumor radiation treatment plans when using a surrogate patient's CT scan?

In retrospective radiation treatment (RT) dosimetry, a surrogate anatomy is often used for patients without 3D CT. To gain insight in what the crucial aspects in a surrogate anatomy are to enable accurate dose reconstruction, we investigated the relation of patient characteristics and internal anatomical features with deviations in reconstructed organ dose using surrogate patient's CT scans. Abdominal CT scans of 35 childhood cancer patients (age: 2.1-5.6 yr; 17 boys, 18 girls) undergoing RT during 2004-2016 were included. Based on whether an intact right or left kidney is present in the CT scan, two groups were formed each containing 24 patients. From each group, four CTs associated with Wilms' tumor RT plans with an anterior-posterior-posterior-anterior field setup were selected as references. For each reference, a 2D digitally reconstructed radiograph was computed from the reference CT to simulate a 2D radiographic image and dose reconstruction was performed on the other CTs in the respective group. Deviations in organ mean dose (DEmean) of the reconstructions versus the references were calculated, as were deviations in patient characteristics (i.e. age, height, weight) and in anatomical features including organ volume, location (in 3D), and spatial overlaps. Per reference, the Pearson's correlation coefficient between deviations in DEmean and patient characteristics/features were studied. Deviation in organ locations and DEmean for the liver, spleen, and right kidney were moderately correlated (R2 > 0.5) for 8/8, 5/8, and 3/4 reference plans, respectively. Deviations in organ volume or spatial overlap and DEmean for the right and left kidney were weakly correlated (0.3 < R2 < 0.5) in 4/4 and 1/4 reference plans. No correlations (R2 < 0.3) were found between deviations in age or height and DEmean. Therefore, the performance of organ dose reconstruction using surrogate patients' CT scans is primarily related to deviation in organ location, followed by volume and spatial overlap. Further, results were plan dependent.

GAMBIT: A Parameterless Model-Based Evolutionary Algorithm for Mixed-Integer Problems.

Learning and exploiting problem structure is one of the key challenges in optimization. This is especially important for black-box optimization (BBO) where prior structural knowledge of a problem is not available. Existing model-based Evolutionary Algorithms (EAs) are very efficient at learning structure in both the discrete, and in the continuous domain. In this article, discrete and continuous model-building mechanisms are integrated for the Mixed-Integer (MI) domain, comprising discrete and continuous variables. We revisit a recently introduced model-based evolutionary algorithm for the MI domain, the Genetic Algorithm for Model-Based mixed-Integer opTimization (GAMBIT). We extend GAMBIT with a parameterless scheme that allows for practical use of the algorithm without the need to explicitly specify any parameters. We furthermore contrast GAMBIT with other model-based alternatives. The ultimate goal of processing mixed dependences explicitly in GAMBIT is also addressed by introducing a new mechanism for the explicit exploitation of mixed dependences. We find that processing mixed dependences with this novel mechanism allows for more efficient optimization. We further contrast the parameterless GAMBIT with Mixed-Integer Evolution Strategies (MIES) and other state-of-the-art MI optimization algorithms from the General Algebraic Modeling System (GAMS) commercial algorithm suite on problems with and without constraints, and show that GAMBIT is capable of solving problems where variable dependences prevent many algorithms from successfully optimizing them.

Exploiting Linkage Information and Problem-Specific Knowledge in Evolutionary Distribution Network Expansion Planning.

This article tackles the Distribution Network Expansion Planning (DNEP) problem that has to be solved by distribution network operators to decide which, where, and/or when enhancements to electricity networks should be introduced to satisfy the future power demands. Because of many real-world details involved, the structure of the problem is not exploited easily using mathematical programming techniques, for which reason we consider solving this problem with evolutionary algorithms (EAs). We compare three types of EAs for optimizing expansion plans: the classic genetic algorithm (GA), the estimation-of-distribution algorithm (EDA), and the Gene-pool Optimal Mixing Evolutionary Algorithm (GOMEA). Not fully knowing the structure of the problem, we study the effect of linkage learning through the use of three linkage models: univariate, marginal product, and linkage tree. We furthermore experiment with the impact of incorporating different levels of problem-specific knowledge in the variation operators. Experiments show that the use of problem-specific variation operators is far more important for the classic GA to find high-quality solutions. In all EAs, the marginal product model and its linkage learning procedure have difficulty in capturing and exploiting the DNEP problem structure. GOMEA, especially when combined with the linkage tree structure, is found to have the most robust performance by far, even when an out-of-the-box variant is used that does not exploit problem-specific knowledge. Based on experiments, we suggest that when selecting optimization algorithms for power system expansion planning problems, EAs that have the ability to effectively model and efficiently exploit problem structures, such as GOMEA, should be given priority, especially in the case of black-box or grey-box optimization.

Keeping your best options open with AI-based treatment planning in prostate and cervix brachytherapy.

PURPOSE: Without a clear definition of an optimal treatment plan, no optimization model can be perfect. Therefore, instead of automatically finding a single "optimal" plan, finding multiple, yet different near-optimal plans, can be an insightful approach to support radiation oncologists in finding the plan they are looking for. METHODS AND MATERIALS: BRIGHT is a flexible AI-based optimization method for brachytherapy treatment planning that has already been shown capable of finding high-quality plans that trade-off target volume coverage and healthy tissue sparing. We leverage the flexibility of BRIGHT to find plans with similar dose-volume criteria, yet different dose distributions. We further describe extensions that facilitate fast plan adaptation should planning aims need to be adjusted, and straightforwardly allow incorporating hospital-specific aims besides standard protocols. RESULTS: Results are obtained for prostate (n = 12) and cervix brachytherapy (n = 36). We demonstrate the possible differences in dose distribution for optimized plans with equal dose-volume criteria. We furthermore demonstrate that adding hospital-specific aims enables adhering to hospital-specific practice while still being able to automatically create cervix plans that more often satisfy the EMBRACE-II protocol than clinical practice. Finally, we illustrate the feasibility of fast plan adaptation. CONCLUSIONS: Methods such as BRIGHT enable new ways to construct high-quality treatment plans for brachytherapy while offering new insights by making explicit the options one has. In particular, it becomes possible to present to radiation oncologists a manageable set of alternative plans that, from an optimization perspective are equally good, yet differ in terms of coverage-sparing trade-offs and shape of the dose distribution.

Benchmarking parameter-free AMaLGaM on functions with and without noise.

We describe a parameter-free estimation-of-distribution algorithm (EDA) called the adapted maximum-likelihood Gaussian model iterated density-estimation evolutionary algorithm (AMaLGaM-ID[Formula: see text]A, or AMaLGaM for short) for numerical optimization. AMaLGaM is benchmarked within the 2009 black box optimization benchmarking (BBOB) framework and compared to a variant with incremental model building (iAMaLGaM). We study the implications of factorizing the covariance matrix in the Gaussian distribution, to use only a few or no covariances. Further, AMaLGaM and iAMaLGaM are also evaluated on the noisy BBOB problems and we assess how well multiple evaluations per solution can average out noise. Experimental evidence suggests that parameter-free AMaLGaM can solve a wide range of problems efficiently with perceived polynomial scalability, including multimodal problems, obtaining the best or near-best results among all algorithms tested in 2009 on functions such as the step-ellipsoid and Katsuuras, but failing to locate the optimum within the time limit on skew Rastrigin-Bueche separable and Lunacek bi-Rastrigin in higher dimensions. AMaLGaM is found to be more robust to noise than iAMaLGaM due to the larger required population size. Using few or no covariances hinders the EDA from dealing with rotations of the search space. Finally, the use of noise averaging is found to be less efficient than the direct application of the EDA unless the noise is uniformly distributed. AMaLGaM was among the best performing algorithms submitted to the BBOB workshop in 2009.

Automatic landmark correspondence detection in medical images with an application to deformable image registration.

Purpose: Deformable image registration (DIR) can benefit from additional guidance using corresponding landmarks in the images. However, the benefits thereof are largely understudied, especially due to the lack of automatic landmark detection methods for three-dimensional (3D) medical images. Approach: We present a deep convolutional neural network (DCNN), called DCNN-Match, that learns to predict landmark correspondences in 3D images in a self-supervised manner. We trained DCNN-Match on pairs of computed tomography (CT) scans containing simulated deformations. We explored five variants of DCNN-Match that use different loss functions and assessed their effect on the spatial density of predicted landmarks and the associated matching errors. We also tested DCNN-Match variants in combination with the open-source registration software Elastix to assess the impact of predicted landmarks in providing additional guidance to DIR. Results: We tested our approach on lower abdominal CT scans from cervical cancer patients: 121 pairs containing simulated deformations and 11 pairs demonstrating clinical deformations. The results showed significant improvement in DIR performance when landmark correspondences predicted by DCNN-Match were used in the case of simulated ( p = 0 e 0 ) as well as clinical deformations ( p = 0.030 ). We also observed that the spatial density of the automatic landmarks with respect to the underlying deformation affect the extent of improvement in DIR. Finally, DCNN-Match was found to generalize to magnetic resonance imaging scans without requiring retraining, indicating easy applicability to other datasets. Conclusions: DCNN-match learns to predict landmark correspondences in 3D medical images in a self-supervised manner, which can improve DIR performance.

Towards artificial intelligence-based automated treatment planning in clinical practice: A prospective study of the first clinical experiences in high-dose-rate prostate brachytherapy.

PURPOSE: This prospective study evaluates our first clinical experiences with the novel ``BRachytherapy via artificial Intelligent GOMEA-Heuristic based Treatment planning'' (BRIGHT) applied to high-dose-rate prostate brachytherapy. METHODS AND MATERIALS: Between March 2020 and October 2021, 14 prostate cancer patients were treated in our center with a 15Gy HDR-brachytherapy boost. BRIGHT was used for bi-objective treatment plan optimization and selection of the most desirable plans from a coverage-sparing trade-off curve. Selected BRIGHT plans were imported into the commercial treatment planning system Oncentra Brachy . In Oncentra Brachy a dose distribution comparison was performed for clinical plan choice, followed by manual fine-tuning of the preferred BRIGHT plan when deemed necessary. The reasons for plan selection, clinical plan choice, and fine-tuning, as well as process speed were monitored. For each patient, the dose-volume parameters of the (fine-tuned) clinical plan were evaluated. RESULTS: In all patients, BRIGHT provided solutions satisfying all protocol values for coverage and sparing. In four patients not all dose-volume criteria of the clinical plan were satisfied after manual fine-tuning. Detailed information on tumour coverage, dose-distribution, dwell time pattern, and insight provided by the patient-specific trade-off curve, were used for clinical plan choice. Median time spent on treatment planning was 42 min, consisting of 16 min plan optimization and selection, and 26 min undesirable process steps. CONCLUSIONS: BRIGHT is implemented in our clinic and provides automated prostate high-dose-rate brachytherapy planning with trade-off based plan selection. Based on our experience, additional optimization aims need to be implemented to further improve direct clinical applicability of treatment plans and process efficiency.

Interactive Education on Sleep Hygiene with a Social Robot at a Pediatric Oncology Outpatient Clinic: Feasibility, Experiences, and Preliminary Effectiveness.

OBJECTIVES: Children with cancer often experience sleep problems, which are associated with many negative physical and psychological health outcomes, as well as with a lower quality of life. Therefore, interventions are strongly required to improve sleep in this population. We evaluated interactive education with respect to sleep hygiene with a social robot at a pediatric oncology outpatient clinic regarding the feasibility, experiences, and preliminary effectiveness. METHODS: Researchers approached children (8 to 12 years old) who were receiving anticancer treatment and who were visiting the outpatient clinic with their parents during the two-week study period. The researchers completed observation forms regarding feasibility, and parents completed the Children's Sleep Hygiene Scale before and two weeks after the educational regimen. The experiences of children and parents were evaluated in semi-structured interviews. We analyzed open answers by labeling each answer with a topic reflecting the content and collapsed these topics into categories. We used descriptive statistics to describe the feasibility and experiences, and a dependent-samples t-test to evaluate the preliminary effectiveness. RESULTS: Twenty-eight families participated (58% response rate) and all interactions with the robot were completed. The children and parents reported that they learned something new (75% and 50%, respectively), that they wanted to learn from the robot more often (83% and 75%, respectively), and that they applied the sleeping tips from the robot afterwards at home (54%). Regarding the preliminary effectiveness, children showed a statistically significant improvement in their sleep hygiene (p = 0.047, d = 0.39). CONCLUSIONS: Providing an educational regimen on sleep hygiene in a novel, interactive way by using a social robot at the outpatient clinic seemed feasible, and the children and parents mostly exhibited positive reactions. We found preliminary evidence that the sleep hygiene of children with cancer improved.

Validation and Comparison of Radiograph-Based Organ Dose Reconstruction Approaches for Wilms Tumor Radiation Treatment Plans.

Purpose: Our purpose was to validate and compare the performance of 4 organ dose reconstruction approaches for historical radiation treatment planning based on 2-dimensional radiographs. Methods and Materials: We considered 10 patients with Wilms tumor with planning computed tomography images for whom we developed typical historic Wilms tumor radiation treatment plans, using anteroposterior and posteroanterior parallel-opposed 6 MV flank fields, normalized to 14.4 Gy. Two plans were created for each patient, with and without corner blocking. Regions of interest (lungs, heart, nipples, liver, spleen, contralateral kidney, and spinal cord) were delineated, and dose-volume metrics including organ mean and minimum dose (Dmean and Dmin) were computed as the reference baseline for comparison. Dosimetry for the 20 plans was then independently reconstructed using 4 different approaches. Three approaches involved surrogate anatomy, among which 2 used demographic-matching criteria for phantom selection/building, and 1 used machine learning. The fourth approach was also machine learning-based, but used no surrogate anatomies. Absolute differences in organ dose-volume metrics between the reconstructed and the reference values were calculated. Results: For Dmean and Dmin (average and minimum point dose) all 4 dose reconstruction approaches performed within 10% of the prescribed dose (≤1.4 Gy). The machine learning-based approaches showed a slight advantage for several of the considered regions of interest. For Dmax (maximum point dose), the absolute differences were much higher, that is, exceeding 14% (2 Gy), with the poorest agreement observed for near-beam and out-of-beam organs for all approaches. Conclusions: The studied approaches give comparable dose reconstruction results, and the choice of approach for cohort dosimetry for late effects studies should still be largely driven by the available resources (data, time, expertise, and funding).

Robust optimization for HDR prostate brachytherapy applied to organ reconstruction uncertainty.

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 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.

Machine learning for the prediction of pseudorealistic pediatric abdominal phantoms for radiation dose reconstruction.

Purpose: Current phantoms used for the dose reconstruction of long-term childhood cancer survivors lack individualization. We design a method to predict highly individualized abdominal three-dimensional (3-D) phantoms automatically. Approach: We train machine learning (ML) models to map (2-D) patient features to 3-D organ-at-risk (OAR) metrics upon a database of 60 pediatric abdominal computed tomographies with liver and spleen segmentations. Next, we use the models in an automatic pipeline that outputs a personalized phantom given the patient's features, by assembling 3-D imaging from the database. A step to improve phantom realism (i.e., avoid OAR overlap) is included. We compare five ML algorithms, in terms of predicting OAR left-right (LR), anterior-posterior (AP), inferior-superior (IS) positions, and surface Dice-Sørensen coefficient (sDSC). Furthermore, two existing human-designed phantom construction criteria and two additional control methods are investigated for comparison. Results: Different ML algorithms result in similar test mean absolute errors: ∼ 8 mm for liver LR, IS, and spleen AP, IS; ∼ 5 mm for liver AP and spleen LR; ∼ 80 % for abdomen sDSC; and ∼ 60 % to 65% for liver and spleen sDSC. One ML algorithm (GP-GOMEA) significantly performs the best for 6/9 metrics. The control methods and the human-designed criteria in particular perform generally worse, sometimes substantially ( + 5 - mm error for spleen IS, - 10 % sDSC for liver). The automatic step to improve realism generally results in limited metric accuracy loss, but fails in one case (out of 60). Conclusion: Our ML-based pipeline leads to phantoms that are significantly and substantially more individualized than currently used human-designed criteria.

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

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.

Automatic generation of three-dimensional dose reconstruction data for two-dimensional radiotherapy plans for historically treated patients.

Performing large-scale three-dimensional radiation dose reconstruction for patients requires a large amount of manual work. We present an image processing-based pipeline to automatically reconstruct radiation dose. The pipeline was designed for childhood cancer survivors that received abdominal radiotherapy with anterior-to-posterior and posterior-to-anterior field set-up. First, anatomical landmarks are automatically identified on two-dimensional radiographs. Second, these landmarks are used to derive parameters to emulate the geometry of the plan on a surrogate computed tomography. Finally, the plan is emulated and used as input for dose calculation. For qualitative evaluation, 100 cases of automatic and manual plan emulations were assessed by two experienced radiation dosimetrists in a blinded comparison. The two radiation dosimetrists approved 100%/100% and 92%/91% of the automatic/manual plan emulations, respectively. Similar approval rates of 100% and 94% hold when the automatic pipeline is applied on another 50 cases. Further, quantitative comparisons resulted in on average < 5 mm difference in plan isocenter/borders, and < 0.9 Gy in organ mean dose (prescribed dose: 14.4 Gy) calculated from the automatic and manual plan emulations. No statistically significant difference in terms of dose reconstruction accuracy was found for most organs at risk. Ultimately, our automatic pipeline results are of sufficient quality to enable effortless scaling of dose reconstruction data generation.

GPU-accelerated bi-objective treatment planning for prostate high-dose-rate brachytherapy.

PURPOSE: The purpose of this study is to improve upon a recently introduced bi-objective treatment planning method for prostate high-dose-rate (HDR) brachytherapy (BT), both in terms of resulting plan quality and runtime requirements, to the extent that its execution time is clinically acceptable. METHODS: Bi-objective treatment planning is done using a state-of-the-art multiobjective evolutionary algorithm, which produces a large number of potential treatment plans with different trade-offs between coverage of the target volumes and sparing organs at risk. A graphics processing unit (GPU) is used for large-scale parallelization of dose calculations and the calculation of the dose-volume (DV) indices of potential treatment plans. Moreover, the objectives of the previously used bi-objective optimization model are modified to produce better results. RESULTS: We applied the GPU-accelerated bi-objective treatment planning method to a set of 18 patients, resulting in a set containing a few hundred potential treatment plans with different trade-offs for each of these patients. Due to accelerations introduced in this article, results previously achieved after 1 hour are now achieved within 30 seconds of optimization. We found plans satisfying the clinical protocol for 15 of 18 patients, whereas this was the case for only 4 of 18 clinical plans. Higher quality treatment plans are obtained when the accuracy of DV index calculation is increased using more dose calculation points, requiring still no more than 3 minutes of optimization for 100 000 points. CONCLUSIONS: Large sets of high-quality treatment plans that trade-off coverage and sparing are now achievable within 30 seconds, due to the GPU-acceleration of a previously introduced bi-objective treatment planning method for prostate HDR brachytherapy. Higher quality plans can be achieved when optimizing for 3 minutes, which we still consider to be clinically acceptable. This allows for more insightful treatment plan selection in a clinical setting.

Evaluation of bi-objective treatment planning for high-dose-rate prostate brachytherapy-A retrospective observer study.

PURPOSE: Bi-objective treatment planning for high-dose-rate prostate brachytherapy is a novel treatment planning method with two separate objectives that represent target coverage and organ-at-risk sparing. In this study, we investigated the feasibility and plan quality of this method by means of a retrospective observer study. METHODS AND MATERIALS: Current planning sessions were recorded to configure a bi-objective optimization model and to assess its applicability to our clinical practice. Optimization software, GOMEA, was then used to automatically generate a large set of plans with different trade-offs in the two objectives for each of 18 patients treated with high-dose-rate prostate brachytherapy. From this set, five plans per patient were selected for comparison to the clinical plan in terms of satisfaction of planning criteria and in a retrospective observer study. Three brachytherapists were asked to evaluate the blinded plans and select the preferred one. RESULTS: Recordings demonstrated applicability of the bi-objective optimization model to our clinical practice. For 14/18 patients, GOMEA plans satisfied all planning criteria, compared with 4/18 clinical plans. In the observer study, in 53/54 cases, a GOMEA plan was preferred over the clinical plan. When asked for consensus among observers, this ratio was 17/18 patients. Observers highly appreciated the insight gained from comparing multiple plans with different trade-offs simultaneously. CONCLUSIONS: The bi-objective optimization model adapted well to our clinical practice. GOMEA plans were considered equal or superior to the clinical plans. In addition, presenting multiple high-quality plans provided novel insight into patient-specific trade-offs.

Simplex-based navigation tool for a posteriori selection of the preferred deformable image registration outcome from a set of trade-off solutions obtained with multiobjective optimization for the case of breast MRI.

Multiobjective optimization approaches for deformable image registration (DIR) remove the need for manual adjustment of key parameters and provide a set of solutions that represent high-quality trade-offs between objectives of interest. Choosing a desired outcome a posteriori is potentially far more insightful as differences between solutions can be immediately visualized. The purpose of this work is to investigate whether such an approach allows clinical experts to intuitively select their preferred DIR outcome. To this end, we developed a simplex-based tool for solution navigation and asked 10 clinical experts to use it to choose their preferred DIR outcome from sets of trade-off solutions obtained for 10 breast magnetic resonance DIR cases of low (prone-prone DIR; n = 5 ) and high (prone-supine DIR; n = 5 ) difficulty, of patients and volunteers, respectively. The usability of the software is subsequently evaluated by the observers using the system usability scale. Further, the quality of the selected DIR outcomes is evaluated using the mean target registration error. Results show that the users are able to identify and select high-quality DIR outcomes, and attested to high learnability and usability of our software, supporting the validity of the presumed added value of taking a multiobjective perspective on DIR in clinical practice.

On the feasibility of automatically selecting similar patients in highly individualized radiotherapy dose reconstruction for historic data of pediatric cancer survivors.

PURPOSE: The aim of this study is to establish the first step toward a novel and highly individualized three-dimensional (3D) dose distribution reconstruction method, based on CT scans and organ delineations of recently treated patients. Specifically, the feasibility of automatically selecting the CT scan of a recently treated childhood cancer patient who is similar to a given historically treated child who suffered from Wilms' tumor is assessed. METHODS: A cohort of 37 recently treated children between 2- and 6-yr old are considered. Five potential notions of ground-truth similarity are proposed, each focusing on different anatomical aspects. These notions are automatically computed from CT scans of the abdomen and 3D organ delineations (liver, spleen, spinal cord, external body contour). The first is based on deformable image registration, the second on the Dice similarity coefficient, the third on the Hausdorff distance, the fourth on pairwise organ distances, and the last is computed by means of the overlap volume histogram. The relationship between typically available features of historically treated patients and the proposed ground-truth notions of similarity is studied by adopting state-of-the-art machine learning techniques, including random forest. Also, the feasibility of automatically selecting the most similar patient is assessed by comparing ground-truth rankings of similarity with predicted rankings. RESULTS: Similarities (mainly) based on the external abdomen shape and on the pairwise organ distances are highly correlated (Pearson rp ≥ 0.70) and are successfully modeled with random forests based on historically recorded features (pseudo-R2 ≥ 0.69). In contrast, similarities based on the shape of internal organs cannot be modeled. For the similarities that random forest can reliably model, an estimation of feature relevance indicates that abdominal diameters and weight are the most important. Experiments on automatically selecting similar patients lead to coarse, yet quite robust results: the most similar patient is retrieved only 22% of the times, however, the error in worst-case scenarios is limited, with the fourth most similar patient being retrieved. CONCLUSIONS: Results demonstrate that automatically selecting similar patients is feasible when focusing on the shape of the external abdomen and on the position of internal organs. Moreover, whereas the common practice in phantom-based dose reconstruction is to select a representative phantom using age, height, and weight as discriminant factors for any treatment scenario, our analysis on abdominal tumor treatment for children shows that the most relevant features are weight and the anterior-posterior and left-right abdominal diameters.