Angelos Barmpoutis

Regularized positive-definite fourth order tensor field estimation from DW-MRI.

In Diffusion Weighted Magnetic Resonance Image (DW-MRI) processing, a 2nd order tensor has been commonly used to approximate the diffusivity function at each lattice point of the DW-MRI data. From this tensor approximation, one can compute useful scalar quantities (e.g. anisotropy, mean diffusivity) which have been clinically used for monitoring encephalopathy, sclerosis, ischemia and other brain disorders. It is now well known that this 2nd-order tensor approximation fails to capture complex local tissue structures, e.g. crossing fibers, and as a result, the scalar quantities derived from these tensors are grossly inaccurate at such locations. In this paper we employ a 4th order symmetric positive-definite (SPD) tensor approximation to represent the diffusivity function and present a novel technique to estimate these tensors from the DW-MRI data guaranteeing the SPD property. Several articles have been reported in literature on higher order tensor approximations of the diffusivity function but none of them guarantee the positivity of the estimates, which is a fundamental constraint since negative values of the diffusivity are not meaningful. In this paper we represent the 4th-order tensors as ternary quartics and then apply Hilberts theorem on ternary quartics along with the Iwasawa parametrization to guarantee an SPD 4th-order tensor approximation from the DW-MRI data. The performance of this model is depicted on synthetic data as well as real DW-MRIs from a set of excised control and injured rat spinal cords, showing accurate estimation of scalar quantities such as generalized anisotropy and trace as well as fiber orientations.

Tensor Splines for Interpolation and Approximation of DT-MRI With Applications to Segmentation of Isolated Rat Hippocampi

In this paper, we present novel algorithms for statistically robust interpolation and approximation of diffusion tensors-which are symmetric positive definite (SPD) matrices-and use them in developing a significant extension to an existing probabilistic algorithm for scalar field segmentation, in order to segment diffusion tensor magnetic resonance imaging (DT-MRI) datasets. Using the Riemannian metric on the space of SPD matrices, we present a novel and robust higher order (cubic) continuous tensor product of -splines algorithm to approximate the SPD diffusion tensor fields. The resulting approximations are appropriately dubbed tensor splines. Next, we segment the diffusion tensor field by jointly estimating the label (assigned to each voxel) field, which is modeled by a Gauss Markov measure field (GMMF) and the parameters of each smooth tensor spline model representing the labeled regions. Results of interpolation, approximation, and segmentation are presented for synthetic data and real diffusion tensor fields from an isolated rat hippocampus, along with validation. We also present comparisons of our algorithms with existing methods and show significantly improved results in the presence of noise as well as outliers.

Symmetric positive 4th order tensors & their estimation from diffusion weighted MRI

In Diffusion Weighted Magnetic Resonance Image (DW-MRI) processing a 2nd order tensor has been commonly used to approximate the diffusivity function at each lattice point of the DW-MRI data. It is now well known that this 2nd-order approximation fails to approximate complex local tissue structures, such as fibers crossings. In this paper we employ a 4th order symmetric positive semi-definite (PSD) tensor approximation to represent the diffusivity function and present a novel technique to estimate these tensors from the DW-MRI data guaranteeing the PSD property. There have been several published articles in literature on higher order tensor approximations of the diffusivity function but none of them guarantee the positive semi-definite constraint, which is a fundamental constraint since negative values of the diffusivity coefficients are not meaningful. In our methods, we parameterize the 4th order tensors as a sum of squares of quadratic forms by using the so called Gram matrix method from linear algebra and its relation to the Hilberts theorem on ternary quartics. This parametric representation is then used in a nonlinear-least squares formulation to estimate the PSD tensors of order 4 from the data. We define a metric for the higher-order tensors and employ it for regularization across the lattice. Finally, performance of this model is depicted on synthetic data as well as real DW-MRI from an isolated rat hippocampus.

A unified framework for estimating diffusion tensors of any order with symmetric positive-definite constraints

Cartesian tensors of various orders have been employed for either modeling the diffusivity or the orientation distribution function in Diffusion-Weighted MRI datasets. In both cases, the estimated tensors have to be positive-definite since they model positive-valued functions. In this paper we present a novel unified framework for estimating positive-definite tensors of any order, in contrast to the existing methods in literature, which are either order-specific or fail to handle the positive-definite property. The proposed framework employs a homogeneous polynomial parametrization that covers the full space of any order positive-definite tensors and explicitly imposes the positive-definite constraint on the estimated tensors. We show that this parametrization leads to a linear system that is solved using the non-negative least squares technique. The framework is demonstrated using synthetic and real data from an excised rat hippocampus.

Tensor Body: Real-Time Reconstruction of the Human Body and Avatar Synthesis From RGB-D

Real-time 3-D reconstruction of the human body has many applications in anthropometry, telecommunications, gaming, fashion, and other areas of human-computer interaction. In this paper, a novel framework is presented for reconstructing the 3-D model of the human body from a sequence of RGB-D frames. The reconstruction is performed in real time while the human subject moves arbitrarily in front of the camera. The method employs a novel parameterization of cylindrical-type objects using Cartesian tensor and b-spline bases along the radial and longitudinal dimension respectively. The proposed model, dubbed tensor body, is fitted to the input data using a multistep framework that involves segmentation of the different body regions, robust filtering of the data via a dynamic histogram, and energy-based optimization with positive-definite constraints. A Riemannian metric on the space of positive-definite tensor splines is analytically defined and employed in this framework. The efficacy of the presented methods is demonstrated in several real-data experiments using the Microsoft Kinect sensor.

Registration of high angular resolution diffusion MRI images using 4th order tensors

Registration of Diffusion Weighted (DW)-MRI datasets has been commonly achieved to date in literature by using either scalar or 2nd-order tensorial information. However, scalar or 2nd-order tensors fail to capture complex local tissue structures, such as fiber crossings, and therefore, datasets containing fiber-crossings cannot be registered accurately by using these techniques. In this paper we present a novel method for non-rigidly registering DW-MRI datasets that are represented by a field of 4th-order tensors. We use the Hellinger distance between the normalized 4th-order tensors represented as distributions, in order to achieve this registration. Hellinger distance is easy to compute, is scale and rotation invariant and hence allows for comparison of the true shape of distributions. Furthermore, we propose a novel 4th-order tensor re-transformation operator, which plays an essential role in the registration procedure and shows significantly better performance compared to the re-orientation operator used in literature for DTI registration. We validate and compare our technique with other existing scalar image and DTI registration methods using simulated diffusion MR data and real HARDI datasets.

Extracting Tractosemas from a Displacement Probability Field for Tractography in DW-MRI

In this paper we present a novel method for estimating a field of asymmetric spherical functions, dubbed tractosemas, given the intra-voxel displacement probability information. The peaks of tractosemas correspond to directions of distinct fibers, which can have either symmetric or asymmetric local fiber structure. This is in contrast to the existing methods that estimate fiber orientation distributions which are naturally symmetric and therefore cannot model asymmetries such as splaying fibers. We propose a method for extracting tractosemas from a given field of displacement probability iso-surfaces via a diffusion process. The diffusion is performed by minimizing a kernel convolution integral, which leads to an update formula expressed in the convenient form of a discrete kernel convolution. The kernel expresses the probability of diffusion between two neighboring spherical functions and we model it by the product of Gaussian and von Mises distributions. The model is validated via experiments on synthetic and real diffusion-weighted magnetic resonance (DW-MRI) datasets from a rat hippocampus and spinal cord.

Fast displacement probability profile approximation from HARDI using 4th-order tensors

Cartesian tensor basis have been widely used to approximate spherical functions. In Medical Imaging, tensors of various orders have been used to model the diffusivity function in Diffusion-weighted MRI data sets. However, it is known that the peaks of the diffusivity do not correspond to orientations of the underlying fibers and hence the displacement probability profiles should be employed instead. In this paper, we present a novel representation of the probability profile by a 4th order tensor, which is a smooth spherical function that can approximate single-fibers as well as multiple-fiber structures. We also present a method for efficiently estimating the unknown tensor coefficients of the probability profile directly from a given high-angular resolution diffusion-weighted (HARDI) data set. The accuracy of our model is validated by experiments on synthetic and real HARDI datasets from a fixed rat spinal cord.

Symmetric positive semi-definite Cartesian Tensor fiber orientation distributions (CT-FOD)

A novel method for estimating a field of fiber orientation distribution (FOD) based on signal de-convolution from a given set of diffusion weighted magnetic resonance (DW-MR) images is presented. We model the FOD by higher order Cartesian tensor basis using a parametrization that explicitly enforces the positive semi-definite property to the computed FOD. The computed Cartesian tensors, dubbed Cartesian Tensor-FOD (CT-FOD), are symmetric positive semi-definite tensors whose coefficients can be efficiently estimated by solving a linear system with non-negative constraints. Next, we show how to use our method for converting higher-order diffusion tensors to CT-FODs, which is an essential task since the maxima of higher-order tensors do not correspond to the underlying fiber orientations. Finally, we propose a diffusion anisotropy index computed directly from CT-FODs using higher order tensor distance measures thus consolidating the whole analysis pipeline of diffusion imaging solely using CT-FODs. We evaluate our method qualitatively and quantitatively using simulated DW-MR images, phantom images, and human brain real dataset. The results conclusively demonstrate the superiority of the proposed technique over several existing multi-fiber reconstruction methods.

Multi-fiber reconstruction from DW-MRI using a continuous mixture of von Mises-Fisher distributions

In this paper we propose a method for reconstructing the Diffusion Weighted Magnetic Resonance (DW-MR) signal at each lattice point using a novel continuous mixture of von Mises-Fisher distribution functions. Unlike most existing methods, neither does this model assume a fixed functional form for the MR signal attenuation (e.g. 2nd or 4th order tensor) nor does it arbitrarily fix important mixture parameters like the number of components. We show that this continuous mixture has a closed form expression and leads to a linear system which can be easily solved. Through extensive experimentation with synthetic data we show that this technique outperforms various other state-of-the-art techniques in resolving fiber crossings. Finally, we demonstrate the effectiveness of this method using real DW-MRI data from rat brain and optic chiasm.