Thèse soutenue

Optimisation quadratique et géométrique de problèmes de dosimétrie inverse
FR  |  
EN
Accès à la thèse
Auteur / Autrice : Yann Menguy
Direction : Philippe Cinquin
Type : Thèse de doctorat
Discipline(s) : Mathématiques appliquées
Date : Soutenance en 1996
Etablissement(s) : Grenoble 1

Résumé

FR

La radiothérapie utilise les densités massiques d'énergie (appelées doses) que délivrent les rayonnements ionisants, afin de détruire des tumeurs : cette technique présente l'énorme avantage d'être non invasive. Elle nécessite cependant une grande prudence car lors d'une irradiation, toutes les parties du corps reçoivent une certaine dose. Le travail de dosimétrie est confié aux radiophysiciens qui, par essais successifs, déterminent une balistique de traitement permettant de concentrer la dose sur la tumeur en épargnant autant que possible les tissus sains. Le but de notre travail, qui se limite à la radiothérapie de la prostate, est d'automatiser cette étape. Le problème de dosimétrie directe consiste, pour une balistique donnée, à calculer la dose en un point du corps. Nous rappelons les principales méthodes déjà existantes, et détaillons l'une d'elles, la méthode de Clarkson. Nous construisons ensuite, pour un faisceau circulaire, un modèle basé sur le lissage par fonctions spline de données expérimentales tridimensionnelles et le comparons à une interpolation multilinéaire sur les données. Nous traitons, afin de cibler au mieux la prostate, le problème posé par la recherche de la plus petite boule englobant un ensemble de points. Nous développons ainsi un algorithme basé sur des notions géométriques, généralisation au cas multidimensionnel de l'algorithme de Chrystal-Peirce. Nous montrons que cet algorithme converge en un nombre fini d'itérations et effectuons de nombreux tests numériques afin de prouver son efficacité. Nous exposons enfin son côté analytique et montrons qu'il s'agit d'un algorithme de sous-gradient. La dernière partie de cette thèse est consacrée à l'optimisation dosimétrique. En nous plaçant dans le cadre d'un protocole multi-faisceaux, nous définissons ce problème comme minimisation d'une fonctionnelle sur un ensemble de contraintes et comparons les résultats pour différentes géométries de faisceaux. De manière générale, nous notons une réelle amélioration par rapport aux traitements actuels