Finsler Geometry on Higher Order Tensor Fields and Applications to High Angular Resolution Diffusion Imaging

Automatic tractography and segmentation using finsler geometry based on higher-order tensor fields

BACKGROUND AND OBJECTIVE

We focus on three-dimensional higher-order tensorial (HOT) images using Finsler geometry. In biomedical image analysis, these images are widely used, and they are based on the diffusion profiles inside the voxels. The diffusion information is stored in the so-called diffusion tensor D. Our objective is to present new methods revealing the architecture of neural fibers in presence of crossings and high curvatures. After tracking the fibers, we achieve direct 3D image segmentation to analyse the brain's white matter structures.

METHODS

To deal with the construction of the underlying fibers, the inverse of the second-order diffusion tensor D, understood as the metric tensor D-1, is commonly used in DTI modality. For crossing and highly curved fibers, higher order tensors are more relevant, but it is challenging to find an analogue of such an inverse in the HOT case. We employ an innovative approach to metrics based on higher order tensors to track the fibers properly. We propose to feed the tracked fibers as the internal initial contours in an efficient version of 3D segmentation.

RESULTS

We propose a brand-new approach to the inversion of a diffusion HOT, and an effective way of fiber tracking in the Finsler setting, based on innovative classification of the individual voxels. Thus, we can handle complex structures with high curvatures and crossings, even in the presence of noise. Based on our novel tractography approach, we also introduce a new segmentation method. We feed the detected fibers as the initial position of the contour surfaces to segment the image using a relevant active contour method (i.e., initiating the segmentation from inside the structures).

CONCLUSIONS

This is a pilot work, enhancing methods for fiber tracking and segmentation. The implemented algorithms were successfully tested on both synthetic and real data. The new features make our algorithms robust and fast, and they allow distinguishing individual objects in complex structures, even under noise.

Geodesic Uncertainty in Diffusion MRI

We study theoretical and operational issues of geodesic tractography, a geometric methodology for retrieving biologically plausible neural fibers in the brain from diffusion weighted magnetic resonance imaging. The premise is that true positives are geodesics in a suitably constructed metric space, but unlike traditional first order methods these are not a priori constrained to connect nongeneric points on subdimensional manifolds, such as the characteristics in traditional streamline methods. By virtue of the Hopf-Rinow theorem geodesic tractography furnishes a huge amount of redundancy, ensuring the a priori existence of at least one tentative fiber between any two points and permitting additional tractometric and data-extrinsic constraints for (fuzzy or crisp) classification of true and false positives. In our feasibility study we consider a hybrid paradigm that unifies existing ideas on tractography, combining deterministic and probabilistic elements in a way naturally supported by metric geometry. Particular attention is paid to an analytical prediction of geodesic deviation on numerically computed geodesics, a ‘tidal’ effect induced by small perturbations resulting from data noise. Taking these effects into account clarifies the inherent uncertainty of geodesics, while simultaneosuly offering a dimensionality reduction of the tractography problem.

Optimal Paths for Variants of the 2D and 3D Reeds-Shepp Car with Applications in Image Analysis

We present a PDE-based approach for finding optimal paths for the Reeds-Shepp car. In our model we minimize a (data-driven) functional involving both curvature and length penalization, with several generalizations. Our approach encompasses the two- and three-dimensional variants of this model, state-dependent costs, and moreover, the possibility of removing the reverse gear of the vehicle. We prove both global and local controllability results of the models. Via eikonal equations on the manifoldR d × S d - 1 we compute distance maps w.r.t. highly anisotropic Finsler metrics, which approximate the singular (quasi)-distances underlying the model. This is achieved using a fast-marching (FM) method, building on Mirebeau (Numer Math 126(3):515-557, 2013; SIAM J Numer Anal 52(4):1573-1599, 2014). The FM method is based on specific discretization stencils which are adapted to the preferred directions of the Finsler metric and obey a generalized acuteness property. The shortest paths can be found with a gradient descent method on the distance map, which we formalize in a theorem. We justify the use of our approximating metrics by proving convergence results. Our curve optimization model inR d × S d - 1 with data-driven cost allows to extract complex tubular structures from medical images, e.g., crossings, and incomplete data due to occlusions or low contrast. Our work extends the results of Sanguinetti et al. (Progress in Pattern Recognition, Image Analysis, Computer Vision, and Applications LNCS 9423, 2015) on numerical sub-Riemannian eikonal equations and the Reeds-Shepp car to 3D, with comparisons to exact solutions by Duits et al. (J Dyn Control Syst 22(4):771-805, 2016). Numerical experiments show the high potential of our method in two applications: vessel tracking in retinal images for the case d = 2 and brain connectivity measures from diffusion-weighted MRI data for the case d = 3 , extending the work of Bekkers et al. (SIAM J Imaging Sci 8(4):2740-2770, 2015). We demonstrate how the new model without reverse gear better handles bifurcations.

APPROXIMATING SYMMETRIC POSITIVE SEMIDEFINITE TENSORS OF EVEN ORDER()

Tensors of various orders can be used for modeling physical quantities such as strain and diffusion as well as curvature and other quantities of geometric origin. Depending on the physical properties of the modeled quantity, the estimated tensors are often required to satisfy the positivity constraint, which can be satisfied only with tensors of even order. Although the space [Formula: see text] of 2m(th)-order symmetric positive semi-definite tensors is known to be a convex cone, enforcing positivity constraint directly on [Formula: see text] is usually not straightforward computationally because there is no known analytic description of [Formula: see text] for m > 1. In this paper, we propose a novel approach for enforcing the positivity constraint on even-order tensors by approximating the cone [Formula: see text] for the cases 0 < m < 3, and presenting an explicit characterization of the approximation Σ(2) (m) ⊂ Ω(2) (m) for m ≥ 1, using the subset [Formula: see text] of semi-definite tensors that can be written as a sum of squares of tensors of order m. Furthermore, we show that this approximation leads to a non-negative linear least-squares (NNLS) optimization problem with the complexity that equals the number of generators in Σ(2) (m). Finally, we experimentally validate the proposed approach and we present an application for computing 2m(th)-order diffusion tensors from Diffusion Weighted Magnetic Resonance Images.

Medians and means in Finsler geometry

We investigate existence and uniqueness of p-means ep and the median e1 of a probability measure μ on a Finsler manifold, in relation with the convexity of the support of μ. We prove that ep is the limit point of a continuous time gradient flow. Under some additional condition which is always satisfied for p≥2, a discretization of this path converges to ep. This provides an algorithm for determining the Finsler center points.