# Jumping Manifolds: Geometry Aware Dense Non-Rigid Structure from Motion

Given dense image feature correspondences of a non-rigidly moving object across multiple frames, this paper proposes an algorithm to estimate its 3D shape for each frame. To solve this problem accurately, the recent state-of-the-art algorithm reduces this task to set of local linear subspace reconstruction and clustering problem using Grassmann manifold representation kumar2018scalable. Unfortunately, their method missed on some of the critical issues associated with the modeling of surface deformations, for e.g., the dependence of a local surface deformation on its neighbors. Furthermore, their representation to group high dimensional data points inevitably introduce the drawbacks of categorizing samples on the high-dimensional Grassmann manifold huang2015projection, harandi2014manifold. Hence, to deal with such limitations with kumar2018scalable, we propose an algorithm that jointly exploits the benefit of high-dimensional Grassmann manifold to perform reconstruction, and its equivalent lower-dimensional representation to infer suitable clusters. To accomplish this, we project each Grassmannians onto a lower-dimensional Grassmann manifold which preserves and respects the deformation of the structure w.r.t its neighbors. These Grassmann points in the lower-dimension then act as a representative for the selection of high-dimensional Grassmann samples to perform each local reconstruction. In practice, our algorithm provides a geometrically efficient way to solve dense NRSfM by switching between manifolds based on its benefit and usage. Experimental results show that the proposed algorithm is very effective in handling noise with reconstruction accuracy as good as or better than the competing methods.

Comments

There are no comments yet.

## Authors

• 7 publications
• ### Scalable Dense Non-rigid Structure-from-Motion: A Grassmannian Perspective

This paper addresses the task of dense non-rigid structure from motion (...
03/01/2018 ∙ by Suryansh Kumar, et al. ∙ 0

read it

• ### Manifold Approximation by Moving Least-Squares Projection (MMLS)

In order to avoid the curse of dimensionality, frequently encountered in...
06/22/2016 ∙ by Barak Sober, et al. ∙ 0

read it

• ### Hierarchic Neighbors Embedding

Manifold learning now plays a very important role in machine learning an...
09/16/2019 ∙ by Shenglan Liu, et al. ∙ 14

read it

• ### Efficient moving point handling for incremental 3D manifold reconstruction

As incremental Structure from Motion algorithms become effective, a good...
07/20/2015 ∙ by Andrea Romanoni, et al. ∙ 0

read it

• ### Superpixel Soup: Monocular Dense 3D Reconstruction of a Complex Dynamic Scene

This work addresses the task of dense 3D reconstruction of a complex dyn...
11/19/2019 ∙ by Suryansh Kumar, et al. ∙ 19

read it

• ### A Theoretical Analysis of Joint Manifolds

The emergence of low-cost sensor architectures for diverse modalities ha...
01/07/2009 ∙ by Mark A. Davenport, et al. ∙ 0

read it

• ### 3D Morphology Prediction of Progressive Spinal Deformities from Probabilistic Modeling of Discriminant Manifolds

We introduce a novel approach for predicting the progression of adolesce...
01/17/2017 ∙ by Samuel Kadoury, et al. ∙ 0

read it

##### This week in AI

Get the week's most popular data science and artificial intelligence research sent straight to your inbox every Saturday.

## 1 Introduction

Non-rigid Structure-from-Motion (NRSfM), a problem where the task is to recover the three-dimensional structure of a deforming object from a set of image feature correspondences across frames. Any solution to this problem depends on the proper modeling of structure and an efficient estimation of motion , where denotes some structure manifold and denotes special Euclidean group which is a differentiable manifold [20]. Though, after Bregler et al. factorization framework to NRSfM [48], motion estimations are mostly relaxed to rotation estimation . Even after such relaxation, the problem still remains unsolved for any arbitrary motion. The main difficulty in NRSfM comes from the fact that both the camera and the object are moving and, along with it the object themselves are deforming, hence, it becomes difficult to distinguish camera motion from object motion using only image data. Despite such difficulties, many efficient and reliable solutions based on the priors are proposed to solve NRSfM. A reliable solution to this problem is important as it covers a wide range of applications from medical industry to the entertainment industry and many more.

To solve NRSfM, the algorithms proposed in the past can broadly be divided into two major classes 1) sparse NRSfM and 2) dense NRSfM. This classification is based on the number of feature points that the algorithm can efficiently process to model the deformation of the object. Although many reliable solution to this problem exists for sparse NRSfM [18, 37, 35, 3, 47, 42, 40, 26, 29], very few work have been done towards solving the dense NRSfM reliably and efficiently [23, 17, 34, 4]. Also, the existing solutions to dense NRSfM are computationally expensive and are mostly constrained to analyze the global deformation of the non-rigid shape [23, 17]. The basis for this gradual progress in dense NRSfM is perhaps due to its dependence on per pixel reliable correspondences across frames, and the absence of a resilient structure modeling framework to capture the local non-linearities. One may argue on the efficient motion estimation, however, from image correspondences, we can only estimate relative motion and reliable algorithms with solid theory exists to perform this task well [18, 40]

. Also, with the recent progress in deep learning algorithms, per pixel correspondences can be achieved with a remarkable accuracy

[46], which leaves structure modeling as a potential gray area in dense NRSfM to focus.

Very recently, Kumar et al. [34] has exploited the Grassmann manifold to model non-rigid surfaces in dense NRSfM. The key insight in their work is; even though the overall complexity of the deforming shape is high, each local deformation may be less complex [13, 14, 15, 16]. Using this idea, they proposed a union of local linear subspace approach to solve dense NRSfM problem. Nevertheless, their work overlooked on some of the intrinsic issues associated with the modeling of non-rigidly deforming surface. Firstly, their method represents each local linear subspace independently via a high-dimensional Grassmannian representation. Now, such representation may help reconstruct complex 3D deformation but can lead to wrong clustering, and it’s very important in joint reconstruction and clustering framework to have suitable clustering of subspaces, else reconstruction may suffer. Secondly, their approach to represent local non-linear deformation completely ignored the neighboring surfaces, which may result in an inefficient representation of the Grassmannians in the trajectory space. Thirdly, the representation of Grassmannians in the shape space adopted by [34] results in irredeemable discontinuity of the trajectories (see Fig.(2)). Hence, temporal representation of the set of shapes using Grassmannians seems not an extremely beneficial choice for modeling dense NRSfM on Grassmannian manifold111Purpose behind NRSfM is not the same as activity/action recognition. See supplementary material for a detail discussion on this.. Lastly, although the dense NRSfM algorithm proposed in [34] works better and faster than the previous methods, it depends on several manual parameters which are inadmissible for practical applications.

This paper introduces an algorithm that overcomes the aforementioned limitations with Kumar et al. method [34]. The main point we are trying to make is that; reconstruction and grouping of subspace on the same high dimensional Grassmann manifold seem an unreasonable choice. Even recent research in the Riemannian geometry has shown that the low-dimensional representation of the corresponding high dimensional Grassmann manifold is more favorable for grouping Grassmannians [32, 31]. So, we formulate dense NRSfM in a way that it takes advantage of both high and low dimensional representation of Grassmannians i.e., perform reconstruction in the original high-dimension and cluster subspace in its lower-dimension representation.

We devise an unsupervised approach to efficiently represent the high-dimensional non-rigid surface on a lower dimensional Grassmann manifold. These low-dimensional Grassmannians are represented in such a way that it preserves the local structure of the surface deformation in accordance with its neighboring surfaces when projected. Now, these low-dimensional Grassmannians serves as a potential representative for its high-dimensional Grassmannians for suitable grouping, which subsequently help improve the reconstruction and representation of the Grassmannians on the high-dimensional Grassmann manifold, hence, the term Jumping Manifolds (MoJu). Further, we drop the temporal grouping of shapes using Grassmannians to discourage the discontinuity of trajectories (see Fig.(2)).

In essence, our work is inspired by [34] and is oriented towards settling its important limitations. Moreover, in contrast to [34], we capture the notion of dependent local subspace in a union of subspace algorithm [39] via Grassmannian modeling. The algorithm we proposed is an attempt to supply a more efficient, reliable and practical solution to this problem. Our formulation gives an efficient framework for modeling dense NRSfM on the Grassmann manifold than [34]. Experimental results show that our method is as accurate as other algorithms and is numerically more efficient in handling noise. The main contributions of our work are as follows:

• [noitemsep]

• An efficient framework for modeling non-rigidly deforming surface that exploits the advantage of Grassmann manifold representation of different dimensions based on its geometry.

• A formulation that encapsulates the local non-linearity of the deforming surface w.r.t its neighbors to enable the proper inference and representation of local linear subspaces.

• An iterative solution to the proposed cost function based on ADMM [8], which is simple to implement and provide results as good as the best available methods. Additionally, it helps improve the 3D reconstruction substantially, in the case of noisy trajectories.

## 2 Related Work

A working solution to NRSfM was first introduced in the seminal work by Bregler et al[9] which is an extension to Tomsai-Kannade rigid factorization method[48]. Although the problem still remains unsolved for arbitrary deformations, many profound works have been done to achieve a reliable solution to this problem under some or the other prior assumptions about the object [18, 3, 43, 40, 50, 37, 35, 23]. Since the literature on this topic is very extensive, we review the works that are of close relevance to the dense NRSfM methods under classical setting222By classical setting, we mean without using RGB-D or 3D template..

Earlier attempts to solve dense NRSfM used piecewise reconstruction of the shape parts which were further processed via a stitching step to get a global 3D shape [12, 45]. To our knowledge, Garg et al. variational approach [24] was the first to propose and demonstrate per pixel dense NRSfM algorithm without any 3D template prior. This method introduced a discrete total variational constraint with trace norm constraint on the global shape, which resulted in a biconvex optimization problem. Despite the algorithm outstanding results, it’s computationally very expensive and needs a GPU to provide the solution.

In contrast, Dai et al. extended his simple prior free approach to solve dense NRSfM problem [18, 17]. The algorithm proposed a spatial-temporal formulation to solve the problem. The author revisits the temporal smoothness term from [18] and integrate it with a spatial smoothness term using the Laplacian of the non-rigid shape. The resultant optimization leads to a series of least squares to be minimized which makes it extremely slow to process. Recently, Kumar et al. modeled this problem on the Grassmann manifold [34]. Their work extended the spatiotemporal multi-body framework to solve dense NRSfM [37]. The algorithm demonstrated that such an approach is more efficient, faster and accurate than all the other recent approaches to solve dense NRSfM task [24, 17, 4].

Consecutive frame-based approach has recently shown some promising results to solve dense 3D reconstruction of a general dynamic scene including non-rigid object [36, 44]. Nevertheless, motion segmentation, triangulation, as rigid as possible constraint and scale consistency quite often breaks down for the deforming object over frames. Therefore, dense NRSfM becomes extremely challenging for such algorithms. Not long ago, Gallardo et al. combined shading, motion and generic physical deformation to model dense NRSfM [21].

## 3 Preliminaries

In this paper, , denotes the Frobenius norm and nuclear norm respectively. represent the notion of norm on the Grassmann manifold. Single angle bracket denotes the Euclidean inner product. For ease of understanding and completeness, in this section, we briefly review few important definitions related to the Grassmann manifold. Firstly, a manifold is a topological space that is locally similar to the Euclidean space. Out of several manifolds, the Grassmann manifold is a topologically rich non-linear manifold, each point of which represent the set of all right invariant subspace of the Euclidean space [19, 1, 34].

###### Definition 1.

The Grassmann manifold, denoted by , consists of all the linear -dimensional subspace embedded in a ‘’ dimensional Euclidean space such that [Absil et al., 2009] [1].

A point ‘’ on the Grassmann manifold can be represented by matrix whose columns are composed of orthonormal basis. The space of such matrices with orthonormal columns is a Riemannian manifold such that , where is a identity matrix.

###### Definition 2.

Grassmann manifold can be embedded into the space of symmetric matrices via mapping , where is a Grassmann point [28, 30]. Given two Grassmann points and , then the distance between them can be measured using the projection metric [28].

In the past, these two properties of Grassmann manifold has been used in many computer vision applications

[28, 10, 34]. Second definition is very important as it allows to measure the distance on the Grassmann manifold, hence, ) forms a metric space. We used these properties in the construction of our formulation. For comprehensive details on this topic readers may refer to [28].

## 4 Problem Formulation

Let ‘’ be the total number of feature points tracked across ‘’ frames. Concatenating these 2D coordinates of each feature points for all frames across the columns of a matrix gives ‘ matrix. This matrix is popularly known as measurement matrix [48]. Our goal is, given the image measurement matrix, estimate the camera motion and 3D coordinates of every 2D feature points across all frames.

We start our formulation with the classical representation to NRSfM i.e. , where, is a block diagonal rotation matrix with each block as a orthographic rotation matrix, and as the 3D structure matrix. With such a representation, the entire problem simplifies to the estimation of correct rotation matrix ‘’ and structure matrix ‘’ such that the above relation holds. Following the assumption of the previous work [34], we estimate the rotation using Intersection method [18]. As a result, the task reduces to composing of an efficient algorithm that correctly models the surface deformations and provide better reconstruction results. Recent algorithms in NRSfM have demonstrated that clustering benefits reconstruction and vice-versa, however, the existing framework to employ this idea is not scalable to millions of points. To establish this idea for dense NRSfM, Kumar et al. [34] used LRR on Grassmannian manifold. Using the similar notions, we model dense deforming surface using Grassmannian representation to provide more reliable and accurate solution.

In the following subsection, we first introduce the Grassmannian representation of the surface and how to project these Grassmannians onto the lower dimension Grassmann manifold by preserving the neighboring information. In the later subsection, we use these representations to formulate the overall cost function for solving dense NRSfM problem.

### 4.1 Grassmannian representation

Let ‘ be a Grassmann point representing the local linear subspace spanned by set of columns of ‘’. Using this notion, we decompose the entire trajectories of the structure into a set of ‘’ Grassmannians . Now, such a representation treats each subspace independently and therefore, its low-dimensional linear representation may not be suitable to capture the surface dependent non-linearity. To properly represent Grassmannian which respects the neighboring non-linearity in low-dimension, we introduce a different strategy to model non-rigid surface in low-dimension. For now, let be a matrix that maps ‘ to ‘ such that . Mathematically,

 ϕ\bfi=Δ%$T$Φ%$i$ (1)

Its quite easy to examine that

is not a orthogonal matrix and, therefore, does not qualifies as a potential point on a Grassmann manifold. However, by performing a orthogonal-triangular (QR) decomposition of

, we estimate the new representative of on the Grassmann manifold of ‘’ dimension.

 Θ\bfi\bfU\bfi=qr(ϕ%$i$)=Δ%$T$Φ%$i$ (2)

Here, is a function that returns the QR decomposition of the matrix. The is an orthogonal matrix and is the upper triangular matrix333Note: The value of , Use in MATLAB to get a square matrix (). Using Eq.(2), we represent the equivalence of in low dimension as

 Θ\bfi=Δ\bfT(Φ\bfi\bfU−1%$i$) (3) Θ\bfi=Δ\bfTΩ\bfi

where, . The key-point to note is that both and has the same column space. In principle such a representation is useful however, it does not serve the purpose of preserving the non-linearity w.r.t its neighbors. In order to encapsulate the local dependencies (see Fig.(3), Fig.(4)), we further constrain our representation as:

 \bfE(Δ)=minimizeΔ\bfK∑(\bfi,\bfj)\bfw\bfi\bfj12∥Π(Θ\bfi)−Π(Θ\bfj)∥\bfF2 (4)

The parameter ‘’ accommodate the similarity knowledge between the two Grassmannians. Using the Definition(2) and Eq.(3), we further simplify Eq.(4) as

 \bfE(Δ)≡minimizeΔ\bfK∑(\bfi,\bfj)\bfw\bfi\bfj12∥Δ\bfTΩ\bfiΩ\bfT\bfiΔ−Δ\bfTΩ\bfjΩ\bfT\bfj\bfjΔ∥2\bfF (5) \bfE(Δ)≡minimizeΔ\bfK∑(\bfi,\bfj)\bfw\bfi\bfj12∥Δ\bfT(Ω\bfiΩ\bfT\bfi−Ω\bfjΩ\bfT\bfj\bfj)Δ∥2\bfF \bfE(Δ)≡minimizeΔ\bfK∑(\bfi,\bfj)\bfw\bfi\bfj12∥Δ\bfT(Λ\bfi\bfj)Δ∥2\bfF

where, . The parameter ‘’ (similarity graph) is set as with as the projection metric (see Definition (2)). Eq.(5) is an unconstrained optimization problem and its solution may provide a trivial solution. To estimate the useful solution, we further constrain the problem. Using Grassmann point ‘’ and its neighbors, expand Eq.(5). By performing some simple algebraic manipulation, Eq.(5) reduces to

 trace(Δ\bfT(\bfK∑\bfi=1λ\bfi\bfiΩ\bfi\bfiΩ\bfT\bfi\bfi)Δ) (6)

where, . Constraining the value of Eq.(6) to 1 provides the overall optimization for an efficient representation of the local non-rigid surface on the Grassmann manifold.

 \bfE(Δ)≡minimizeΔ\bfK∑(\bfi,\bfj)\bfw\bfi\bfj12∥Δ\bfT(Λ\bfi\bfj)Δ∥2\bfF (7) subject to: trace(Δ\bfT(\bfK∑\bfi=1λ\bfi\bfiΩ\bfi\bfiΩ\bfT\bfi\bfi)Δ)=1

Its easy to verify that the matrix and are symmetric and positive semi-definite, and therefore, the above optimization can be solved as a generalized eigen value problem —refer supplementary material for details.

### 4.2 Dense NRSfM formulation

To solve the dense non-rigid structure from motion with the representation formulated in the previous sub-section §4.1, we propose to jointly optimize the objective function over the 3D structure and its local group representation. In order to build the overall objective function, we introduce each constraint equation one by one for clear understanding of our overall cost function.

 E\bfp(\bfS)=minimize% %$S$ 12∥\bfW−\bfR\bfS∥2\bfF (8)

The first term constrain the 3D structure such that it satisfies the re-projection error.

 E\bfs(\bfS♯)=minimize\bfS♯ ∥\bfS♯∥∗ (9)

The second term caters the global assumption about the non-rigid object; that is the overall shape matrix is low-rank. To establish this assumption, we perform rank minimization of the shape matrix. Although the rank minimization of a matrix is NP-hard, it’s relaxed to nuclear norm minimization to find an approximate solution. This term mainly penalizes the total number of independent shape required to represent the shape. The choice of minimizing instead to is inspired from Dai et al.’s work [18]. Since the dense deforming shape is composed of several local linear low-dimensional subspace, the global constraint (Eq.(9)) may not reflect their local dependency. Therefore, in order to introduce the local subspace constraint on the shape, we use the notion of self-expressiveness on the non-linear Grassmann manifold space.

 %minimize\bfE,\bfC,%$S$♯ ∥\bfE∥2G+β2∥\bfS♯∥∗+β3∥\bfC∥∗ (10) subject to: \bfS♯=f(\bfS),\bfS=\bfS\bfC+\bfE

Here, we define and as the coefficient matrix. We know from the literature that the Grassmann manifold is isometrically equivalent to the symmetric idempotent matrix [11]. So, we embed the Grassmann manifold into symmetric matrix manifold to define the self-expressiveness. Let be the set of Grassmannians on a low-dimensional Grassmann manifold. The elements of are the projection of high-dimensional Grassmannian representation of the columns of ‘’ matrix. Let be its embedding onto symmetric matrix manifold. Using such embedding techniques we re-write Eq.(10) as

 minimize\bfE,~\bfC,\bfS♯ ∥\bfE∥2\bfF+β2∥\bfS♯∥∗+β3∥~%$C$∥∗ (11) subject to: \bfS♯=f(\bfS),\raisebox0.0pt$χ$=\raisebox0.0pt$χ$\raisebox0.0pt$χ$~\bfC+\bfE

where, and

denotes the coefficient matrix of Grassmannians and structure tensor respectively, with

as the total number of Grassmannians. Generally, , which makes such representation scalable.

The third term we introduce is composed of few constraint functions that provides a way to group Grassmannians and recover 3D shape simultaneously. Let

be an ordering vector that contains the index of columns of

. Our function definition is of the form . Using it, we define the function , , and as follows:

 {(ξ,fg(P,\bfS)):order {\bfS\bfi}\bfK\bfi=1 columns of \bfS of% using P, (12) ξ:={Φ\bfi}\bfK\bfi=1where,[ Φ\bfi,Σ\bfi,ξvi]=svds(\bfS\bfi,% \bfp)}
 {(~ξ,fh(Δ,ξ)):~ξ={Θ\bfi}\bfK\bfK%$i$=1,Θ\bfi\bfi=Δ\bfT\bfT(Φ\bfi% \bfU−1\bfi), (13) where, Δ=solution to the minimization of Eq% .(???)}
 {(P,fp(~ξ,~% \bfC,P\bfo):P=spectral\_% clustering(~ξ,~\bfC,P\bfo)} (14)
 {(\bfS,fs(ξ,Σ,ξ\bfv)):\bfS\bfi=[ξ\bfi Σ\bfi ξ\bfv\bfi],where Σ\bfi∈R\bfp×\bfp} (15)

Intuitively, the first function () uses the ordering vector to refine the grouping of the trajectories for suitable Grassmannian representation. The second function () projects the Grassmannians to a lower dimension in accordance with the neighbors using Eq.(7). The third function () uses the projected Grassmannians to assign proper labeling to the Grassmann points and update the given ordering vector

using spectral clustering. The fourth function (

) uses the group of trajectories to reconstruct back the set of local surface. ,

are the singular values and right singular vector matrices in the high-dimension.

Objective Function: Combining all the above terms and constraints provides our overall cost function.

 %minimize\bfE,~\bfC,\bfS,%$S$♯12∥\bfW−\bfR\bfS∥2\bfF+β1∥\bfE∥2\bfF+β2∥\bfS♯∥∗+β3∥~\bfC∥∗ (16) subject to: \bfS♯=f(\bfS),\raisebox0.0pt$χ$=\raisebox0.0pt$χ$\raisebox0.0pt$χ$~\bfC+\bfE, ξ=fg(P,\bfS),~ξ=fh(Δ,ξ), \bfS=fs(ξ,Σ,ξ\bfv),P=fp(~ξ,~\bfC,P%$o$)

where vector contains the initial ordering of the columns of ‘’ and ‘’. The function provides the ordering index to rearrange the columns of ‘’ matrix to be consistent with ‘’ matrix. This is important because, grouping the set of columns of ‘’ over iteration, disturbs its initial arrangements.

## 5 Solution

The optimization proposed in Eq.(16) is a coupled optimization problem. Several methods of Bi-level optimization can be used to solve such minimization problem [6, 27]. Nevertheless, we propose ADMM [8] based solution due to its application in many non-convex optimization problems. The key point to note is that one of our constraint is composed of separate optimization problem i.e., the solution to Eq.(7), and therefore, we cannot directly embed the constraint to the main objective function. Instead, we only introduce two Lagrange multiplier to concatenate a couple of constraints back to the original objective function. The remaining constraints are enforced over iteration. To decouple the variable from , we introduce auxiliary variable . We apply these operations to our optimization problem to get the following Augmented Lagrangian form:

 %minimize\bfZ,~\bfC,\bfS,%$S$♯12∥\bfW−\bfR\bfS∥2\bfF+β1∥\raisebox0.0pt$χ$−\raisebox0.0pt$χ$\raisebox0.0pt$χ$~\bfC∥2\bfF+β2∥\bf% S♯∥∗+ (17) ρ2∥\bfS♯−f(\bfS)∥2\bfF\bfF++β3∥\bfZ∥∗+ ρ2∥~\bfC−\bfZ∥2\bfF+ subject to: ξ=fg(P,% \bfS),~ξ=fh(Δ,ξ), \bfS=fs(ξ,Σ,ξ\bfv),P=fp(~ξ,~\bfC,P%$o$)

Note that provides the information about the subspace, not the vectorial points. However, we have the chart of the trajectories and its corresponding subspace. Once, we group the trajectories based on , provides new Grassmann sample corresponding to each group. The definition of and is provided in Eq.(7) and Eq.(14) respectively. More generally, the solution to the optimization in Eq.(7

) is obtained by solving it as a generalized eigenvalue problem. To keep the order of columns of ‘

’ matrix consistent with ‘’ matrix provides the ordering index. We provide the implementation details of our method with suitable MATLAB commands in the Algorithm Table (1). For details on the derivation to each sub-problem, kindly refer to the supplementary material.

## 6 Initialization and Evaluation

We performed experiments and evaluation on the available standard benchmark datasets [23, 49, 7]. To keep our evaluations consistent with the previous methods, we compute the mean normalized 3D reconstruction error of the estimated shape ‘’ after convergence as

 \bfe3\bfD=1\bfF\bfF\bfF∑%$i$=1∥\bfS%$i$ est−\bfS\bfi\bfG\bfT∥\bfF∥\bfS\bfi\bfG\bfT∥\bfF (18)

here ‘’ denotes the ground-truth 3D shape matrix.

#### Initialization:

We used Intersection method [18] to estimate the rotation matrix and initialize . The initial grouping of the trajectories or columns of

is done using k-means++ algorithm

[5]. These groups are then used to initialize , and the Grassmann points via subset of singular vectors. To represent the Grassmannians in the lower-dimension, we solve Eq.(7) to initialize and store corresponding singular values. The similarity matrix or graph in Eq.(7) is build using the distance measure between the Grassmannians in the embedding space §4.1.
1. Results on synthetic Face dataset: The synthetic face dataset is composed of four distinct sequence [23] with 28,880 feature points tracked over multiple frames. Each sequence captures the human facial expression with a different range of deformations and camera motion. Sequence 1 and Sequence 2 are 10 frame long video with rotation in the range and respectively. Sequence 3 and Sequence 4 are 99 frame long video that contains high frequencies and low frequencies rotation respectively which captures real human facial deformations. Table (1) shows the statistical results obtained on these sequences using our algorithm. For qualitative results on these sequences kindly refer to the supplementary material.

2. Results on Paper and T-shirt dataset: Varol et al. introduced ‘kinect_paper’ and ‘kinect_tshirt’ datasets to test the performance of NRSfM algorithm under real conditions [49]. This dataset provides sparse SIFT [41] feature tracks and noisy depth information captured from Microsoft Kinect for all the frames. As a result, to get dense 2D feature correspondences of the non-rigid object for all the frames becomes difficult. To circumvent this issue, we used Garg et al. [22] algorithm to estimate the measurement matrix. To keep the numerical comparison consistent with the previous work in dense NRSfM [34], we used the same coordinate range for tracking the features. Numerically, its = , = rectangular window across 193 frames for kinect_paper sequence. For kinect_tshirt sequence, we considered rectangular window of = , = across 313 frames, same as used in Kumar et al. work [34]. Fig.(5) shows the reconstruction results on these sequence with comparative results provided in Table (1).

3. Results on Actor dataset: Beeler et al. [7] introduced Actor dataset for high-quality facial performance capture. This dataset is composed of 346 frames captured from seven cameras with 1,180,232 vertices. The dataset captures the fine details of facial expressions which is extremely useful in the testing of NRSfM algorithms. Nevertheless, for our experiment, we require dense 2D image feature correspondences across all images as input, which we synthesized using ground-truth 3D points and synthetically generated orthographic camera rotations. To maintain the consistency with the previous works in dense NRSfM for performance evaluations, we synthesized two different datasets namely Actor Sequence1 and Actor Sequence2 based on the head movement as described in Ansari et al. work [4]. Fig.(6) shows the dense detailed reconstruction that is achieved using our algorithm. Table (1) clearly indicates the benefit of our approach to reconstruct such complex deformations.

4. Results on Face, Heart, Back dataset: To evaluate the variational approach to dense NRSfM [23] Garg et al. introduced these datasets. Its sequences are composed of monocular video’s captured in a natural environment with varying lighting condition and large displacements. It consists of three different videos with 120, 150 and 80 frames for face sequence, back sequence and heart sequence respectively. Additionally, it provides dense 2D feature track for the same with 28332, 20561, and 68295 features track over the frames for face, back and heart sequence. No ground-truth 3D is available with this dataset for evaluation. Fig.(5) show reconstruction results on back and heart sequence. For more qualitative results on these sequences, kindly refer to the supplementary material.

### 6.1 Algorithmic Analysis

We performed some more experiments to understand the behavior of our algorithm under different input parameters and evaluation setups. In practice these experiment help analyze the practical applicability of our algorithm.
1. Performance over noisy trajectories: We utilized the standard experimental procedure to analyze the behavior of our algorithm under different noise levels. Similar to the work of Kumar et al. [34]

, we added the Gaussian noise to the input trajectories. The standard deviation of the noise are adjusted as

with varying from 0.01 to 0.055. Fig.(7(a)) show the quantitative comparison of our approach with recent algorithms namely DS [17] and SDG [34]. The graph is plotted by taking the average reconstruction error of all the four synthetic face dataset [23]. The procured statistics indicate that our algorithm is more resilient to noise than other competing methods.
2. Performance with change in the number of singular values: The selection of ‘’ in i.e. the number of top singular vectors for Grassmannian representation and its corresponding singular values to perform reconstruction can directly affect the performance of our algorithm. However, it has been observed over several experiments that we need very few singular value and singular vectors to recover dense detailed 3D reconstruction of the deforming object. Fig.(7(b)) show the variation in average 3D reconstruction with the values of ‘’ for synthetic face dataset [23].
3. Processing Time and Convergence: Our algorithm execution time is almost at par or a bit slower than SDG [34]. Fig.(7(c)) show the processing time taken by our method on different datasets. Fig.(7(d)) show a typical convergence curve of our algorithm. Ideally, it takes 120-150 iteration to provide an optimal solution to the problem.

## 7 Conclusion

Our Grassmannian representation of a non-rigidly deforming surface exploits the advantage of Grassmannians of different dimensions to jointly estimate better grouping of subspaces and their corresponding 3D geometry. Our approach explicitly leverages the geometric structure of the non-rigidly moving object w.r.t its neighbors on manifold via similarity graph and, it’s embedding in the lower dimension. We empirically demonstrated that our method is able to achieve 3D reconstruction accuracy which is better or as good as the state-of-the-art, with significant improvement in handling noisy trajectories.

Acknowledgement. The author was supported in part by the ARC Centre of Excellence for Robotic Vision (CE140100016), ARC Discovery project on 3D computer vision for geo-spatial localisation (DP190102261) and ARC DECRA project (DE140100180). The author thank his elder brother Aditya Kumar for his relentless support and advise. Thank you Roya Safaei for her constant help and constructive suggestions. Also, the author thank Dr. J. Varela and Dr. C. Punch, ACT for the medical treatment during the course of this work.

## References

• [1] P-A Absil, Robert Mahony, and Rodolphe Sepulchre. Optimization algorithms on matrix manifolds. Princeton University Press, 2009.
• [2] Ijaz Akhter, Yaser Sheikh, Sohaib Khan, and Takeo Kanade. Nonrigid structure from motion in trajectory space. In Advances in neural information processing systems, pages 41–48, 2009.
• [3] Ijaz Akhter, Yaser Sheikh, Sohaib Khan, and Takeo Kanade. Trajectory space: A dual representation for nonrigid structure from motion. IEEE Transactions on Pattern Analysis and Machine Intelligence, 33(7):1442–1456, 2011.
• [4] Mohammad Dawud Ansari, Vladislav Golyanik, and Didier Stricker. Scalable dense monocular surface reconstruction. arXiv preprint arXiv:1710.06130, 2017.
• [5] David Arthur and Sergei Vassilvitskii. k-means++: The advantages of careful seeding. In Proceedings of the eighteenth annual ACM-SIAM symposium on Discrete algorithms, pages 1027–1035. Society for Industrial and Applied Mathematics, 2007.
• [6] Jonathan F Bard. Practical bilevel optimization: algorithms and applications, volume 30. Springer Science & Business Media, 2013.
• [7] Thabo Beeler, Fabian Hahn, Derek Bradley, Bernd Bickel, Paul Beardsley, Craig Gotsman, Robert W Sumner, and Markus Gross. High-quality passive facial performance capture using anchor frames. In ACM Transactions on Graphics (TOG), volume 30, page 75. ACM, 2011.
• [8] Stephen Boyd, Neal Parikh, Eric Chu, Borja Peleato, and Jonathan Eckstein. Distributed optimization and statistical learning via the alternating direction method of multipliers.

Foundations and Trends® in Machine Learning

, 3(1):1–122, 2011.
• [9] Christoph Bregler, Aaron Hertzmann, and Henning Biermann. Recovering non-rigid 3d shape from image streams. In

IEEE Conference on Computer Vision and Pattern Recognition

, volume 2, pages 690–696. IEEE, 2000.
• [10] Hasan Ertan Cetingul and René Vidal. Intrinsic mean shift for clustering on stiefel and grassmann manifolds. In IEEE Conference on Computer Vision and Pattern Recognition, pages 1896–1902, 2009.
• [11] Yasuko Chikuse. Statistics on special manifolds, volume 174. Springer Science & Business Media, 2012.
• [12] T Collins and A Bartoli. Locally affine and planar deformable surface reconstruction from video. In International Workshop on Vision, Modeling and Visualization, pages 339–346, 2010.
• [13] Keenan Crane. Conformal geometry processing. California Institute of Technology, 2013.
• [14] Keenan Crane. Discrete differential geometry: An applied introduction. 2015.
• [15] Keenan Crane, Fernando De Goes, Mathieu Desbrun, and Peter Schröder. Digital geometry processing with discrete exterior calculus. In ACM SIGGRAPH 2013 Courses, page 7. ACM, 2013.
• [16] Keenan Crane and Max Wardetzky. A glimpse into discrete differential geometry. Notices of the American Mathematical Society, 64(10):1153–1159, November 2017.
• [17] Yuchao Dai, Huizhong Deng, and Mingyi He. Dense non-rigid structure-from-motion made easy-a spatial-temporal smoothness based solution. arXiv preprint arXiv:1706.08629, 2017.
• [18] Yuchao Dai, Hongdong Li, and Mingyi He. A simple prior-free method for non-rigid structure-from-motion factorization. International Journal of Computer Vision, 107(2):101–122, 2014.
• [19] Piotr Dollár, Vincent Rabaud, and Serge Belongie. Non-isometric manifold learning: Analysis and an algorithm. In International Conference on Machine Learning, pages 241–248, 2007.
• [20] Marián Fecko. Differential geometry and Lie groups for physicists. Cambridge University Press, 2006.
• [21] Mathias Gallardo, Toby Collins, Adrien Bartoli, and France Mathias. Dense non-rigid structure-from-motion and shading with unknown albedos. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pages 3884–3892, 2017.
• [22] Ravi Garg, Anastasios Roussos, and Lourdes Agapito. Robust trajectory-space tv-l1 optical flow for non-rigid sequences. In Energy Minimization Methods in Computer Vision and Pattern Recognition, pages 300–314. Springer, 2011.
• [23] Ravi Garg, Anastasios Roussos, and Lourdes Agapito. Dense variational reconstruction of non-rigid surfaces from monocular video. In IEEE Conference on Computer Vision and Pattern Recognition, pages 1272–1279, 2013.
• [24] Ravi Garg, Anastasios Roussos, and Lourdes Agapito. A variational approach to video registration with subspace constraints. International journal of computer vision, 104(3):286–314, 2013.
• [25] Paulo FU Gotardo and Aleix M Martinez. Computing smooth time trajectories for camera and deformable shape in structure from motion with occlusion. IEEE Transactions on Pattern Analysis and Machine Intelligence, 33(10):2051–2065, 2011.
• [26] Paulo FU Gotardo and Aleix M Martinez. Non-rigid structure from motion with complementary rank-3 spaces. In Computer Vision and Pattern Recognition (CVPR), 2011 IEEE Conference on, pages 3065–3072. IEEE, 2011.
• [27] Stephen Gould, Basura Fernando, Anoop Cherian, Peter Anderson, Rodrigo Santa Cruz, and Edison Guo. On differentiating parameterized argmin and argmax problems with application to bi-level optimization. arXiv preprint arXiv:1607.05447, 2016.
• [28] Jihun Hamm and Daniel D Lee. Grassmann discriminant analysis: a unifying view on subspace-based learning. In International conference on Machine learning, pages 376–383. ACM, 2008.
• [29] Onur C Hamsici, Paulo FU Gotardo, and Aleix M Martinez. Learning spatially-smooth mappings in non-rigid structure from motion. In European Conference on Computer Vision, pages 260–273. Springer, 2012.
• [30] Mehrtash Harandi, Conrad Sanderson, Chunhua Shen, and Brian Lovell. Dictionary learning and sparse coding on gr