Variable length genetic algorithm with continuous parameters optimization of beam layout in proton therapy

Proton therapy is a modality in fast development. Characterized by a maximum dose deposition at the end of the proton trajectory followed by a sharp fall-off, proton beams can deliver a highly conformal dose to the tumor while sparing organs at risk and surrounding healthy tissues. New treatment planning systems based on spot scanning techniques can now propose multi-field optimization. However, in most cases, this optimization only processes the field fluences whereas the choice of ballistics (field geometry) is left to the oncologist and medical physicist. In this work, we investigate a new optimization framework based on a genetic approach. This tool is intended to explore new irradiation schemes and to evaluate the potential of actual or future irradiation systems. We propose to optimize simultaneously the target points and beam incidence angles in a continuous manner and with a variable number of beams. No \textit{a priori} technological constraints are taken into account, \textit{i.e.}~the beam energy values, incidence directions and target points are free parameters. The proposed algorithm is based on a modified version of classical genetic operators: mutation, crossover and selection. We use the real coding associated with random perturbations of the parameters to obtain a continuous variation of the potential solutions. We also introduce a perturbation in the exchange points of the crossover to allow variations of the number of beams. These variations are controlled by introducing a beam fluence lower limit. In this paper, we present a complete description of the algorithm and of its behaviour in an elementary test case. The proposed method is finally assessed in a clinically-realistic test case.

Towards the optimization of ballistics in proton therapy using genetic algorithms: implementation issues

The dose delivered to the planning target volume by proton beams is highly conformal, sparing organs at risk and normal tissues. New treatment planning systems adapted to spot scanning techniques have been recently proposed to simultaneously optimize several fields and thus improve dose delivery. In this paper, we investigate a new optimization framework based on a genetic algorithm approach. This tool is intended to make it possible to explore new schemes of treatment delivery, possibly with future enhanced technologies. The optimization framework is designed to be versatile and to account for many degrees of freedom, without any {\it a priori} technological constraint. To test the behavior of our algorithm, we propose in this paper, as an example, to optimize beam fluences, target points and irradiation directions at the same time. The proposed optimization routine takes typically into account several thousands of spots of fixed size. The evolution is carried out by the three standard genetic operators: mutation, crossover and selection. The figure-of-merit (or fitness) is based on an objective function relative to the dose prescription to the tumor and to the limits set for organs at risk and normal tissues. Fluence optimization is carried out via a specific scheme based on a plain gradient with analytical solution. Several specific genetic algorithm issues are addressed: (i) the mutation rate is tuned to balance the search and selection forces, (ii) the initial population is selected using a bootstrap technique and (iii) to scale down the computation time, dose calculations are carried out with a fast analytical ray tracing method and are multi-threaded. In this paper implementation issues of the optimization framework are thoroughly described. The behavior of the proposed genetic algorithm is illustrated in both elementary and clinically-realistic test cases.