Multiscale Representations for Manifold-Valued Data

We describe multiscale representations for data observed on equispaced grids and taking values in manifolds such as the sphere $S^2$, the special orthogonal group $SO(3)$, the positive definite matrices $SPD(n)$, and the Grassmann manifolds $G(n,k)$. The representations are based on the deployment of Deslauriers--Dubuc and average-interpolating pyramids "in the tangent plane" of such manifolds, using the $Exp$ and $Log$ maps of those manifolds. The representations provide "wavelet coefficients" which can be thresholded, quantized, and scaled in much the same way as traditional wavelet coefficients. Tasks such as compression, noise removal, contrast enhancement, and stochastic simulation are facilitated by this representation. The approach applies to general manifolds but is particularly suited to the manifolds we consider, i.e., Riemannian symmetric spaces, such as $S^{n-1}$, $SO(n)$, $G(n,k)$, where the $Exp$ and $Log$ maps are effectively computable. Applications to manifold-valued data sources of a geometric nature (motion, orientation, diffusion) seem particularly immediate. A software toolbox, SymmLab, can reproduce the results discussed in this paper

A global approach to the refinement of manifold data

A refinement of manifold data is a computational process, which produces a denser set of discrete data from a given one. Such refinements are closely related to multiresolution representations of manifold data by pyramid transforms, and approximation of manifold-valued functions by repeated refinements schemes. Most refinement methods compute each refined element separately, independently of the computations of the other elements. Here we propose a global method which computes all the refined elements simultaneously, using geodesic averages. We analyse repeated refinements schemes based on this global approach, and derive conditions guaranteeing strong convergence.Comment: arXiv admin note: text overlap with arXiv:1407.836

H\"older Regularity of Geometric Subdivision Schemes

We present a framework for analyzing non-linear $\mathbb{R}^d$-valued subdivision schemes which are geometric in the sense that they commute with similarities in $\mathbb{R}^d$. It admits to establish $C^{1,\alpha}$-regularity for arbitrary schemes of this type, and $C^{2,\alpha}$-regularity for an important subset thereof, which includes all real-valued schemes. Our results are constructive in the sense that they can be verified explicitly for any scheme and any given set of initial data by a universal procedure. This procedure can be executed automatically and rigorously by a computer when using interval arithmetics.Comment: 31 pages, 1 figur

Quasi-interpolation in Riemannian manifolds

We consider quasi-interpolation operators for functions assuming their values in a Riemannian manifold. We construct such operators from corresponding linear quasi-interpolation operators by replacing affine averages with the Riemannian centre of mass. As a main result, we show that the approximation rate of such a nonlinear operator is the same as for the linear operator it has been derived from. In order to formulate this result in an intrinsic way, we use the Sasaki metric to compare the derivatives of the function to be approximated with the derivatives of the nonlinear approximant. Numerical experiments confirm our theoretical finding

Subdivision schemes with general dilation in the geometric and nonlinear setting

AbstractWe establish results on convergence and smoothness of subdivision rules operating on manifold-valued data which are based on a general dilation matrix. In particular we cover irregular combinatorics. For the regular grid case results are not restricted to isotropic dilation matrices. The nature of the results is that intrinsic subdivision rules which operate on geometric data inherit smoothness properties of their linear counterparts

Point-Normal Subdivision Curves and Surfaces

This paper proposes to generalize linear subdivision schemes to nonlinear subdivision schemes for curve and surface modeling by refining vertex positions together with refinement of unit control normals at the vertices. For each round of subdivision, new control normals are obtained by projections of linearly subdivided normals onto unit circle or sphere while new vertex positions are obtained by updating linearly subdivided vertices along the directions of the newly subdivided normals. Particularly, the new position of each linearly subdivided vertex is computed by weighted averages of end points of circular or helical arcs that interpolate the positions and normals at the old vertices at one ends and the newly subdivided normal at the other ends. The main features of the proposed subdivision schemes are three folds: (1) The point-normal (PN) subdivision schemes can reproduce circles, circular cylinders and spheres using control points and control normals; (2) PN subdivision schemes generalized from convergent linear subdivision schemes converge and can have the same smoothness orders as the linear schemes; (3) PN $C^2$ subdivision schemes generalizing linear subdivision schemes that generate $C^2$ subdivision surfaces with flat extraordinary points can generate visually $C^2$ subdivision surfaces with non-flat extraordinary points. Experimental examples have been given to show the effectiveness of the proposed techniques for curve and surface modeling.Comment: 30 pages, 17 figures, 22.5M

Single basepoint subdivision schemes for manifold-valued data: Timesymmetry without space-symmetry.

Abstract. This paper establishes smoothness results for a class of nonlinear subdivision schemes, known as the single basepoint manifold-valued subdivision schemes, which shows up in the construction of waveletlike transform for manifold-valued data. This class includes the (single basepoint) Log-Exp subdivision scheme as a special case. In these schemes, the exponential map is replaced by a so-called retraction map f from the tangent bundle of a manifold to the manifold. It is known that any choice of retraction map yields a C 2 scheme, provided the underlying linear scheme is C 2 (this is called "C 2 equivalence"). But when the underlying linear scheme is C 3 , Navayazdani and Yu have shown that to guarantee C 3 equivalence, a certain tensor P f associated to f must vanish. They also show that P f vanishes when the underlying manifold is a symmetric space and f is the exponential map. In the present paper, a geometric interpretation of the tensor P f is given. Associated to the retraction map f is a torsion-free affine connection, which in turn defines an exponential map. The condition P f = 0 is shown to be equivalent to the condition that f agrees with the exponential map of the connection up to the 3rd order. In particular, when f is the exponential map of a connection, one recovers the original connection and P f vanishes. It then follows that the condition P f = 0 is satisfied by a wider class of manifolds than was previously known. Under the additional assumption that the subdivision rule satisfies a time-symmetry, it is shown that the vanishing of P f also guarantees C 4 equivalence. Finally, the analysis in the paper strongly indicates that vanishing curvature of the connection associated to f is a necessary condition for C k equivalence for k ≥ 5