A Geometric Framework for Statistical Analysis of Trajectories with Distinct Temporal Spans

Analyzing data representing multifarious trajectories is central to the many fields in Science and Engineering; for example, trajectories representing a tennis serve, a gymnast’s parallel bar routine, progression/remission of disease and so on. We present a novel geometric algorithm for performing statistical analysis of trajectories with distinct number of samples representing longitudinal (or temporal) data. A key feature of our proposal is that unlike existing schemes, our model is deployable in regimes where each participant provides a different number of acquisitions (trajectories have different number of sample points or temporal span). To achieve this, we develop a novel method involving the parallel transport of the tangent vectors along each given trajectory to the starting point of the respective trajectories and then use the span of the matrix whose columns consist of these vectors, to construct a linear subspace in Rm. We then map these linear subspaces (possibly of distinct dimensions) of Rm on to a single high dimensional hypersphere. This enables computing group statistics over trajectories by instead performing statistics on the hypersphere (equipped with a simpler geometry). Given a point on the hypersphere representing a trajectory, we also provide a “reverse mapping” algorithm to uniquely (under certain assumptions) reconstruct the subspace that corresponds to this point. Finally, by using existing algorithms for recursive Fr´echet mean and exact principal geodesic analysis on the hypersphere, we present several experiments on synthetic and real (vision and medical) data sets showing how group testing on such diversely sampled longitudinal data is possible by analyzing the reconstructed data in the subspace spanned by the first few principal components.

[1]  Gunnar Rätsch,et al.  Kernel PCA and De-Noising in Feature Spaces , 1998, NIPS.

[2]  Anuj Srivastava,et al.  Statistical analysis of trajectories on Riemannian manifolds: Bird migration, hurricane tracking and video surveillance , 2014, 1405.0803.

[3]  Fatih Murat Porikli,et al.  Covariance Tracking using Model Update Based on Lie Algebra , 2006, 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'06).

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

[5]  Fatih Murat Porikli,et al.  Region Covariance: A Fast Descriptor for Detection and Classification , 2006, ECCV.

[6]  Bill Triggs,et al.  Histograms of oriented gradients for human detection , 2005, 2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'05).

[7]  Baba C. Vemuri,et al.  Recursive Karcher Expectation Estimators And Geometric Law of Large Numbers , 2013, AISTATS.

[8]  René Vidal,et al.  Distances on Spaces of High-Dimensional Linear Stochastic Processes: A Survey , 2014 .

[9]  P. Cochat,et al.  Et al , 2008, Archives de pediatrie : organe officiel de la Societe francaise de pediatrie.

[10]  Daniel D. Lee,et al.  Grassmann discriminant analysis: a unifying view on subspace-based learning , 2008, ICML '08.

[11]  Aarnout Brombacher,et al.  Probability... , 2009, Qual. Reliab. Eng. Int..

[12]  Xavier Pennec,et al.  Intrinsic Statistics on Riemannian Manifolds: Basic Tools for Geometric Measurements , 2006, Journal of Mathematical Imaging and Vision.

[13]  Stefano Soatto,et al.  Recognition of human gaits , 2001, Proceedings of the 2001 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. CVPR 2001.

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

[15]  Hamid Laga,et al.  Covariance Descriptors for 3D Shape Matching and Retrieval , 2014, 2014 IEEE Conference on Computer Vision and Pattern Recognition.

[16]  Takeo Kanade,et al.  An Iterative Image Registration Technique with an Application to Stereo Vision , 1981, IJCAI.

[17]  Rudrasis Chakraborty,et al.  A Nonlinear Regression Technique for Manifold Valued Data with Applications to Medical Image Analysis , 2016, 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR).

[18]  Vikas Singh,et al.  Riemannian Nonlinear Mixed Effects Models: Analyzing Longitudinal Deformations in Neuroimaging , 2017, 2017 IEEE Conference on Computer Vision and Pattern Recognition (CVPR).

[19]  Rama Chellappa,et al.  Statistical Computations on Grassmann and Stiefel Manifolds for Image and Video-Based Recognition , 2011, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[20]  Paul A. Yushkevich,et al.  Deformable registration of diffusion tensor MR images with explicit orientation optimization , 2006, Medical Image Anal..

[21]  Rudrasis Chakraborty,et al.  An Efficient Exact-PGA Algorithm for Constant Curvature Manifolds , 2016, 2016 IEEE Conference on Computer Vision and Pattern Recognition (CVPR).

[22]  Shuicheng Yan,et al.  Pursuing Informative Projection on Grassmann Manifold , 2006, 2006 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR'06).

[23]  Essa Yacoub,et al.  The WU-Minn Human Connectome Project: An overview , 2013, NeuroImage.

[24]  Anand Rangarajan,et al.  A Riemannian Framework for Matching Point Clouds Represented by the Schrödinger Distance Transform , 2014, 2014 IEEE Conference on Computer Vision and Pattern Recognition.

[25]  J. Ross Beveridge,et al.  Grassmann Registration Manifolds for Face Recognition , 2008, ECCV.

[26]  B. Afsari Riemannian Lp center of mass: existence, uniqueness, and convexity , 2011 .

[27]  U. Grenander Elements of Pattern Theory , 1996 .

[28]  N. J. A. Sloane,et al.  Packing Lines, Planes, etc.: Packings in Grassmannian Spaces , 1996, Exp. Math..

[29]  Anna Varentsova,et al.  Development of a high angular resolution diffusion imaging human brain template , 2014, NeuroImage.

[30]  St. Thomas,et al.  Linear Subspace and Manifold Learning via Extrinsic Geometry , 2015 .

[31]  Marc Niethammer,et al.  Group Testing for Longitudinal Data , 2015, IPMI.

[32]  Guido Gerig,et al.  Toward a Comprehensive Framework for the Spatiotemporal Statistical Analysis of Longitudinal Shape Data , 2012, International Journal of Computer Vision.

[33]  Baba C. Vemuri,et al.  A Novel Dynamic System in the Space of SPD Matrices with Applications to Appearance Tracking , 2013, SIAM J. Imaging Sci..

[34]  Sigurdur Helgason Elementary differential geometry , 2001 .

[35]  P. Thomas Fletcher,et al.  Geodesic Regression and the Theory of Least Squares on Riemannian Manifolds , 2012, International Journal of Computer Vision.

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

[37]  N. Ayache,et al.  Log‐Euclidean metrics for fast and simple calculus on diffusion tensors , 2006, Magnetic resonance in medicine.

[38]  P. Thomas Fletcher,et al.  Sasaki metrics for analysis of longitudinal data on manifolds , 2012, 2012 IEEE Conference on Computer Vision and Pattern Recognition.

[39]  Anuj Srivastava,et al.  Bayesian and geometric subspace tracking , 2004, Advances in Applied Probability.

[40]  R. Christensen Plane Answers to Complex Questions: The Theory of Linear Models. , 1997 .

[41]  Moo K. Chung,et al.  Multivariate General Linear Models (MGLM) on Riemannian Manifolds with Applications to Statistical Analysis of Diffusion Weighted Images , 2014, 2014 IEEE Conference on Computer Vision and Pattern Recognition.

[42]  James W. Davis,et al.  A recursive filter for linear systems on Riemannian manifolds , 2008, 2008 IEEE Conference on Computer Vision and Pattern Recognition.

[43]  Rudrasis Chakraborty,et al.  An efficient recursive estimator of the Fréchet mean on a hypersphere with applications to Medical Image Analysis , 2015 .

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

[45]  Rudrasis Chakraborty,et al.  Statistics on the (compact) Stiefel manifold: Theory and Applications , 2017, ArXiv.

[46]  Ivan Laptev,et al.  On Space-Time Interest Points , 2005, International Journal of Computer Vision.