Comparison of Distances for Supervised Segmentation of White Matter Tractography

Tractograms are mathematical representations of the main paths of axons within the white matter of the brain, from diffusion MRI data. Such representations are in the form of polylines, called streamlines, and one streamline approximates the common path of tens of thousands of axons. The analysis of tractograms is a task of interest in multiple fields, like neurosurgery and neurology. A basic building block of many pipelines of analysis is the definition of a distance function between streamlines. Multiple distance functions have been proposed in the literature, and different authors use different distances, usually without a specific reason other than invoking the "common practice". To this end, in this work we want to test such common practices, in order to obtain factual reasons for choosing one distance over another. For these reasons, in this work we compare many streamline distance functions available in the literature. We focus on the common task of automatic bundle segmentation and we adopt the recent approach of supervised segmentation from expert-based examples. Using the HCP dataset, we compare several distances obtaining guidelines on the choice of which distance function one should use for supervised bundle segmentation

Fuzzy Fibers: Uncertainty in dMRI Tractography

Fiber tracking based on diffusion weighted Magnetic Resonance Imaging (dMRI) allows for noninvasive reconstruction of fiber bundles in the human brain. In this chapter, we discuss sources of error and uncertainty in this technique, and review strategies that afford a more reliable interpretation of the results. This includes methods for computing and rendering probabilistic tractograms, which estimate precision in the face of measurement noise and artifacts. However, we also address aspects that have received less attention so far, such as model selection, partial voluming, and the impact of parameters, both in preprocessing and in fiber tracking itself. We conclude by giving impulses for future research

White matter multi-resolution segmentation using fuzzy set theory

International audienceThe neural architecture of the white matter of the brain, obtained using tractography algorithms, can be divided into different tracts. Their function is, in many cases, still an object of study and might be affected in some syndromes or conditions. Obtaining a reproducible and correct segmentation is therefore crucial both in clinics and in research. However, it is difficult to obtain due to the huge number of fibers and high inter-subject variability. In this paper, we propose to segment and recognize tracts by directly modeling their anatomical definitions, which are usually based on relationships between structures. Since these definitions are mainly qualitative, we propose to model their intrinsic vagueness using fuzzy spatial relations and combine them into a single quantitative score mapped to each fiber. To cope with the high redundancy of tractograms and ease interpretation , we also take advantage of a simplification scheme based on a multi-resolution representation. This allows for an interactive and real-time navigation through different levels of detail. We illustrate our method using the Human Connectome Project dataset and compare it to other well-known white matter segmentation techniques

Alignment of Tractograms As Graph Matching

The white matter pathways of the brain can be reconstructed as 3D polylines, called streamlines, through the analysis of diffusion magnetic resonance imaging (dMRI) data. The whole set of streamlines is called tractogram and represents the structural connectome of the brain. In multiple applications, like group-analysis, segmentation, or atlasing, tractograms of different subjects need to be aligned. Typically, this is done with registration methods, that transform the tractograms in order to increase their similarity. In contrast with transformation-based registration methods, in this work we propose the concept of tractogram correspondence, whose aim is to find which streamline of one tractogram corresponds to which streamline in another tractogram, i.e., a map from one tractogram to another. As a further contribution, we propose to use the relational information of each streamline, i.e., its distances from the other streamlines in its own tractogram, as the building block to define the optimal correspondence. We provide an operational procedure to find the optimal correspondence through a combinatorial optimization problem and we discuss its similarity to the graph matching problem. In this work, we propose to represent tractograms as graphs and we adopt a recent inexact sub-graph matching algorithm to approximate the solution of the tractogram correspondence problem. On tractograms generated from the Human Connectome Project dataset, we report experimental evidence that tractogram correspondence, implemented as graph matching, provides much better alignment than affine registration and comparable if not better results than non-linear registration of volumes

Doctor of Philosophy

dissertationDiffusion magnetic resonance imaging (dMRI) has become a popular technique to detect brain white matter structure. However, imaging noise, imaging artifacts, and modeling techniques, etc., create many uncertainties, which may generate misleading information for further analysis or applications, such as surgical planning. Therefore, how to analyze, effectively visualize, and reduce these uncertainties become very important research questions. In this dissertation, we present both rank-k decomposition and direct decomposition approaches based on spherical deconvolution to decompose the fiber directions more accurately for high angular resolution diffusion imaging (HARDI) data, which will reduce the uncertainties of the fiber directions. By applying volume rendering techniques to an ensemble of 3D orientation distribution function (ODF) glyphs, which we call SIP functions of diffusion shapes, one can elucidate the complex heteroscedastic structural variation in these local diffusion shapes. Furthermore, we quantify the extent of this variation by measuring the fraction of the volume of these shapes, which is consistent across all noise levels, the certain volume ratio. To better understand the uncertainties in white matter fiber tracks, we propose three metrics to quantify the differences between the results of diffusion tensor magnetic resonance imaging (DT-MRI) fiber tracking algorithms: the area between corresponding fibers of each bundle, the Earth Mover's Distance (EMD) between two fiber bundle volumes, and the current distance between two fiber bundle volumes. Based on these metrics, we discuss an interactive fiber track comparison visualization toolkit we have developed to visualize these uncertainties more efficiently. Physical phantoms, with high repeatability and reproducibility, are also designed with the hope of validating the dMRI techniques. In summary, this dissertation provides a better understanding about uncertainties in diffusion magnetic resonance imaging: where and how much are the uncertainties? How do we reduce these uncertainties? How can we possibly validate our algorithms

New tractography methods based on parametric models of white matter fibre dispersion

Diffusion weighted magnetic resonance imaging (DW-MRI) is a powerful imaging technique that can probe the complex structure of the body, revealing structural trends which exist at scales far below the voxel resolution. Tractography utilises the information derived from DW-MRI to examine the structure of white matter. Using information derived from DW-MRI, tractography can estimate connectivity between distinct, functional cortical and sub-cortical regions of grey matter. Understanding how seperate functional regions of the brain are connected as part of a network is key to understanding how the brain works. Tractography has been used to deliniate many known white matter structures and has also revealed structures not fully understood from anatomy due to limitations of histological examination. However, there still remain many shortcomings of tractography, many anatomical features for which tractography algorithms are known to fail, which leads to discrepancies between known anatomy and tractography results. With the aim of approaching a complete picture of the human connectome via tractography, we seek to address the shortcomings in current tractography techniques by exploiting new advances in modelling techniques used in DW-MRI, which provide more accurate representation of underlying white matter anatomy. This thesis introduces a methodology for fully utilising new tissue models in DWMRI to improve tractography. It is known from histology that there are regions of white matter where fibres disperse or curve rapidly at length scales below the DW-MRI voxel resolution. One area where dispersion is particularly prominent is the corona radiata. New DW-MRI models capture dispersion utilising specialised parametric probability distributions. We present novel tractography algorithms utilising these parametric models of dispersion in tractography to improve connectivity estimation in areas of dispersing fibres. We first present an algorithm utilising the the new parametric models of dispersion for tractography in a simple Bayesian framework. We then present an extension to this algorithm which introduces a framework to pool neighbourhood information from multiple voxels in the neighbournhood surrounding the tract in order to better estimate connectivity, introducing the new concept of the neighbourhood-informed orientation distribution function (NI-ODF). Specifically, using neighbourhood exploration we address the ambiguity arising in ’fanning polarity’. In regions of dispersing fibres, the antipodal symmetry inherent in DW-MRI makes it impossible to resolve the polarity of a dispersing fibre configuration from a local voxel-wise model in isolation, by pooling information from neighbouring voxels, we show that this issue can be addressed. We evaluate the newly proposed tractography methods using synthetic phantoms simulating canonical fibre configurations and validate the ability to effectively navigate regions of dispersing fibres and resolve fanning polarity. We then validate that the algorithms perform effectively in real in vivo data, using DW-MRI data from 5 healthy subjects. We show that by utilising models of dispersion, we recover a wider range of connectivity compared to other standard algorithms when tracking through an area of the brain known to have significant white fibre dispersion - the corona radiata. We then examine the impact of the new algorithm on global connectivity estimates in the brain. We find that whole brain connectivity networks derived using the new tractography method feature strong connectivity between frontal lobe regions. This is in contrast to networks derived using competing tractography methods which do not account for sub-voxel fibre dispersion. We also compare thalamo-cortical connectivity estimated using the newly proposed tractography method and compare with a compteing tractography method, finding that the recovered connectivity profiles are largely similar, with some differences in thalamo-cortical connections to regions of the frontal lobe. The results suggest that fibre dispersion is an important structural feature to model in the basis of a tractography algorithm, as it has a strong effect on connectivity estimation