Background: Dynamic contrast-enhanced MRI (DCE-MRI) has the potential to produce images of physiological quantities such as blood flow, blood vessel volume fraction, and blood vessel permeability. Such information is highly valuable, e.g., in oncology. The focus of this work was to improve the quantitative aspects of DCE-MRI in terms of better understanding of error sources and their effect on estimated physiological quantities. Methods: Firstly, a novel parameter estimation algorithm was developed to overcome a problem with sensitivity to the initial guess in parameter estimation with a specific pharmacokinetic model. Secondly, the accuracy of the arterial input function (AIF), i.e., the estimated arterial blood contrast agent concentration, was evaluated in a phantom environment for a standard magnitude-based AIF method commonly used in vivo. The accuracy was also evaluated in vivo for a phase-based method that has previously shown very promising results in phantoms and in animal studies. Finally, a method was developed for estimation of uncertainties in the estimated physiological quantities. Results: The new parameter estimation algorithm enabled significantly faster parameter estimation, thus making it more feasible to obtain blood flow and permeability maps from a DCE-MRI study. The evaluation of the AIF measurements revealed that inflow effects and non-ideal radiofrequency spoiling seriously degrade magnitude-based AIFs and that proper slice placement and improved signal models can reduce this effect. It was also shown that phase-based AIFs can be a feasible alternative provided that the observed difficulties in quantifying low concentrations can be resolved. The uncertainty estimation method was able to accurately quantify how a variety of different errors propagate to uncertainty in the estimated physiological quantities. Conclusion: This work contributes to a better understanding of parameter estimation and AIF quantification in DCE-MRI. The proposed uncertainty estimation method can be used to efficiently calculate uncertainties in the parametric maps obtained in DCE-MRI.

This thesis deals with the magnetic resonance perfusion data analysis especially DCEMRI. In its introduction the thesis describes the problem of DCE-MRI data aquisition, the necessity of appropriate contrast agent and basic principles of perfusion analysis. The dynamic behavior of contrast agent vascular distribution can be described by arterial input function (AIF). The shape of the curves close to the area of interest is affected by dispersion which is called vascular transport function (VTF) due to the distribution of the contrast agent to the region of interest. Finally the tissue residual function describes system behavior of tissue. The practical part of the diploma thesis is aimed at implementation of model curves AIF, VTF and TRF. Furthermore, a simulation program was created for easy manipulation with introduced models moreover the program is used to perform an estimation of perfusion parameters based on nonblind deconvolution. The method is validated on synthetic data and illustrated on clinical data of the renal cell carcinoma patient.

La tomographie par émission de positons (TEP) est une méthode d’imagerie fonctionnelle, utilisée en particulier lors du développement de nouveaux médicaments et pour imager les tumeurs. En TEP, l’estimation de la concentration plasmatique artérielle d’activité du traceur non métabolisé (nommée « fonction d’entrée ») est nécessaire pour l’extraction des paramètres pharmacocinétiques. Ceux-ci permettent de quantifier le comportement du traceur dans les tissus, ou plus précisément le traitement du traceur par les tissus. Cette thèse constitue une contribution à l’étude de la fonction d’entrée, par l’élaboration d’une méthode d’estimation de la fonction d’entrée peu invasive à partir des images TEP et de prélèvements veineux. L’exemple du traceur FDG (analogue du glucose) dans le cerveau humain a été choisi. La méthode proposée repose sur la modélisation compartimentale de l’organisme : elle déconvolue le modèle à trois compartiments utilisé pour le FDG. L’originalité de la méthode repose sur trois points : l’utilisation d’un grand nombre de régions d’intérêt ; l’utilisation d’un grand nombre de jeux de trois régions d’intérêt différentes; une estimation itérative. Pour la validation de la méthode, un soin particulier a été porté à la simulation d’images TEP (simulation d’acquisition, reconstruction, corrections) de plus en plus réalistes, depuis une image simple simulée avec un simulateur analytique jusqu’à une image la plus proche possible de la réalité, simulée avec simulateur Monte-Carlo. Une chaîne de pré-traitement (segmentation des IRM associés, recalage entre images TEP et IRM et correction de l’effet de volume partiel par une variante de la méthode de Rousset) a ensuite été appliquée à ces images afin d’extraire les cinétiques des régions d’intérêt, données d’entrée de la méthode d’estimation de la fonction d’entrée. L’évaluation de la méthode sur différentes données, simulées et réelles, est présentée, ainsi que l’étude de la sensibilité de la méthode à différents facteurs tels que les erreurs de segmentation, de recalage, de mesure de l’activité des prélèvements sanguins. / Positron Emission Tomography (PET) is a method of functional imaging, used in particular for drug development and tumor imaging. In PET, the estimation of the arterial plasmatic activity concentration of the non-metabolized compound (the "input function") is necessary for the extraction of the pharmacokinetic parameters. These parameters enable the quantification of the compound dynamics in the tissues. This PhD thesis contributes to the study of the input function by the development of a minimally invasive method to estimate the input function. This method uses the PET image and a few blood samples. In this work, the example of the FDG tracer is chosen. The proposed method relies on compartmental modeling: it deconvoluates the three-compartment-model. The originality of the method consists in using a large number of regions of interest (ROIs), a large number of sets of three ROIs, and an iterative process. To validate the method, simulations of PET images of increasing complexity have been performed, from a simple image simulated with an analytic simulator to a complex image simulated with a Monte-Carlo simulator. After simulation of the acquisition, reconstruction and corrections, the images were segmented (through segmentation of an IRM image and registration between PET and IRM images) and corrected for partial volume effect by a variant of Rousset’s method, to obtain the kinetics in the ROIs, which are the input data of the estimation method. The evaluation of the method on simulated and real data is presented, as well as a study of the method robustness to different error sources, for example in the segmentation, in the registration or in the activity of the used blood samples.

<p> </p><p>Acquired data from dynamic contrast enhanced MRI measurements can be used to non-invasively assess tumour vascular characteristics through pharmacokinetic modelling. The modelling requires an arterial input function which is the concentration of contrast agent in the blood reaching the volume of interest as a function of time. The aim of this work is testing and optimizing a turboFLASH sequence to appraise its suitability for measuring the arterial input function by measuring phase.</p><p>Contrast concentration measurements in a phantom were done with both phase and relaxivity techniques. The results were compared to simulations of the experiment conditions to compare the conformance. The results using the phase technique were promising, and the method was carried on to in-vivo testing. The in-vivo data displayed a large signal loss which motivated a new phantom experiment to examine the cause of this signal reduction. Dynamic measurements were made in a phantom with pulsatile flow to mimic a blood vessel with a somewhat modified turboFLASH sequence. The conclusions drawn from analyzing the data were used to further improve the sequence and this modified turboFLASH sequence was tested in an in-vivo experiment. The obtained concentration curve showed significant improvement and was deemed to be a good representation of the true blood concentration.</p><p>The conclusion is that phase measurements can be recommended over relaxivity based measurements. This recommendation holds for using a slice selective saturation recovery turboFLASH sequence and measuring the arterial input function in the neck. Other areas of application need more thorough testing.</p><p> </p>

Acquired data from dynamic contrast enhanced MRI measurements can be used to non-invasively assess tumour vascular characteristics through pharmacokinetic modelling. The modelling requires an arterial input function which is the concentration of contrast agent in the blood reaching the volume of interest as a function of time. The aim of this work is testing and optimizing a turboFLASH sequence to appraise its suitability for measuring the arterial input function by measuring phase. Contrast concentration measurements in a phantom were done with both phase and relaxivity techniques. The results were compared to simulations of the experiment conditions to compare the conformance. The results using the phase technique were promising, and the method was carried on to in-vivo testing. The in-vivo data displayed a large signal loss which motivated a new phantom experiment to examine the cause of this signal reduction. Dynamic measurements were made in a phantom with pulsatile flow to mimic a blood vessel with a somewhat modified turboFLASH sequence. The conclusions drawn from analyzing the data were used to further improve the sequence and this modified turboFLASH sequence was tested in an in-vivo experiment. The obtained concentration curve showed significant improvement and was deemed to be a good representation of the true blood concentration. The conclusion is that phase measurements can be recommended over relaxivity based measurements. This recommendation holds for using a slice selective saturation recovery turboFLASH sequence and measuring the arterial input function in the neck. Other areas of application need more thorough testing.

The theoretical part of this work describes perfusion analysis of dynamic contrast enhanced magnetic resonance imaging from data acquisition to estimation of perfusion parameters. The main application fields are oncology, cardiology and neurology. The thesis is focused on quantitative perfusion analysis, specifically it contributes to solving of the the main challenge of this method – correct estimation of the contrast-agent concentration sequence in the arterial input of the region of interest (arterial input function). The goals of the thesis are stated based on literature review and on the expertise of our group. Blind deconvolution is selected as the method of choice. In the practical part of this thesis, a new method for arterial input function identification based on blind deconvolution is proposed. The method is designed for both preclinical and clinical applications. It was validated on synthetic, preclinical and clinical data. Furthermore, possibilities of the longer temporal sampling provided by blind deconvolution were analyzed. This can be used for improved spatial resolution and possibly for higher SNR. For easier deployment of the proposed methods into clinical and preclinical use, a software tool for perfusion data processing was designed.

No description available.

La séquence d’IRM de perfusion pondérée T1 après injection de Gadolinium (Gd), appelée dynamic contrast enhanced magnetic resonance imaging (DCE-MRI) fait partie du protocole d’IRM multiparamétrique (IRM-mp) réalisée pour le bilan d’extension du cancer prostatique (CaP). Le rationnel pour l’utilisation de cette séquence est d’une part le rôle capital de la néoangiogénèse dans le développement et la dissémination du CaP et d’autre part la possibilité d’imager l’angiogénèse in vivo et de façon non invasive. L’analyse quantitative nécessite un post-traitement multi-étapes complexe, dont le principe repose sur la modélisation pharmacocinétique (PC) de la biodistrubtion du Gd. Elle permet de calculer des paramètres PC reflétant la perméabilité capillaire et/ou la perfusion. Dans le CaP, ces paramètres PC ont montré leur potentiel pour évaluer l’agressivité tumorale, le pronostic, l’efficacité d’un traitement et/ou pour déterminer la dose efficace d’une nouvelle molécule anti-angiogéniques ou antivasculaires en cours de développement. Néanmoins, ils manquent de reproductibilité, notamment du fait des différentes techniques de quantifications utilisées par les logiciels de post-traitement.Nous avons développé au sein du laboratoire un outil de quantification capable de calculer une cartographie T1(0) à partir de la méthode des angles de bascule variables, nécessaire pour convertir les courbes du signal en courbe de concentration du Gd (Ct); de déterminer la fonction d’entrée artérielle (AIF – arterial input function) dans l’artère fémorale (Indivuduelle – Ind) ou lorsque cela n’était pas possible, d’utiliser une AIF issue de la littérature, telle que celle de Weinmann (W) ou de Fritz-Hansen (FH) ; et d’utiliser deux modèles PC, celui de Tofts et celui de Tofts modifié. Le logiciel a été validé sur des données simulées et sur une petite série clinique.Nous avons ensuite étudié l’impact du choix de la fonction d’entrée artériel sur les paramètres PC et notamment sur leur capacité à distinguer le CaP du tissu sain. 38 patients avec un CaP (>0,5cc) de la zone périphérique (ZP) ont été rétrospectivement inclus. Chaque patient avait bénéficié d’une IRM-mp sur laquelle deux régions d’intérêt (RI) : une tumorale et une bénigne ont été sélectionnées en utilisant une corrélation avec des cartes histo-morphométriques obtenues après prostatectomie radicale. En utilisant trois logiciels d’analyse quantitative différents, les valeurs moyennes de Ktrans (constante de transfert), ve (fraction du volume interstitiel) and vp (fraction du volume plasmatique) dans les RI ont été calculées avec trois AIF différentes (AIF Ind, AIF de W et AIF de FH). Ktrans était le paramètre PC qui permettait de mieux distinguer le CaP du tissu sain et ses valeurs étaient significativement supérieures dans le CaP, quelque soit l’AIF ou le logiciel. L’AIF de W donnait des aires sous les courbes ROC (Receiver Operating Characteristic) significativement plus grandes que l’AIF de FH (0.002≤p≤0.045) et plus grandes ou égales à l’AIF Ind (0.014≤p≤0.9). L’AIF Ind et de FH avaient des aires sous les courbes ROC comparables (0.34≤p≤0.81). Nous avons donc montré que les valeurs de Ktrans et sa capacité à distinguer CaP du tissu sain variaient significativement avec le choix de l’AIF et que les meilleures performances étaient obtenues avec l’AIF de W. / Dynamic contrast enhanced (DCE)-MRI is a T1 weighted sequence performed before, during and after a bolus injection of a contrast agent (CA). It is included in the multi-parametric prostate MRI (mp-MRI) protocol using to assess the extent of prostate cancer (PCa). The rationale for using DCE-MRI in PCa is that on one hand angiogenesis has been shown to play a central role in the PCa development and metastasis and on the other hand that DCE-MRI is a non invasive method able to depict this angiogenesis in vivo. The quantitative analysis of DCE-MRI data is a complex and multi-step process. The principle is to use a pharmacokinetic (PK) model reflecting the theoretical distribution of the CA in a tissue to extract PK parameters that describe the perfusion and capillary permeability. These parameters are of growing interest, especially in the field of oncology, for their use in assessing the aggressiveness, the prognosis and the efficacy of anti-angiogenic or anti-vascular treatments. The potential utility of these parameters is significant; however, the parameters often lack reproducibility, particularly between different quantitative analysis software programs.Firstly, we developed a quantitative analysis software solution using the variable flip angle method to estimate the T1 mapping which is needed to convert the signal-time curves to CA concentration-time curves; using three different arterial input functions (AIF): an individual AIF (Ind) measured manually in a large artery, and two literature population average AIFs of Weinmann (W) and of Fritz-Hansen (FH); and using two PK models (Tofts and modified Tofts). The robustness of the software programs was assessed on synthetic DCE-MRI data set and on a clinical DCE-MRI data set. Secondly, we assessed the impact of the AIF selection on the PK parameters to distinguish PCa from benign tissue. 38 patients with clinically important peripheral PCa (≥0.5cc) were retrospectively included. These patients underwent 1.5T multiparametric prostate MR with PCa and benign regions of interest (ROI) selected using a visual registration with morphometric reconstruction obtained from radical prostatectomy. Using three pharmacokinetic (PK) analysis software programs, the mean Ktrans, ve and vp of ROIs were computed using three AIFs: Ind-AIF, W-AIF and FH-AIF. The Ktrans provided higher area under the receiver operating characteristic curves (AUROCC) than ve and vp. The Ktrans was significantly higher in the PCa ROIs than in the benign ROIs. AUROCCs obtained with W-AIF were significantly higher than FH-AIF (0.002≤p≤0.045) and similar to or higher than Ind-AIF (0.014≤p≤0.9). Ind-AIF and FH-AIF provided similar AUROCC (0.34≤p≤0.81).We have then demonstrated that the selection of AIF can modify the capacity of the PK parameter Ktrans to distinguish PCa from benign tissue and that W-AIF yielded a similar or higher performance than Ind-AIF and a higher performance than FH-AIF.

This thesis deals with perfusion analysis using DCE-MRI (Dynamic contrast-enhanced magnetic resonance imaging). DCE-MRI is commonly used for microcirculation evaluation mainly in oncology and in recent years also in cardiology. The theoretical overview focuses on the issue of pharmacokinetic modeling and the estimation of perfusion parameters using selected models. The experimental part describes research software PerfLab and then it is aimed at the proposed program for synthetic data generation. Simulated data obtained under ideal conditions and in the presence of noise were used to compare models for the accuracy and reliability of DCE-MRI analysis.

Perfusion magnetic resonance imaging is modern diagnostic method used mainly in oncology. In this method, contrast agent is injected to the subject and then is continuously monitored the progress of its concentration in the affected area in time. Correct determination of the arterial input function (AIF) is very important for perfusion analysis. One possibility is to model AIF by multichannel blind deconvolution but the estimated AIF is necessary to be scaled. This master´s thesis is focused on description of scaling methods and their influence on perfussion parameters in dependence on used model of AIF in different tissues.