Segmentation of corpus callosum using diffusion tensor imaging: validation in patients with glioblastoma

Abstract Background This paper presents a three-dimensional (3D) method for segmenting corpus callosum in normal subjects and brain cancer patients with glioblastoma. Methods Nineteen patients with histologically confirmed treatment naïve glioblastoma and eleven normal control subjects underwent DTI on a 3T scanner. Based on the information inherent in diffusion tensors, a similarity measure was proposed and used in the proposed algorithm. In this algorithm, diffusion pattern of corpus callosum was used as prior information. Subsequently, corpus callosum was automatically divided into Witelson subdivisions. We simulated the potential rotation of corpus callosum under tumor pressure and studied the reproducibility of the proposed segmentation method in such cases. Results Dice coefficients, estimated to compare automatic and manual segmentation results for Witelson subdivisions, ranged from 94% to 98% for control subjects and from 81% to 95% for tumor patients, illustrating closeness of automatic and manual segmentations. Studying the effect of corpus callosum rotation by different Euler angles showed that although segmentation results were more sensitive to azimuth and elevation than skew, rotations caused by brain tumors do not have major effects on the segmentation results. Conclusions The proposed method and similarity measure segment corpus callosum by propagating a hyper-surface inside the structure (resulting in high sensitivity), without penetrating into neighboring fiber bundles (resulting in high specificity)

Segmentation of diffusion weighted MRI using the level set framework

Medical imaging is a rapidly growing field in which diffusion imaging is a recently developed modality. This novel imaging contrast permits in-vivo measurement of the diffusion of water molecules. This is particularly interesting in brain imaging where the diffusion reveals an amazing insight into the neuronal organization and cerebral cytoarchitecture. Diffusion images contain from six up to hundreds of values in each voxel and are represented as tensor fields (Diffusion Tensor Imaging (DTI)) or as fields of functions (High Angular Resolution Diffusion (HARD) imaging). To fully extract the large amount of data contained within these images new processing and analysis tools are needed. The aim of this thesis is the development of such tools. The method we have been mainly focusing on for this purpose is the level set method. The level set method is a numerical and theoretical tool for propagating interfaces. In image processing it is used for propagating curves in 2D or surfaces in 3D for delineation of objects or for regularization purposes. In this thesis we have explored some of the numerous aspects of the level set frame work to see how the diffusion properties can be used for segmentation. For segmentation of tensor fields we have considered similarity measures for comparison of tensors. From these similarity measures several applications of the level set method have been developed for the segmentation of different structures. Different measures of similarity have been used dependent on the application. When segmenting white matter regions in DTI, the similarity measure emphasizes anisotropic regions. The segmentation algorithm itself has a very local dependence since white matter, in general fiber tracts, experiences different diffusion in different parts of the structure. The presented results show segmentations of the major fiber tracts in the brain. Other structures, such as the deep cerebral nuclei, that are mainly composed of gray matter, have more homogenous diffusion properties than white matter structures. Therefore, in these structures we maximize the internal coherence within the entire structure by using a region based approach to the segmentation problem. Segmentations of the thalamus and its nuclei as well as on tensor fields from fluid mechanics are presented. For High Angular Resolution Diffusion (HARD) images, two methods for fiber tract segmentation are presented based on different types of coherence. The coherence is either measured as the similarity between fibers obtained from a tractography algorithm, or the similarity of scalar values in a five-dimensional non-Euclidean space. The similarity between two fibers is determined by a counting strategy and is equal to the number of voxels they have in common. A spectral clustering algorithm is then used for grouping fibers with a high inter-resemblance. When segmenting white matter with the level set method, we propose to expand the space we are working in from a three-dimensional space of Orientation Distribution Functions (ODF) to a five-dimensional space of position and orientation. By a careful definition of this space and an adaptation of the level set to five dimensions the fibers tracts can be segmented as separated structures. We show some preliminary results from segmentations in this 5D space. The approaches proposed in this thesis permit a consideration of the fiber tracts and gray matter structures as an entity, allowing quantitative measures of the diffusion without losing information by simplifying the images to scalar representations

Improving the Tractography Pipeline: on Evaluation, Segmentation, and Visualization

Recent advances in tractography allow for connectomes to be constructed in vivo. These have applications for example in brain tumor surgery and understanding of brain development and diseases. The large size of the data produced by these methods lead to a variety problems, including how to evaluate tractography outputs, development of faster processing algorithms for tractography and clustering, and the development of advanced visualization methods for verification and exploration. This thesis presents several advances in these fields. First, an evaluation is presented for the robustness to noise of multiple commonly used tractography algorithms. It employs a Monte–Carlo simulation of measurement noise on a constructed ground truth dataset. As a result of this evaluation, evidence for obustness of global tractography is found, and algorithmic sources of uncertainty are identified. The second contribution is a fast clustering algorithm for tractography data based on k–means and vector fields for representing the flow of each cluster. It is demonstrated that this algorithm can handle large tractography datasets due to its linear time and memory complexity, and that it can effectively integrate interrupted fibers that would be rejected as outliers by other algorithms. Furthermore, a visualization for the exploration of structural connectomes is presented. It uses illustrative rendering techniques for efficient presentation of connecting fiber bundles in context in anatomical space. Visual hints are employed to improve the perception of spatial relations. Finally, a visualization method with application to exploration and verification of probabilistic tractography is presented, which improves on the previously presented Fiber Stippling technique. It is demonstrated that the method is able to show multiple overlapping tracts in context, and correctly present crossing fiber configurations

Geodesic tractography segmentation for directional medical image analysis

Acknowledgements page removed per author's request, 01/06/2014.Geodesic Tractography Segmentation is the two component approach presented in this thesis for the analysis of imagery in oriented domains, with emphasis on the application to diffusion-weighted magnetic resonance imagery (DW-MRI). The computeraided analysis of DW-MRI data presents a new set of problems and opportunities for the application of mathematical and computer vision techniques. The goal is to develop a set of tools that enable clinicians to better understand DW-MRI data and ultimately shed new light on biological processes. This thesis presents a few techniques and tools which may be used to automatically find and segment major neural fiber bundles from DW-MRI data. For each technique, we provide a brief overview of the advantages and limitations of our approach relative to other available approaches.Ph.D.Committee Chair: Tannenbaum, Allen; Committee Member: Barnes, Christopher F.; Committee Member: Niethammer, Marc; Committee Member: Shamma, Jeff; Committee Member: Vela, Patrici

Image processing methods for human brain connectivity analysis from in-vivo diffusion MRI

In this PhD Thesis proposal, the principles of diffusion MRI (dMRI) in its application to the human brain mapping of connectivity are reviewed. The background section covers the fundamentals of dMRI, with special focus on those related to the distortions caused by susceptibility inhomogeneity across tissues. Also, a deep survey of available correction methodologies for this common artifact of dMRI is presented. Two methodological approaches to improved correction are introduced. Finally, the PhD proposal describes its objectives, the research plan, and the necessary resources

Extraction of Structural Metrics from Crossing Fiber Models

Diffusion MRI (dMRI) measurements allow us to infer the microstructural properties of white matter and to reconstruct fiber pathways in-vivo. High angular diffusion imaging (HARDI) allows for the creation of more and more complex local models connecting the microstructure to the measured signal. One of the challenges is the derivation of meaningful metrics describing the underlying structure from the local models. The aim hereby is to increase the specificity of the widely used metric fractional anisotropy (FA) by using the additional information contained within the HARDI data. A local model which is connected directly to the underlying microstructure through the model of a single fiber population is spherical deconvolution. It produces a fiber orientation density function (fODF), which can often be interpreted as superposition of multiple peaks, each associated to one relatively coherent fiber population (bundle). Parameterizing these peaks one is able to disentangle and characterize these bundles. In this work, the fODF peaks are approximated by Bingham distributions, capturing first and second order statistics of the fiber orientations, from which metrics for the parametric quantification of fiber bundles are derived. Meaningful relationships between these measures and the underlying microstructural properties are proposed. The focus lies on metrics derived directly from properties of the Bingham distribution, such as peak length, peak direction, peak spread, integral over the peak, as well as a metric derived from the comparison of the largest peaks, which probes the complexity of the underlying microstructure. These metrics are compared to the conventionally used fractional anisotropy (FA) and it is shown how they may help to increase the specificity of the characterization of microstructural properties. Visualization of the micro-structural arrangement is another application of dMRI. This is done by using tractography to propagate the fiber layout, extracted from the local model, in each voxel. In practice most tractography algorithms use little of the additional information gained from HARDI based local models aside from the reconstructed fiber bundle directions. In this work an approach to tractography based on the Bingham parameterization of the fODF is introduced. For each of the fiber populations present in a voxel the diffusion signal and tensor are computed. Then tensor deflection tractography is performed. This allows incorporating the complete bundle information, performing local interpolation as well as using multiple directions per voxel for generating tracts. Another aspect of this work is the investigation of the spherical harmonic representation which is used most commonly for the fODF by means of the parameters derived from the Bingham distribution fit. Here a strong connection between the approximation errors in the spherical representation of the Dirac delta function and the distribution of crossing angles recovered from the fODF was discovered. The final aspect of this work is the application of the metrics derived from the Bingham fit to a number of fetal datasets for quantifying the brain’s development. This is done by introducing the Gini-coefficient as a metric describing the brain’s age

Higher-Order Tensors and Differential Topology in Diffusion MRI Modeling and Visualization

Diffusion Weighted Magnetic Resonance Imaging (DW-MRI) is a noninvasive method for creating three-dimensional scans of the human brain. It originated mostly in the 1970s and started its use in clinical applications in the 1980s. Due to its low risk and relatively high image quality it proved to be an indispensable tool for studying medical conditions as well as for general scientific research. For example, it allows to map fiber bundles, the major neuronal pathways through the brain. But all evaluation of scanned data depends on mathematical signal models that describe the raw signal output and map it to biologically more meaningful values. And here we find the most potential for improvement. In this thesis we first present a new multi-tensor kurtosis signal model for DW-MRI. That means it can detect multiple overlapping fiber bundles and map them to a set of tensors. Compared to other already widely used multi-tensor models, we also add higher order kurtosis terms to each fiber. This gives a more detailed quantification of fibers. These additional values can also be estimated by the Diffusion Kurtosis Imaging (DKI) method, but we show that these values are drastically affected by fiber crossings in DKI, whereas our model handles them as intrinsic properties of fiber bundles. This reduces the effects of fiber crossings and allows a more direct examination of fibers. Next, we take a closer look at spherical deconvolution. It can be seen as a generalization of multi-fiber signal models to a continuous distribution of fiber directions. To this approach we introduce a novel mathematical constraint. We show, that state-of-the-art methods for estimating the fiber distribution become more robust and gain accuracy when enforcing our constraint. Additionally, in the context of our own deconvolution scheme, it is algebraically equivalent to enforcing that the signal can be decomposed into fibers. This means, tractography and other methods that depend on identifying a discrete set of fiber directions greatly benefit from our constraint. Our third major contribution to DW-MRI deals with macroscopic structures of fiber bundle geometry. In recent years the question emerged, whether or not, crossing bundles form two-dimensional surfaces inside the brain. Although not completely obvious, there is a mathematical obstacle coming from differential topology, that prevents general tangential planes spanned by fiber directions at each point to be connected into consistent surfaces. Research into how well this constraint is fulfilled in our brain is hindered by the high precision and complexity needed by previous evaluation methods. This is why we present a drastically simpler method that negates the need for precisely finding fiber directions and instead only depends on the simple diffusion tensor method (DTI). We then use our new method to explore and improve streamsurface visualization.<br /