Intensity Segmentation of the Human Brain with Tissue dependent Homogenization

High-precision segmentation of the human cerebral cortex based on T1-weighted MRI is still a challenging task. When opting to use an intensity based approach, careful data processing is mandatory to overcome inaccuracies. They are caused by noise, partial volume effects and systematic signal intensity variations imposed by limited homogeneity of the acquisition hardware. We propose an intensity segmentation which is free from any shape prior. It uses for the first time alternatively grey (GM) or white matter (WM) based homogenization. This new tissue dependency was introduced as the analysis of 60 high resolution MRI datasets revealed appreciable differences in the axial bias field corrections, depending if they are based on GM or WM. Homogenization starts with axial bias correction, a spatially irregular distortion correction follows and finally a noise reduction is applied. The construction of the axial bias correction is based on partitions of a depth histogram. The irregular bias is modelled by Moody Darken radial basis functions. Noise is eliminated by nonlinear edge preserving and homogenizing filters. A critical point is the estimation of the training set for the irregular bias correction in the GM approach. Because of intensity edges between CSF (cerebro spinal fluid surrounding the brain and within the ventricles), GM and WM this estimate shows an acceptable stability. By this supervised approach a high flexibility and precision for the segmentation of normal and pathologic brains is gained. The precision of this approach is shown using the Montreal brain phantom. Real data applications exemplify the advantage of the GM based approach, compared to the usual WM homogenization, allowing improved cortex segmentation

Improving the clinico-radiological association in neurological diseases

Despite the key role of magnetic resonance imaging (MRI) in the diagnosis and monitoring of multiple sclerosis (MS) and cerebral small vessel disease (SVD), the association between clinical and radiological disease manifestations is often only moderate, limiting the use of MRI-derived markers in the clinical routine or as endpoints in clinical trials. In the projects conducted as part of this thesis, we addressed this clinico-radiological gap using two different approaches. Lesion-symptom association: In two voxel-based lesion-symptom mapping studies, we aimed at strengthening lesion-symptom associations by identifying strategic lesion locations. Lesion mapping was performed in two large cohorts: a dataset of 2348 relapsing-remitting MS patients, and a population-based cohort of 1017 elderly subjects. T2-weighted lesion masks were anatomically aligned and a voxel-based statistical approach to relate lesion location to different clinical rating scales was implemented. In the MS lesion mapping, significant associations between white matter (WM) lesion location and several clinical scores were found in periventricular areas. Such lesion clusters appear to be associated with impairment of different physical and cognitive abilities, probably because they affect commissural and long projection fibers. In the SVD lesion mapping, the same WM fibers and the caudate nucleus were identified to significantly relate to the subjects’ cerebrovascular risk profiles, while no other locations were found to be associated with cognitive impairment. Atrophy-symptom association: With the construction of an anatomical physical phantom, we aimed at addressing reliability and robustness of atrophy-symptom associations through the provision of a “ground truth” for atrophy quantification. The built phantom prototype is composed of agar gels doped with MRI and computed tomography (CT) contrast agents, which realistically mimic T1 relaxation times of WM and grey matter (GM) and showing distinguishable attenuation coefficients using CT. Moreover, due to the design of anatomically simulated molds, both WM and GM are characterized by shapes comparable to the human counterpart. In a proof-of-principle study, the designed phantom was used to validate automatic brain tissue quantification by two popular software tools, where “ground truth” volumes were derived from high-resolution CT scans. In general, results from the same software yielded reliable and robust results across scans, while results across software were highly variable reaching volume differences of up to 8%

Brain-Blood Partition Coefficient and Cerebral Blood Flow in Canines Using Calibrated Short TR Recovery (CaSTRR) Correction Method.

The brain-blood partition coefficient (BBPC) is necessary for quantifying cerebral blood flow (CBF) when using tracer based techniques like arterial spin labeling (ASL). A recent improvement to traditional MRI measurements of BBPC, called calibrated short TR recovery (CaSTRR), has demonstrated a significant reduction in acquisition time for BBPC maps in mice. In this study CaSTRR is applied to a cohort of healthy canines (n = 17, age = 5.0 - 8.0 years) using a protocol suited for application in humans at 3T. The imaging protocol included CaSTRR for BBPC maps, pseudo-continuous ASL for CBF maps, and high resolution anatomical images. The standard CaSTRR method of normalizing BBPC to gadolinium-doped deuterium oxide phantoms was also compared to normalization using hematocrit (Hct) as a proxy value for blood water content. The results show that CaSTRR is able to produce high quality BBPC maps with a 4 min acquisition time. The BBPC maps demonstrate significantly higher BBPC in gray matter (0.83 ± 0.05 mL/g) than in white matter (0.78 ± 0.04 mL/g, p = 0.006). Maps of CBF acquired with pCASL demonstrate a negative correlation between gray matter perfusion and age (p = 0.003). Voxel-wise correction for BBPC is also shown to improve contrast to noise ratio between gray and white matter in CBF maps. A novel aspect of the study was to show that that BBPC measurements can be calculated based on the known Hct of the blood sample placed in scanner. We found a strong correlation (R 2 = 0.81 in gray matter, R 2 = 0.59 in white matter) established between BBPC maps normalized to the doped phantoms and BBPC maps normalized using Hct. This obviates the need for doped water phantoms which simplifies both the acquisition protocol and the post-processing methods. Together this suggests that CaSTRR represents a feasible, rapid method to account for BBPC variability when quantifying CBF. As canines have been used widely for aging and Alzheimer's disease studies, the CaSTRR method established in the animals may further improve CBF measurements and advance our understanding of cerebrovascular changes in aging and neurodegeneration

Estimation of Fiber Orientations Using Neighborhood Information

Data from diffusion magnetic resonance imaging (dMRI) can be used to reconstruct fiber tracts, for example, in muscle and white matter. Estimation of fiber orientations (FOs) is a crucial step in the reconstruction process and these estimates can be corrupted by noise. In this paper, a new method called Fiber Orientation Reconstruction using Neighborhood Information (FORNI) is described and shown to reduce the effects of noise and improve FO estimation performance by incorporating spatial consistency. FORNI uses a fixed tensor basis to model the diffusion weighted signals, which has the advantage of providing an explicit relationship between the basis vectors and the FOs. FO spatial coherence is encouraged using weighted l1-norm regularization terms, which contain the interaction of directional information between neighbor voxels. Data fidelity is encouraged using a squared error between the observed and reconstructed diffusion weighted signals. After appropriate weighting of these competing objectives, the resulting objective function is minimized using a block coordinate descent algorithm, and a straightforward parallelization strategy is used to speed up processing. Experiments were performed on a digital crossing phantom, ex vivo tongue dMRI data, and in vivo brain dMRI data for both qualitative and quantitative evaluation. The results demonstrate that FORNI improves the quality of FO estimation over other state of the art algorithms.Comment: Journal paper accepted in Medical Image Analysis. 35 pages and 16 figure

3D medical volume segmentation using hybrid multiresolution statistical approaches

This article is available through the Brunel Open Access Publishing Fund. Copyright © 2010 S AlZu’bi and A Amira.3D volume segmentation is the process of partitioning voxels into 3D regions (subvolumes) that represent meaningful physical entities which are more meaningful and easier to analyze and usable in future applications. Multiresolution Analysis (MRA) enables the preservation of an image according to certain levels of resolution or blurring. Because of multiresolution quality, wavelets have been deployed in image compression, denoising, and classification. This paper focuses on the implementation of efficient medical volume segmentation techniques. Multiresolution analysis including 3D wavelet and ridgelet has been used for feature extraction which can be modeled using Hidden Markov Models (HMMs) to segment the volume slices. A comparison study has been carried out to evaluate 2D and 3D techniques which reveals that 3D methodologies can accurately detect the Region Of Interest (ROI). Automatic segmentation has been achieved using HMMs where the ROI is detected accurately but suffers a long computation time for its calculations

Physical and digital phantoms for validating tractography and assessing artifacts

Fiber tractography is widely used to non-invasively map white-matter bundles in vivo using diffusion-weighted magnetic resonance imaging (dMRI). As it is the case for all scientific methods, proper validation is a key prerequisite for the successful application of fiber tractography, be it in the area of basic neuroscience or in a clinical setting. It is well-known that the indirect estimation of the fiber tracts from the local diffusion signal is highly ambiguous and extremely challenging. Furthermore, the validation of fiber tractography methods is hampered by the lack of a real ground truth, which is caused by the extremely complex brain microstructure that is not directly observable non-invasively and that is the basis of the huge network of long-range fiber connections in the brain that are the actual target of fiber tractography methods. As a substitute for in vivo data with a real ground truth that could be used for validation, a widely and successfully employed approach is the use of synthetic phantoms. In this work, we are providing an overview of the state-of-the-art in the area of physical and digital phantoms, answering the following guiding questions: “What are dMRI phantoms and what are they good for?”, “What would the ideal phantom for validation fiber tractography look like?” and “What phantoms, phantom datasets and tools used for their creation are available to the research community?”. We will further discuss the limitations and opportunities that come with the use of dMRI phantoms, and what future direction this field of research might take

IMAGE PROCESSING, SEGMENTATION AND MACHINE LEARNING MODELS TO CLASSIFY AND DELINEATE TUMOR VOLUMES TO SUPPORT MEDICAL DECISION

Techniques for processing and analysing images and medical data have become the main’s translational applications and researches in clinical and pre-clinical environments. The advantages of these techniques are the improvement of diagnosis accuracy and the assessment of treatment response by means of quantitative biomarkers in an efficient way. In the era of the personalized medicine, an early and efficacy prediction of therapy response in patients is still a critical issue. In radiation therapy planning, Magnetic Resonance Imaging (MRI) provides high quality detailed images and excellent soft-tissue contrast, while Computerized Tomography (CT) images provides attenuation maps and very good hard-tissue contrast. In this context, Positron Emission Tomography (PET) is a non-invasive imaging technique which has the advantage, over morphological imaging techniques, of providing functional information about the patient’s disease. In the last few years, several criteria to assess therapy response in oncological patients have been proposed, ranging from anatomical to functional assessments. Changes in tumour size are not necessarily correlated with changes in tumour viability and outcome. In addition, morphological changes resulting from therapy occur slower than functional changes. Inclusion of PET images in radiotherapy protocols is desirable because it is predictive of treatment response and provides crucial information to accurately target the oncological lesion and to escalate the radiation dose without increasing normal tissue injury. For this reason, PET may be used for improving the Planning Treatment Volume (PTV). Nevertheless, due to the nature of PET images (low spatial resolution, high noise and weak boundary), metabolic image processing is a critical task. The aim of this Ph.D thesis is to develope smart methodologies applied to the medical imaging field to analyse different kind of problematic related to medical images and data analysis, working closely to radiologist physicians. Various issues in clinical environment have been addressed and a certain amount of improvements has been produced in various fields, such as organs and tissues segmentation and classification to delineate tumors volume using meshing learning techniques to support medical decision. In particular, the following topics have been object of this study: • Technique for Crohn’s Disease Classification using Kernel Support Vector Machine Based; • Automatic Multi-Seed Detection For MR Breast Image Segmentation; • Tissue Classification in PET Oncological Studies; • KSVM-Based System for the Definition, Validation and Identification of the Incisinal Hernia Reccurence Risk Factors; • A smart and operator independent system to delineate tumours in Positron Emission Tomography scans; 3 • Active Contour Algorithm with Discriminant Analysis for Delineating Tumors in Positron Emission Tomography; • K-Nearest Neighbor driving Active Contours to Delineate Biological Tumor Volumes; • Tissue Classification to Support Local Active Delineation of Brain Tumors; • A fully automatic system of Positron Emission Tomography Study segmentation. This work has been developed in collaboration with the medical staff and colleagues at the: • Dipartimento di Biopatologia e Biotecnologie Mediche e Forensi (DIBIMED), University of Palermo • Cannizzaro Hospital of Catania • Istituto di Bioimmagini e Fisiologia Molecolare (IBFM) Centro Nazionale delle Ricerche (CNR) of Cefalù • School of Electrical and Computer Engineering at Georgia Institute of Technology The proposed contributions have produced scientific publications in indexed computer science and medical journals and conferences. They are very useful in terms of PET and MRI image segmentation and may be used daily as a Medical Decision Support Systems to enhance the current methodology performed by healthcare operators in radiotherapy treatments. The future developments of this research concern the integration of data acquired by image analysis with the managing and processing of big data coming from a wide kind of heterogeneous sources