A polynomial approach for extracting the extrema of a spherical function and its application in diffusion MRI

Aurobrata Ghosh 1 Elias Tsigaridas 2 Bernard Mourrain 3 Rachid Deriche 1
1 ATHENA - Computational Imaging of the Central Nervous System
CRISAM - Inria Sophia Antipolis - Méditerranée
2 PolSys - Polynomial Systems
LIP6 - Laboratoire d'Informatique de Paris 6, Inria Paris-Rocquencourt
3 GALAAD - Geometry, algebra, algorithms
CRISAM - Inria Sophia Antipolis - Méditerranée , UNS - Université Nice Sophia Antipolis, CNRS - Centre National de la Recherche Scientifique : UMR6621
Abstract : Antipodally symmetric spherical functions play a pivotal role in diffusion MRI in representing sub-voxel-resolution microstructural information of the underlying tissue. This information is described by the geometry of the spherical function. In this paper we propose a method to automatically compute all the extrema of a spherical function. We then classify the extrema as maxima, minima and saddle-points to identify the maxima. We take advantage of the fact that a spherical function can be described equivalently in the spherical harmonic (SH) basis, in the symmetric tensor (ST) basis constrained to the sphere, and in the homogeneous polynomial (HP) basis constrained to the sphere. We extract the extrema of the spherical function by computing the stationary points of its constrained HP representation. Instead of using traditional optimization approaches, which are inherently local and require exhaustive search or re-initializations to locate multiple extrema, we use a novel polynomial system solver which analytically brackets all the extrema and refines them numerically, thus missing none and achieving high precision. To illustrate our approach we consider the Orientation Distribution Function (ODF). In diffusion MRI the ODF is a spherical function which represents a state-of-the-art reconstruction algorithm whose maxima are aligned with the dominant fiber bundles. It is, therefore, vital to correctly compute these maxima to detect the fiber bundle directions. To demonstrate the potential of the proposed polynomial approach we compute the extrema of the ODF to extract all its maxima. This polynomial approach is, however, not dependent on the ODF and the framework presented in this paper can be applied to any spherical function described in either the SH basis, ST basis or the HP basis.
Type de document :
Article dans une revue
Medical Image Analysis, Elsevier, 2013, 17 (5), pp.503-514. <10.1016/j.media.2013.03.004>
Liste complète des métadonnées

Contributeur : Rachid Deriche <>
Soumis le : jeudi 18 avril 2013 - 11:13:42
Dernière modification le : samedi 7 novembre 2015 - 01:06:06




Aurobrata Ghosh, Elias Tsigaridas, Bernard Mourrain, Rachid Deriche. A polynomial approach for extracting the extrema of a spherical function and its application in diffusion MRI. Medical Image Analysis, Elsevier, 2013, 17 (5), pp.503-514. <10.1016/j.media.2013.03.004>. <hal-00815120>



Consultations de la notice