Directional functions for orientation distribution estimation

Asymmetric fiber trajectory distribution estimation using streamline differential equation

Fiber orientation distribution estimation with diffusion magnetic resonance imaging (dMRI) is critical in white matter fiber tractography which is most commonly based on symmetric crossing structures. Ambiguous spatial correspondence between estimated diffusion directions and fiber geometry, such as asymmetric crossing, bending, fanning, or kissing, makes tractography challenging. Consequently, numerous tracts suggest intertwined connections in unexpected regions of the white matter or actually stop prematurely in the white matter. In this work, we propose a novel asymmetric fiber trajectory distribution (FTD) function defined on neighboring voxels based on a streamline differential equation from fluid mechanics. The spatial consistency with intra- and inter-voxel constraints is derived for FTD estimation by introducing the concept of divergence. At a local level, the FTD is a series of curve flows that minimize the energy function, characterizes the relations between fibers and joint fiber fragments within the same fiber bundle. Experiments are performed on FiberCup phantom, ISMRM 2015 Tractography challenge data, and in vivo brain dMRI data for qualitative and quantitative evaluations. Results show that our approach can reveal continuous asymmetric FTD details that are potentially useful for robust tractography.

White Matter Fiber Tractography Using Nonuniform Rational B-Splines Curve Fitting

The study of neural connectivity has grown rapidly in the past decade. Revealing brain anatomical connection improves not only clinical measures but also cognition understanding. In order to achieve this goal, we have to track neural fiber pathways first. Aiming to estimate 3D fiber pathways more accurately from orientation distribution function (ODF) fields, we presented a novel tracking method based on nonuniform rational B-splines (NURBS) curve fitting. First, we constructed ODF fields from high angular resolution diffusion imaging (HARDI) datasets using diffusion orientation transform (DOT) method. Second, under the angular and length constraints, the consecutive diffusion directions were extracted along each fiber pathway starting from a seed voxel. Finally, after the coordinates of the control points and their corresponding weights were determined, NURBS curve fitting was employed to track fiber pathways. The performance of the proposal has been evaluated on the tractometer phantom and real brain datasets. Based on several measure metrics, the resulting fiber pathways show promising anatomic consistency.

Parcellation of Human Amygdala Subfields Using Orientation Distribution Function and Spectral K-means Clustering

Amygdala plays an important role in fear and emotional learning, which are critical for human survival. Despite the functional relevance and unique circuitry of each human amygdaloid subnuclei, there has yet to be an efficient imaging method for identifying these regions in vivo. A data-driven approach without prior knowledge provides advantages of efficient and objective assessments. The present study uses high angular and high spatial resolution diffusion magnetic resonance imaging to generate orientation distribution function, which bears distinctive microstructural features. The features were extracted using spherical harmonic decomposition to assess microstructural similarity within amygdala subfields are identified via similarity matrices using spectral k-mean clustering. The approach was tested on 32 healthy volunteers and three distinct amygdala subfields were identified including medial, posterior-superior lateral, and anterior-inferior lateral.

Joint Multi-Fiber NODDI Parameter Estimation and Tractography Using the Unscented Information Filter

Tracing white matter fiber bundles is an integral part of analyzing brain connectivity. An accurate estimate of the underlying tissue parameters is also paramount in several neuroscience applications. In this work, we propose to use a joint fiber model estimation and tractography algorithm that uses the NODDI (neurite orientation dispersion diffusion imaging) model to estimate fiber orientation dispersion consistently and smoothly along the fiber tracts along with estimating the intracellular and extracellular volume fractions from the diffusion signal. While the NODDI model has been used in earlier works to estimate the microstructural parameters at each voxel independently, for the first time, we propose to integrate it into a tractography framework. We extend this framework to estimate the NODDI parameters for two crossing fibers, which is imperative to trace fiber bundles through crossings as well as to estimate the microstructural parameters for each fiber bundle separately. We propose to use the unscented information filter (UIF) to accurately estimate the model parameters and perform tractography. The proposed approach has significant computational performance improvements as well as numerical robustness over the unscented Kalman filter (UKF). Our method not only estimates the confidence in the estimated parameters via the covariance matrix, but also provides the Fisher-information matrix of the state variables (model parameters), which can be quite useful to measure model complexity. Results from in-vivo human brain data sets demonstrate the ability of our algorithm to trace through crossing fiber regions, while estimating orientation dispersion and other biophysical model parameters in a consistent manner along the tracts.

Estimating Fiber Orientation Distribution Functions in 3D-Polarized Light Imaging

Research of the human brain connectome requires multiscale approaches derived from independent imaging methods ideally applied to the same object. Hence, comprehensible strategies for data integration across modalities and across scales are essential. We have successfully established a concept to bridge the spatial scales from microscopic fiber orientation measurements based on 3D-Polarized Light Imaging (3D-PLI) to meso- or macroscopic dimensions. By creating orientation distribution functions (pliODFs) from high-resolution vector data via series expansion with spherical harmonics utilizing high performance computing and supercomputing technologies, data fusion with Diffusion Magnetic Resonance Imaging has become feasible, even for a large-scale dataset such as the human brain. Validation of our approach was done effectively by means of two types of datasets that were transferred from fiber orientation maps into pliODFs: simulated 3D-PLI data showing artificial, but clearly defined fiber patterns and real 3D-PLI data derived from sections through the human brain and the brain of a hooded seal.

A localized Richardson-Lucy algorithm for fiber orientation estimation in high angular resolution diffusion imaging

PURPOSE

Diffusion tensor imaging is widely used for studying neural fiber trajectories in white matter and for quantifying changes in tissue using diffusion properties at each voxel in the brain. To better model the nature of crossing fibers within complex architectures, rather than using a simplified tensor model that assumes only a single fiber direction at each image voxel, a model mixing multiple diffusion tensors is used to profile diffusion signals from high angular resolution diffusion imaging (HARDI) data. Based on the HARDI signal and a multiple tensors model, spherical deconvolution methods have been developed to overcome the limitations of the diffusion tensor model when resolving crossing fibers. The Richardson-Lucy algorithm is a popular spherical deconvolution method used in previous work. However, it is based on a Gaussian distribution, while HARDI data are always very noisy, and the distribution of HARDI data follows a Rician distribution. This current work aims to present a novel solution to address these issues.

METHODS

By simultaneously considering both the Rician bias and neighbor correlation in HARDI data, the authors propose a localized Richardson-Lucy (LRL) algorithm to estimate fiber orientations for HARDI data. The proposed method can simultaneously reduce noise and correct the Rician bias.

RESULTS

Mean angular error (MAE) between the estimated Fiber orientation distribution (FOD) field and the reference FOD field was computed to examine whether the proposed LRL algorithm offered any advantage over the conventional RL algorithm at various levels of noise. Normalized mean squared error (NMSE) was also computed to measure the similarity between the true FOD field and the estimated FOD filed. For MAE comparisons, the proposed LRL approach obtained the best results in most of the cases at different levels of SNR and b-values. For NMSE comparisons, the proposed LRL approach obtained the best results in most of the cases at b-value = 3000 s/mm(2), which is the recommended schema for HARDI data acquisition. In addition, the FOD fields estimated by the proposed LRL approach in regions of fiber crossing regions using real data sets also showed similar fiber structures which agreed with common acknowledge in these regions.

CONCLUSIONS

The novel spherical deconvolution method for improved accuracy in investigating crossing fibers can simultaneously reduce noise and correct Rician bias. With the noise smoothed and bias corrected, this algorithm is especially suitable for estimation of fiber orientations in HARDI data. Experimental results using both synthetic and real imaging data demonstrated the success and effectiveness of the proposed LRL algorithm.

Non-parametric representation and prediction of single- and multi-shell diffusion-weighted MRI data using Gaussian processes

Diffusion MRI offers great potential in studying the human brain microstructure and connectivity. However, diffusion images are marred by technical problems, such as image distortions and spurious signal loss. Correcting for these problems is non-trivial and relies on having a mechanism that predicts what to expect. In this paper we describe a novel way to represent and make predictions about diffusion MRI data. It is based on a Gaussian process on one or several spheres similar to the Geostatistical method of "Kriging". We present a choice of covariance function that allows us to accurately predict the signal even from voxels with complex fibre patterns. For multi-shell data (multiple non-zero b-values) the covariance function extends across the shells which means that data from one shell is used when making predictions for another shell.

Strengths and weaknesses of state of the art fiber tractography pipelines--A comprehensive in-vivo and phantom evaluation study using Tractometer

Many different tractography approaches and corresponding isolated evaluation attempts have been presented over the last years, but a comparative and quantitative evaluation of tractography algorithms still remains a challenge, particularly in-vivo. The recently presented evaluation framework Tractometer is the first attempt to approach this challenge in a quantitative, comparative, persistent and open-access way. Tractometer is currently based on the evaluation of several global connectivity and tract-overlap metrics on hardware phantom data. The work presented in this paper focuses on extending Tractometer with a metric that enables the assessment of the local consistency of tractograms with the underlying image data that is not only applicable to phantom dataset but allows the quantitative and purely data-driven evaluation of in-vivo tractography. We furthermore present an extensive reference-based evaluation study of 25,000 tractograms obtained on phantom and in-vivo datasets using the presented local metric as well as all the methods already established in Tractometer. The experiments showed that the presented local metric successfully reflects the behavior of in-vivo tractography under different conditions and that it is consistent with the results of previous studies. Additionally our experiments enabled a multitude of conclusions with implications for fiber tractography in general, including recommendations regarding optimal choice of a local modeling technique, tractography algorithm, and parameterization, confirming and complementing the results of earlier studies.

Sparse deconvolution of higher order tensor for fiber orientation distribution estimation

PURPOSE

Higher order tensor (HOT) imaging approaches based on the spherical deconvolution framework have attracted much interest for their effectiveness in estimating fiber orientation distribution (FOD). However, sparse regularization techniques are still needed to obtain stable FOD in solving the deconvolution problem, particularly in very high orders. Our goal is to adequately characterize the actual sparsity lying in the FOD domain to develop accurate estimation approach for fiber orientation in HOT framework.

MATERIALS AND METHODS

We propose a sparse HOT regularization model by enforcing the sparse constraint directly on the representation of FOD instead of imposing it on coefficients of basis function. Then, we incorporate both the stabilizing effect of the l2 penalty and the sparsity encouraging effect of the l1 penalty in the sparse model to adequately characterize the actual sparsity lying in the FOD domain. Furthermore, a weighted regularization scheme is developed to iteratively solve the deconvolution problem. The deconvolution technique is compared against existing methods using l2 or l1 regularizer and tested on synthetic data and real human brain.

RESULTS

Experiments were conducted on synthetic data and real human brain data. The synthetic experimental results indicate that crossing fibers are more easily detected and the angular resolution limit is improved by our method by approximately 20°-30° compared to existing HOT method. The detection accuracy is considerably improved compared with that of spherical deconvolution approaches using the l2 regularizer and the reweighted l1 scheme.

CONCLUSIONS

Results of testing the deconvolution technique demonstrate that it allows HOTs to obtain increasingly clean and sharp FOD, which in turn significantly increases the angular resolution of current HOT methods. With sparsity on FOD domain, this method efficiently improves the ability of HOT in resolving crossing fibers.

Estimating constrained multi-fiber diffusion MR volumes by orientation clustering

Diffusion MRI is a valuable tool for mapping tissue microstructure; however, multi-fiber models present challenges to image analysis operations. In this paper, we present a method for estimating models for such operations by clustering fiber orientations. Our approach is applied to ball-and-stick diffusion models, which include an isotropic tensor and multiple sticks encoding fiber volume and orientation. We consider operations which can be generalized to a weighted combination of fibers and present a method for representing such combinations with a mixture-of-Watsons model, learning its parameters by Expectation Maximization. We evaluate this approach with two experiments. First, we show it is effective for filtering in the presence of synthetic noise. Second, we demonstrate interpolation and averaging by construction of a tractography atlas, showing improved reconstruction of white matter pathways. These experiments indicate that our method is useful in estimating multi-fiber ball-and-stick diffusion volumes resulting from a range of image analysis operations.

Linear transforms for Fourier data on the sphere: application to high angular resolution diffusion MRI of the brain

This paper presents a novel family of linear transforms that can be applied to data collected from the surface of a 2-sphere in three-dimensional Fourier space. This family of transforms generalizes the previously-proposed Funk-Radon Transform (FRT), which was originally developed for estimating the orientations of white matter fibers in the central nervous system from diffusion magnetic resonance imaging data. The new family of transforms is characterized theoretically, and efficient numerical implementations of the transforms are presented for the case when the measured data is represented in a basis of spherical harmonics. After these general discussions, attention is focused on a particular new transform from this family that we name the Funk-Radon and Cosine Transform (FRACT). Based on theoretical arguments, it is expected that FRACT-based analysis should yield significantly better orientation information (e.g., improved accuracy and higher angular resolution) than FRT-based analysis, while maintaining the strong characterizability and computational efficiency of the FRT. Simulations are used to confirm these theoretical characteristics, and the practical significance of the proposed approach is illustrated with real diffusion weighted MRI brain data. These experiments demonstrate that, in addition to having strong theoretical characteristics, the proposed approach can outperform existing state-of-the-art orientation estimation methods with respect to measures such as angular resolution and robustness to noise and modeling errors.

Computational Representation of White Matter Fiber Orientations

We present a new methodology based on directional data clustering to represent white matter fiber orientations in magnetic resonance analyses for high angular resolution diffusion imaging. A probabilistic methodology is proposed for estimating intravoxel principal fiber directions, based on clustering directional data arising from orientation distribution function (ODF) profiles. ODF reconstructions are used to estimate intravoxel fiber directions using mixtures of von Mises-Fisher distributions. The method focuses on clustering data on the unit sphere, where complexity arises from representing ODF profiles as directional data. The proposed method is validated on synthetic simulations, as well as on a real data experiment. Based on experiments, we show that by clustering profile data using mixtures of von Mises-Fisher distributions it is possible to estimate multiple fiber configurations in a more robust manner than currently used approaches, without recourse to regularization or sharpening procedures. The method holds promise to support robust tractographic methodologies and to build realistic models of white matter tracts in the human brain.

Spatial warping of DWI data using sparse representation

Registration of DWI data, unlike scalar image data, is complicated by the need of reorientation algorithms for keeping the orientation architecture of each voxel aligned with the rest of the image. This paper presents an algorithm for effective and efficient warping and reconstruction of diffusion-weighted imaging (DWI) signals for the purpose of spatial transformation. The key idea is to decompose the DWI signal profile, a function defined on a unit sphere, into a series of weighted fiber basis functions (FBFs), reorient these FBFs independently based on the local affine transformation, and then recompose the reoriented FBFs to obtain the final transformed DWI signal profile. We enforce a sparsity constraint on the weights of the FBFs during the decomposition to reflect the fact that the DWI signal profile typically gains its information from a limited number of fiber populations. A non-negative constraint is further imposed so that noise-induced negative lobes in the profile can be avoided. The proposed framework also explicitly models the isotropic component of the diffusion signals to avoid undesirable reorientation artifacts in signal reconstruction. In contrast to existing methods, the current algorithm is executed directly in the DWI signal space, thus allowing any diffusion models to be fitted to the data after transformation.

A nonparametric Riemannian framework for processing high angular resolution diffusion images and its applications to ODF-based morphometry

High angular resolution diffusion imaging (HARDI) has become an important technique for imaging complex oriented structures in the brain and other anatomical tissues. This has motivated the recent development of several methods for computing the orientation probability density function (PDF) at each voxel. However, much less work has been done on developing techniques for filtering, interpolation, averaging and principal geodesic analysis of orientation PDF fields. In this paper, we present a Riemannian framework for performing such operations. The proposed framework does not require that the orientation PDFs be represented by any fixed parameterization, such as a mixture of von Mises-Fisher distributions or a spherical harmonic expansion. Instead, we use a nonparametric representation of the orientation PDF. We exploit the fact that under the square-root re-parameterization, the space of orientation PDFs forms a Riemannian manifold: the positive orthant of the unit Hilbert sphere. We show that various orientation PDF processing operations, such as filtering, interpolation, averaging and principal geodesic analysis, may be posed as optimization problems on the Hilbert sphere, and can be solved using Riemannian gradient descent. We illustrate these concepts with numerous experiments on synthetic, phantom and real datasets. We show their application to studying left/right brain asymmetries.

Spatially regularized compressed sensing for high angular resolution diffusion imaging

Despite the relative recency of its inception, the theory of compressive sampling (aka compressed sensing) (CS) has already revolutionized multiple areas of applied sciences, a particularly important instance of which is medical imaging. Specifically, the theory has provided a different perspective on the important problem of optimal sampling in magnetic resonance imaging (MRI), with an ever-increasing body of works reporting stable and accurate reconstruction of MRI scans from the number of spectral measurements which would have been deemed unacceptably small as recently as five years ago. In this paper, the theory of CS is employed to palliate the problem of long acquisition times, which is known to be a major impediment to the clinical application of high angular resolution diffusion imaging (HARDI). Specifically, we demonstrate that a substantial reduction in data acquisition times is possible through minimization of the number of diffusion encoding gradients required for reliable reconstruction of HARDI scans. The success of such a minimization is primarily due to the availability of spherical ridgelet transformation, which excels in sparsifying HARDI signals. What makes the resulting reconstruction procedure even more accurate is a combination of the sparsity constraints in the diffusion domain with additional constraints imposed on the estimated diffusion field in the spatial domain. Accordingly, the present paper describes an original way to combine the diffusion- and spatial-domain constraints to achieve a maximal reduction in the number of diffusion measurements, while sacrificing little in terms of reconstruction accuracy. Finally, details are provided on an efficient numerical scheme which can be used to solve the aforementioned reconstruction problem by means of standard and readily available estimation tools. The paper is concluded with experimental results which support the practical value of the proposed reconstruction methodology.

Recent advances in diffusion MRI modeling: Angular and radial reconstruction

Recent advances in diffusion magnetic resonance image (dMRI) modeling have led to the development of several state of the art methods for reconstructing the diffusion signal. These methods allow for distinct features to be computed, which in turn reflect properties of fibrous tissue in the brain and in other organs. A practical consideration is that to choose among these approaches requires very specialized knowledge. In order to bridge the gap between theory and practice in dMRI reconstruction and analysis we present a detailed review of the dMRI modeling literature. We place an emphasis on the mathematical and algorithmic underpinnings of the subject, categorizing existing methods according to how they treat the angular and radial sampling of the diffusion signal. We describe the features that can be computed with each method and discuss its advantages and limitations. We also provide a detailed bibliography to guide the reader.

Reconstruction of the orientation distribution function in single- and multiple-shell q-ball imaging within constant solid angle

q-Ball imaging is a high-angular-resolution diffusion imaging technique that has been proven very successful in resolving multiple intravoxel fiber orientations in MR images. The standard computation of the orientation distribution function (the probability of diffusion in a given direction) from q-ball data uses linear radial projection, neglecting the change in the volume element along each direction. This results in spherical distributions that are different from the true orientation distribution functions. For instance, they are neither normalized nor as sharp as expected and generally require postprocessing, such as artificial sharpening. In this paper, a new technique is proposed that, by considering the solid angle factor, uses the mathematically correct definition of the orientation distribution function and results in a dimensionless and normalized orientation distribution function expression. Our model is flexible enough so that orientation distribution functions can be estimated either from single q-shell datasets or by exploiting the greater information available from multiple q-shell acquisitions. We show that the latter can be achieved by using a more accurate multiexponential model for the diffusion signal. The improved performance of the proposed method is demonstrated on artificial examples and high-angular-resolution diffusion imaging data acquired on a 7-T magnet.

A full bi-tensor neural tractography algorithm using the unscented Kalman filter

We describe a technique that uses tractography to visualize neural pathways in human brains by extending an existing framework that uses overlapping Gaussian tensors to model the signal. At each point on the fiber, an unscented Kalman filter is used to find the most consistent direction as a mixture of previous estimates and of the local model. In our previous framework, the diffusion ellipsoid had a cylindrical shape, i.e., the diffusion tensor's second and third eigenvalues were identical. In this paper, we extend the tensor representation so that the diffusion tensor is represented by an arbitrary ellipsoid. Experiments on synthetic data show a reduction in the angular error at fiber crossings and branchings. Tests on in vivo data demonstrate the ability to trace fibers in areas containing crossings or branchings, and the tests also confirm the superiority of using a full tensor representation over the simplified model.

Filtered multitensor tractography

We describe a technique that uses tractography to drive the local fiber model estimation. Existing techniques use independent estimation at each voxel so there is no running knowledge of confidence in the estimated model fit. We formulate fiber tracking as recursive estimation: at each step of tracing the fiber, the current estimate is guided by those previous. To do this we perform tractography within a filter framework and use a discrete mixture of Gaussian tensors to model the signal. Starting from a seed point, each fiber is traced to its termination using an unscented Kalman filter to simultaneously fit the local model to the signal and propagate in the most consistent direction. Despite the presence of noise and uncertainty, this provides a causal estimate of the local structure at each point along the fiber. Using two- and three-fiber models we demonstrate in synthetic experiments that this approach significantly improves the angular resolution at crossings and branchings. In vivo experiments confirm the ability to trace through regions known to contain such crossing and branching while providing inherent path regularization.

Tensor kernels for simultaneous fiber model estimation and tractography

This paper proposes a novel framework for joint orientation distribution function estimation and tractography based on a new class of tensor kernels. Existing techniques estimate the local fiber orientation at each voxel independently so there is no running knowledge of confidence in the measured signal or estimated fiber orientation. In this work, fiber tracking is formulated as recursive estimation: at each step of tracing the fiber, the current estimate of the orientation distribution function is guided by the previous. To do this, second-and higher-order tensor-based kernels are employed. A weighted mixture of these tensor kernels is used for representing crossing and branching fiber structures. While tracing a fiber, the parameters of the mixture model are estimated based on the orientation distribution function at that location and a smoothness term that penalizes deviation from the previous estimate along the fiber direction. This ensures smooth estimation along the direction of propagation of the fiber. In synthetic experiments, using a mixture of two and three components it is shown that this approach improves the angular resolution at crossings. In vivo experiments using two and three components examine the corpus callosum and corticospinal tract and confirm the ability to trace through regions known to contain such crossing and branching.

A Riemannian framework for orientation distribution function computing

Compared with Diffusion Tensor Imaging (DTI), High Angular Resolution Imaging (HARDI) can better explore the complex microstructure of white matter. Orientation Distribution Function (ODF) is used to describe the probability of the fiber direction. Fisher information metric has been constructed for probability density family in Information Geometry theory and it has been successfully applied for tensor computing in DTI. In this paper, we present a state of the art Riemannian framework for ODF computing based on Information Geometry and sparse representation of orthonormal bases. In this Riemannian framework, the exponential map, logarithmic map and geodesic have closed forms. And the weighted Frechet mean exists uniquely on this manifold. We also propose a novel scalar measurement, named Geometric Anisotropy (GA), which is the Riemannian geodesic distance between the ODF and the isotropic ODF. The Renyi entropy H1/2 of the ODF can be computed from the GA. Moreover, we present an Affine-Euclidean framework and a Log-Euclidean framework so that we can work in an Euclidean space. As an application, Lagrange interpolation on ODF field is proposed based on weighted Frechet mean. We validate our methods on synthetic and real data experiments. Compared with existing Riemannian frameworks on ODF, our framework is model-free. The estimation of the parameters, i.e. Riemannian coordinates, is robust and linear. Moreover it should be noted that our theoretical results can be used for any probability density function (PDF) under an orthonormal basis representation.

On approximation of orientation distributions by means of spherical ridgelets

Visualization and analysis of the micro-architecture of brain parenchyma by means of magnetic resonance imaging is nowadays believed to be one of the most powerful tools used for the assessment of various cerebral conditions as well as for understanding the intracerebral connectivity. Unfortunately, the conventional diffusion tensor imaging (DTI) used for estimating the local orientations of neural fibers is incapable of performing reliably in the situations when a voxel of interest accommodates multiple fiber tracts. In this case, a much more accurate analysis is possible using the high angular resolution diffusion imaging (HARDI) that represents local diffusion by its apparent coefficients measured as a discrete function of spatial orientations. In this note, a novel approach to enhancing and modeling the HARDI signals using multiresolution bases of spherical ridgelets is presented. In addition to its desirable properties of being adaptive, sparsifying, and efficiently computable, the proposed modeling leads to analytical computation of the orientation distribution functions associated with the measured diffusion, thereby providing a fast and robust analytical solution for q-ball imaging.

A filtered approach to neural tractography using the Watson directional function

We propose a technique to simultaneously estimate the local fiber orientations and perform multi-fiber tractography. Existing techniques estimate the local fiber orientation at each voxel independently so there is no running knowledge of confidence in the measured signal or estimated fiber orientation. Further, to overcome noise, many algorithms use a filter as a post-processing step to obtain a smooth trajectory. We formulate fiber tracking as causal estimation: at each step of tracing the fiber, the current estimate of the signal is guided by the previous. To do this, we model the signal as a discrete mixture of Watson directional functions and perform tractography within a filtering framework. Starting from a seed point, each fiber is traced to its termination using an unscented Kalman filter to simultaneously fit the signal and propagate in the most consistent direction. Despite the presence of noise and uncertainty, this provides an accurate estimate of the local structure at each point along the fiber. We choose the Watson function since it provides a compact representation of the signal parameterized by the principal diffusion direction and a scaling parameter describing anisotropy, and also allows analytic reconstruction of the oriented diffusion function from those parameters. Using a mixture of two and three components (corresponding to two-fiber and three-fiber models) we demonstrate in synthetic experiments that this approach reduces signal reconstruction error and significantly improves the angular resolution at crossings and branchings. In vivo experiments examine the corpus callosum and internal capsule and confirm the ability to trace through regions known to contain such crossing and branching while providing inherent path regularization.

Neural tractography using an unscented Kalman filter

We describe a technique to simultaneously estimate a local neural fiber model and trace out its path. Existing techniques estimate the local fiber orientation at each voxel independently so there is no running knowledge of confidence in the estimated fiber model. We formulate fiber tracking as recursive estimation: at each step of tracing the fiber, the current estimate is guided by the previous. To do this we model the signal as a mixture of Gaussian tensors and perform tractography within a filter framework. Starting from a seed point, each fiber is traced to its termination using an unscented Kalman filter to simultaneously fit the local model and propagate in the most consistent direction. Despite the presence of noise and uncertainty, this provides a causal estimate of the local structure at each point along the fiber. Synthetic experiments demonstrate that this approach reduces signal reconstruction error and significantly improves the angular resolution at crossings and branchings. In vivo experiments confirm the ability to trace out fibers in areas known to contain such crossing and branching while providing inherent path regularization.

Adaptive kernels for multi-fiber reconstruction

In this paper we present a novel method for multi-fiber reconstruction given a diffusion-weighted MRI dataset. There are several existing methods that employ various spherical deconvolution kernels for achieving this task. However the kernels in all of the existing methods rely on certain assumptions regarding the properties of the underlying fibers, which introduce inaccuracies and unnatural limitations in them. Our model is a non trivial generalization of the spherical deconvolution model, which unlike the existing methods does not make use of a fix-shaped kernel. Instead, the shape of the kernel is estimated simultaneously with the rest of the unknown parameters by employing a general adaptive model that can theoretically approximate any spherical deconvolution kernel. The performance of our model is demonstrated using simulated and real diffusion-weighed MR datasets and compared quantitatively with several existing techniques in literature. The results obtained indicate that our model has superior performance that is close to the theoretic limit of the best possible achievable result.