Principal geodesic analysis for the study of nonlinear statistics of shape

A primary goal of statistical shape analysis is to describe the variability of a population of geometric objects. A standard technique for computing such descriptions is principal component analysis. However, principal component analysis is limited in that it only works for data lying in a Euclidean vector space. While this is certainly sufficient for geometric models that are parameterized by a set of landmarks or a dense collection of boundary points, it does not handle more complex representations of shape. We have been developing representations of geometry based on the medial axis description or m-rep. While the medial representation provides a rich language for variability in terms of bending, twisting, and widening, the medial parameters are not elements of a Euclidean vector space. They are in fact elements of a nonlinear Riemannian symmetric space. In this paper, we develop the method of principal geodesic analysis, a generalization of principal component analysis to the manifold setting. We demonstrate its use in describing the variability of medially-defined anatomical objects. Results of applying this framework on a population of hippocampi in a schizophrenia study are presented.

[1]  HARRY BLUM,et al.  Shape description using weighted symmetric axis features , 1978, Pattern Recognit..

[2]  W. Boothby An introduction to differentiable manifolds and Riemannian geometry , 1975 .

[3]  D. Kendall SHAPE MANIFOLDS, PROCRUSTEAN METRICS, AND COMPLEX PROJECTIVE SPACES , 1984 .

[4]  Samuel R. Buss,et al.  Spherical averages and applications to spherical splines and interpolation , 2001, TOGS.

[5]  Timothy F. Cootes,et al.  Active Appearance Models , 1998, ECCV.

[6]  W. Kendall Probability, Convexity, and Harmonic Maps with Small Image I: Uniqueness and Fine Existence , 1990 .

[7]  Maher Moakher,et al.  To appear in: SIAM J. MATRIX ANAL. APPL. MEANS AND AVERAGING IN THE GROUP OF ROTATIONS∗ , 2002 .

[8]  H. Karcher Riemannian center of mass and mollifier smoothing , 1977 .

[9]  P. Thomas Fletcher,et al.  Statistics of shape via principal geodesic analysis on Lie groups , 2003, 2003 IEEE Computer Society Conference on Computer Vision and Pattern Recognition, 2003. Proceedings..

[10]  Martin Styner,et al.  Medial Models Incorporating Object Variability for 3D Shape Analysis , 2001, IPMI.

[11]  D. Kendall A Survey of the Statistical Theory of Shape , 1989 .

[12]  M. Fréchet Les éléments aléatoires de nature quelconque dans un espace distancié , 1948 .

[13]  Paul A. Yushkevich,et al.  Multiscale deformable model segmentation and statistical shape analysis using medial descriptions , 2002, IEEE Transactions on Medical Imaging.

[14]  Xavier Pennec,et al.  Probabilities and statistics on Riemannian manifolds: Basic tools for geometric measurements , 1999, NSIP.

[15]  F. Bookstein Size and Shape Spaces for Landmark Data in Two Dimensions , 1986 .

[16]  Gábor Székely,et al.  Multiscale Medial Loci and Their Properties , 2003, International Journal of Computer Vision.

[17]  C. Goodall Procrustes methods in the statistical analysis of shape , 1991 .

[18]  Anuj Srivastava,et al.  Monte Carlo extrinsic estimators of manifold-valued parameters , 2002, IEEE Trans. Signal Process..

[19]  S. Helgason Differential Geometry, Lie Groups, and Symmetric Spaces , 1978 .

[20]  Michael I. Miller,et al.  Hilbert-Schmidt Lower Bounds for Estimators on Matrix Lie Groups for ATR , 1998, IEEE Trans. Pattern Anal. Mach. Intell..

[21]  Andrew Thall,et al.  Fast C 2 Interpolating Subdivision Surfaces using Iterative Inversion of Stationary Subdivision Rules , 2003 .

[22]  Guido Gerig,et al.  Elastic model-based segmentation of 3-D neuroradiological data sets , 1999, IEEE Transactions on Medical Imaging.

[23]  Rikard Berthilsson,et al.  A Statistical Theory of Shape , 1998, SSPR/SPR.

[24]  Paul A. Yushkevich,et al.  Deformable M-Reps for 3D Medical Image Segmentation , 2003, International Journal of Computer Vision.

[25]  P. Thomas Fletcher,et al.  Gaussian Distributions on Lie Groups and Their Application to Statistical Shape Analysis , 2003, IPMI.

[26]  Timothy F. Cootes,et al.  Active Shape Models-Their Training and Application , 1995, Comput. Vis. Image Underst..

[27]  U. Grenander,et al.  Hippocampal morphometry in schizophrenia by high dimensional brain mapping. , 1998, Proceedings of the National Academy of Sciences of the United States of America.