Les calculs de dose en curiethérapie à bas débit reposent depuis 1995 sur un formalisme qui considère le milieu irradié comme étant homogène et constitué d’eau afin de mieux accommoder les temps de calcul clinique. Ce travail présente une amélioration de ce formalisme dosimétrique issu du protocole TG-43 de l’AAPM (American Association of Physicists in Medicine), utilisé actuellement dans les systèmes de planification de traitement (SPT) clinique. Avec l’ajout d’un tracé radiologique pour tenir compte des hétérogénéités et en séparant la dose primaire de la dose diffusée, il est possible de raffiner les calculs dosimétriques. Or, cette modification au calcul le rendant plus complexe était coûteuse en temps d’exécution jusqu’à très récemment. Elle se voit matérialisée dans ce travail avec l’arrivée récente de la technologie GPGPU, les calculs scientifiques sur périphériques graphiques. En exploitant le parallélisme des calculs de dose en curiethérapie à bas débit pour les implants permanents, l’algorithme baptisé TG-43-RT pour TG-43 avec tracé de rayons (Ray-Tracing en anglais) permet d’obtenir des facteurs d’accélération de l’ordre de 103 par rapport au CPU. Seulement 0.5 s par source est requis dans une géométrie de 1003 voxels, et les dosimétries d’anatomies hétérogènes sont sensiblement améliorées. Le TG-43-RT corrige les dépôts de dose en aval des hétérogénéités et réduit les effets indésirables de l’atténuation inter-sources. Des écarts de dose de plus de 80% au-delà d’une calcification avaient de quoi influencer les HDV, or, le TG-43-RT ramène ces déviations en-deçà de l’incertitude de 8.7% (2σ) concédée par les SPT, tel qu’établi dans le TG-138. Les simulations Monte Carlo ont été employées pour servir de référence absolue à la méthode développée et quantifier les améliorations dosimétriques relatives à l’actuel calcul de dose. L’algorithme a été testé d’un point de vue clinique dans un fantôme de prostate avec des sources d’125I et dans un fantôme de sein avec des sources de 103Pd pour unifier les corrections apportées. / Brachytherapy dose calculations have been relying since 1995 on a formalism that considers the whole geometry as a homogeneous water tank. This gives the opportunity to compute dose distributions within a reasonable clinical timeframe, but with considerable approximations that can influence or even change the treatment. This work presents an upgraded version of this formalism derived from the TG-43 protocol (AAPM), currently used in actual treatment planning systems. The new algorithm includes a primary and scatter dose separation using ray-tracing operation to account for heterogeneities through the medium. These modifications to the dose calculation are very power consuming and too long for the clinical needs when executed on modern CPUs. With the GPGPU technology, a GPU-driven algorithm allows a complex handling of the anatomic heterogeneities in the dose calculation and keeps execution times below 0.5 s/seed. The algorithm was named TG-43-RT for TG-43 with Ray-Tracing and accelerations factors of three order of magnitude were obtained over a regular CPU implementation. The TG-43-RT algorithm adequatly corrects dose deposition along heterogeneities and it reduces interseed attenuation effects. Dose is scored in the medium instead of water, which also applies a severe correction for high heterogeneous medium like the breast. Deviations of more than 80% in dose deposition were obtained over calcification inside the prostate when compared to Monte Carlo simulations. This kind of deviation influences DVH shape. However, the TG-43-RT was keeping those deviations within the 8.7% uncertainty range (2σ) associated with a regular TG-43 evaluation as stated in the TG-138. Monte Carlo simulations were used as a gold standard to get an absolute dose calculation reference. The algorithm was also tested in a prostate phantom with 125I seeds and in a breast phantom with 103Pd seeds to mimic real anatomic geometries and unify the applied corrections.

Les cellules cancéreuses du sein sont souvent réfractaires aux morts radio-induites, diminuant ainsi l'efficacité des traitements de radiothérapie. L'irradiation d'atomes radiosensibilisants, comme le brome, l'iode et le platine, peut augmenter la concentration locale de radicaux libres dans les cellules. Cet effet peut être obtenu en amplifiant la création d'électrons Auger générés par un mécanisme résonant, ou non résonant impliquant des électrons secondaires créés suite à l'ionisation de l'eau par la radiation. Le but de cette étude était d'évaluer et de comparer des dérivés de porphyrines cationiques bromées comme nouveaux radiosensibilisateurs potentiels pour le traitement du cancer du sein. La capacité radiosensibilisante de trois porphyrines bromées a été mesurée et comparée avec celle de leur homologue non bromé par test de formation de colonies sur les cellules cancéreuses mammaires humaines MDA-MB-231. Seulement les deux porphyrines portant des atomes de bromes en périphérie ont démontré un pouvoir radiosensibilisant, soit le CuPBr[indice inférieur 8](PyPrBr[indice inférieur 2])[indice inférieur 4] et le CuP(PyPreBr[indice inférieur 2])[indice inférieur 4]. Le CupBr[indice inférieur 8](PyPrBr[indice inférieur 2)[indice inférieur 4] est un dérivé de porphyrine avec 16 bromes répartis équitablement sur les chaînes propyles périphériques et sur le centre de la molécule. Ce composé à une concentration de 10 [micro]M, en combinaison avec une dose de radiation de 2 Gy, a conduit à une radiosensibilisation de 76,6% avec des rayons-x de 14,8 keV, de 53,8% avec des rayons-x de 70,1 keV et de 48,4% avec des rayons gamma ([indice supérieur 60]Co). Son effet radiosensibilisateur a été dépendant de la dose de radiation pour les 3 sources entre 1 et 3,5 Gy et de la concentration entre 5 et 10 [micro]M au [indice supérieur 60]Co. L'effet radiosensibilisateur du CuPBr[indice inférieur 8](PyPrBr[indice inférieur 2])[indice inférieur 4] mesuré suite à l'irradiation de 2 Gy au [indice supérieur 60]Co a été indépendant de la présence du métal au centre de la porphyrine bromée et semblable à celui du CuP(PyPreBr[indice inférieur 2])[indice inférieur 4], un dérivé portant 8 bromes uniquement sur ses groupements périphériques. Alors qu'aucun effet radiosensibilisateur n'a été observé avec la troisième porphyrine, le CuPBr[indice inférieur 8](PyM)[indice inférieur 4] dont les huit bromes ont été ajoutés uniquement sur le centre de la molécule. D'autre part l'ajout du CuPBr[indice inférieur 8](PyPrBr[indice inférieur 2)[indice inférieur 4] sur les cellules cancéreuses mammaires de souris MC7-L1, combiné avec 2 Gy de rayons-x de 14,8 keV, a amené une radiosensibilisation de 41,4%. Dans le cas des cellules MDA-MB-231, la radiosensibilisation des porphyrines bromées en périphérie a corrélé avec une accumulation cellulaire de drogue supérieure à leur homologue non bromé de 2,5 fois pour le CuPBr[indice inférieur 8](PyPrBr[indice inférieur 2])[indice inférieur 4], la porphyrine portant 16 bromes et de 4 fois pour le CuP(PyPreBr[indice inférieur 2)[indice inférieur 4], la porphyrine avec 8 bromes uniquement en périphérie. Cette augmentation de la radiosensibilité n'a cependant pas été associée à une formation d'espèces toxiques issues de la dégradation des porphyrines lors de l'irradiation, démontré par analyse des porphyrines irradiées en solution et traitement des cellules avec les radiosensibilisateurs pré-irradiés seuls. La radiosensibilisation par ces composés bromés n'a pas non plus conduit à une externalisation des phosphatidylsérines, ni à une perméabilisation de la membrane cellulaire et aucune variation du potentiel de la membrane mitochondriale n'a été mesurée. Ces observations ont permis de conclure que la mort cellulaire principalement associée à la radio sensibilisation n'était pas l'apoptose. En conclusion, l'étude de ces porphyrines bromées en périphérie a démontré que cette classe de radiosensibilisateurs pourrait avoir un bon potentiel d'application pour traiter les cellules cancéreuses mammaires.

Des métastases cérébrales vont se développer chez 10 à 30% des patients atteints de cancer. La radiothérapie fait partie des possibilités de traitement, et ceci même si les dommages induits au cerveau par des rayonnements ionisants sont potentiellement importants. Nous proposons l’utilisation de l’Imagerie par Résonance Magnétique fonctionnelle (IRMf) sur le rat pour mieux comprendre ces effets. Ce mémoire traite de la mise en place d’un tel protocole d’IRMf. Les principaux points abordés sont la préparation de l’animal, les différentes insultes et stimulations sensorielles possibles ainsi que la méthode d’acquisition. Notre protocole d’insulte hyperoxique permet de déceler des dommages physiques d’origine vasculaire suite à une intense irradiation dans le cerveau du rat. Toutefois, la même procédure associée à une stimulation mécanique de la patte arrière de l’animal n’amène pas de changement observable à l’IRMf sur un sujet sain. Malgré tout, ce type de stimulation induit une réponse respiratoire, même sous anesthésie d’isoflurane. Une telle méthode n’est donc pas adéquate à l’étude d’animaux anesthésiés, surtout ceux dont la réponse cérébrale pourra avoir été réduite par une irradiation. Quelques améliorations et modifications du protocole seraient possiblement à même de permettre une mesure reproductible de la réponse d’IRMf à une stimulation sensorielle. Le présent mémoire décrit les tentatives de mise en place d’une stimulation sensorielle donnant lieu à une activation IRMf reproductible et localisée. De plus, un protocole de traitement d’image adapté au petit animal ainsi qu’une implémentation de la méthode keyhole ont été mis en place. L’insulte hyperoxique et ses effets sur le cerveau de rat ont été explorés plus en détail.

Mémoire numérisé par la Direction des bibliothèques de l'Université de Montréal.

Dans une démarche d'évaluation rigoureuse des innovations technologiques de la radiothérapie, notre travail de thèse a été construit autours de 2 thématiques : 1) Evaluation clinique des résultats de la radiothérapie conformationnelle avec modulation d'intensité (RCMI) appliquée au traitement des cancers ORL. La première étude rapporte les conclusions d'un enregistrement prospectif monocentrique du devenir au long cours des patients ORL traités par RCMI. Le contrôle locorégional est excellent et les niveaux de toxicité faibles. Des recommandations dosimétriques (dose aux parotides) adaptées à nos pratiques ont été établies. La deuxième étude rapporte les résultats d'une comparaison nationale multicentrique de la qualité de vie des patients traités pour un cancer ORL, soit par radiothérapie conventionnelle, soit par RCMI. Le bénéfice de la RCMI est remarquable dans les domaines du confort de la sphère buccale et de la prise alimentaire orale. 2) Evaluation de l'apport de l'imagerie embarquée (Cone Beam CT ou CBCT) dans la prise en charge des patients traités par RCMI pour un cancer ORL. Un premier projet propose, sur un panel de patients ORL traités par RCMI, le contrôle qualité de la distribution de la dose réellement délivrée. Le retentissement dosimétrique des variations anatomiques est discuté. L'apport d'outils dosimétriques spécifiques est étudié. Le second projet a pour objectif de quantifier la mobilité relative des régions anatomiques qui composent la sphère ORL et d'étudier l'impact dosimétrique lié au choix de différentes procédures de recalage. Une recommandation sur les modalités d'utilisation des images CBCT a pu être proposée. / Numerous and exciting technological innovations were recently developed in radiotherapy. We aimed to assess benefits in two specific fields. 1) Clinical results of Intensity Modulated Radiotherapy (IMRT) applied to the treatment of Head and Neck (H&N) patients. The first study was a long-term monocentric prospective registration of all H&N patients treated with IMRT in our institution. Locoregional control was excellent and toxicities limited. Recurrences were in-field. Dosimetric recommendations (parotids mean dose) were established. The second study assessed the impact of IMRT on health-related quality of life for H&N patients through a multicentric matched-pair comparison with conventional radiotherapy. Outstanding benefits were observed particularly in the fields of salivary dysfunction and oral discomfort. 2) Contribution of Image Guided Radiotherapy (IGRT) in the management of H&N patients treated with IMRT. The first study was a monitoring of delivered dose, using 3D dose recalculation from Megavoltage Cone-Beam CT (CBCT), as a quality assurance measure of a panel of H&N IMRT patients aligned with IGRT. Dosimetric consequences of anatomical changes were assessed. Contribution of color-coded MVCBCT dose-difference maps was studied. The aim of the second study was to quantify the inherent relative mobility between anatomic regions of the H&N area and to assess the dosimetric impact of several different matching procedures. Recommendations for the use of CBCT images in a daily practice were established.

Le but de cette étude exploratrice est d'investiguer le rôle potentiellement paradoxal de la radiothérapie sur le caractère invasif des glioblastomes multiformes (GBM). Utilisant le modèle murin de tumeur astrocytaire F98/Fisher, nous avons dans un premier temps irradié soit le cerveau (BrainlR) de l'animal, les cellules cancéreuses implantées (F98IR), ou ces deux composantes (Brain+F98IR), afin d'étudier l'effet des radiations sur le développement tumoral. De plus, un regard particulier a été porté sur la cinétique de moléculaires inflammatoire et pro-migratoires induite par les radiations. Finalement, une approche pharmacologique a été utilisée dans le but de contrôler la réaction inflammatoire radio-induite pour diminuer les effets potentiellement néfastes des radiations sur l'augmentation de la migration des cellules F98 dans le cerveau de rat. Les analyses histologiques des cerveaux portant une tumeur F98 ont démontré une augmentation de l'infiltration pour les rats ayant été irradiés avant l'implantation. Ce même groupe a démontré une augmentation de létalité et une immunoréactivité accrue pour la nestine, illustrant une différence phénotypique. De façon opposée, le groupe où les cellules F98 sont irradiées avant l'implantation n'a démontré aucune différence flagrante dans leur profil de migration comparativement au groupe contrôle. Les analyses moléculaires ont illustré la présence d'une réaction inflammatoire suite à une exposition du cerveau aux radiations. Une stimulation de la production d'IL-1[beta] et de TGF-[beta] aptes l'exposition, supporte cette observation. De plus, une augmentation des niveaux de collagénase (MMP-2), permettant l'infiltration cellulaire, a été démontrée pour une période allant jusqu'à 20 jours post-irradiation. L'analyse de prostanoïdes pro-inflammatoires a révélé des niveaux augmentés des prostaglandines (PGE [indice inférieur 2] /PGD [indice inférieur 2] ) jusqu'à une période de 20 jours également. A la suite de ces résultats d'analyse moléculaire, nous avons choisi d'utiliser un anti-inflammatoire non stéroïdien (AINS) pour diminuer cette réaction due aux radiations pour ainsi réprimer la formation de molécules potentiellement néfastes. Le traitement au Meloxicam des rats avec des doses de 1 mg/kg quotidiennement, s'est traduit par une diminution du caractère invasif des tumeurs implantées dans un cerveau irradié. Ces données sur l'utilisation du Meloxicam comme approche pharmacologique sont encourageantes dans l'optique de diminuer l'invasion radio-induite. Notre étude illustre donc l'importance de se soucier de l'irradiation des tissus sains lors des traitements de radiothérapie pouvant jouer un rôle lors de la récurrence rapide observée en clinque chez les patients atteint de GBM. De plus, l'utilisation d'AINS pourrait s'avérer un outil efficace pour augmenter l'espérance de vie des patients atteints de glioblastomes.

Thèse ou mémoire avec insertion d'articles / Le contrôle tumoral en radiothérapie nécessite une exposition des cellules cancéreuses au rayonnement ionisant tout en limitant l'exposition des cellules saines pour minimiser les effets secondaires indésirables. Pour atteindre cet objectif, des techniques sophistiquées telles que l'IMRT, le VMAT, le SRS et le SBRT ont été développées, utilisant plusieurs champs d'irradiation de forme variable et de petites tailles. Le succès de ces techniques repose sur la caractérisation des petits champs, pour lesquels il n'existe pas encore de norme établie quant au type de dosimètre à utiliser. En effet, aucun des dosimètres couramment utilisés en radiothérapie ne permet de mesurer la dose avec une précision inférieure à 2 % pour ces champs de petite taille. De façon générale, ce sont les détecteurs ayant une équivalence dosimétrique à l'eau ainsi qu'un petit volume de détection qui performent le mieux pour la dosimétrie des petits champs. Cette thèse vise donc à développer une nouvelle génération de calorimètres multi-point spécifiquement adapté à la dosimétrie des petits champs en radiothérapie possédant à la fois une équivalence dosimétrique à l'eau ainsi qu'un petit volume de détection. Ces calorimètres utilisent une fibre optique à réseaux de Bragg comme thermomètre, permettant des mesures de haute résolution pour les petits champs de traitement. Cette thèse se divise en 4 volets distincts, soit la caractérisation et la maximisation du changement de longueur d'onde de Bragg produit par la radiation ; la caractérisation des mécanismes menant à un changement de longueur d'onde de Bragg induit par la radiation ; la correction de la dépendance aux variations de température ambiante et finalement, la conception d'un dosimètre multi-point utilisant une fibre optique à réseaux de Bragg ayant une résolution spatiale adaptée à la dosimétrie standard et à la dosimétrie des petits champs. D'abord, la caractérisation a permis de déterminer que, pour une fibre à réseau de Bragg recouverte de polymère, la longueur d'onde de Bragg varie linéairement avec la dose et cette variation est indépendante du débit de dose (2.8-11.6 Gy/min) et de l'énergie du faisceau de photons ou d'électrons (6 MV, 18 MV, 6 MeV, 9 MeV, 12 MeV, 18 MeV). Les facteurs permettant de maximiser cette variation de longueur d'onde et de diminuer le bruit de mesure ont ensuite été identifiés. Le type de revêtement a été identifiée comme étant un paramètre important pour maximiser le signal, plus particulièrement ses propriété thermiques (coefficient d'expansion thermique, capacité thermique massique), son module de Young et sa solidité d'adhésion ainsi que le type de colle maximisant cette dernière. La taille du revêtement et de la fibre optique contribuent également à l'optimisation du détecteur. L'optimisation des paramètres de réseaux ont également permis de minimiser le bruit de mesure. En ce qui concerne les mécanismes menant à un changement de longueur d'onde de Bragg, il a d'abord été déterminé que la réponse du détecteur est indépendante du dopage du cœur de la fibre et provient du revêtement de polymère puisque sans ce dernier, aucun signal n'est mesuré. Il a ensuite été démontré que le détecteur développé est un calorimètre et non un détecteur de défauts radio-induits, contrairement aux dosimètres à réseaux de Bragg généralement utilisés dans le secteur de l'énergie nucléaire. Cette démonstration est réalisée en validant un modèle théorique basé sur la théorie des calorimètres et en confirmant le comportement thermique suite à une irradiation à l'aide de simulations thermiques. Quatre techniques ont été comparées afin de corriger la dépendance aux variations de température ambiante, dont la technique standard utilisée en calorimétrie ainsi qu'une nouvelle technique de gradient de température interpolé pour effectuer une dosimétrie multi-point. Cette dernière technique se révèle la plus performante tout en permettant une correction de température en temps réel, réduisant les erreurs de mesure de 200% à environ 10% pour une irradiation de 20 Gy. Deux prototypes ont été développés : un ayant une résolution spatiale optimisée pour la dosimétrie standard et un ayant une résolution spatiale optimisée pour la dosimétrie des petits champs. Le premier détecteur est composé de vingt (20) réseaux d'environ 4 mm de long répartis sur 20 cm fixé à une plaque de PMMA de 5.5 x 107.5 x 205 mm³ permettant d'obtenir une erreur statistique de 0.03 pm sur les points de données limitant ainsi la dose détectable à 0.4 Gy. Le second détecteur se compose d'un fibre optique de 80 µm contenant trente (30) réseaux d'environ 1 mm de long fixés à l'intérieur d'un cylindre de PMMA de 0.63 cm de diamètre sur une longueur de 30 cm permettant d'obtenir la même erreur statistique et la même dose détectable que précédemment. Ce dernier prototype a également permis de mesurer un profil de dose d'un faisceau de 2 x 2 cm² ayant une énergie de 6 MV avec une différence relative moyenne de 1.8% (en excluant la région de pénombre) pour une irradiation de 12.37 Gy. Ainsi la thèse présente le premier calorimètre multi-point à réseaux de Bragg spécifiquement adapté à la dosimétrie des petits champs en radiothérapie. Ce type de détecteur pourrait s'avérer extrêmement bénéfique pour la dosimétrie des petits champs, et présenter également un fort potentiel pour des applications en IRM-LINAC et en radiothérapie FLASH, compte tenu de la grande résistance des fibres de silice aux radiations. / Radiotherapy relies on precise control of tumor exposure to ionizing radiation while limiting the dose to healthy tissues to minimize undesirable side effects. Sophisticated techniques such as IMRT, VMAT, SRS, and SBRT have been developed, utilizing multiple small and variableshaped radiation fields. The success of these techniques depends on the characterization of small fields, for which no standard dosimeter has been established. Commonly used dosimeters in radiotherapy do not provide dose measurements with accuracy below 2% for small fields. Generally, dosimeters with near water-equivalence and small detection volumes perform best for small-field dosimetry. This thesis focuses on developing a new generation of multi-point calorimeters specifically adapted to the dosimetry of small fields in radiotherapy, possessing both water dosimetric equivalence and a small detection volume. These calorimeters use fiber Bragg gratings as thermometers, enabling high-resolution dose measurements for small fields radiotherapy treatments. The thesis comprises four distinct aspects : characterizing and maximizing the Bragg wavelength change induced by radiation ; understanding the mechanisms causing the Bragg wavelength shift ; compensating for ambient temperature variations ; and designing a multipoint dosimeter using gratings with a suitable spatial resolution for standard and small field dosimetry. The characterization revealed that the Bragg wavelength of a polymer-embedded fiber Bragg grating varies linearly with the dose and this variation is independent of the dose rate (2.8-11.6 Gy/min) and the beam energy (6 MV, 18 MV, 6 MeV, 9 MeV, 12 MeV, 18 MeV). Factors optimizing this wavelength variation and reducing measurement noise were identified. The type of coating, particularly its thermal properties (thermal expansion coefficient, specific heat capacity), Young's modulus, and adhesive strength play a critical role in maximizing the signal. Additionally, optimizing the fiber and coating sizes contributes to detector optimization. The optimization of gratings parameters also allowed the minimization of measurement noise. Regarding the mechanisms leading to a Bragg wavelength shift, it was first determined that the detector's response is independent of the fiber core doping and originates from the polymer coating since no signal is measured without it. The developed detector was confirmed to be a calorimeter, rather than a radiation-induced defect detector like conventional fiber Bragg gratings dosimeters used in the nuclear industry. This was demonstrated by validating a theoretical model based on calorimeter theory and confirming the thermal behavior through thermal simulations after irradiation. Four temperature compensation techniques were compared to correct ambient temperature variations, including standard calorimetric techniques and a new interpolated temperature gradient technique for multi-point dosimetry. The interpolated technique proved most effective, providing real-time temperature correction and reducing measurement errors from 200% to approximately 10% for a 20 Gy irradiation. Two prototypes were developed : one with spatial resolution optimized for standard dosimetry and the other with spatial resolution optimized for small-field dosimetry. The first detector consists of twenty (20) gratings, each approximately 4 mm long, distributed over 20 cm and attached to a PMMA plate measuring 5.5 x 107.5 x 205 mm³ , achieving a statistical error of 0.03 pm on data points, limiting the detectable dose to 0.4 Gy. The second detector comprises a 80 µm optical fiber containing thirty (30) gratings, each approximately 1 mm long, fixed inside a PMMA cylinder measuring 0.63 cm in diameter and 30 cm in length, allowing for the same statistical error and detectable dose as the previous one. This latter prototype also enabled the measurement of a dose profile for a 2 x 2 cm² beam with a 6 MV energy, showing a mean relative difference of 1.8% (excluding the penumbra region) for an irradiation of 12.37 Gy.

Thèse ou mémoire avec insertion d'articles / L'objectif principal de cette thèse est d'explorer les applications de l'imagerie par résonance magnétique multiparamétrique (IRMmp) en radio-oncologie. Les images IRM sont utilisées à différentes étapes de la radiothérapie, notamment la planification, le guidage et le suivi. L'étape de planification utilise généralement le contraste des tissus mous de l'IRM ou des images physiologiques pour délimiter les cibles. L'étape de guidage utilise le contraste des tissus mous de l'IRM pour guider les traitements de radiothérapie. L'IRM quantitative (qMRI) et l'IRMmp sont utilisées pour évaluer une réponse au traitement. Cette thèse vise à utiliser les cartes qMRI pour quantifier la réponse des tissus de la tête et du cou à l'oxygénation, réduire les artefacts et les distorsions de l'IRM, et enfin fusionner l'IRM et la tomodensitométrie (TDM). Un des objectif est d'optimiser les cartes qMRI incluant T₁ et R₂* pour quantifier la préservation de l'oxygène après avoir administré de l'oxygène pur (pO₂) aux volontaires sains. Cette étude était la première in vivo quantifiant la préservation de l'oxygène dans les tissus de la tête et du cou après l'administration de pO₂. Deux expériences distinctes ont été réalisées : (1) délivrer une oxygénation normobare (NBO) et comparer les cartes qMRI en série T₁ et R₂* acquises pendant et après la NBO avec la ligne de base (avant de délivrer la NBO), et (2 ) délivrant une oxygénation hyperbare (HBO) et comparant les cartes qMRI en série T₁ et R₂* acquises après HBO avec la ligne de base (avant de délivrer HBO). Nous avons réalisé que le niveau d'oxygène de la glande parotide et du muscle masséter augmentait de manière significative (*p* < 0,05) après l'administration du NBO pendant 15 minutes et revenait à la valeur de base environ sept minutes après le NBO. Dans l'expérience HBO, une tendance similaire a été observée. Nous n'avons pas trouvé de différence entre les images qMRI de base et post-HBO. Ensuite, le travail s'est concentrée sur les distorsions et les artefacts de l'IRM. Dans cette étude, nous avons proposé un modèle d'apprentissage profond pour prédire les distorsions géométriques spécifiques aux patients (PSgD) causées par les différences de susceptibilité entre deux tissus. Ensuite, le T₁w (ceT₁w) amélioré en contraste a été corrigé pour le PSgD. Nous avons formé notre modèle d'apprentissage profond (DL) à l'aide d'un ensemble de données de 380 sujets. Nous avons comparé les distributions de doses de radiochirurgie stéréotaxique Gamma Knife (GK-SRS) pour les patients atteints de schwannome vestibulaire, calculées à l'aide d'images corrigées et corrompues par le PSgD. Nous avons constaté qu'en moyenne D$_μ$, D$_\textup{95\%}$ et D$_\textup{min}$ en Gy ont augmenté de manière significative après correction PSgD de 16,85 à 17,25, 12,30 à 12,77 et de 8,98 à 9,92. La couverture moyenne des objectifs et l'indice de conformité de Paddick ont considérablement augmenté, passant de 0,97 à 0,98 et de 0,94 à 0,96. L'indice de gradient de Paddick moyen a diminué de manière significative, passant de 2,24 à 2,15 après correction du PSgD. Les troisième et quatrième études visaient à réduire les artefacts de mouvement de l'IRM à l'aide de nouveaux modèles génératifs. Les deux études ont été formées et testées auprès de 230 patients postopératoires. De plus, un ensemble de données de 148 patients présentant des artefacts de mouvement réels a été utilisé pour valider notre modèle afin de supprimer les artefacts de mouvement réels. Cette troisième étude a exploité la T₂ fluid-attenuated inversion recovery pour former un modèle génératif non supervisé permettant de supprimer les artefacts de mouvement ceT₁w. Nous avons comparé notre modèle avec deux modèles génératifs bien connus. Nous avons également réalisé une étude d'ablation pour évaluer différentes variantes de notre modèle proposé. Notre modèle a mieux performé quantitativement et qualitativement par rapport aux autres modèles. Dans la quatrième étude, nous avons utilisé un modèle probabiliste de diffusion conditionnelle stable pour supprimer les artefacts de mouvement ceT₁w. Cette étude a proposé une nouvelle méthode pour simuler le modèle de diffusion directe dans l'espace *k*, permettant la génération d'ensembles de données in silico avec différents niveaux d'artefacts de mouvement. Notre méthode a surpassé quantitativement et qualitativement deux modèles génératifs et un modèle supervisé. Le dernier volet a proposé une méthode non supervisée pour fusionner la tomodensitométrie et l'IRM, car elles fournissent des informations complémentaires pour la planification de la radiothérapie. Les images TDM sont géométriquement corrigées et fournissent des informations sur la densité électronique pour un calcul précis de la dose et l'IRM fournit des images avec un excellent contraste des tissus mous. Cette étude a proposé un modèle génératif de fusion médicale non supervisé pour fusionner la tomodensitométrie et le ceT₁w haute résolution afin de générer une image avec la structure osseuse de la tomodensitométrie et le contraste des tissus mous de l'IRM. Ainsi, la délimitation des tumeurs et des organes à risque sera améliorée et le temps de planification de la radiothérapie sera réduit. Les résultats de notre méthode ont été comparés quantitativement et qualitativement à sept méthodes traditionnelles et à huit méthodes de pointe DL, où notre méthode les a toutes surpassées. / The main goal of this thesis is to explore multi-parametric magnetic resonance imaging (mpMRI) applications in radiation oncology. MRI images are used in different stages of the treatments including planning, guiding, and following up. The planning step usually utilizes MRI soft-tissue contrast or physiological images to delineate targets. The guiding step uses MRI soft-tissue contrast to guide radiation treatments. Quantitative MRI (qMRI) and mpMRI are used to evaluate a response to the treatment. This thesis aimed to use qMRI maps to quantify head and neck tissues' response to oxygenation, reduce MRI artifacts and distortions, and finally fuse MRI and computed tomography (CT). This thesis was started to optimize qMRI maps including T₁ and R₂* to quantify oxygen preservation after delivering pure oxygen (pO₂) to the healthy volunteers. This study was the first in vivo quantifying head and neck tissue oxygen preservation after delivering pO₂ . Two separate experiments were carried out : (1) delivering normobaric oxygenation (NBO) and comparing serial qMRI T₁ and R₂* maps acquired during and after NBO with the baseline (before delivering NBO), and (2) delivering hyperbaric oxygenation (HBO) and comparing the serial qMRI T₁ and R₂* maps acquired after HBO with the baseline (before delivering HBO). We realized that the parotid gland and masseter muscle oxygen level increased significantly (*p* < 0.05) after delivering NBO for 15 minutes and returned to the baseline around seven minutes after the NBO. In the HBO experiment, a similar trend was observed. We did not find a difference between baseline and post-HBO qMRI images. Then, our thesis concentrated on the MRI distortions and artifacts. In the second study, we proposed a deep learning model to predict MRI patient-specific geometry distortions (PSgD) caused by susceptibility differences between two tissues. Then, the contrast-enhanced T₁W (ceT₁W) was corrected for the PSgD. We trained our deep learning (DL) model using a dataset of 380 subjects. We compared the Gamma Knife stereotactic radiosurgery (GK-SRS) dose distributions for the patients with vestibular schwannoma calculated using PSgD-corrected and PSgD-corrupted images. We found that on average D$_μ$, D$_\textup{95\%}$ et D$_\textup{min}$ in Gy significantly increased after PSgD correction from 16.85 to 17.25, 12.30 to 12.77, and from 8.98 to 9.92. Average target coverage and Paddick's conformity index significantly increased from 0.97 to 0.98, and 0.94 to 0.96. Average Paddick's gradient index decreased significantly from 2.24 to 2.15 after PSgD correction. The third and fourth studies aimed to reduce MRI motion artifacts using novel generative models. Both studies were trained and tested using post-surgery 230 patients. In addition, a dataset of 148 patients with real motion artifacts was used to validate our model to remove real motion artifacts. The third study leveraged T₂-fluid-attenuated inversion recovery to train an unsupervised generative model for removing ceT₁W motion artifacts. We compared our model with two well-known generative models. Also, we performed an ablation study to evaluate different variations of our proposed model. Our model outperformed quantitatively and qualitatively other models. In the fourth study, we used a stable conditional diffusion probabilistic model to remove ceT₁W motion artifacts. This study proposed a novel method for simulating the forward diffusion model in *k*-space, enabling the generation of in-silico datasets with varying levels of motion artifacts. Our method quantitatively and qualitatively outperformed two generative models and one supervised model. The fifth study proposed an unsupervised method to fuse CT and MRI because they provide complementary information for radiation therapy planning. CT images are geometrically corrected and provide electron density information for accurate dose calculation and MRI provides images with excellent soft-tissue contrast. This study proposed an end-to-end unsupervised medical fusion generative model to fuse CT and high-resolution ceT₁W to generate an image with CT bone structure and MRI soft-tissue contrast. Thus, tumor and organs at risk delineation will be improved and the radiotherapy planning time will be reduced. The results of our method were quantitatively and qualitatively compared with seven traditional and eight DL state-of-the-art methods where our method outperformed the traditional and DL approaches.

Cette thèse de doctorat porte sur le traitement personnalisé des tumeurs neuroendocrines (TNE) en utilisant la thérapie radiométabolique par ¹⁷⁷Lu-octréotate. Comme la plupart des TNE surexpriment des récepteurs de la somatostatine (RSST)s, les analogues de la somatostatine radiomarqués sont utilisés dans le diagnostic des TNE par imagerie moléculaire, ainsi que pour la radiothérapie interne en utilisant des radionucléides émetteurs bêta (PRRT : peptide receptor radionuclide therapy). Les protocoles actuels de la PRRT impliquent l’administration de 4 cycles d’une activité fixe de ¹⁷⁷Lu-octréotate (p. ex. 7,4 GBq), ce qui donne des doses de radiation absorbée très variables pour les organes critiques et les lésions tumorales. Les faibles taux de toxicité sérieuse rénale et hématopoïétique (< 1 %) après traitement au ¹⁷⁷Lu-octréotate suggèrent que la plupart des patients reçoivent un traitement sous-optimal. La mise en place d’une personnalisation de l’activité administrée en fonction de la dosimétrie individualisée par imagerie moléculaire quantitative pourrait donc permettre de maximiser la réponse chez chaque patient, comparativement au régime empirique, et ce, en augmentant la dose tumorale absorbée sans dépasser les seuils de toxicité aux organes critiques. Les objectifs de cette thèse étaient donc, (i) de simuler un régime de PRRT personnalisée (P-PRRT) à partir d’une cohorte de patients atteints de TNE traités par PRRT empirique, (ii) de rapporter les premiers résultats dosimétriques, d’efficacité et de toxicité de notre essai prospectif de PRRT personnalisée (NCT02754297) et finalement (iii) d’évaluer l’exactitude et la reproductibilité interobservateurs des protocoles simplifiés de dosimétrie à deux temps (Jour 1 et 3) ou à un temps unique (Jour 3) basés sur la tomographie quantitative par émission de photon (QSPECT : quantitative single-photon emission computed tomography) effectuée après chaque administration de thérapie. Un autre objectif de la dernière étude était la mise à jour de notre protocole personnalisé de prescription de l’activité injectée (AI). (i) Pour répondre au premier objectif, nous nous sommes basés sur une étude rétrospective afin de concevoir un protocole personnalisé de PRRT basé sur les calculs dosimétriques. Par la suite, une simulation de ce protocole a été effectuée rétrospectivement chez des patients atteints de TNE traités par PRRT depuis 2012 à l’Hôpital de l’Hôtel-Dieu de Québec. Selon notre régime personnalisé de PRRT, la dose rénale absorbée par AI était prédite en fonction de la surface corporelle (SC) et le taux de filtration glomérulaire (TFG) pour le premier cycle de traitement, et par la dosimétrie rénale du cycle précédent (ou la moyenne de la dosimétrie rénale des deux derniers cycles) pour les cycles suivants. En normalisant la dose rénale absorbée administrée pendant les cycles d’induction, notre approche personnalisée a démontré une augmentation moyenne significative de 1,48 fois la dose cumulative maximale absorbée par la tumeur par rapport au régime empirique, augmentant potentiellement le bénéfice thérapeutique sans dépasser le seuil de toxicité rénale. / Cette thèse de doctorat porte sur le traitement personnalisé des tumeurs neuroendocrines (TNE) en utilisant la thérapie radiométabolique par 177Lu-octréotate. Comme la plupart des TNE surexpriment des récepteurs de la somatostatine (RSST)s, les analogues de la somatostatine radiomarqués sont utilisés dans le diagnostic des TNE par imagerie moléculaire, ainsi que pour la radiothérapie interne en utilisant des radionucléides émetteurs bêta (PRRT : peptide receptor radionuclide therapy). Les protocoles actuels de la PRRT impliquent l’administration de 4 cycles d’une activité fixe de 177Lu-octréotate (p. ex. 7,4 GBq), ce qui donne des doses de radiation absorbée très variables pour les organes critiques et les lésions tumorales. Les faibles taux de toxicité sérieuse rénale et hématopoïétique (< 1 %) après traitement au 177Lu-octréotate suggèrent que la plupart des patients reçoivent un traitement sous-optimal. La mise en place d’une personnalisation de l’activité administrée en fonction de la dosimétrie individualisée par imagerie moléculaire quantitative pourrait donc permettre de maximiser la réponse chez chaque patient, comparativement au régime empirique, et ce, en augmentant la dose tumorale absorbée sans dépasser les seuils de toxicité aux organes critiques. Les objectifs de cette thèse étaient donc, (i) de simuler un régime de PRRT personnalisée (P-PRRT) à partir d’une cohorte de patients atteints de TNE traités par PRRT empirique, (ii) de rapporter les premiers résultats dosimétriques, d’efficacité et de toxicité de notre essai prospectif de PRRT personnalisée (NCT02754297) et finalement (iii) d’évaluer l’exactitude et la reproductibilité interobservateurs des protocoles simplifiés de dosimétrie à deux temps (Jour 1 et 3) ou à un temps unique (Jour 3) basés sur la tomographie quantitative par émission de photon (QSPECT : quantitative single-photon emission computed tomography) effectuée après chaque administration de thérapie. Un autre objectif de la dernière étude était la mise à jour de notre protocole personnalisé de prescription de l’activité injectée (AI). (i) Pour répondre au premier objectif, nous nous sommes basés sur une étude rétrospective afin de concevoir un protocole personnalisé de PRRT basé sur les calculs dosimétriques. Par la suite, une simulation de ce protocole a été effectuée rétrospectivement chez des patients atteints de TNE traités par PRRT depuis 2012 à l’Hôpital de l’Hôtel-Dieu de Québec. Selon notre régime personnalisé de PRRT, la dose rénale absorbée par AI était prédite en fonction de la surface corporelle (SC) et le taux de filtration glomérulaire (TFG) pour le premier cycle de traitement, et par la dosimétrie rénale du cycle précédent (ou la moyenne de la dosimétrie rénale des deux derniers cycles) pour les cycles suivants. En normalisant la dose rénale absorbée administrée pendant les cycles d’induction, notre approche personnalisée a démontré une augmentation moyenne significative de 1,48 fois la dose cumulative maximale absorbée par la tumeur par rapport au régime empirique, augmentant potentiellement le bénéfice thérapeutique sans dépasser le seuil de toxicité rénale. (ii) Sur la base de notre étude rétrospective, et en appliquant le protocole de PRRT personnalisée développé aux patients atteints de TNE métastatique traités depuis avril 2016 à l’Hôpital de l’Hôtel-Dieu de Québec, nous avons rapporté les premiers résultats dosimétriques sur l’efficacité et la toxicité de notre essai prospectif de PRRT personnalisée. Notre approche personnalisée a démontré une augmentation significative de la dose tumorale absorbée médiane de 1,26 fois par rapport au régime de PRRT empirique simulé sans augmenter le taux de toxicité sérieuse des organes critiques liée au traitement. Le traitement personnalisé a donné un taux de réponse partielle ou mineure de 59 % avec une stabilité dans 33,3 % des cas. En particulier, la PRRT personnalisée a induit une réponse partielle ou mineure chez 84,6 % des patients atteints de TNE d’origine pancréatique (pTNE). Finalement, le profil de sécurité est similaire à celui du régime empirique alors qu’aucun patient n’a présenté d'insuffisance rénale sévère. (iii) La dosimétrie, obtenue à l’aide de petits volumes d’intérêt pour l'échantillonnage de concentration d’activité pour le rein, la moelle osseuse et la tumeur dominante, est basée sur trois temps (Jour 0, 1 et 3) selon la pratique courante. Toutefois, en raison de la complexité de cette procédure, des méthodes de dosimétrie plus pratiques mais robustes sont nécessaires pour une adoption clinique étendue. Nous avons démontré qu’un protocole simplifié de dosimétrie basé sur des images QSPECT à deux temps fournit des estimations de dose reproductibles et plus précises que les techniques basées sur un seul temps. Nous avons aussi observé une excellente reproductibilité interobservateurs, en particulier pour la dose rénale, en utilisant les deux méthodes les plus précises basées sur trois ou deux temps avec des erreurs moyennes < 1 % et des écarts-types ≤ 5 %. Finalement, nous avons mis à jour notre protocole de P-PRRT : la dose rénale absorbée peut être normalisée en personnalisant l’AI en fonction du produit du TFG avec la SC pour le premier cycle et de la dosimétrie rénale précédente pour les cycles subséquents. En conclusion, nous avons conçu et validé un protocole de personnalisation de la PRRT au 177Lu-octréotate basé sur la dosimétrie, dans lequel la dose absorbée au rein est contrôlée et normalisée. La personnalisation de la PRRT permet d’augmenter de façon significative la dose absorbée à la tumeur, et ainsi la probabilité d’obtenir une meilleure réponse thérapeutique, tout en limitant le risque de toxicité aux organes critiques. Nous proposons également un protocole de dosimétrie simplifiée, basé sur l'imagerie QSPECT en deux temps, en utilisant des petits volumes d’intérêt. Cette méthode donne des résultats de dosimétrie précis, notamment pour les reins et la tumeur, et une reproductibilité interobservateur élevée. / (ii) Sur la base de notre étude rétrospective, et en appliquant le protocole de PRRT personnalisée développé aux patients atteints de TNE métastatique traités depuis avril 2016 à l’Hôpital de l’Hôtel-Dieu de Québec, nous avons rapporté les premiers résultats dosimétriques sur l’efficacité et la toxicité de notre essai prospectif de PRRT personnalisée. Notre approche personnalisée a démontré une augmentation significative de la dose tumorale absorbée médiane de 1,26 fois par rapport au régime de PRRT empirique simulé sans augmenter le taux de toxicité sérieuse des organes critiques liée au traitement. Le traitement personnalisé a donné un taux de réponse partielle ou mineure de 59 % avec une stabilité dans 33,3 % des cas. En particulier, la PRRT personnalisée a induit une réponse partielle ou mineure chez 84,6 % des patients atteints de TNE d’origine pancréatique (pTNE). Finalement, le profil de sécurité est similaire à celui du régime empirique alors qu’aucun patient n’a présenté d'insuffisance rénale sévère. (iii) La dosimétrie, obtenue à l’aide de petits volumes d’intérêt pour l'échantillonnage de concentration d’activité pour le rein, la moelle osseuse et la tumeur dominante, est basée sur trois temps (Jour 0, 1 et 3) selon la pratique courante. Toutefois, en raison de la complexité de cette procédure, des méthodes de dosimétrie plus pratiques mais robustes sont nécessaires pour une adoption clinique étendue. Nous avons démontré qu’un protocole simplifié de dosimétrie basé sur des images QSPECT à deux temps fournit des estimations de dose reproductibles et plus précises que les techniques basées sur un seul temps. Nous avons aussi observé une excellente reproductibilité interobservateurs, en particulier pour la dose rénale, en utilisant les deux méthodes les plus précises basées sur trois ou deux temps avec des erreurs moyennes < 1 % et des écarts-types ≤ 5 %. Finalement, nous avons mis à jour notre protocole de P-PRRT : la dose rénale absorbée peut être normalisée en personnalisant l’AI en fonction du produit du TFG avec la SC pour le premier cycle et de la dosimétrie rénale précédente pour les cycles subséquents. En conclusion, nous avons conçu et validé un protocole de personnalisation de la PRRT au ¹⁷⁷Lu-octréotate basé sur la dosimétrie, dans lequel la dose absorbée au rein est contrôlée et normalisée. La personnalisation de la PRRT permet d’augmenter de façon significative la dose absorbée à la tumeur, et ainsi la probabilité d’obtenir une meilleure réponse thérapeutique, tout en limitant le risque de toxicité aux organes critiques. Nous proposons également un protocole de dosimétrie simplifiée, basé sur l'imagerie QSPECT en deux temps, en utilisant des petits volumes d’intérêt. Cette méthode donne des résultats de dosimétrie précis, notamment pour les reins et la tumeur, et une reproductibilité interobservateur élevée. / (ii) Based on our retrospective study and applying our developed personalized PRRT protocol to the patients with metastatic NET treated since April 2016 at l’Hotel-Dieu de Québec hospital, we reported the preliminary data on dosimetry, efficacy and safety from our prospective PRRT trial. Our personalized approach yielded a significant increase in the median absorbed tumor dose (1.26 fold) compared with the simulated empirical PRRT regimen, without increasing the rate of serious treatment-related toxicity for critical organs. The personalized treatment gave a partial or minor response rate of 59 % and stability in 33.3 % of cases. In particular, the personalized PRRT regime resulted in a partial or minor response in 84.6 % of patients with NET of pancreatic origin (pNET). Finally, the safety profile was similar to that of the empirical regimen whereas, no patient has had severe renal impairment. (iii) According to the current practice, the dosimetry, obtained using small volumes of interest activity for concentration sampling in the kidney, bone marrow and the dominant tumor, is based on three time points (day 0, 1 and 3). However, due to the complexity of this procedure, more practical yet robust dosimetry methods are required for extended clinical use. We have observed that a simplified dosimetry protocol based on two-time point QSPECT scanning on Day 1 and 3 after PRRT treatment provides reproducible and more accurate dose estimates than the techniques relying on a single time point for noninitial or all cycles. We also observed an excellent inter-observer reproducibility, in particular for the renal absorbed dose using the two most precise dosimetry methods based on three or two time points with mean relative errors < 1 % and standard deviations ≤ 5 %. In addition, we have updated our P-PRRT protocol. In particular, renal absorbed radiation dose over a four-cycle induction PRRT course can be standardized by personalizing IA based on the product of estimated GFR with BSA for the first cycle, and on prior renal dosimetry for the subsequent cycles. In conclusion, we have designed and validated a personalized PRRT protocol with ¹⁷⁷Lu-octreotate based on dosimetry, in which the absorbed dose to the kidney is controlled and standardized. The personalization of PRRT can significantly increase the dose absorbed to the tumor, and thus the probability of obtaining a better therapeutic response, while limiting the risk of toxicity to critical organs. We also proposed a simplified dosimetry protocol, based on two time-points QSPECT scans, using small volumes of interest. This method gives accurate dosimetry results, especially for the kidneys and the tumors, and high interobserver reproducibility. / This doctoral thesis is focused on the personalized treatment of neuroendocrine tumours (NET) using ¹⁷⁷Lu-octreotate radiometabolic therapy. Since NETs overexpress somatostatin receptors (SSRT)s, radiolabeled somatostatin analogues are used in the diagnosis of NETs by molecular imaging, and for internal radiotherapy using beta-emitting radionuclides (PRRT: peptide receptor radionuclide therapy). The existing protocols for PRRT involve the administration of 4 cycles of a fixed activity of ¹⁷⁷Lu-octreotate (i.e., 7.4 GBq), which gives highly variable absorbed radiation doses to critical organs and tumor lesions. The low incidence of serious renal and hematopoietic toxicities (< 1 %) of ¹⁷⁷Luoctreotate suggests that the majority of patients are undertreated. The personalization of the administered activity according to the individualized dosimetry obtained from quantitative molecular imaging could make possible to maximize the response in each patient, compared to the empirical regime, without exceeding the thresholds of toxicity to the critical organs. The objectives of this thesis were, therefore (i) to simulate a personalized PRRT (P-PRRT) regime in a cohort of patients with NET treated with empiric PRRT, (ii) to report the preliminary data on dosimetry, efficacy and safety from our prospective PRRT trial and (iii) to assess the accuracy and the inter-observer reproducibility of simplified two time points (Day 1 and 3) or one time point (Day 3) dosimetry protocols based on quantitative single-photon emission computed tomography (QSPECT) performed after each PRRT cycle. Another objective of the last study concerned the update of our P-PRRT protocol. (i) We first reviewed retrospective data on patients with NETs treated with PRRT since 2012 at l’Hôtel-Dieu de Québec hospital. According to our P-PRRT regimen, the renal absorbed dose per injected activity (IA) was predicted by the body surface (BSA) and the glomerular filtration rate (GFR) for the first induction cycle, and by the renal dosimetry of the previous cycle (or the average renal dosimetry of the last two cycles) for the following cycles. By normalizing the absorbed renal dose administered during the induction cycles, our personalized approach resulted in an average 1.48-fold increase in the maximum cumulative tumour absorbed dose compared with the empiric regimen, potentially increasing the therapeutic benefit without exceeding the threshold of renal toxicity.

Thèse ou mémoire avec insertion d'articles / Cette thèse porte sur le développement de modèles mathématiques qui permettent de prédire des métriques quantifiant la dose pour la planification de radiothérapie externe. La planification de traitement en radiothérapie externe est un processus complexe et itératif. Le planificateur spécifie des objectifs à atteindre et l'optimisation du plan de traitement est faite par ordinateur. Chaque patient ayant une anatomie propre, il peut être difficile de savoir si le plan de traitement est optimal c.à.d. où les objectifs de couverture en dose du volume cible sont remplis tout en minimisant la dose aux organes sains pour diminuer la possibilité d'effets secondaires néfastes. Afin de pallier cette problématique, plusieurs types de contrôles de qualité (CQ) en planification de traitement ont été développés. Plus précisément, le projet proposé utilise l'analyse de frontières stochastiques, une méthode mathématique conçue pour l'économétrie, afin de déterminer la dose minimale atteignable aux organes à risque. La planification de traitement en radiothérapie externe, les techniques utilisées ainsi que les principaux défis sont d'abord présentés dans le cadre de cette thèse. Une revue des principaux types de contrôles de qualité des plans de traitement est par la suite présentée afin de mettre en contexte le projet doctoral. Les principaux principes mathématiques employés tels que l'analyse de frontières stochastiques et la gestion de données manquantes sont par la suite abordés. Finalement, l'analyse de frontières stochastiques appliquée à la prédiction de paramètres dosimétriques en planification de traitement est présentée. Son approche fréquentiste est utilisée pour la prédiction dosimétrique pour le rectum et la vessie pour le cancer de la prostate traité par VMAT. Sept paramètres géométriques ont été extraits afin de caractériser la relation entre les organes à risque et le volume cible. Au total, 37 paramètres dosimétriques ont été prédits pour les deux organes à risque. Le modèle développé a été testé sur une cohorte de 30 patients avec une dose de prescription de 60-70 Gy où 77% (23 sur 30) des DVH prédits pour le rectum et la vessie présentent une déviation de 5% ou moins avec le DVH planifié. De plus, l'analyse de frontières stochastiques bayésienne, en plus d'un modèle de gestion de données manquantes, sont utilisés pour la prédiction dosimétrique de six organes à risque pour le cancer du poumon traité par SBRT. Au total, 16 indices dosimétriques ont été prédits pour l'arbre bronchique principal, le cœur, l'œsophage, le PRV de la moelle, les vaisseaux ainsi que la paroi thoracique. Le modèle prédictif est testé sur une cohorte de 50 patients. La différence moyenne entre la dose prédite et planifiée pour le se situe à 1.5 ± 1.9 Gy et 4.9 ± 5.3 Gy pour le D0.35cc du PRV de la moelle et le D0.035cc de l'arbre bronchique principal respectivement. L'analyse de frontières stochastiques est ainsi démontrée comme une solide base pour un CQ en planification de traitement en radiothérapie externe. / This thesis is on the development of mathematical models predicting dosimetric metrics for treatment planning in external beam radiation therapy (EBRT). EBRT treatment planning is complex and iterative. Planner sets plan objectives to attain during the computerized plan optimization process. Because of specific patient anatomical variation, it is difficult to determine if a plan is optimal regarding the target coverage and the sparing of the organs-at-risk. Several quality control such as knowledge-based planning and multicriteria optimization were developed in order to address this problematic. This thesis focus on the use of an econometric method, Stochastic Frontier Analysis, to predict minimum achievable dose to organs-at-risk for several treatment sites. EBRT treatment planning, techniques and challenges are first discussed. A review of the most common quality controls is made in order to assess the importance of the project presented in this thesis. Mathematical considerations such as the theoretical foundation of stochastic frontier analysis and missing data management are also discussed. Stochastic frontier analysis is used to develop predictive models of dose distribution in treatment planning. Frequentist stochastic frontier analysis is used to predict dosimetric metrics for the rectum and the bladder for prostate cancer treated by VMAT. Seven geometric parameters are extracted to characterize the relationship between the organs-at-risk with the planning volume. In total, 37 dosimetric parameters are tested. The developped model is tested using validation cohort (30 patients with prescribed dose between 60 and 70 Gy) where 77% (23 out of 30) of the predicted DVHs present a 5% or less dose deterioration with the planned DVH. Bayesian stochastic frontier analysis with a missing data management is then presented for the prediction of dosimetric parameters for 6 organs-at-risk for lung cancer treated by SBRT. 16 DVH metrics were predicted for the main bronchus, heart, esophagus, spinal cord PRV, great vessels and chest wall. The predictive model is tested on a test group of 50 patients. The mean difference between the observed and predicted values ranges between 1.5 ± 1.9 Gy and 4.9 ± 5.3 Gy for the spinal cord PRV D0.35cc and the main bronchus D0.035cc respectively. Stochastic frontier analysis is considered to be a new and valid method used in order to develop a quality control for EBRT treatment planning.