In translational medicine, complementary functional and morphological imaging techniques are used extensively to observe physiological processes in vivo and to assess structural changes as a result of disease progression. The combination of magnetic resonance imaging (MRI) and positron emission tomography (PET) provides excellent soft tissue contrast from MRI with exceptional sensitivity and specificity from PET. This thesis explores the use of sequentially acquired PET and MR images to improve the quantification of small animal PET data. The primary focus was to improve image-based estimates of the arterial input function (AIF), which defines the amount of PET tracer within blood plasma over time. The AIF is required to produce physiological parameters quantifying key processes such as metabolism or perfusion from dynamic PET images. The gold standard for AIF measurement, however, requires serial blood sampling over the course of a PET scan, which is invasive in rat studies but prohibitive in mice due to small total blood volumes. To address this issue, the geometric transfer matrix (GTM) and recovery coefficient (RC) techniques were applied using anatomical MR images to enable the extraction of partial volume corrected image based AIFs from mouse PET images. A non-invasive AIF extraction method was also developed for rats, beginning with the optimization of an automated voxel selection algorithm to assist in extracting MR contrast agent signal time courses from dynamic susceptibility contrast (DSC) MRI data. This procedure was then combined with dynamic contrast enhanced (DCE) MRI to track a combined injection of Gadolinium-based contrast agent and PET tracer through the rat brain. By comparison with gold standard tracer blood sample data, it was found that normalized MRI-based AIFs could be successfully converted into PET tracer AIFs in the first pass phase when fitted with gamma variate functions. Finally, a MR image segmentation method used to provide PET attenuation correction in mice was validated using the Cambridge split magnet PET/MR scanner?s transmission scanning capabilities. This work recommends that contributions from MR hardware in the PET field of view must be accounted forto gain accurate estimates of tracer uptake and standard uptake values (SUVs). This thesis concludes that small animal MR data taken in the same imaging session can provide non-invasive methods to improve PET image quantification, giving added value to combined PET/MR studies over those conducted using PET alone.

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.

This master thesis deals with the contrast agents and their application in the ultrasound perfusion analysis. It is focused on Bolus & Burst method which, as a combination of two approaches that have been used so far, allows an absolute quantification of perfusion parameters in the region of interest. Contrast agent concentration time sequence is modeled as a convolution of the parametrically defined arterial input function and the tissue residual funkction. Thesis discusses different mathematical models of these functions as well as the methods of the parameters estimation. The methods functionality is validated on simulated and also preclinical data.

La modélisation pharmacocinétique en tomographie par émission par positrons (TEP) permet d’estimer les paramètres physiologiques liés à l’accumulation dynamique d’un radiotraceur. Les paramètres estimés sont biaisés par le bruit dans les images TEP dynamiques durant l’ajustement des courbes d’activité des tissus, plus communément appelées TAC de l’anglais Time Activity Curve. La qualité des images TEP dynamiques est limitée par la statistique de comptage et influencée par les paramètres de reconstruction choisis en termes de résolution spatiale et temporelle. Il n’existe pas de recommandations claires pour les paramètres de reconstruction à utiliser pour les images dynamiques TEP. L’objectif de ce projet de maitrise est d’évaluer le biais dans l’estimation des paramètres pharmacocinétiques afin de trouver les paramètres de reconstruction TEP les plus optimaux en termes de résolution spatiale et de niveau de bruit. Plus précisément, ce projet cherche à déterminer quel modèle d’AIF offre les meilleurs ajustements, mais aussi quel modèle de poids permet la meilleure estimation des paramètres pharmacocinétiques pour le modèle à deux compartiments. Ce faisant, il serait possible de mieux planifier la reconstruction d’images TEP dynamique et potentiellement améliorer leur résolution spatiale. Afin de tester les biais dans les paramètres pharmacocinétiques sous différents niveaux de bruit, un objet de référence numérique (DRO) avec les informations trouvées dans la littérature sera construit. Ensuite, des simulations numériques seront effectuées avec ce DRO afin de trouver les paramètres de reconstruction et le niveau de bruit le plus optimal. Un biais réduit des paramètres pharmacocinétiques et une meilleure résolution spatiale des images TEP dynamique permettrait de détecter des cancers ou tumeurs à des stades moins avancés de la maladie, permettant potentiellement un traitement plus efficace et avec moins de séquelles et d’effets secondaires pour les patients. En outre, cela permettrait aussi de visualiser l’hétérogénéité des tumeurs. / Pharmacokinetic models in positron emission tomography (PET) allow for the estimation of physiological parameters linked to the dynamic accumulation of a radiotracer. Estimated parameters are biased by noise in dynamic PET images during the fitting of Time Activity Curves (TAC). Image quality in dynamic PET is limited by counting statistics and influenced by the chosen reconstruction parameters in terms of spatial and temporal resolution. Clear recommendations and guidelines for the reconstruction parameters that should be used do not exist at the moment for dynamic PET. The goal of this masters project is to evaluate the bias in the pharmacokinetic parameters estimation to find the optimal PET reconstruction parameters in terms of spatial resolution and noise levels. More precisely, this project aims to determine which AIF model produces the best fits, but also which weight noise model allows for the best parameters estimation with the two compartment model. It would then be possible to plan the PET image reconstruction more finely and potentially improve spatial resolution. To test the pharmacokinetic parameters’ biases under different noise levels, a Digital Reference Object (DRO) with information and specifications found from the litterature will be built. Then, numerical simulations will be done with that DRO to find the optimal noise level and value for the pharmacokinetic parameter. A reduced bias in these parameters and an improved spatial resolution would allow the detection of tumors or lesions at earlier stages, which could potentially allow for a more potent treatment with less short and long term side effects. It would also allow the visualization and quantification of lesion heterogeneity.

Cerebral blood flow quantification is a vital diagnostic tool for disease monitoring and used for diagnosing a variation of pathological conditions. The human brain requires roughly about 20 % of the total cardiac output to sustain normal functioning, hence the perfusion of blood is an important factor to deliver oxygenated blood. The golden standard for quantifying the cerebral blood flow follows by measurement with dynamic positron emission tomography of 15O-labelled water modelled by tracer kinetic compartments. For implementation, knowledge of an input function must exist which is in general being sampled through arterial cannulation of the radial artery with a continuous sampling instrument. The core of this thesis is to establish if two sampling instruments contradicts in comparison to each other when sampling the data to the input function.  In total 22 subjects underwent a 10-minute dynamic  15O-labeled water brain PET scan on two imaging modalities PET/CT and PET/MR. Continuous arterial sampling was performed either by a Veenstra on PET/CT or a Swisstrace on PET/MR during a baseline scan. In two subjects the two sampling instruments were coupled in series and imaged solely on the PET/CT.  Cerebral blood flow analysis was done comparing varying dispersion times, the two imaging modalities compared each other and comparing the calculated and measured blood flows obtained through this study with the values obtained prior. To be able to compare the values showing inconsistency to the values obtained through this thesis, a comparison between two different iterative reconstruction methods was done. Here the method of ordered subsets expectation maximum was compared to a Bayesian penalized-likelihood method. To further compare the two sampling instruments an image derived input function was constructed and compared with the blood sampled input function. The results showed that there was no significant difference between measured cerebral blood flow between the two imaging modalities with the currently used reconstruction method based on Bayesian penalized likelihood but presented in the earlier data there was an inconsistency. A dispersion analysis with variation on the external dispersion time shows that if the time was chosen to low or to high compared to the standard time used it introduced distorted fitted models of the activity curves. This distortion creates further errors in the calculation of the cerebral blood flow, however with the analysis the standard dispersion time could be confirmed as an accurate fit. Subjects imaged with the two sampling instruments in series showed no significant difference except for the measured values on Veenstra to be slightly higher. Lastly the correlation between the image derived input function and the blood sampled input function showed poorly performance. Only a R2 value of 0.42 was achieved on the PET/CT while a meagre R2 value of 0.18 was achieved on the PET/MR. Although the correlation was poorly, the plotted activity curves from the two functions showed a representable appearance between each other. / Kvantifiering av det cerebrala blodflödet är ett nödvändigt diagnostiskt verktyg som används för att kontrollera och diagnostisera en variation av patologiska sjukdomstillstånd. Den mänskliga hjärnan kräver kring 20 %av den totala produktionen från hjärtat för att upprätthålla normal funktion, följaktligen är perfusion av blod en viktig faktor för att distribuera syrerikt blod runt om i kroppen. Den gyllene standarden för kvantifiering av det cerebrala blodflödet följer som undersökning med dynamisk positron emission tomografi av 15O-märkt vatten, modellerat  med kinetisk kompartment teori. För att kunna implementera detta måste information om en input-funktion erhållas, generellt erhålls detta genom att blod tags genom arteriell kanylering av antingen den radiella artären med ett kontinuerligt samplings instrument. Målet med detta arbete är att fastställa om två samplings instrument motsäger varandra vid mätning av data till input-funktionen. Totalt deltagande är 22 patienter som genomgick en 10-minuters dynamisk 15O-märkt vatten PET undersökning av hjärnan på två bildtagningsmodaliteter PET/CT och PET/MR. Kontinuerlig blodtagning genomfördes antingen med en Veenstra sampler instrument på PET/CT eller en Swisstrace sampler instrument på PET/MR tillsammans med en baseline undersökning. Vid två undersökningar seriekopplades de två instrumenten och patienterna blev endast undersökta vid PET/CT. För ytterligare kunna utvärdera de två instrumenten, konstruerades en bild framtagen input-funktion som sedan kunde jämföras med den blod samplade input-funktionen. Cerebrala blodflödes analyser gjordes med olika dispersions tider, även för att kunna jämföra de två bildtagningsmodaliteterna mot varandra och jämföra erhållna värden framtagna under denna studie med en tidigare studie. För att kunna jämföra avvikelserna i de uppmätta värdena har även två olika rekonstruerings metoder studerats. Resultaten visar ingen signifikant skillnad mellan de uppmätta cerebrala blodflödena mellan de två bildtagningsmodaliteterna rekonstruerade med den nuvarande standarden. Dispersions analysen med varierande extern dispersions tid visar att om tiden är för kort eller för lång jämfört med standardtiden, introduceras en osann anpassning av aktivitets kurvorna. Denna förvrängning av datat resulterar till fler avvikelser i beräkningarna av blodflödet, likväl var det möjligt att bekräfta standardtiderna som används. Patienter som undersöktes med instrumenten i seriekoppling visade ingen signifikant skillnad förutom att det uppmättes en aningens högre värden hos patienter med Veenstra som blod sampler. Slutligen, korrelationen mellan den bild framtagna input-funktionen och den blod samplade input-funktionen visade ett dåligt resultat. Endast ett R2 värde på 0.42 erhölls för PET/CT medan endast ett R2 värde på 0.18 på PET/MR erhölls. Trotts att korrelationen var dålig, visade de plottade aktivitets kurvorna ett representativt utseende mellan de två typerna av input funktion.