Publications by Year: 2006

The hippocampus is a critical structure for learning and memory formation injured by diverse neuropathologies such as epilepsy or Alzheimer’s disease. Recently, clinical investigations have attempted to use diffusion tensor MRI as a more specific surrogate marker for hippocampal damage. To first better understand the tissue architecture of healthy hippocampal regions, this study characterized 10 rat hippocampi with diffusion tensor imaging (DTI) at 50-microm in-plane image resolution using a 14.1-T magnet. Chemical fixation of the dissected and straightened rat hippocampus provided a simple, effective way to reduce partial volume effects when segmenting hippocampal regions and improved mean signal-to-noise per unit time (e.g. 50.6+/-4.4 at b=1250 s/mm2 in 27 min). Contrary to previous reports that water diffusion is homogeneous throughout the nervous system, statistically different mean diffusivities were observed (e.g. 0.238+/-0.054 and 0.318+/-0.084 microm2/ms for the molecular and granule cell layers respectively) (ANOVA, P

We present a two-step process including white matter atlas generation and automatic segmentation. Our atlas generation method is based on population fiber clustering. We produce an atlas which contains high-dimensional descriptors of fiber bundles as well as anatomical label information. We use the atlas to automatically segment tractography in the white matter of novel subjects and we present quantitative results (FA measurements) in segmented white matter regions from a small population. We demonstrate reproducibility of these measurements across scans. In addition, we introduce the idea of using clustering for automatic matching of anatomical structures across hemispheres.

This article describes an accurate and fast method for fiber orientation mapping using multidirectional diffusion-weighted magnetic resonance (MR) data. This novel approach utilizes the Fourier transform relationship between the water displacement probabilities and diffusion-attenuated MR signal expressed in spherical coordinates. The radial part of the Fourier integral is evaluated analytically under the assumption that MR signal attenuates exponentially. The values of the resulting functions are evaluated at a fixed distance away from the origin. The spherical harmonic transform of these functions yields the Laplace series coefficients of the probabilities on a sphere of fixed radius. Alternatively, probability values can be computed nonparametrically using Legendre polynomials. Orientation maps calculated from excised rat nervous tissue data demonstrate this technique’s ability to accurately resolve crossing fibers in anatomical regions such as the optic chiasm. This proposed methodology has a trivial extension to multiexponential diffusion-weighted signal decay. The developed methods will improve the reliability of tractography schemes and may make it possible to correctly identify the neural connections between functionally connected regions of the nervous system.

Observation of translational self-diffusion of water molecules using magnetic resonance (MR) techniques has proven to be a powerful means to probe tissue microstructure. The collected MR signal depends on experimentally controllable parameters as well as the descriptors of tissue geometry. In order to obtain the latter, one needs to employ accurate models to characterize the dependence of the signal on the varied experimental parameters. In this work, a simple model describing diffusion in disordered media and fractal spaces is shown to describe the diffusion-time dependence of the diffusion attenuated MR signal obtained from biological specimens successfully. The model enables one to quantify the evolution of the average water displacement probabilities in terms of two exponents—dw and ds. The experiments performed on excised human neural tissue samples and human red blood cell ghosts indicate that these two parameters are sensitive to tissue microstructure. Therefore, it may be possible to use the proposed scheme to generate novel contrast mechanism for classifying and segmenting tissue.

Segmentation involves separating an object from the background. In this work, we propose a novel segmentation method combining image information with prior shape knowledge, within the level-set framework. Following the work of Leventon et al., we revisit the use of principal component analysis (PCA) to introduce prior knowledge about shapes in a more robust manner. To this end, we utilize Kernel PCA and show that this method of learning shapes outperforms linear PCA, by allowing only shapes that are close enough to the training data. In the proposed segmentation algorithm, shape knowledge and image information are encoded into two energy functionals entirely described in terms of shapes. This consistent description allows to fully take advantage of the Kernel PCA methodology and leads to promising segmentation results. In particular, our shape-driven segmentation technique allows for the simultaneous encoding of multiple types of shapes, and offers a convincing level of robustness with respect to noise, clutter, partial occlusions, or smearing.

A classical neural tract tracer, WGA-HRP, was injected at multiple sites within the brain of a macaque monkey. Histological sections of the labeled fiber tracts were reconstructed in 3D, and the fibers were segmented and registered with the anatomical post-mortem MRI from the same animal. Fiber tracing along the same pathways was performed on the DTI data using a classical diffusion tracing technique. The fibers derived from the DTI were compared with those segmented from the histology in order to evaluate the performance of DTI fiber tracing. While there was generally good agreement between the two methods, our results reveal certain limitations of DTI tractography, particularly at regions of fiber tract crossing or bifurcation.

In vivo white matter tractography by diffusion tensor imaging (DTI) has become a popular tool for investigation of white matter architecture in the normal brain. Despite some unresolved issues regarding the accuracy of DTI, recent studies applied DTI for delineating white matter organization in the vicinity of brain lesions and especially brain tumors. Apart from the intrinsic limitations of DTI, the tracking of fibers in the vicinity or within lesions is further complicated due to changes in diseased tissue such as elevated water content (edema), tissue compression and degeneration. These changes deform the architecture of the white matter and in some cases prevent definite selection of the seed region of interest (ROI) from which fiber tracking begins. We show here that for displaced fiber systems, the use of anatomical approach for seed ROI selection yields insufficient results. Alternatively, we propose to select the seed points based on functional MRI activations which constrain the subjective seed ROI selection. The results are demonstrated on two major fiber systems: the pyramidal tract and the superior longitudinal fasciculus that connect critical motor and language areas, respectively. The fMRI based seed ROI selection approach enabled a more comprehensive mapping of these fiber systems. Furthermore, this procedure enabled the characterization of displaced white matter using the eigenvalue decomposition of DTI. We show that along the compressed fiber system, the diffusivity parallel to the fiber increases, while that perpendicular to the fibers decreases, leading to an overall increase in the fractional anisotropy index reflecting the compression of the fiber bundle. We conclude that definition of the functional network of a subject with deformed white matter should be done carefully. With fMRI, one can more accurately define the seed ROI for DTI based tractography and to provide a more comprehensive, functionally related, white matter mapping, a very important tool used in pre-surgical mapping.

MR diffusion tensor imaging (DTI) of the brain and spine provides a unique tool for both visualizing directionality and assessing intactness of white matter fiber tracts in vivo. At the spatial resolution of clinical MRI, much of primate white matter is composed of interdigitating fibers. Analyses based on an assumed single diffusion tensor per voxel yield important information about the average diffusion in the voxel but fail to reveal structure in the presence of crossing tracts. Until today, all clinical scans assume only one tensor, causing potential serious errors in tractography. Since high angular resolution imaging remains, so far, untenable for routine clinical use, a method is proposed whereby the single-tensor field is augmented with additional information gleaned from standard clinical DTI. The method effectively resolves two distinct tract directions within voxels, in which only two tracts are assumed to exist. The underlying constrained two-tensor model is fitted in two stages, utilizing the information present in the single-tensor fit. As a result, the necessary MRI time can be drastically reduced when compared with other approaches, enabling widespread clinical use. Upon evaluation in simulations and application to in vivo human brain DTI data, the method appears to be robust and practical and, if correctly applied, could elucidate tract directions at critical points of uncertainty.

A new and complementary method to assess image quality is presented. It is based on the comparison of the local variance distribution of two images. This new quality index is better suited to assess the non-stationarity of images, therefore it explicitly focuses on the image structure. We show that this new index outperforms other methods for the assessment of image quality in medical images.

This paper describes a novel user-guided method for grouping fibers from diffusion tensor MRI tractography into bundles. The method finds fibers, that passing through user-defined ROIs, still fit to the underlying data model given by the diffusion tensor. This is achieved by filtering the data and the ROIs with a kernel derived from a geodesic metric between tensors. A standard approach using binary decisions defining tracts passing through ROIs is critically dependent on ROIs that includes all trace lines of interest. The method described in this paper uses a softer decision mechanism through a kernel which enables grouping of bundles driven less exact, or even single point, ROIs. The method analyzes the responses obtained from the convolution with a kernel function along the fiber with the ROI data. Results in real data shows the feasibility of the approach to fiber bundling.