Files:
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.
BibTeX:
@article{gtmd-mia-13,
hal_id = {hal-00815120},
url = {http://hal.archives-ouvertes.fr/hal-00815120},
title = {{A polynomial approach for extracting the extrema of
a spherical function and its application in
diffusion MRI}},
author = {Ghosh, Aurobrata and Tsigaridas, Elias and Mourrain,
Bernard and Deriche, Rachid},
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.}},
keywords = {* HARDI; * ODF; * Maxima; * Tensors; * Polynomials},
language = {English},
affiliation = {Athena - INRIA Sophia Antipolis , Laboratoire
d'Informatique de Paris 6 - LIP6 , POLSYS - INRIA
Paris-Rocquencourt , GALAAD - INRIA Sophia
Antipolis},
pages = {503--514},
journal = {Medical Image Analysis},
volume = 17,
number = 5,
audience = {international},
year = 2013,
month = Jul,
}
Generated by bib2html.pl (written by Patrick Riley , modified by Elias ) on Wed Oct 23, 2019 21:41:02