Brain Connectivity Mapping Using Riemannian Geometry, Control Theory, and PDEs

We introduce an original approach for the cerebral white matter connectivity mapping from diffusion tensor imaging (DTI). Our method relies on a global modeling of the acquired magnetic resonance imaging volume as a Riemannian manifold whose metric directly derives from the diffusion tensor. These tensors will be used to measure physical three-dimensional distances between different locations of a brain diffusion tensor image. The key concept is the notion of geodesic distance that will allow us to find optimal paths in the white matter. We claim that such optimal paths are reasonable approximations of neural fiber bundles. The geodesic distance function can be seen as the solution of two theoretically equivalent but, in practice, significantly different problems in the partial differential equation framework: an initial value problem which is intrinsically dynamic, and a boundary value problem which is, on the contrary, intrinsically stationary. The two approaches have very different properties which make them more or less adequate for our problem and more or less computationally efficient. The dynamic formulation is quite easy to implement but has several practical drawbacks. On the contrary, the stationary formulation is much more tedious to implement; we will show, however, that it has many virtues which make it more suitable for our connectivity mapping problem. Finally, we will present different possible measures of connectivity, reflecting the degree of connectivity between different regions of the brain. We will illustrate these notions on synthetic and real DTI datasets.

[1]  D. Doman,et al.  Optimal Control Problems on Parallelizable Riemannian Manifolds: Theory and Applications , 2006 .

[2]  J. Tsitsiklis,et al.  Efficient algorithms for globally optimal trajectories , 1994, Proceedings of 1994 33rd IEEE Conference on Decision and Control.

[3]  J A Sethian,et al.  A fast marching level set method for monotonically advancing fronts. , 1996, Proceedings of the National Academy of Sciences of the United States of America.

[4]  Stanley Osher,et al.  Level Set Methods , 2003 .

[5]  Rachid Deriche,et al.  Inferring White Matter Geometry from Di.usion Tensor MRI: Application to Connectivity Mapping , 2004, ECCV.

[6]  Hongkai Zhao,et al.  A Fast Sweeping Method for Static Convex Hamilton–Jacobi Equations , 2007, J. Sci. Comput..

[7]  Qingfen Lin,et al.  Enhancement, Extraction, and Visualization of 3D Volume Data , 2001 .

[8]  R. Passingham,et al.  Initial Demonstration of in Vivo Tracing of Axonal Projections in the Macaque Brain and Comparison with the Human Brain Using Diffusion Tensor Imaging and Fast Marching Tractography , 2002, NeuroImage.

[9]  Ross T. Whitaker A FAST EIKONAL EQUATION SOLVER FOR PARALLEL SYSTEMS , 2007 .

[10]  Jean Daunizeau,et al.  Accurate Anisotropic Fast Marching for Diffusion-Based Geodesic Tractography , 2007, Int. J. Biomed. Imaging.

[11]  Arthur W. Toga,et al.  A Diffusion Tensor Imaging Tractography Algorithm Based on Navier–Stokes Fluid Mechanics , 2009, IEEE Transactions on Medical Imaging.

[12]  Alexander Vladimirsky,et al.  Ordered Upwind Methods for Static Hamilton-Jacobi Equations: Theory and Algorithms , 2003, SIAM J. Numer. Anal..

[13]  Carl-Fredrik Westin,et al.  Uncertainty in White Matter Fiber Tractography , 2005, MICCAI.

[14]  Jean-Philippe Thiran,et al.  Mapping Brain Connectivity with Statistical Fibre Tracking and Virtual Dissection , 2003 .

[15]  A. Alexander,et al.  White matter tractography using diffusion tensor deflection , 2003, Human brain mapping.

[16]  Tim McGraw,et al.  Stochastic DT-MRI Connectivity Mapping on the GPU , 2007, IEEE Transactions on Visualization and Computer Graphics.

[17]  S. Osher A level set formulation for the solution of the Dirichlet problem for Hamilton-Jacobi equations , 1993 .

[18]  B. Vemuri,et al.  Fiber tract mapping from diffusion tensor MRI , 2001, Proceedings IEEE Workshop on Variational and Level Set Methods in Computer Vision.

[19]  Roland Bammer,et al.  A Physical Model for DT-MRI Based Connectivity Map Computation , 2005, MICCAI.

[20]  Emmanuel Prados eprados,et al.  Control Theory and Fast Marching Methods for Brain Connectivity Mapping , 2022 .

[21]  Alper Yilmaz,et al.  Level Set Methods , 2007, Wiley Encyclopedia of Computer Science and Engineering.

[22]  Yoshikazu Giga,et al.  A level set approach for computing discontinuous solutions of Hamilton-Jacobi equations , 2003, Math. Comput..

[23]  Stefano Soatto,et al.  Fast Marching Method for Generic Shape from Shading , 2005, VLSM.

[24]  D. Le Bihan,et al.  A framework based on spin glass models for the inference of anatomical connectivity from diffusion‐weighted MR data – a technical review , 2002, NMR in biomedicine.

[25]  E. Rouy,et al.  A viscosity solutions approach to shape-from-shading , 1992 .

[26]  O. Faugeras,et al.  A VISCOSITY SOLUTION METHOD FOR SHAPE-FROM-SHADING WITHOUT IMAGE BOUNDARY DATA , 2006 .

[27]  R. Davidson,et al.  Probabilistic Connectivity Measure in Diffusion Tensor Imaging via Anisotropic Kernel Smoothing , 2003 .

[28]  C. Büchel,et al.  Estimation of anatomical connectivity from diffusion tensor data , 2001, NeuroImage.

[29]  R. W. R. Darling Intrinsic Location Parameter of a Diffusion Process , 1998 .

[30]  Olivier D. Faugeras,et al.  Shape from Shading and Viscosity Solutions , 2002, ECCV.

[31]  S. Osher,et al.  Regular ArticleUniformly High Order Accurate Essentially Non-oscillatory Schemes, III , 1997 .

[32]  R. Kimmel,et al.  Finding shortest paths on surfaces , 1994 .

[33]  David B. Doman,et al.  Optimal Control Problems on Riemannian Manifolds: Theory and Applications , 2002 .

[34]  Jean-Francois Mangin,et al.  Fiber Tracking in q-Ball Fields Using Regularized Particle Trajectories , 2005, IPMI.

[35]  Christopher I. Byrnes,et al.  Directions in Mathematical Systems : Theory and Optimization , 2003 .

[36]  Simon J. Doran,et al.  NMR imaging of fluids in porous solids , 1990, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[37]  R. LeVeque Numerical methods for conservation laws , 1990 .

[38]  Daniel Gembris,et al.  White matter fiber tractography via anisotropic diffusion simulation in the human brain , 2005, IEEE Transactions on Medical Imaging.

[39]  H. L. Dryden,et al.  Investigations on the Theory of the Brownian Movement , 1957 .

[40]  Carl-Fredrik Westin,et al.  Regularized Stochastic White Matter Tractography Using Diffusion Tensor MRI , 2002, MICCAI.

[41]  I. Holopainen Riemannian Geometry , 1927, Nature.

[42]  M. Cohen de Lara,et al.  GEOMETRIC AND SYMMETRY PROPERTIES OF A NONDEGENERATE DIFFUSION PROCESS , 1995 .

[43]  S. Osher,et al.  Weighted essentially non-oscillatory schemes , 1994 .

[44]  Guido Gerig,et al.  Probabilistic Fiber Tracking Using Particle Filtering , 2007, MICCAI.

[45]  L. Cohen Minimal Paths and Fast Marching Methods for Image Analysis , 2006, Handbook of Mathematical Models in Computer Vision.

[46]  H. Soner Optimal control with state-space constraint I , 1986 .

[47]  Lawrence H. Staib,et al.  White matter tractography by anisotropic wavefront evolution and diffusion tensor imaging , 2005, Medical Image Anal..

[48]  Emmanuel Prados Application of the theory of the viscosity solutions to the Shape From Shading problem , 2004 .

[49]  Gareth J. Barker,et al.  Estimating distributed anatomical connectivity using fast marching methods and diffusion tensor imaging , 2002, IEEE Transactions on Medical Imaging.

[50]  Timothy Edward John Behrens,et al.  Characterization and propagation of uncertainty in diffusion‐weighted MR imaging , 2003, Magnetic resonance in medicine.

[51]  D. Tuch Diffusion MRI of complex tissue structure , 2002 .

[52]  Stanley Osher,et al.  Fast Sweeping Algorithms for a Class of Hamilton-Jacobi Equations , 2003, SIAM J. Numer. Anal..

[53]  Rachid Deriche,et al.  Control Theory and Fast Marching Techniques for Brain Connectivity Mapping , 2006, 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'06).

[54]  Alan Connelly,et al.  Diffusion-weighted magnetic resonance imaging fibre tracking using a front evolution algorithm , 2003, NeuroImage.

[55]  Christophe Lenglet,et al.  Geometric and variational methods for diffusion tensor MRI processing. (Méthodes géométriques et variationnelles pour le traitement d'IRM du tenseur de diffusion) , 2006 .

[56]  J. Tsuruda,et al.  Diffusion-weighted MR imaging of anisotropic water diffusion in cat central nervous system. , 1990, Radiology.

[57]  Hervé Delingette,et al.  A Recursive Anisotropic Fast Marching Approach to Reaction Diffusion Equation: Application to Tumor Growth Modeling , 2007, IPMI.

[58]  Ross T. Whitaker,et al.  A Volumetric Approach to Quantifying Region-to-Region White Matter Connectivity in Diffusion Tensor MRI , 2007, IPMI.

[59]  David Atkinson,et al.  Study of Connectivity in the Brain Using the Full Diffusion Tensor from MRI , 2001, IPMI.

[60]  D. Chopp Computing Minimal Surfaces via Level Set Curvature Flow , 1993 .

[61]  James A. Sethian,et al.  The Fast Construction of Extension Velocities in Level Set Methods , 1999 .

[62]  Yun-Gang Chen,et al.  Uniqueness and existence of viscosity solutions of generalized mean curvature flow equations , 1989 .

[63]  P. V. van Zijl,et al.  Three‐dimensional tracking of axonal projections in the brain by magnetic resonance imaging , 1999, Annals of neurology.

[64]  R. Temam,et al.  Analyse convexe et problèmes variationnels , 1974 .

[65]  J A Sethian,et al.  Computing geodesic paths on manifolds. , 1998, Proceedings of the National Academy of Sciences of the United States of America.

[66]  Rachid Deriche,et al.  Statistics on the Manifold of Multivariate Normal Distributions: Theory and Application to Diffusion Tensor MRI Processing , 2006, Journal of Mathematical Imaging and Vision.

[67]  S. Osher,et al.  Algorithms Based on Hamilton-Jacobi Formulations , 1988 .

[68]  Jens Frahm,et al.  Self-diffusion NMR imaging using stimulated echoes , 1985 .

[69]  Daniel C. Alexander,et al.  Probabilistic Monte Carlo Based Mapping of Cerebral Connections Utilising Whole-Brain Crossing Fibre Information , 2003, IPMI.

[70]  Andrew L. Alexander,et al.  Bootstrap white matter tractography (BOOT-TRAC) , 2005, NeuroImage.

[71]  S. Osher,et al.  Uniformly high order accurate essentially non-oscillatory schemes, 111 , 1987 .

[72]  M. Bardi,et al.  Optimal Control and Viscosity Solutions of Hamilton-Jacobi-Bellman Equations , 1997 .

[73]  J. Sethian,et al.  Fronts propagating with curvature-dependent speed: algorithms based on Hamilton-Jacobi formulations , 1988 .

[74]  P. Lions Generalized Solutions of Hamilton-Jacobi Equations , 1982 .

[75]  Susumu Mori,et al.  Fiber tracking: principles and strategies – a technical review , 2002, NMR in biomedicine.

[76]  Abbas F. Sadikot,et al.  Flow-based fiber tracking with diffusion tensor and q-ball data: Validation and comparison to principal diffusion direction techniques , 2005, NeuroImage.

[77]  J. Sethian,et al.  FRONTS PROPAGATING WITH CURVATURE DEPENDENT SPEED: ALGORITHMS BASED ON HAMILTON-JACOB1 FORMULATIONS , 2003 .

[78]  M C Bushell,et al.  PRELIMINARY COMMUNICATION: The spatial mapping of translational diffusion coefficients by the NMR imaging technique , 1985 .

[79]  Geoffrey J M Parker,et al.  A framework for a streamline‐based probabilistic index of connectivity (PICo) using a structural interpretation of MRI diffusion measurements , 2003, Journal of magnetic resonance imaging : JMRI.

[80]  P. Basser,et al.  Estimation of the effective self-diffusion tensor from the NMR spin echo. , 1994, Journal of magnetic resonance. Series B.

[81]  Marc Niethammer,et al.  Finsler Tractography for White Matter Connectivity Analysis of the Cingulum Bundle , 2007, MICCAI.

[82]  S C Williams,et al.  Non‐invasive assessment of axonal fiber connectivity in the human brain via diffusion tensor MRI , 1999, Magnetic resonance in medicine.

[83]  P. Basser,et al.  In vivo fiber tractography using DT‐MRI data , 2000, Magnetic resonance in medicine.

[84]  P. Dupuis,et al.  An Optimal Control Formulation and Related Numerical Methods for a Problem in Shape Reconstruction , 1994 .

[85]  Carl-Fredrik Westin,et al.  New Approaches to Estimation of White Matter Connectivity in Diffusion Tensor MRI: Elliptic PDEs and Geodesics in a Tensor-Warped Space , 2002, MICCAI.

[86]  D. Parker,et al.  Analysis of partial volume effects in diffusion‐tensor MRI , 2001, Magnetic resonance in medicine.

[87]  Rachid Deriche,et al.  Diffusion tensor sharpening improves white matter tractography , 2007, SPIE Medical Imaging.

[88]  A. Mennucci,et al.  Hamilton—Jacobi Equations and Distance Functions on Riemannian Manifolds , 2002, math/0201296.

[89]  Guillermo Sapiro,et al.  O(N) implementation of the fast marching algorithm , 2006, Journal of Computational Physics.

[90]  Carl-Fredrik Westin,et al.  A Hamilton-Jacobi-Bellman Approach to High Angular Resolution Diffusion Tractography , 2005, MICCAI.

[91]  Denis Le Bihan,et al.  Imagerie de diffusion in-vivo par résonance magnétique nucléaire , 1985 .

[92]  Carl-Fredrik Westin,et al.  A Bayesian approach for stochastic white matter tractography , 2006, IEEE Transactions on Medical Imaging.

[93]  Arthur W. Toga,et al.  A diffusion tensor imaging tractography algorithm based on Navier-Stokes fluid mechanics , 2006, ISBI.

[94]  Alfred M. Bruckstein,et al.  Finding Shortest Paths on Surfaces Using Level Sets Propagation , 1995, IEEE Trans. Pattern Anal. Mach. Intell..