A new optimization method using a compressed sensing inspired solver for real-time LDR-brachytherapy treatment planning

Automatic planning for head and neck seed implant brachytherapy based on deep convolutional neural network dose engine

BACKGROUND

Seed implant brachytherapy (SIBT) is an effective treatment modality for head and neck (H&N) cancers; however, current clinical planning requires manual setting of needle paths and utilizes inaccurate dose calculation algorithms.

PURPOSE

This study aims to develop an accurate and efficient deep convolutional neural network dose engine (DCNN-DE) and an automatic SIBT planning method for H&N SIBT.

METHODS

A cohort of 25 H&N patients who received SIBT was utilized to develop and validate the methods. The DCNN-DE was developed based on 3D-unet model. It takes single seed dose distribution from a modified TG-43 method, the CT image and a novel inter-seed shadow map (ISSM) as inputs, and predicts the dose map of accuracy close to the one from Monte Carlo simulations (MCS). The ISSM was proposed to better handle inter-seed attenuation. The accuracy and efficacy of the DCNN-DE were validated by comparing with other methods taking MCS dose as reference. For SIBT planning, a novel strategy inspired by clinical practice was proposed to automatically generate parallel or non-parallel potential needle paths that avoid puncturing bone and critical organs. A heuristic-based optimization method was developed to optimize the seed positions to meet clinical prescription requirements. The proposed planning method was validated by re-planning the 25 cases and comparing with clinical plans.

RESULTS

The absolute percentage error in the TG-43 calculation for CTV V100 and D90 was reduced from 5.4% and 13.2% to 0.4% and 1.1% with DCNN-DE, an accuracy improvement of 93% and 92%, respectively. The proposed planning method could automatically obtain a plan in 2.5 ± 1.5 min. The generated plans were judged clinically acceptable with dose distribution comparable with those of the clinical plans.

CONCLUSIONS

The proposed method can generate clinically acceptable plans quickly with high accuracy in dose evaluation, and thus has a high potential for clinical use in SIBT.

Centre-specific autonomous treatment plans for prostate brachytherapy using cGANs

PURPOSE

In low-dose-rate prostate brachytherapy (LDR-PB), treatment planning is the process of determining the arrangement of implantable radioactive sources that radiates the prostate while sparing healthy surrounding tissues. Currently, these plans are prepared manually by experts incorporating the centre's planning style and guidelines. In this article, we develop a novel framework that can learn a centre's planning strategy and automatically reproduce rapid clinically acceptable plans.

METHODS

The proposed framework is based on conditional generative adversarial networks that learn our centre's planning style using a pool of 931 historical LDR-PB planning data. Two additional losses that help constrain prohibited needle patterns and produce similar-looking plans are also proposed. Once trained, this model generates an initial distribution of needles which is passed to a planner. The planner then initializes the sources based on the predicted needles and uses a simulated annealing algorithm to optimize their locations further.

RESULTS

Quantitative analysis was carried out on 170 cases which showed the generated plans having similar dosimetry to that of the manual plans but with significantly lower planning durations. Indeed, on the test cases, the clinical target volumes achieving [Formula: see text] of the prescribed dose for the generated plans was on average [Formula: see text] ([Formula: see text] for manual plans) with an average planning time of [Formula: see text] min ([Formula: see text] min for manual plans). Further qualitative analysis was conducted by an expert planner who accepted [Formula: see text] of the plans with some changes ([Formula: see text] requiring minor changes & [Formula: see text] requiring major changes).

CONCLUSION

The proposed framework demonstrated the ability to rapidly generate quality treatment plans that not only fulfil the dosimetric requirements but also takes into account the centre's planning style. Adoption of such a framework would save significant amount of time and resources spent on every patient; boosting the overall operational efficiency of this treatment.

Optimization in treatment planning of high dose-rate brachytherapy - Review and analysis of mathematical models

Treatment planning in high dose-rate brachytherapy has traditionally been conducted with manual forward planning, but inverse planning is today increasingly used in clinical practice. There is a large variety of proposed optimization models and algorithms to model and solve the treatment planning problem. Two major parts of inverse treatment planning for which mathematical optimization can be used are the decisions about catheter placement and dwell time distributions. Both these problems as well as integrated approaches are included in this review. The proposed models include linear penalty models, dose-volume models, mean-tail dose models, quadratic penalty models, radiobiological models, and multiobjective models. The aim of this survey is twofold: (i) to give a broad overview over mathematical optimization models used for treatment planning of brachytherapy and (ii) to provide mathematical analyses and comparisons between models. New technologies for brachytherapy treatments and methods for treatment planning are also discussed. Of particular interest for future research is a thorough comparison between optimization models and algorithms on the same dataset, and clinical validation of proposed optimization approaches with respect to patient outcome.

Knowledge-based inverse treatment planning for low-dose-rate prostate brachytherapy

PURPOSE

Permanent low-dose-rate brachytherapy is a widely used treatment modality for managing prostate cancer. In such interventions, treatment planning can be a challenging task and requires experience and skills of the planner. We developed a novel knowledge-based (KB) optimization method based on previous treatment plans. The purpose of this method was to generate clinically acceptable plans that do not require extensive manual adjustments in clinical scenarios.

METHODS

Objective functions used in current inverse planning methods are preferably based on spatial invariant dose objectives rather than spatial dose distributions. Therefore, they are prone to return suboptimal plans resulting in time consuming plan adjustments. To overcome this limitation, a KB approach is introduced. The KB model uses the dose distributions of previous clinical plans projected onto a standardized geometry. From those standardized distributions a template plan is generated. The treatment plans were optimized with an in-house developed planning system by solving a constraint inverse optimization problem that mimics the projected template dose plan constraint to DVH metrics. The method is benchmarked under an IRB-approved retrospective study by comparing optimization time, dosimetric performance, and clinical acceptability against current clinical practice. The quality of the KB model is evaluated with a Turing test.

RESULTS

The KB model consists of five high-quality treatment plans. Those plans were selected by one of our experts and showed all desired dosimetric features. After generating the model treatment plans were created with one run of the optimizer for the remaining 20 patients. The optimization time including needle optimization ranged from 6 to 29 s. Based on a Wilcoxon signed rank test the new plans are dosimetrically equivalent to current clinical practice. The Turing test showed that the proposed method generates plans that are equivalent to current clinical practice and that the dose prediction drives the optimization to achieve high-quality treatment plans.

CONCLUSIONS

This study demonstrated that the proposed KB model was able to capture user-specific features in isodose lines which can be used to generate acceptable treatment plans with a single run of the optimization engine in under a minute. This could potentially reduce the time in the operating room and the time a patient is under anesthesia.

DVH-Based Inverse Planning Using Monte Carlo Dosimetry for LDR Prostate Brachytherapy

PURPOSE

Inverse planning is an integral part of modern low-dose-rate brachytherapy. Current clinical planning systems do not exploit the total dose information and largely use the American Association of Physicists in Medicine TG-43 dosimetry formalism to ensure clinically acceptable planning times. Thus, suboptimal plans may be derived as a result of TG-43-related dose overestimation and nonconformity with dose distribution requirements. The purpose of this study was to propose an inverse planning approach that can improve planning quality by combining dose-volume information and precision without compromising the overall execution times.

METHODS AND MATERIALS

The dose map was generated by accumulating precomputed Monte Carlo (MC) dose kernels for each candidate source implantation site. The MC computational burden was reduced by using graphics processing unit acceleration, allowing accurate dosimetry calculations to be performed in the intraoperative environment. The proposed dose-volume histogram (DVH) fast simulated annealing optimization algorithm was evaluated using clinical plans that were delivered to 18 patients who underwent low-dose-rate prostate brachytherapy.

RESULTS

Our method generated plans in 37.5 ± 3.2 seconds with similar prostate dose coverage, improved prostate dose homogeneity of up to 6.1%, and lower dose to the urethra of up to 4.0%.

CONCLUSIONS

A DVH-based optimization algorithm using MC dosimetry was developed. The inclusion of the DVH requirements allowed for increased control over the optimization outcome. The optimal plan's quality was further improved by considering tissue heterogeneity.

Determining optimal eluter design by modeling physical dose enhancement in brachytherapy

PURPOSE

In situ drug release concurrent with radiation therapy has been proposed to enhance the therapeutic ratio of permanent prostate brachytherapy. Both brachytherapy sources and brachytherapy spacers have been proposed as potential eluters to release compounds, such as nanoparticles or chemotherapeutic agents. The relative effectiveness of the approaches has not been compared yet. This work models the physical dose enhancement of implantable eluters in conjunction with brachytherapy to determine which delivery mechanism provides greatest opportunity to enhance the therapeutic ratio.

MATERIALS AND METHODS

The combined effect of implanted eluters and radioactive sources were modeled in a manner that allowed the comparison of the relative effectiveness of different types of implantable eluters over a range of parameters. Prostate geometry, source, and spacer positions were extracted from treatment plans used for 125 I permanent prostate implants. Compound concentrations were calculated using steady-state solution to the diffusion equation including an elimination term characterized by the diffusion-elimination modulus (ϕb ). Does enhancement was assumed to be dependent on compound concentration up to a saturation concentration (csat ). Equivalent uniform dose (EUD) was used as an objective to determine the optimal configuration of eluters for a range of diffusion-elimination moduli, concentrations, and number of eluters. The compound delivery vehicle that produced the greatest enhanced dose was tallied for points in parameter space mentioned to determine the conditions under whether there are situations where one approach is preferable to the other.

RESULTS

The enhanced effect of implanted eluters was calculated for prostate volumes from 14 to 45 cm3 , ϕb from 0.01 to 4 mm-1 , csat from 0.05 to 7.5 times the steady-state compound concentration released from the surface of the eluter. The number of used eluters (ne ) was simulated from 10 to 60 eluters. For the region of (csat , Φ)-space that results in a large fraction of the gland being maximally sensitized, compound eluting spacers or sources produce equal increase in EUD. In the majority of the remaining (csat , Φ)-space, eluting spacers result in a greater EUD than sources even where sources often produce greater maximal physical dose enhancement. Placing eluting implants in planned locations throughout the prostate results in even greater enhancement than using only source or spacer locations.

CONCLUSIONS

Eluting brachytherapy spacers offer an opportunity to increase EUD during the routine brachytherapy process. Incorporating additional needle placements permits compound eluting spacer placement independent of source placement and thereby allowing a further increase in the therapeutic ratio. Additional work is needed to understand the in vivo spatial distribution of compound around eluters, and to incorporate time dependence of both compound release and radiation dose.

A radiobiology-based inverse treatment planning method for optimisation of permanent l-125 prostate implants in focal brachytherapy

Treatment plans for ten patients, initially treated with a conventional approach to low dose-rate brachytherapy (LDR, 145 Gy to entire prostate), were compared with plans for the same patients created with an inverse-optimisation planning process utilising a biologically-based objective. The 'biological optimisation' considered a non-uniform distribution of tumour cell density through the prostate based on known and expected locations of the tumour. Using dose planning-objectives derived from our previous biological-model validation study, the volume of the urethra receiving 125% of the conventional prescription (145 Gy) was reduced from a median value of 64% to less than 8% whilst maintaining high values of TCP. On average, the number of planned seeds was reduced from 85 to less than 75. The robustness of plans to random seed displacements needs to be carefully considered when using contemporary seed placement techniques. We conclude that an inverse planning approach to LDR treatments, based on a biological objective, has the potential to maintain high rates of tumour control whilst minimising dose to healthy tissue. In future, the radiobiological model will be informed using multi-parametric MRI to provide a personalised medicine approach.