4D hyperspherical harmonic (HyperSPHARM) representation of surface anatomy: A holistic treatment of multiple disconnected anatomical structures

Image-based parcellation of the brain often leads to multiple disconnected anatomical structures, which pose significant challenges for analyses of morphological shapes. Existing shape models, such as the widely used spherical harmonic (SPHARM) representation, assume topological invariance, so are unable to simultaneously parameterize multiple disjoint structures. In such a situation, SPHARM has to be applied separately to each individual structure. We present a novel surface parameterization technique using 4D hyperspherical harmonics in representing multiple disjoint objects as a single analytic function, terming it HyperSPHARM. The underlying idea behind HyperSPHARM is to stereographically project an entire collection of disjoint 3D objects onto the 4D hypersphere and subsequently simultaneously parameterize them with the 4D hyperspherical harmonics. Hence, HyperSPHARM allows for a holistic treatment of multiple disjoint objects, unlike SPHARM. In an imaging dataset of healthy adult human brains, we apply HyperSPHARM to the hippocampi and amygdalae. The HyperSPHARM representations are employed as a data smoothing technique, while the HyperSPHARM coefficients are utilized in a support vector machine setting for object classification. HyperSPHARM yields nearly identical results as SPHARM, as will be shown in the paper. Its key advantage over SPHARM lies computationally; HyperSPHARM possess greater computational efficiency than SPHARM because it can parameterize multiple disjoint structures using much fewer basis functions and stereographic projection obviates SPHARM's burdensome surface flattening. In addition, HyperSPHARM can handle any type of topology, unlike SPHARM, whose analysis is confined to topologically invariant structures.

[1]  Richard M. Leahy,et al.  Optimization method for creating semi-isometric flat maps of the cerebral cortex , 2000, Medical Imaging: Image Processing.

[2]  G. Domokos,et al.  Four-Dimensional Symmetry , 1967 .

[3]  Corinna Cortes,et al.  Support-Vector Networks , 1995, Machine Learning.

[4]  Cheng Guan Koay,et al.  A signal transformational framework for breaking the noise floor and its applications in MRI. , 2009, Journal of magnetic resonance.

[5]  James S. Duncan,et al.  Medical Image Analysis , 1999, IEEE Pulse.

[6]  A M Dale,et al.  Measuring the thickness of the human cerebral cortex from magnetic resonance images. , 2000, Proceedings of the National Academy of Sciences of the United States of America.

[7]  A 3D shape descriptor: 4D hyperspherical harmonics an exploration into the fourth dimension , 2007 .

[8]  E. P. Ferreira,et al.  The interactions of π - -mesons with complex nuclei in the energy range (100–800) MeV. III. The interaction lengths and elastic scattering of 300 MeV π - -mesons in G5 emulsion , 1959 .

[9]  Fillia Makedon,et al.  A surface-based approach for classification of 3D neuroanatomic structures , 2004, Intell. Data Anal..

[10]  Farid G. Mitri,et al.  Improving the use of vibro-acoustography for brachytherapy metal seed imaging: A feasibility study , 2004, IEEE Transactions on Medical Imaging.

[11]  Cheng Guan Koay,et al.  A simple scheme for generating nearly uniform distribution of antipodally symmetric points on the unit sphere , 2011, J. Comput. Sci..

[12]  Brian B. Avants,et al.  Symmetric diffeomorphic image registration with cross-correlation: Evaluating automated labeling of elderly and neurodegenerative brain , 2008, Medical Image Anal..

[13]  Cheng Guan Koay,et al.  Distributing points uniformly on the unit sphere under a mirror reflection symmetry constraint , 2014, J. Comput. Sci..

[14]  Paul M. Thompson,et al.  Genus zero surface conformal mapping and its application to brain surface mapping , 2004, IEEE Transactions on Medical Imaging.

[15]  Douglas W. Jones,et al.  Shape analysis of brain ventricles using SPHARM , 2001, Proceedings IEEE Workshop on Mathematical Methods in Biomedical Image Analysis (MMBIA 2001).

[16]  Stephen M Smith,et al.  Fast robust automated brain extraction , 2002, Human brain mapping.

[17]  Christopher A. Schuh,et al.  Hyperspherical harmonics for the representation of crystallographic texture , 2008 .

[18]  V. Bargmann,et al.  Zur Theorie des Wasserstoffatoms , 1936 .

[19]  Moo K. Chung,et al.  General multivariate linear modeling of surface shapes using SurfStat , 2010, NeuroImage.

[20]  Yalin Wang,et al.  Disease classification with hippocampal shape invariants , 2009, Hippocampus.

[21]  Moo K. Chung,et al.  Weighted Fourier Series Representation and Its Application to Quantifying the Amount of Gray Matter , 2007, IEEE Transactions on Medical Imaging.

[22]  Vincenzo Aquilanti,et al.  The d-dimensional hydrogen atom: hyperspherical harmonics as momentum space orbitals and alternative Sturmian basis sets , 1997 .

[23]  William E. Lorensen,et al.  Marching cubes: A high resolution 3D surface construction algorithm , 1987, SIGGRAPH.

[24]  Ron Kikinis,et al.  On the Laplace-Beltrami operator and brain surface flattening , 1999, IEEE Transactions on Medical Imaging.

[25]  Moo K. Chung,et al.  4D Hyperspherical Harmonic (HyperSPHARM) Representation of Multiple Disconnected Brain Subcortical Structures , 2013, MICCAI.

[26]  Catherine J. Norris,et al.  Aging is associated with positive responding to neutral information but reduced recovery from negative information. , 2011, Social cognitive and affective neuroscience.

[27]  Y. Benjamini,et al.  Controlling the false discovery rate: a practical and powerful approach to multiple testing , 1995 .

[28]  N. Given 4 D Hyperspherical Harmonic ( HyperSPHARM ) Representation of Multiple Disconnected Brain Subcortical Structures , 2013 .

[29]  Kenneth Stephenson,et al.  Cortical cartography using the discrete conformal approach of circle packings , 2004, NeuroImage.

[30]  Martin Styner,et al.  Framework for the Statistical Shape Analysis of Brain Structures using SPHARM-PDM. , 2006, The insight journal.

[31]  Cheng Guan Koay,et al.  Pseudometrically constrained centroidal voronoi tessellations: Generating uniform antipodally symmetric points on the unit sphere with a novel acceleration strategy and its applications to diffusion and three‐dimensional radial MRI , 2012, Magnetic resonance in medicine.

[32]  Guido Gerig,et al.  Parametrization of Closed Surfaces for 3-D Shape Description , 1995, Comput. Vis. Image Underst..