Treatment planning optimization by conjugate gradients and simulated annealing methods in stereotactic radiosurgery

This paper presents a new optimization method of treatment planning in linac stereotactic radiosurgery.

On a workstation integrating x-rays, computed tomography (CT), magnetic resonance imaging (MRI), and digital subtracted angiography (DSA) images, we first determine the outlines of the target volume and surrounding healthy tissues to spare. To achieve complete optimization of the treatment plans, this method decomposes the optimization process in two steps. The position of the isocenters and the diameter of the collimators are first deduced by a conjugate gradients method, from the position and size of ellipsoids or spheres modeling the target volume. The other irradiation parameters, such as the isocenter dose, the aperture, and the weight of each irradiation plane and of their irradiation sectors are finally deduced by a simulated annealing optimization algorithm.

The system can perform multitarget/multisector treatment plans that are automatically obtained in a satisfactory time (as a rule, 20 min for a two-target irradiation), much faster than the time needed for a manual treatment planning. We present the results in two cases: the simulation of a single-target treatment and a two-target real treatment with constraints. In these two cases, we can control the dose received by target and sensitive volumes.

This method achieves an excellent conformation of the estimated isodose curves with the outlines of the target volume, which allows us to avoid the surrounding healthy tissues, thanks to the different weighting factors given on each volume concerned according to the importance we grant to each of them.