Radiothérapie guidée par l'image du cancer de la prostate (vers l'intégration des déformations anatomiques)

Ce travail de thèse porte sur la quantification et la prise en compte des variations anatomiques en cours de radiothérapie guidée par l'image pour le cancer de la prostate. Nous proposons tout d'abord une approche basée population pour quantifier et analyser les incertitudes géométriques, notamment à travers des matrices de probabilité de présence de la cible en cours de traitement. Nous proposons ensuite une méthode d'optimisation des marges suivant des critères de couverture géométrique de la cible tumorale. Cette méthode permet d'obtenir des marges objectives associées aux différents types d'incertitudes géométriques et aux différentes modalités de repositionnement du patient. Dans un second temps, nous proposons une méthode d'estimation de la dose cumulée reçue localement par les tissus pendant un traitement de radiothérapie de la prostate. Cette méthode repose notamment sur une étape de recalage d'images de façon à estimer les déformations des organes entre les séances de traitement et la planification. Différentes méthodes de recalage sont proposées, suivant les informations disponibles (délinéations ou points homologues) pour contraindre la déformation estimée. De façon à évaluer les méthodes proposées au regard de l'objectif de cumul de dose, nous proposons ensuite la génération et l'utilisation d'un fantôme numérique reposant sur un modèle biomécanique des organes considérés. Les résultats de l'approche sont présentés sur ce fantôme numérique et sur données réelles. Nous montrons ainsi que l'apport de contraintes géométriques permet d'améliorer significativement la précision du cumul et que la méthode reposant sur la sélection de contraintes ponctuelles présente un bon compromis entre niveau d'interaction et précision du résultat. Enfin, nous abordons la question de l'analyse de données de populations de patients dans le but de mieux comprendre les relations entre dose délivrée localement et effets cliniques. Grâce au recalage déformable d'une population de patients sur une référence anatomique, les régions dont la dose est significativement liée aux événements de récidive sont identifiées. Il s'agit d'une étude exploratoire visant à terme à mieux exploiter l'information portée par l'intégralité de la distribution de dose, et ce en fonction du profil du cancer.This work deals with the quantification and the compensation of anatomical deformations during image-guided radiotherapy of prostate cancer. Firstly, we propose a population-based approach to quantify the geometrical uncertainties by means of coverage probability matrices of the target tumor during the treatment. We then propose a margins optimization method based on geometrical coverage criteria of the tumor target. This method provides rationnal margins models associated to the different geometrical uncertainties and patient repositioning protocols. Secondly, we propose a method to estimate the locally accumulated dose during the treatment. This method relies on a deformable image registration process in order to estimate the organ deformations between each treatment fraction and the planning. Different registration methods are proposed, using different level of user interactions (landmarks specification or delineations) to constrain the deformation estimation. In order to evaluate the performance of the proposed methods, we then describe the generation of a numerical phantom based on a biomechanical model. The results are presented for the numerical phantom and real clinical cases. We show that the benefit brought by the manual placement of some landmarks to constrain the registration represents a good compromise between the required interaction level and the dose estimation accuracy. Finally, we address the issue of the analysis of population data in order to better understand the relationship between the locally delivered dose and clinical effects. With deformable image registration of a population of patients on an anatomical template, regions whose dose is significantly associated with recurrence events are identified. This last part is an exploratory study aiming to better use the information carried by the entire distribution dose, and according to the cancer profile.RENNES1-Bibl. électronique (352382106) / SudocSudocFranceF

Analysis of contrast-enhanced medical images.

Early detection of human organ diseases is of great importance for the accurate diagnosis and institution of appropriate therapies. This can potentially prevent progression to end-stage disease by detecting precursors that evaluate organ functionality. In addition, it also assists the clinicians for therapy evaluation, tracking diseases progression, and surgery operations. Advances in functional and contrast-enhanced (CE) medical images enabled accurate noninvasive evaluation of organ functionality due to their ability to provide superior anatomical and functional information about the tissue-of-interest. The main objective of this dissertation is to develop a computer-aided diagnostic (CAD) system for analyzing complex data from CE magnetic resonance imaging (MRI). The developed CAD system has been tested in three case studies: (i) early detection of acute renal transplant rejection, (ii) evaluation of myocardial perfusion in patients with ischemic heart disease after heart attack; and (iii), early detection of prostate cancer. However, developing a noninvasive CAD system for the analysis of CE medical images is subject to multiple challenges, including, but are not limited to, image noise and inhomogeneity, nonlinear signal intensity changes of the images over the time course of data acquisition, appearances and shape changes (deformations) of the organ-of-interest during data acquisition, determination of the best features (indexes) that describe the perfusion of a contrast agent (CA) into the tissue. To address these challenges, this dissertation focuses on building new mathematical models and learning techniques that facilitate accurate analysis of CAs perfusion in living organs and include: (i) accurate mathematical models for the segmentation of the object-of-interest, which integrate object shape and appearance features in terms of pixel/voxel-wise image intensities and their spatial interactions; (ii) motion correction techniques that combine both global and local models, which exploit geometric features, rather than image intensities to avoid problems associated with nonlinear intensity variations of the CE images; (iii) fusion of multiple features using the genetic algorithm. The proposed techniques have been integrated into CAD systems that have been tested in, but not limited to, three clinical studies. First, a noninvasive CAD system is proposed for the early and accurate diagnosis of acute renal transplant rejection using dynamic contrast-enhanced MRI (DCE-MRI). Acute rejection–the immunological response of the human immune system to a foreign kidney–is the most sever cause of renal dysfunction among other diagnostic possibilities, including acute tubular necrosis and immune drug toxicity. In the U.S., approximately 17,736 renal transplants are performed annually, and given the limited number of donors, transplanted kidney salvage is an important medical concern. Thus far, biopsy remains the gold standard for the assessment of renal transplant dysfunction, but only as the last resort because of its invasive nature, high cost, and potential morbidity rates. The diagnostic results of the proposed CAD system, based on the analysis of 50 independent in-vivo cases were 96% with a 95% confidence interval. These results clearly demonstrate the promise of the proposed image-based diagnostic CAD system as a supplement to the current technologies, such as nuclear imaging and ultrasonography, to determine the type of kidney dysfunction. Second, a comprehensive CAD system is developed for the characterization of myocardial perfusion and clinical status in heart failure and novel myoregeneration therapy using cardiac first-pass MRI (FP-MRI). Heart failure is considered the most important cause of morbidity and mortality in cardiovascular disease, which affects approximately 6 million U.S. patients annually. Ischemic heart disease is considered the most common underlying cause of heart failure. Therefore, the detection of the heart failure in its earliest forms is essential to prevent its relentless progression to premature death. While current medical studies focus on detecting pathological tissue and assessing contractile function of the diseased heart, this dissertation address the key issue of the effects of the myoregeneration therapy on the associated blood nutrient supply. Quantitative and qualitative assessment in a cohort of 24 perfusion data sets demonstrated the ability of the proposed framework to reveal regional perfusion improvements with therapy, and transmural perfusion differences across the myocardial wall; thus, it can aid in follow-up on treatment for patients undergoing the myoregeneration therapy. Finally, an image-based CAD system for early detection of prostate cancer using DCE-MRI is introduced. Prostate cancer is the most frequently diagnosed malignancy among men and remains the second leading cause of cancer-related death in the USA with more than 238,000 new cases and a mortality rate of about 30,000 in 2013. Therefore, early diagnosis of prostate cancer can improve the effectiveness of treatment and increase the patient’s chance of survival. Currently, needle biopsy is the gold standard for the diagnosis of prostate cancer. However, it is an invasive procedure with high costs and potential morbidity rates. Additionally, it has a higher possibility of producing false positive diagnosis due to relatively small needle biopsy samples. Application of the proposed CAD yield promising results in a cohort of 30 patients that would, in the near future, represent a supplement of the current technologies to determine prostate cancer type. The developed techniques have been compared to the state-of-the-art methods and demonstrated higher accuracy as shown in this dissertation. The proposed models (higher-order spatial interaction models, shape models, motion correction models, and perfusion analysis models) can be used in many of today’s CAD applications for early detection of a variety of diseases and medical conditions, and are expected to notably amplify the accuracy of CAD decisions based on the automated analysis of CE images

Quantifying, Understanding and Predicting Differences Between Planned and Delivered Dose to Organs at Risk in Head & Neck Cancer Patients Undergoing Radical Radiotherapy to Promote Intelligently Targeted Adaptive Radiotherapy

Introduction: Radical radiotherapy (RT) is an effective but toxic treatment for head and neck cancer (HNC). Contemporary radiotherapy techniques sculpt dose to target disease and avoid organs at risk (OARs), but anatomical change during treatment mean that the radiation dose delivered to the patient – delivered dose (DA), is different to that anticipated at planning – planned dose (DP). Modifying the RT plan during treatment – Adaptive Radiotherapy (ART) – could mitigate these risks by reducing dose to OARs. However, clinical data to guide patient selection for, and timing of ART, are for lacking. Methods: 337 patients with HNC were recruited to the Cancer Research UK VoxTox study. Demographic, disease and treatment data were collated, and both DP and DA to organs at risk (OARs) were computed from daily megavoltage CT image guidance scans, using an open-source deformable image registration package (Elastix). Toxicity data were prospectively collected. Relationships between DP, DA and late toxicities were investigated with univariate, and logistic regression normal tissue complication probability (NTCP) modelling approaches. A sub-study of VoxTox recruited 18 patients who had MRI scans before RT fractions 1, 6, 16, and 26. Changes in salivary gland volumes and relative apparent diffusion coefficient (ADC) values were measured and related to toxicity events. Results: Spinal cord dose differences were small, and not predicted by weight loss or shape change. Mean DA to all other OARs was higher than DP; factors predicting higher DA included primary disease site, concomitant therapy, shape change and advanced neck disease. Nine patients (3.7%) saw DA>DP by 2Gy to more than half of the OARs assessed. These patients all had received bilateral neck RT for N-stage 2b oropharyngeal cancer. Strong uni- and multivariate relationships between OAR dose and toxicity were seen. Differences between DA and DP-based dose-toxicity models were minimal, and not statistically significant. On MRI, both parotid and submandibular glands shrank during treatment, whilst relative ADC rose. Relationships with toxicity were inconclusive. Conclusions: Small differences between OAR DP and DA mean that DA-based toxicity prediction models confer negligible additional benefit at the population level. Factors such as primary disease sub-site, concomitant systemic therapy, staging and shape change may help to select the patients that do develop clinically significant dose differences, and would benefit most from ART for toxicity reduction

Computational methods for the analysis of functional 4D-CT chest images.

Medical imaging is an important emerging technology that has been intensively used in the last few decades for disease diagnosis and monitoring as well as for the assessment of treatment effectiveness. Medical images provide a very large amount of valuable information that is too huge to be exploited by radiologists and physicians. Therefore, the design of computer-aided diagnostic (CAD) system, which can be used as an assistive tool for the medical community, is of a great importance. This dissertation deals with the development of a complete CAD system for lung cancer patients, which remains the leading cause of cancer-related death in the USA. In 2014, there were approximately 224,210 new cases of lung cancer and 159,260 related deaths. The process begins with the detection of lung cancer which is detected through the diagnosis of lung nodules (a manifestation of lung cancer). These nodules are approximately spherical regions of primarily high density tissue that are visible in computed tomography (CT) images of the lung. The treatment of these lung cancer nodules is complex, nearly 70% of lung cancer patients require radiation therapy as part of their treatment. Radiation-induced lung injury is a limiting toxicity that may decrease cure rates and increase morbidity and mortality treatment. By finding ways to accurately detect, at early stage, and hence prevent lung injury, it will have significant positive consequences for lung cancer patients. The ultimate goal of this dissertation is to develop a clinically usable CAD system that can improve the sensitivity and specificity of early detection of radiation-induced lung injury based on the hypotheses that radiated lung tissues may get affected and suffer decrease of their functionality as a side effect of radiation therapy treatment. These hypotheses have been validated by demonstrating that automatic segmentation of the lung regions and registration of consecutive respiratory phases to estimate their elasticity, ventilation, and texture features to provide discriminatory descriptors that can be used for early detection of radiation-induced lung injury. The proposed methodologies will lead to novel indexes for distinguishing normal/healthy and injured lung tissues in clinical decision-making. To achieve this goal, a CAD system for accurate detection of radiation-induced lung injury that requires three basic components has been developed. These components are the lung fields segmentation, lung registration, and features extraction and tissue classification. This dissertation starts with an exploration of the available medical imaging modalities to present the importance of medical imaging in today’s clinical applications. Secondly, the methodologies, challenges, and limitations of recent CAD systems for lung cancer detection are covered. This is followed by introducing an accurate segmentation methodology of the lung parenchyma with the focus of pathological lungs to extract the volume of interest (VOI) to be analyzed for potential existence of lung injuries stemmed from the radiation therapy. After the segmentation of the VOI, a lung registration framework is introduced to perform a crucial and important step that ensures the co-alignment of the intra-patient scans. This step eliminates the effects of orientation differences, motion, breathing, heart beats, and differences in scanning parameters to be able to accurately extract the functionality features for the lung fields. The developed registration framework also helps in the evaluation and gated control of the radiotherapy through the motion estimation analysis before and after the therapy dose. Finally, the radiation-induced lung injury is introduced, which combines the previous two medical image processing and analysis steps with the features estimation and classification step. This framework estimates and combines both texture and functional features. The texture features are modeled using the novel 7th-order Markov Gibbs random field (MGRF) model that has the ability to accurately models the texture of healthy and injured lung tissues through simultaneously accounting for both vertical and horizontal relative dependencies between voxel-wise signals. While the functionality features calculations are based on the calculated deformation fields, obtained from the 4D-CT lung registration, that maps lung voxels between successive CT scans in the respiratory cycle. These functionality features describe the ventilation, the air flow rate, of the lung tissues using the Jacobian of the deformation field and the tissues’ elasticity using the strain components calculated from the gradient of the deformation field. Finally, these features are combined in the classification model to detect the injured parts of the lung at an early stage and enables an earlier intervention

A Method for Predicting Dose Changes for HN Treatment Using Surface Imaging

Head and neck cancer is commonly treated with a six- to seven-week course of radiotherapy, during which a patient’s anatomy may change substantially, due to target reduction or weight loss. Anatomical changes lead to reduction in treatment quality due to decreased setup reproducibility and altered dose deposition compared to the original plan. Few clinics have developed a standard method for triggering resimulation and replan due to anatomic changes. This work investigates a new method for determining when to resimulate and replan HNC patients by utilizing their topographic anatomical changes to predict differences in planned versus delivered dose distributions. The first part of the work presents a method for deformable image registration of CT to CBCT which addresses the challenges of inaccurate Hounsfield units and truncated field of view present in CBCT. The registration method was validated on 10 HN patients using contour comparison, with average DSC of 0.82, 0.74, 0.72, and 0.69 for mandible, cord, and left and right parotid. The registration method was then used to generate dose maps and surface contours for 47 patients for the second part of this work, the development of a U-Net which takes the original dose distribution, the original surface, and the treatment day surface as input and predicts the treatment day dose distribution as output. The average RMSE and MAE between the true and predicted dose distributions for a test set of 6 patients was 4.25 and 2.15. This work proves feasibility of a dose prediction neural network using surface imaging

Adaptive Radiotherapy Enabled by MRI Guidance.

Adaptive radiotherapy (ART) strategies systematically monitor variations in target and neighbouring structures to inform treatment-plan modification during radiotherapy. This is necessary because a single plan designed before treatment is insufficient to capture the actual dose delivered to the target and adjacent critical structures during the course of radiotherapy. Magnetic resonance imaging (MRI) provides superior soft-tissue image contrast over current standard X-ray-based technologies without additional radiation exposure. With integrated MRI and radiotherapy platforms permitting motion monitoring during treatment delivery, it is possible that adaption can be informed by real-time anatomical imaging. This allows greater treatment accuracy in terms of dose delivered to target with smaller, individualised treatment margins. The use of functional MRI sequences would permit ART to be informed by imaging biomarkers, so allowing both personalised geometric and biological adaption. In this review, we discuss ART solutions enabled by MRI guidance and its potential gains for our patients across tumour types

Modelling and verification of doses delivered to deformable moving targets in radiotherapy

During the last two decades, advanced treatment techniques have been developed in radiotherapy to achieve more conformal beam targeting of cancerous lesions. The advent of these techniques, such as intensity modulated radiotherapy (IMRT), volumetric modulated arc radiothreapy (VMAT), Tomotherapy etc., allows more precise localisation of higher doses to complex-shaped target volumes, thereby sparing more healthy tissue. In this context, motion management is a critical issue in contemporary radiotherapy (RT). That anatomic structures move during respiration is well known and much research is presently being devoted to strategies to contend with organ motion. However, moving structures are typically regarded as rigid bodies. The fact that many structures deform as a result of motion makes their resultant dose distributions difficult to measure and calculate, and has not been fully accounted for. The potential for ineffective treatments that do not take into account motion and anatomic deformation is self-evident. This thesis addresses the pressing need to investigate dose distributions in targets that deform during and/or between treatments, to ensure robust calculations for dose accumulation and delivery, thus providing the most positive outcomes for patients. This involves the direct measurement of complex and re-distributed dose in deforming objects (an experimental model), as well as calculations of the deformed dose distribution (a mathematical model). The comparison thereof aims to validate the dose deformation technique, thereby to apply the method to a clinical example such as liver stereotactic body radiotherapy. To facilitate four-dimensional deformable dosimetry for both external beam radiotherapy and brachytherapy, methodologies for three-dimensional deformed dose measurements were developed and employed using radiosensitive polymer gel combined with a cone beam optical computed tomography (CT) scanner. This includes the development of a novel prototype deformable target volume using a tissue-equivalent, deformable gel dosimetric phantom, dubbed “defgel”. This can reproducibly simulate targets subject to a range of mass- and density-conserving deformations representative of those observable in anatomical targets. This novel tool was characterised in terms of its suitability for the measurement of dose in deforming geometries. It was demonstrated that planned doses could be delivered to the deformable gel dosimeter in the presence of different deformations and complex spatial re-distributions of dose in all three dimensions could be quantified. For estimating the cumulative dose in different deformed states, deformable image registration (DIR) algorithms were implemented to ‘morph’ a dose distribution calculated by a treatment planning system. To investigate the performance of DIR and dose-warping technique, two key studies were undertaken. The first was to systematically assess the accuracy of a range of different DIR algorithms available in the public domain and quantitatively examine, in particular, low-contrast regions, where accuracy had not previously been established. This work investigates DIR algorithms in 3D via a systematic evaluation process using defgel suitable for verification of mass- and density-conserving deformations. The second study was a full three-dimensional experimental validation of the dose-warping technique using the evaluated DIR algorithm and comparing it to directly measured deformed dose distributions from defgel. It was shown that the dose-warping can be accurate, i.e. over 95% passing rate of 3D-gamma analysis with 3%/3mm criteria for given extents of deformation up to 20 mm For the application of evaluating patient treatment planning involving tumour motion/deformation, two key studies were undertaken in the context of liver stereotactic body radiotherapy. The first was a 4D evaluation of conventional 3D treatment planning, combined with 4D computed tomography, in order to investigate the extent of dosimetric differences between conventional 3D-static and path-integrated 4D-cumulative dose calculation. This study showed that the 3D planning approach overestimated doses to targets by ≤ 9% and underestimated dose to normal liver by ≤ 8%, compared to the 4D methodology. The second study was to assess a consequent reduction of healthy tissue sparing, which may increase risk for surrounding healthy tissues. Estimates for normal tissue complications probabilities (NTCP) based on the two dose calculation schemes are provided. While all NTCP were low for the employed fractionation scheme, analysis of common alternative schemes suggests potentially larger uncertainties exist in the estimation of NTCP for healthy liver and that substantial differences in these values may exist across the different fractionation schemes. These bodies of work have shown the potential to quantify such issues of under- and/or over-dosages which are quite patient dependent in RT. Studies presented in this work consolidate gel dosimetry, image guidance, DIR, dose-warping and consequent dose accumulation calculation to investigate the dosimetric impact and make more accurate evaluation of conventional 3D treatment plans. While liver stereotactic body radiotherapy (SBRT) was primarily concerned for immediate clinical application, the findings of this thesis are also applicable to other organs with various RT techniques. Most importantly, however, it is hoped that the outcomes of this thesis will help to improve treatment plan accuracy. By considering both computation and measurement, it is also hoped that this work will open new windows for future work and hence provide building blocks to further enhance the benefit of radiotherapy treatment

Optimization of Decision Making in Personalized Radiation Therapy using Deformable Image Registration

Cancer has become one of the dominant diseases worldwide, especially in western countries, and radiation therapy is one of the primary treatment options for 50% of all patients diagnosed. Radiation therapy involves the radiation delivery and planning based on radiobiological models derived primarily from clinical trials. Since 2015 improvements in information technologies and data storage allowed new models to be created using the large volumes of treatment data already available and correlate the actually delivered treatment with outcomes. The goals of this thesis are to 1) construct models of patient outcomes after receiving radiation therapy using available treatment and patient parameters and 2) provide a method to determine real accumulated radiation dose including the impact of registration uncertainties. In Chapter 2, a model was developed predicting overall survival for patients with hepatocellular carcinoma or liver metastasis receiving radiation therapy. These models show which patients benefit from curative radiation therapy based on liver function, and the survival benefit of increased radiation dose on survival. In Chapter 3, a method was developed to routinely evaluate deformable image registration (DIR) with computer-generated landmark pairs using the scale-invariant feature transform. The method presented in this chapter created landmark sets for comparing lung 4DCT images and provided the same evaluation of DIR as manual landmark sets. In Chapter 4, an investigation was performed on the impact of DIR error on dose accumulation using landmarked 4DCT images as the ground truth. The study demonstrated the relationship between dose gradient, DIR error and dose accumulation error, and presented a method to determine error bars on the dose accumulation process. In Chapter 5, a method was presented to determine quantitatively when to update a treatment plan during the course of a multi-fraction radiation treatment of head and neck cancer. This method investigated the ability to use only the planned dose with deformable image registration to predict dose changes caused by anatomical deformations. This thesis presents the fundamental elements of a decision support system including patient pre-treatment parameters and the actual delivered dose using DIR while considering registration uncertainties

Consistent and invertible deformation vector fields for a breathing anthropomorphic phantom: a post-processing framework for the XCAT phantom.

Breathing motion is challenging for radiotherapy planning and delivery. This requires advanced four-dimensional (4D) imaging and motion mitigation strategies and associated validation tools with known deformations. Numerical phantoms such as the XCAT provide reproducible and realistic data for simulation-based validation. However, the XCAT generates partially inconsistent and non-invertible deformations where tumours remain rigid and structures can move through each other. We address these limitations by post-processing the XCAT deformation vector fields (DVF) to generate a breathing phantom with realistic motion and quantifiable deformation. An open-source post-processing framework was developed that corrects and inverts the XCAT-DVFs while preserving sliding motion between organs. Those post-processed DVFs are used to warp the first XCAT-generated image to consecutive time points providing a 4D phantom with a tumour that moves consistently with the anatomy, the ability to scale lung density as well as consistent and invertible DVFs. For a regularly breathing case, the inverse consistency of the DVFs was verified and the tumour motion was compared to the original XCAT. The generated phantom and DVFs were used to validate a motion-including dose reconstruction (MIDR) method using isocenter shifts to emulate rigid motion. Differences between the reconstructed doses with and without lung density scaling were evaluated. The post-processing framework produced DVFs with a maximum [Formula: see text]-percentile inverse-consistency error of 0.02 mm. The generated phantom preserved the dominant sliding motion between the chest wall and inner organs. The tumour of the original XCAT phantom preserved its trajectory while deforming consistently with the underlying tissue. The MIDR was compared to the ground truth dose reconstruction illustrating its limitations. MIDR with and without lung density scaling resulted in small dose differences up to 1 Gy (prescription 54 Gy). The proposed open-source post-processing framework overcomes important limitations of the original XCAT phantom and makes it applicable to a wider range of validation applications within radiotherapy