GROUP ACTION INDUCED AVERAGING FOR HARDI PROCESSING

Flow-based Geometric Interpolation of Fiber Orientation Distribution Functions

The fiber orientation distribution function (FOD) is an advanced model for high angular resolution diffusion MRI representing complex fiber geometry. However, the complicated mathematical structures of the FOD function pose challenges for FOD image processing tasks such as interpolation, which plays a critical role in the propagation of fiber tracts in tractography. In FOD-based tractography, linear interpolation is commonly used for numerical efficiency, but it is prone to generate false artificial information, leading to anatomically incorrect fiber tracts. To overcome this difficulty, we propose a flowbased and geometrically consistent interpolation framework that considers peak-wise rotations of FODs within the neighborhood of each location. Our method decomposes a FOD function into multiple components and uses a smooth vector field to model the flows of each peak in its neighborhood. To generate the interpolated result along the flow of each vector field, we develop a closed-form and efficient method to rotate FOD peaks in neighboring voxels and realize geometrically consistent interpolation of FOD components. By combining the interpolation results from each peak, we obtain the final interpolation of FODs. Experimental results on Human Connectome Project (HCP) data demonstrate that our method produces anatomically more meaningful FOD interpolations and significantly enhances tractography performance.

Advances in computational and statistical diffusion MRI

Computational methods are crucial for the analysis of diffusion magnetic resonance imaging (MRI) of the brain. Computational diffusion MRI can provide rich information at many size scales, including local microstructure measures such as diffusion anisotropies or apparent axon diameters, whole-brain connectivity information that describes the brain's wiring diagram and population-based studies in health and disease. Many of the diffusion MRI analyses performed today were not possible five, ten or twenty years ago, due to the requirements for large amounts of computer memory or processor time. In addition, mathematical frameworks had to be developed or adapted from other fields to create new ways to analyze diffusion MRI data. The purpose of this review is to highlight recent computational and statistical advances in diffusion MRI and to put these advances into context by comparison with the more traditional computational methods that are in popular clinical and scientific use. We aim to provide a high-level overview of interest to diffusion MRI researchers, with a more in-depth treatment to illustrate selected computational advances.

Interpolation on the manifold of K component GMMs

Probability density functions (PDFs) are fundamental objects in mathematics with numerous applications in computer vision, machine learning and medical imaging. The feasibility of basic operations such as computing the distance between two PDFs and estimating a mean of a set of PDFs is a direct function of the representation we choose to work with. In this paper, we study the Gaussian mixture model (GMM) representation of the PDFs motivated by its numerous attractive features. (1) GMMs are arguably more interpretable than, say, square root parameterizations (2) the model complexity can be explicitly controlled by the number of components and (3) they are already widely used in many applications. The main contributions of this paper are numerical algorithms to enable basic operations on such objects that strictly respect their underlying geometry. For instance, when operating with a set of K component GMMs, a first order expectation is that the result of simple operations like interpolation and averaging should provide an object that is also a K component GMM. The literature provides very little guidance on enforcing such requirements systematically. It turns out that these tasks are important internal modules for analysis and processing of a field of ensemble average propagators (EAPs), common in diffusion weighted magnetic resonance imaging. We provide proof of principle experiments showing how the proposed algorithms for interpolation can facilitate statistical analysis of such data, essential to many neuroimaging studies. Separately, we also derive interesting connections of our algorithm with functional spaces of Gaussians, that may be of independent interest.

Nonsmooth ICA contrast minimization using a Riemannian Nelder-Mead method

This brief concerns the design and application of a Riemannian Nelder-Mead algorithm to minimize a Hartley-entropybased contrast function to reliably estimate the sources from their mixtures. Despite its nondifferentiability, the contrast function is endowed with attractive properties such as discriminacy, and hence warrants an effort to be effectively handled by a derivative-free optimizer. Aside from tailoring the Nelder-Mead technique to the constraint set, namely, oblique manifold, the source separation results attained in an empirical study with quasi-correlated synthetic signals and digital images are presented, which favor the proposed method on a comparative basis.

Diffusion of fiber orientation distribution functions with a rotation-induced riemannian metric

Advanced diffusion weighted MR imaging allows non-invasive study on the structural connectivity of human brains. Fiber orientation distributions (FODs) reconstructed from diffusion data are a popular model to represent crossing fibers. For this sophisticated image representation of connectivity, classical image operations such as smoothing must be redefined. In this paper, we propose a novel rotation-induced Riemannian metric for FODs, and introduce a weighted diffusion process for FODs regarding this Riemannian manifold. We show how this Riemannian manifold can be used for smoothing, interpolation and building image-pyramids, yielding more accurate or intuitively more reasonable results than the linear or the unit hyper-sphere manifold.

ROTATIONAL GRADIENT FIELD FOR INTERPOLATION OF FIBER ORIENTATION DISTRIBUTION IN CONNECTIVITY ANALYSIS

Advances in diffusion weighted MR imaging have made it possible to non-invasively study the structral connectivity of human brains at high resolution. To model crossing fibers in white matter, a popular choice is the reconstruction of fiber orientation distributions (FODs) from diffusion data. For this sophiscated image representation of brain connectivity, classical image operations such as differentiation and interpolation must be redefined. In this paper, we introduce rotational gradient fields (RGF) as the spatial differential of FODs' orientations. By taking into accout the rotational effect of traveling in the RGF, an FOD at one location can be transported and aligned with the FOD at the target location. We propose a method for inducing RGFs with FOD metrics. We also show how RGFs can be used for interpolation, yielding intuitively more reasonable results than the Frechet mean on the unit sphere in a Hilbert space [1].

Typical and atypical brain development: a review of neuroimaging studies

In the course of development, the brain undergoes a remarkable process of restructuring as it adapts to the environment and becomes more efficient in processing information. A variety of brain imaging methods can be used to probe how anatomy, connectivity, and function change in the developing brain. Here we review recent discoveries regarding these brain changes in both typically developing individuals and individuals with neurodevelopmental disorders. We begin with typical development, summarizing research on changes in regional brain volume and tissue density, cortical thickness, white matter integrity, and functional connectivity. Space limits preclude the coverage of all neurodevelopmental disorders; instead, we cover a representative selection of studies examining neural correlates of autism, attention deficit/hyperactivity disorder, Fragile X, 22q11.2 deletion syndrome, Williams syndrome, Down syndrome, and Turner syndrome. Where possible, we focus on studies that identify an age by diagnosis interaction, suggesting an altered developmental trajectory. The studies we review generally cover the developmental period from infancy to early adulthood. Great progress has been made over the last 20 years in mapping how the brain matures with MR technology. With ever-improving technology, we expect this progress to accelerate, offering a deeper understanding of brain development, and more effective interventions for neurodevelopmental disorders.

Tracking on the Product Manifold of Shape and Orientation for Tractography from Diffusion MRI

Tractography refers to the process of tracing out the nerve fiber bundles from diffusion Magnetic Resonance Images (dMRI) data acquired either in vivo or ex-vivo. Tractography is a mature research topic within the field of diffusion MRI analysis, nevertheless, several new methods are being proposed on a regular basis thereby justifying the need, as the problem is not fully solved. Tractography is usually applied to the model (used to represent the diffusion MR signal or a derived quantity) reconstructed from the acquired data. Separating shape and orientation of these models was previously shown to approximately preserve diffusion anisotropy (a useful bio-marker) in the ubiquitous problem of interpolation. However, no further intrinsic geometric properties of this framework were exploited to date in literature. In this paper, we propose a new intrinsic recursive filter on the product manifold of shape and orientation. The recursive filter, dubbed IUKFPro, is a generalization of the unscented Kalman filter (UKF) to this product manifold. The salient contributions of this work are: (1) A new intrinsic UKF for the product manifold of shape and orientation. (2) Derivation of the Riemannian geometry of the product manifold. (3) IUKFPro is tested on synthetic and real data sets from various tractography challenge competitions. From the experimental results, it is evident that IUKFPro performs better than several competing schemes in literature with regards to some of the error measures used in the competitions and is competitive with respect to others.