Statistical analysis of trajectories on Riemannian manifolds: Bird migration, hurricane tracking and video surveillance

Jingyong Su,Sebastian Kurtek,Eric Klassen,Anuj Srivastava
DOI: https://doi.org/10.1214/13-AOAS701
2014-05-05
Abstract:We consider the statistical analysis of trajectories on Riemannian manifolds that are observed under arbitrary temporal evolutions. Past methods rely on cross-sectional analysis, with the given temporal registration, and consequently may lose the mean structure and artificially inflate observed variances. We introduce a quantity that provides both a cost function for temporal registration and a proper distance for comparison of trajectories. This distance is used to define statistical summaries, such as sample means and covariances, of synchronized trajectories and "Gaussian-type" models to capture their variability at discrete times. It is invariant to identical time-warpings (or temporal reparameterizations) of trajectories. This is based on a novel mathematical representation of trajectories, termed transported square-root vector field (TSRVF), and the $\mathbb{L}^2$ norm on the space of TSRVFs. We illustrate this framework using three representative manifolds---$\mathbb{S}^2$, $\mathrm {SE}(2)$ and shape space of planar contours---involving both simulated and real data. In particular, we demonstrate: (1) improvements in mean structures and significant reductions in cross-sectional variances using real data sets, (2) statistical modeling for capturing variability in aligned trajectories, and (3) evaluating random trajectories under these models. Experimental results concern bird migration, hurricane tracking and video surveillance.
Applications
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the inaccuracy of statistical summaries when conducting statistical analysis of trajectories on Riemannian manifolds due to the arbitrariness of time evolution (i.e., time re - parameterization or time warping). Specifically: 1. **Problems with cross - sectional mean and variance**: Past methods relied on cross - sectional analysis, which may lead to the loss of the mean structure and artificially inflate the observed variance. For example, in Figure 1(a), if the observations of the trajectory are generated at random discrete observation times, the calculated cross - sectional mean and variance will be distorted. The mean curve is quite different from the original curve, and the variance is also entirely caused by the randomness of the observation times. 2. **The influence of time warping**: Trajectories of different individuals may have the same spatial path but different time evolution patterns. If this time variability is not considered and the cross - sectional mean and variance are directly used, the results will be inaccurate. For example, trajectory data in bird migration, hurricane path tracking, and video surveillance all have such problems. To solve these problems, the author introduced a new mathematical representation method - **Transport Square - Root Vector Field (TSRVF)** and the corresponding \( L^2 \) norm. Through this method, a new distance \( d_s \) can be defined, which can be used as both a cost function for time registration and a distance metric for comparing trajectories, and this distance is invariant to the same time warping. ### Key points of the solution - **TSRVF representation**: Convert the trajectory into the TSRVF form, and then compare these TSRVF through the \( L^2 \) norm. - **New distance \( d_s \)**: Define a new distance \( d_s \), which is the shortest \( d_h \) distance between two equivalence classes, thus ensuring invariance to time warping. - **Optimized time alignment**: Find the optimal time warping function \( \gamma^* \) through the dynamic programming algorithm to make the matching between trajectories more accurate. Through these methods, the paper achieved: - Improving the mean structure and significantly reducing the cross - sectional variance. - Capturing the variability of aligned trajectories. - Evaluating the model performance under random trajectories. These improvements have been verified in the actual data sets of bird migration, hurricane path tracking, and video surveillance.