CT perfusion image processing: analysis of liver tumors

Perfusion CT imaging of the liver has potential to improve evaluation of tumour angiogenesis. Quantitative parameters can be obtained applying mathematical models to Time Attenuation Curve (TAC). However, there are still some difficulties for an accurate quantification of perfusion parameters due, for example, to algorithms employed, to mathematical model, to patient’s weight and cardiac output and to the acquisition system. In this thesis, new parameters and alternative methodologies about liver perfusion CT are presented in order to investigate the cause of variability of this technique. Firstly analysis were made to assess the variability related to the mathematical model used to compute arterial Blood Flow (BFa) values. Results were obtained implementing algorithms based on “ maximum slope method” and “Dual input one compartment model” . Statistical analysis on simulated data demonstrated that the two methods are not interchangeable. Anyway slope method is always applicable in clinical context. Then variability related to TAC processing in the application of slope method is analyzed. Results compared with manual selection allow to identify the best automatic algorithm to compute BFa. The consistency of a Standardized Perfusion Index (SPV) was evaluated and a simplified calibration procedure was proposed. At the end the quantitative value of perfusion map was analyzed. ROI approach and map approach provide related values of BFa and this means that pixel by pixel algorithm give reliable quantitative results. Also in pixel by pixel approach slope method give better results. In conclusion the development of new automatic algorithms for a consistent computation of BFa and the analysis and definition of simplified technique to compute SPV parameter, represent an improvement in the field of liver perfusion CT analysis.

Analysis of morphological and functional heterogeneity in ct perfusion images of lung tumours

Questa tesi si propone di innovare lo stato dell’arte dei metodi di analisi dell’eterogeneità in lesioni polmonari attualmente utilizzati, affiancando l’analisi funzionale (emodinamica) a quella morfologica, grazie allo sviluppo di nuove feature specifiche. Grazie alla collaborazione tra il Computer Vision Group (CVG) dell’Università di Bologna e l’Unità Operativa di Radiologia dell’IRCCS-IRST di Meldola (Istituto di Ricovero e Cura a Carattere Scientifico – Istituto Scientifico Romagnolo per lo Studio e la Cura dei Tumori), è stato possibile analizzare un adeguato numero di casi reali di pazienti affetti da lesioni polmonari primitive, effettuando un’analisi dell’eterogeneità sia su sequenze di immagini TC baseline sia contrast-enhanced, consentendo quindi un confronto tra eterogeneità morfologica e funzionale. I risultati ottenuti sono infine discussi sulla base del confronto con le considerazioni di natura clinica effettuate in cieco da due esperti radiologi dell’IRCCS-IRST.

Parahydrogen induced polarization on a clinical MRI system : polarization transfer of two spin order

Hyperpolarization techniques enhance the nuclear spin polarization and thus allow for new nuclear magnetic resonance applications like in vivo metabolic imaging. One of these techniques is Parahydrogen Induced Polarization (PHIP). It leads to a hyperpolarized 1H spin state which can be transferred to a heteronucleus like 13C by a radiofrequency (RF) pulse sequence. In this work, timing of such a sequence was analyzed and optimized for the molecule hydroxyethyl propionate. The pulse sequence was adapted for the work on a clinical magnetic resonance imaging (MRI) system which is usually equipped only with a single RF transmit channel. Optimal control theory optimizations were performed to achieve an optimized polarization transfer. A drawback of hyperpolarization is its limited lifetime due to relaxation processes. The lifetime can be increased by storing the hyperpolarization in a spin singlet state. The second part of this work therefore addresses the spin singlet state of the Cs-symmetric molecule dimethyl maleate which needs to be converted to the spin triplet state to be detectable. This conversion was realized on a clinical MRI system, both by field cycling and by two RF pulse sequences which were adapted and optimized for this purpose. Using multiple conversions enables the determination of the lifetime of the singlet state as well as the conversion efficiency of the RF pulse sequence. Both, the hyperpolarized 13C spin state and the converted singlet state were utilized for MR imaging. Careful choice of the echo time was shown to be crucial for both molecules.

Multimodal (EEG-fMRI) functional connectivity study of levodopa effect in Parkinson’s disease

Aim: To assess if the intake of levodopa in patients with Parkinson’s Disease (PD) changes cerebral connectivity, as revealed by simultaneous recording of hemodynamic (functional MRI, or fMRI) and electric (electroencephalogram, EEG) signals. Particularly, we hypothesize that the strongest changes in FC will involve the motor network, which is the most impaired in PD. Methods: Eight patients with diagnosis of PD “probable”, therapy with levodopa exclusively, normal cognitive and affective status, were included. Exclusion criteria were: moderate-severe rest tremor, levodopa induced dyskinesia, evidence of gray or white matter abnormalities on structural MRI. Scalp EEG (64 channels) were acquired inside the scanner (1.5 Tesla) before and after the intake of levodopa. fMRI functional connectivity was computed from four regions of interest: right and left supplementary motor area (SMA) and right and left precentral gyrus (primary motor cortex). Weighted partial directed coherence (w-PDC) was computed in the inverse space after the removal of EEG gradient and cardioballistic artifacts. Results and discussion: fMRI group analysis shows that the intake of levodopa increases hemodynamic functional connectivity among the SMAs / primary motor cortex and: sensory-motor network itself, attention network and default mode network. w-PDC analysis shows that EEG connectivity among regions of the motor network has the tendency to decrease after the intake the levodopa; furthermore, regions belonging to the DMN have the tendency to increase their outflow toward the rest of the brain. These findings, even if in a small sample of patients, suggest that other resting state physiological functional networks, beyond the motor one, are affected in patients with PD. The behavioral and cognitive tasks corresponding to the affected networks could benefit from the intake of levodopa.

Weighted geometric mean of large-scale matrices: numerical analysis and algorithms

Computing the weighted geometric mean of large sparse matrices is an operation that tends to become rapidly intractable, when the size of the matrices involved grows. However, if we are not interested in the computation of the matrix function itself, but just in that of its product times a vector, the problem turns simpler and there is a chance to solve it even when the matrix mean would actually be impossible to compute. Our interest is motivated by the fact that this calculation has some practical applications, related to the preconditioning of some operators arising in domain decomposition of elliptic problems. In this thesis, we explore how such a computation can be efficiently performed. First, we exploit the properties of the weighted geometric mean and find several equivalent ways to express it through real powers of a matrix. Hence, we focus our attention on matrix powers and examine how well-known techniques can be adapted to the solution of the problem at hand. In particular, we consider two broad families of approaches for the computation of f(A) v, namely quadrature formulae and Krylov subspace methods, and generalize them to the pencil case f(A\B) v. Finally, we provide an extensive experimental evaluation of the proposed algorithms and also try to assess how convergence speed and execution time are influenced by some characteristics of the input matrices. Our results suggest that a few elements have some bearing on the performance and that, although there is no best choice in general, knowing the conditioning and the sparsity of the arguments beforehand can considerably help in choosing the best strategy to tackle the problem.

The effect of diffusion gradient direction number on tractography of corticospinal tract in human brain: an along-tract analysis

Nel presente lavoro di tesi ho sviluppato un metodo di analisi di dati di DW-MRI (Diffusion-Weighted Magnetic Resonance Imaging)cerebrale, tramite un algoritmo di trattografia, per la ricostruzione del tratto corticospinale, in un campione di 25 volontari sani. Il diffusion tensor imaging (DTI) sfrutta la capacità del tensore di diffusione D di misurare il processo di diffusione dell’acqua, per stimare quantitativamente l’anisotropia dei tessuti. In particolare, nella sostanza bianca cerebrale la diffusione delle molecole di acqua è direzionata preferenzialmente lungo le fibre, mentre è ostacolata perpendicolarmente ad esse. La trattografia utilizza le informazioni ottenute tramite il DW imaging per fornire una misura della connettività strutturale fra diverse regioni del cervello. Nel lavoro si è concentrata l’attenzione sul fascio corticospinale, che è coinvolto nella motricità volontaria, trasmettendo gli impulsi dalla corteccia motoria ai motoneuroni del midollo spinale. Il lavoro si è articolato in 3 fasi. Nella prima ho sviluppato il pre-processing di immagini DW acquisite con un gradiente di diffusione sia 25 che a 64 direzioni in ognuno dei 25 volontari sani. Si è messo a punto un metodo originale ed innovativo, basato su “Regions of Interest” (ROIs), ottenute attraverso la segmentazione automatizzata della sostanza grigia e ROIs definite manualmente su un template comune a tutti i soggetti in esame. Per ricostruire il fascio si è usato un algoritmo di trattografia probabilistica che stima la direzione più probabile delle fibre e, con un numero elevato di direzioni del gradiente, riesce ad individuare, se presente, più di una direzione dominante (seconda fibra). Nella seconda parte del lavoro, ciascun fascio è stato suddiviso in 100 segmenti (percentili). Sono stati stimati anisotropia frazionaria (FA), diffusività media, probabilità di connettività, volume del fascio e della seconda fibra con un’analisi quantitativa “along-tract”, per ottenere un confronto accurato dei rispettivi percentili dei fasci nei diversi soggetti. Nella terza parte dello studio è stato fatto il confronto dei dati ottenuti a 25 e 64 direzioni del gradiente ed il confronto del fascio fra entrambi i lati. Dall’analisi statistica dei dati inter-subject e intra-subject è emersa un’elevata variabilità tra soggetti, dimostrando l’importanza di parametrizzare il tratto. I risultati ottenuti confermano che il metodo di analisi trattografica del fascio cortico-spinale messo a punto è risultato affidabile e riproducibile. Inoltre, è risultato che un’acquisizione con 25 direzioni di DTI, meglio tollerata dal paziente per la minore durata dello scan, assicura risultati attendibili. La principale applicazione clinica riguarda patologie neurodegenerative con sintomi motori sia acquisite, quali sindromi parkinsoniane sia su base genetica o la valutazione di masse endocraniche, per la definizione del grado di contiguità del fascio. Infine, sono state poste le basi per la standardizzazione dell’analisi quantitativa di altri fasci di interesse in ambito clinico o di studi di ricerca fisiopatogenetica.

Modelli compartimentali per l'analisi di segnali ottenuti tramite DCE MRI per studi di perfusione epatica

Il presente lavoro di tesi si inserisce all'interno di uno studio dal titolo: "Strategia di posizionamento multi-step come approccio pragmatico per ridurre il rischio di encefalopatia epatica post-TIPS (shunt trans-giugulare porto-sistemico intraepatico) in pazienti cirrotici con ascite refrattaria". Il progetto di tesi si è concentrato sull'analisi dei segnali ottenuti tramite DCE MRI, con lo scopo di implementare in ambiente MatLab due modelli differenti (Dual input - Mono compartment e Dual input - Dual compartment) che descrivono la cinetica del tracciante all'interno del sistema vascolare epatico e valutare l'efficacia dei parametri di perfusione associati nella descrizione delle variazioni in termini di microcircolazione introdotte dall'inserimento del TIPS. Inizialmente si sono voluti valutare, tramite simulazione, gli effetti in termini di amplificazione del rumore e stima dei parametri perfusionali dell'approssimazione lineare nella conversione da intensità di segnale MR a concentrazione di mezzo di contrasto. Successivamente, sempre attraverso simulazioni, per entrambi i modelli considerati è stato scelto uno schema di model-fitting e quindi testata l'affidabilità in termini di accuratezza e precisione delle stime dei parametri ottenute in funzione del livello di rumore associato alle curve di intensità di segnale. Parallelamente all'implementazione dei modelli per la stima di parametri di perfusione, sono stati realizzati dei phantom con l'obiettivo di simulare il parenchima epatico prima e dopo l'arrivo del mezzo di contrasto e poter testare la sequenza utilizzata durante l'acquisizione dei dati su paziente. Infine sono stati considerati gli esami di DCE MRI effettuati su un campione di nove pazienti pre e post-TIPS, utilizzando per l'analisi dei segnali entrambi i modelli implementati in fase di simulazione e successivamente valutando le variazioni nel valori associati ai parametri di perfusione introdotte dall'inserimento del TIPS.

Geometria frattale: analisi della lacunarita cerebrale in MRI

La geometria frattale descrive la complessità strutturale di oggetti che presentano, entro certi limiti, invarianza a fattori di scala. Obiettivo di questa tesi è l’analisi di indici frattali della morfologia cerebrale e cerebellare da immagini di risonanza magnetica (MRI) pesate T1 e della loro correlazione con l’età. A tale scopo sono state analizzate la dimensione frattale (D0) e la lacunarità (λs), indice di eterogeneità strutturale, della sostanza grigia (GM) e bianca (WM), calcolate mediante algoritmi di box counting e di differential gliding box, implementati in linguaggio C++, e regressione lineare con scelta automatica delle scale spaziali. Gli algoritmi sono stati validati su fantocci 3D ed è stato proposto un metodo per compensare la dipendenza di λs dalle dimensioni dell’immagine e dalla frazione di immagine occupata. L’analisi frattale è stata applicata ad immagini T1 a 3T del dataset ICBM (International Consortium for Brain Mapping) composto da 86 soggetti (età 19-85 anni). D0 e λs sono state rispettivamente 2.35±0.02 (media±deviazione standard) e 0.41±0.05 per la GM corticale, 2.34±0.03 e 0.35±0.05 per la WM cerebrale, 2.19±0.05 e 0.17±0.02 per la GM cerebellare, 1.95±0.06 e 0.30±0.04 per la WM cerebellare. Il coefficiente di correlazione lineare tra età e D0 della GM corticale è r=−0.38 (p=0.003); tra età e λs, r=0.72 (p<0.001) (mostrando che l’eterogeneità strutturale aumenta con l’invecchiamento) e tra età e λs compensata rispetto al volume della GM cerebrale (GMV), r=0.51 (p<0.001), superiore in valore assoluto a quello tra età e GMV (r=−0.45, p<0.001). In un modello di regressione lineare multipla, dove l’età è stata modellata da D0, λ e GMV della GM corticale, λs è risultato l’unico predittore significativo (r parziale=0.62, p<0.001). La lacunarità λs è un indice sensibile alle variazioni strutturali dovute all’invecchiamento cerebrale e si candida come biomarcatore nella valutazione della complessità cerebrale nelle malattie neurodegenerative.

Preoperative MRI findings and functional outcome after selective dorsal rhizotomy in children with bilateral spasticity

PURPOSE: To identify MRI characteristics that may predict the functional effect of selective dorsal rhizotomy (SDR) in children with bilateral spastic paresis. METHODS: We performed SDR in a group of 36 patients. The gross motor functioning measure-66 (GMFM-66) was applied before and after SDR. Available cerebral MRIs were retrospectively classified into three diagnostic groups: periventricular leucomalacia (PVL; n = 10), hydrocephalus (n = 2), and normal (n = 6). In patients with PVL, we scored the severity of the MR abnormalities. We compared the changes in the GMFM-66 after SDR in the diagnostic groups. In patients with PVL, we correlated the severity of the MR abnormalities with the changes in the GMFM-66. RESULTS: The mean follow-up period was 5 years and 4 months (range, 1 year and 1 month to 9 years). The best improvement in gross motor function was observed in patients with normal MRI, and the slightest improvement was observed in patients with hydrocephalus. The severity of the PVL did correlate with the GMFM-66 score before SDR but not with the functional effect of SDR. CONCLUSION: We conclude that with respect to gross motor skills, the improvements after SDR are good in patients with no MRI abnormalities. In the patients with hydrocephalus, the improvements after SDR were insignificant. In patients with PVL, the improvements were intermediate and did not correlate with the degree of PVL.

Intracranial pressure pulse amplitude during changes in head elevation: a new parameter for determining optimum cerebral perfusion pressure?

During short-term postural changes, the factors determining the amplitude of intracranial pulse pressure (ICPPA) remain constant, except for cerebrovascular resistance (CVR). Therefore, it may be possible to draw conclusions from the ICPPA onto the cerebrovascular resistance (CVR) and thus the relative change in cerebral perfusion pressure (CPP).

Preservation of Amputated Extremities by Extracorporeal Blood Perfusion; a Feasibility Study in a Porcine Model

Successful extremity transplantations and replantations have to be performed within 6 h of amputation to avoid irreversible tissue loss. This study investigates ex vivo the technical feasibility and the limb preservation potential of extracorporeal whole blood perfusion in a porcine model.

Comparison of physiological triggering schemes for diffusion-weighted magnetic resonance imaging in kidneys

To determine the potential benefit of combined respiratory-cardiac triggering for diffusion-weighted imaging (DWI) of kidneys compared to respiratory triggering alone (RT).

Cancer of the oral cavity and oropharynx

Tumours in the oral cavity and oropharynx differ in presentation and prognosis and the detection of spread of tumour from one subsite to another is essential for the T-staging. This article reviews the anatomy and describes the pattern of spread of different cancers arising in the oral cavity and oropharynx; the imaging findings on computerized tomography and magnetic resonance imaging are also described. Brief mention is made on the role of newer imaging modalities such as [(18)F]fluorodeoxyglucose-positron emission tomography/computed tomography, perfusion studies and diffusion-weighted magnetic resonance imaging.

Effects of lung recruitment maneuvers on splanchnic organ perfusion during endotoxin-induced pulmonary arterial hypertension

Lung recruitment maneuvers (RMs), used to reopen atelectatic lung units and to improve oxygenation during mechanical ventilation, may result in hemodynamic impairment. We hypothesize that pulmonary arterial hypertension aggravates the consequences of RMs in the splanchnic circulation. Twelve anesthetized pigs underwent laparotomy and prolonged postoperative ventilation. Systemic, regional, and organ blood flows were monitored. After 6 h (= baseline), a recruitment maneuver was performed with sustained inflation of the lungs. Thereafter, the pigs were randomly assigned to group C (control, n = 6) or group E with endotoxin-induced pulmonary arterial hypertension (n = 6). Endotoxemia resulted in a normotensive and hyperdynamic state and a deterioration of the oxygenation index by 33%. The RM was then repeated in both groups. Pulmonary artery pressure increased during lipopolysaccharide infusion from 17 ± 2 mmHg (mean ± SD) to 31 ± 10 mmHg and remained unchanged in controls (P < 0.05). During endotoxemia, RM decreased aortic pulse pressure from 37 ± 14 mmHg to 27 ± 13 mmHg (mean ± SD, P = 0.024). The blood flows of the renal artery, hepatic artery, celiac trunk, superior mesenteric artery, and portal vein decreased to 71% ± 21%, 69% ± 20%, 76% ± 16%, 79% ± 18%, and 81% ± 12%, respectively, of baseline flows before RM (P < 0.05 all). Organ perfusion of kidney cortex, kidney medulla, liver, and jejunal mucosa in group E decreased to 65% ± 19%, 77% ± 13%, 66% ± 26%, and 71% ± 12%, respectively, of baseline flows (P < 0.05 all). The corresponding recovery to at least 90% of baseline regional blood flow and organ perfusion lasted 1 to 5 min. Importantly, the decreases in regional blood flows and organ perfusion and the time to recovery of these flows did not differ from the controls. In conclusion, lipopolysaccharide-induced pulmonary arterial hypertension does not aggravate the RM-induced significant but short-lasting decreases in systemic, regional, and organ blood flows.

Automated Intervertebral Disc Detection from Low Resolution, Sparse MRI Images for the Planning of Scan Geometries

Robust and accurate identification of intervertebral discs from low resolution, sparse MRI scans is essential for the automated scan planning of the MRI spine scan. This paper presents a graphical model based solution for the detection of both the positions and orientations of intervertebral discs from low resolution, sparse MRI scans. Compared with the existing graphical model based methods, the proposed method does not need a training process using training data and it also has the capability to automatically determine the number of vertebrae visible in the image. Experiments on 25 low resolution, sparse spine MRI data sets verified its performance.