A Simple Prior-Free Method for Non-rigid Structure-from-Motion Factorization : Revisited

02/27/2019 ∙ by Suryansh Kumar, et al. ∙ Australian National University 6

A simple prior free factorization algorithmdai2014simple is quite often cited work in the field of Non-Rigid Structure from Motion (NRSfM). The benefit of this work lies in its simplicity of implementation, strong theoretical justification to the motion and structure estimation, and its invincible originality. Despite this, the prevailing view is, that it performs exceedingly inferior to other methods on several benchmark datasetsjensen2018benchmark,akhter2009nonrigid. However, our subtle investigation provides some empirical statistics which made us think against such views. The statistical results we obtained supersedes Dai et. al.dai2014simple originally reported results on the benchmark datasets by a significant margin under some elementary changes in their core algorithmic ideadai2014simple. Now, these results not only exposes some unrevealed areas for research in NRSfM but also give rise to new mathematical challenges for NRSfM researchers. In this paper, we explore some of the hidden intricacies missed by Dai et. al. workdai2014simple and how some elementary measures and modifications can significantly enhance its performance, as high as 18% on the benchmark dataset. The improved performance is justified and empirically verified by extensive experiments on several datasets. We believe our work has both practical and theoretical importance for the development of better NRSfM algorithms. Practically, it can also help improve the recently reported state-of-the-art kumar2017spatio, kumar2016multi, jensen2018benchmark and other similar works in this field which are inspired by Dai et. al. workdai2014simple.



There are no comments yet.


page 1

page 4

page 8

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

Notation: The notation used in this paper is similar to Dai et al. work [9] unless otherwise stated.

(a) Paper Sequence
(b) 3D Reconstruction
(c) Tearing Sequence
(d) 3D Reconstruction
Figure 1: The method recovers 3D dimensional structure of the deforming object over multiple frames. Our elementary but powerful changes provides a substantial improvement in the reconstruction accuracy than the previous results reported for “prior-free" approach. The example images are taken from the recently released NRSfM Challenge Dataset [14]

. Our reconstruction results are nearly as good as the best performing algorithm without using very complex and involved mathematical optimization


Non-rigid Structure from Motion (NRSfM) is a well-known problem in geometric computer vision. The goal of this problem is to reconstruct 3D structure of a deforming object using multiple frames. One of the most popular way to solve NRSfM is the matrix factorization approach. The matrix factorization approach to solve this problem dates back to 2000

[5] with no satisfactory solution in place until 2012. In the year 2012, Dai [8] proposed a ground-breaking approach to solve NRSfM. This method for solving NRSfM is now considered as a classical work in NRSfM [9]. In that paper, the camera motion is estimated by imposing the null space constraint and the rank-3 positive semi-definite matrix cone constraint on the Gram matrix (). Further, nuclear norm minimization of the reshuffled shape matrix () was introduced to proffer stronger rank bound on the shape matrix for non-rigid shape estimation. The striking part of their work is that it not only challenged the myth of the inherent basis ambiguity in NRSfM [33] but also supplied a practical “prior-free" algorithm to solve NRSfM. Nonetheless, over years, it was observed that their remarkable theory performs poorly on benchmark datasets [22, 14]. In this paper, our goal is to make “prior-free idea” work well on real world scenario’s.

Theoretically, the elementary idea of Dai [9] conveniently encapsulates all the basic intuitions which are required to solve a general NRSfM problem. One may immediately argue on its performance when the deforming shape is composed of a union of low-rank subspace [19, 36]. However, in this paper, we restrict our discussion to the classical representation of a NRSfM problem [5], without paying much attention to, how clustering benefits 3D reconstruction of the non-rigid object and other such notions of compact data representation. The reason for this choice is that the improvement in the performance of a classical baseline shall benefit the methods built on top of it [19, 36].

The main purpose of this work is to uncover some of the unexplored mathematical intricacies in the prior free factorization approach to NRSfM, and improve on the idea supplied by Dai [9]. Our exposition leads to the possible reasons for its inferior performance on the benchmark datasets [1, 14, 31]. It is shown in this paper that the rotation estimate using Dai work [9] is not unique under the same model complexity prior (/rank), and they overlooked to utilize full correction matrix space [4]. Our investigation unveil the possibility of procuring motion that satisfies the well-known assumption of smooth non-rigid deformation of the object [27]. A simple search for the proper column-triplet (triads [4]) for the correction matrix () based on the smoothness of camera motion can indeed help improve the accuracy of the algorithm. Further, we argue that the weighted nuclear norm minimization of the shape matrix () is a far better choice than its global trace norm minimization. Lastly, due to our extensive analysis, we are able to posit some unsolved issues in NRSfM under “prior-free” idea which needs attention for further progress in this field.

In this paper, it is not claimed that we achieve state-of-the-art results on the benchmark datasets using our new approach. However, we empirically show that we can get very close to the best performing approaches and the difference is not very great, without the employment of complex and involved mathematical optimization [19, 22]. In this paper, we also argue that the inferior performance of “prior-free" method may not be due the proposed theoretical idea but because they overlooked some of the mathematical construction in their own formulation, and missed on properly utilizing the well-known assumptions about non-rigidly moving object , smooth deformation [27] and low-rank shape [9]. Hence, the conclusion, understanding and use of simple “prior-free” algorithm to NRSfM is not complete and precise. Through this work, we try to amend and nullify the prevailing perception about the “prior-free" approach, and how it can be used to its maximum potential. We feel that our paper touches some critical points which are essential to establish a theoretical closure to some of the elementary problems within the factorization approach to NRSfM.

Contribution: Firstly, our work postulates some rectification to the usage of “Intersection Method" [9] to compute camera motion. With the suitable example, we establish that the generalization made on the rotation matrix estimation by Dai work [9] is not convincing and therefore, the knowledge about the strength of “Intersection theorem” is not completely exploited. Secondly, we provide an analytic solution to estimate suitable rotation using Intersection theorem and conjecture some challenges associated with it. Lastly, we propose a weighted nuclear norm minimization problem to estimate non-rigid 3D shape. Our approach shows a substantial improvement in the 3D reconstruction accuracy (as high as 18%). Moreover, we observed performance improvement in the case of noisy and missing trajectories §4.2 (with a minor adjustment) using our method.

In this work, our attempt is to make the baseline method111By baseline, we mean the methods that solve NRSfM using its classical representation that have withstood the test of time [30, 5]. more accurate, both in terms of understanding and performance, subject to the mathematical simplicity. To achieve this, we attempt to avoid the usage of complex mathematical notions such as union of independent subspace, dependent subspace representation [36, 19, 21]

, procrustean normal distribution

[22], kernelization [10] . Hence, it is simple to understand the theoretical and practical justification of our method. We show that by applying simple but powerful logical and mathematical modifications to the prior free idea [9], we can get close to or even perform better at times than the best available algorithms on the benchmark datasets.

2 Representation and Motion Estimation

1. Classical Representation:

Tomasi and Kanade factorization method to structure-from-motion under orthographic camera projection appropriately summarizes the behavior of the 3D points over frames [30]. The relation between 3D shape, motion and its projection over frames was defined as


where, is the measurement matrix formed by stacking all the image coordinates () for ‘’ points along ‘’ rows , total number of frames. denotes the orthographic camera rotation matrix with each as per frame rotation. represent the shape matrix with each row triplet as a 3D shape. This representation was later extended by Bregler [5] to recover non-rigid 3D shapes. More concretely,


The matrix ‘’ and ‘’ are composed of shape bases and shape coefficients respectively, with ‘’ as the number of shape bases. ‘’ denotes the kronecker product and ‘’ is a identity matrix. It is evident from the above formulation that the rank of and also . However, is not a general rank matrix but own a special structure due to factor [9].

2. Null Space Representation of the Orthonormality Constraint:

An initial step in the factorization approach to NRSfM is to perform a rank decomposition of the measurement matrix

via singular value decomposition (

svd) . This is then followed by the estimation of Euclidean corrective matrix ‘’ to solve rotation and 3D structure. The main reason for such a procedure is due to the fact that the singular value decomposition of ‘’ matrix is not unique as any non-singular matrix in between the two matrices and can form a valid factorization. Mathematically,


Now, once we are able to solve correctly, then rotation and shape can be estimated using the above relations [5]. To solve , orthonormality constraints are imposed . Representing the double row of as and as the column triplet of , then using Eq:(2) and Eq:(3) we can write


Multiplying both sides by from right side gives

This leads to two linear equation constraint

where, . Using the algebraic relation = , Dai transformed these constraints (Eq:5) to a null space representation as follows:


Using the above form and previous work in NRSfM [33], Dai proposed the intersection theorem and supplied a SDP solution to estimate the matrix and the Euclidean corrective matrix using svd().

Theorem 1

Intersection Theorem: Under non-generate and noise-free conditions, any correct solution of must lie in the intersection of the dimensional null-space of and a rank 3 positive semi-definite matrix cone must belong to

Figure 2: (a) The column triplet (1:3) of euclidean corrective matrix () used by Dai work [9] shown in red shade. It is stated with the notion that there is no loss of generality to choose . However, choosing other column triplet may result in better rotation and shape estimate as shown in Figure 4(a) and 4(b)

Dai solution to rotation: They proposed that once the is solved, rather than solving for full Euclidean corrective matrix , use svd() to extract rank 3 . The solved can then be use to find (Eq:4) up to sign (). The method quote “we adopt a simpler approach that directly computes the camera motion from single column-triplet without need to fill in a big and full matrix”. Naturally, this single column-triplet is chosen to be the first column-triplet of the matrix (see Fig:2(a)). Now, such strategy give rise to few legitimate concerns

  1. [nolistsep, itemindent=0em, label=()]

  2. When each column triplet qualifies for a suitable correction matrix, then why has a high preference? Are we loosing useful information by such unwarranted preference?

  3. Will each provide the same solution to the rotation matrix?

  4. Generally, most real world deformations are smooth in nature [27]. Whether such solution to rotation is good enough for the smooth deformation assumption?

Dai overlooked all these intrinsic issues to solve rotation using their proposed intersection theorem.

Plausible Rectification: Our experiments show that Dai [9] solution to rotation estimation actually aborted the useful information present in the . Each of the ‘’ column triplets in () gives a possible rotation matrix which is different from each other (see Fig:(3)). Our empirical evaluations on several datasets show that the first column triplet is not always the best choice to estimate rotation. Hence, the details provided by Dai work [9] is incomplete and there is a loss of generality with such procedure to estimate rotation under the well-known assumption of smooth deformation [27]. Fig:(4(a)) and Fig:(4(b)) provides few statistical results with comparison for both rotation and shape error estimate respectively. For clarity, we also provide the column triplet index that gives the better results for the corresponding data sequence and therefore, provides few counter-examples to such generalization.

Figure 3: The rotation samples on using for Pick-up sequence. Below each manifold is the graph showing the per frame change in the camera motion using Eq:(8) ‘’. A simple observation establishes that all rotation matrix () are not the same. ‘’ graph analysis on this dataset show that the rotation estimate provided by , has a smoother camera motion than other ’s, with being the smoothest. Any one out of these 3 ’s supply better performance than . Note: Each via cross product. (Best viewed on screen)

Theoretically and practically, this result is of significant importance as it helps in inferring that the solution provided by “Intersection Theorem" has a lot of useful information left to be exploited completely and Dai work ignored this. Also, it gives rise to some challenges that finding the best column triplet for is not an easy task. With these results, we conjecture few problems for further research in NRSfM that are: (a) Can we find a best possible column triplet for the corrective matrix with a given rank prior ‘’, or (b) At least can we put an upper bound on the value such that there exists no such ‘’ for which will provide better rotation and structure estimate. The problem seems hard keeping in view that the prior rank in NRSfM factorization methods is an assumed approximation and it changes for different datasets to achieve better results.
A solution: In 2005, Brand. M [4] argued to use full correction matrix to estimate motion which in a way utilizes all the multiple estimates of column-triads of . Recently, Lee [23] briefly mentions on the problems with motion estimates using [9]. In contrast, we use an analytical observation based on the smoothness and regularity222The term <<regularity>> is used in a loose sense (Mathematically). of the camera motion trajectory to filter to infer better ‘’. Let be a function that takes as input and gives ‘’ as output using Intersection theorem. We estimate different for all the column triplets , then compute smoothness of the camera motion for each as:


By examining the smoothness of the camera motion for each , we select the suitable rotation matrix for structure estimation (see Fig. 3). Our strategy to select smooth camera motion over frames based on Eq:(8) consistently supplied us with better performance than the previously proposed approach. We acknowledge that this is not a profound way to infer the best rotation, however, it does provide a possibility to deduce better rotation using “prior-free" approach which respects the well-known assumption of smooth deformation in NRSfM. Further, it helps endorse our claim on the generalization of rotation estimate by [9]. You may use the variable ‘’ Eq:(8) as a smoothness term in the final optimization (Eq:(11)) to further improve rotation, however, to show the competence within the “prior-free” idea [9], we stick to the classical two staged approach.

Figure 4: Counter examples on benchmark dataset [2]. (a) Rotation error in comparison to BMM [9] on synthetic data. (b) 3D reconstruction error using global trace norm minimization of shape matrix as used in BMM with rotation matrix estimate using other column triplet in comparison to (1:3). The column triplets of () for which the method perform better on Drink, Pickup, Yoga and Stretch are (1:3), (19:21), (1:3) and (19:21) respectively. Note that we used the same rank prior value ‘’ used in Dai work [9].

3 Structure Estimation

Once the rotation is estimated based on the smoothness of the camera motion [27], the next step is to solve for 3D structure. The block matrix method (BMM) by Dai [9] proposed the following optimization problem to estimate the non-rigid low-rank shape.


where, is a rearranged shape matrix with each row corresponds to the shape for that frame. The trace norm minimization on ‘’ is enforced instead of ‘’ to provide a stronger rank bound on the shape matrix [9]. The second term in Eq:(9) enforces the re-projection error constraint. The function maps to .

Dai solution to shape: Following the work of Ma [26] on rank minimization problems, Dai [9] proposed a solution to the optimization in Eq:(9). The method enforces low-rank constraint on ‘’ matrix and provide the solution by solving Eq:(9) via ADMM[3] using matrix shrinkage operator = , where = with = and 0 otherwise }.

Plausible Rectification: Despite the trace norm minimization provides a satisfactory solution to non-rigid structure estimation, it has some serious issues. The proposed solution to nuclear norm minimization problem (Eq:(9)) gives equal priority to each singular values, as a result, the shrinkage operator penalizes each singular value with the same quantity (). To estimate 3D structure of a non-rigidly deforming object using matrix factorization approach, we use a prior assumption that the shape lies in a low-rank subspace. Therefore, it’s not a better choice to penalize the major component of the shape data and its very minor component equally. Consequently, nuclear norm minimization of the shape matrix struggles to appropriately conserve the useful component of the non-rigidly deforming shape.

Truncated nuclear norm regularization can be a choice to handle such issues, however, it depends on the binary decision, hence not versatile in nature [35]. To really cater the behavior of the deformations based on its low-rank nature, we propose to use weighted nuclear norm minimization approach to solve for non-rigid structure [28, 12]. In contrast to the previous notation to the nuclear norm of the shape matrix , we introduce a different notation for its weighted nuclear norm


where denotes the singular value of . We assume that the weights ’s are non-negative scalar . Using this representation, we redefine the optimization proposed in the Eq:(9) as follows:


The motivation for such formulation is quite clear, however, the proposed optimization (Eq:11) is generally non-convex, and is more difficult to solve than the nuclear norm minimization. Fortunately, recent results [34, 25, 12] in compressed sensing have shown that we can achieve a global optimal solution to Eq:(11) in the case when §3.1.

3.1 Optimization

This section provides the mathematical derivation to the optimization proposed in Eq:(11). Our solution use the following theorems and proofs as stated and used in [34, 12, 7].

Theorem 2

For all , denoted by , the SVD of it. The solution to

, with non-negative weight vector

, its solution can be written as , where is the solution to the following optimization problem

Theorem 3

If the singular values and the weights satisfy then the weighted nuclear norm minimization problem has a globally optimal solution


where is the SVD of , and is the generalized soft-thresholding operator with weight vector


The readers are encouraged to refer to [34, 12] work for detailed derivations to the lemma’s leading to the proof of the theorems. In conclusion, if the weights satisfies non-descending order, not necessarily with the same value, the weighted nuclear norm minimization problem is still convex and optimal solution can be obtained using a soft-thresholding operator with different weights [34, 12].

3.2 Solution

We propose our solution to the optimization problem defined in Eq:(11) using alternating direction method of multipliers [3] (ADMM), a simple, fast but powerful algorithm used to solve many convex and non-convex problems in computer vision and mathematical optimization. The ADMM algorithm decompose the original problem into several sub-problems, where each of them is solved separately by introducing Lagrange multipliers and penalty parameters to estimate convergence. Using the method of multipliers, the Augmented Lagrangian form for Eq:(11) is written as follows:


here is a Lagrange multiplier and is the penalty parameter. The solution to each variable is obtained by solving the following subproblems over iterations (indexed with the variable ):


The Lagrange multiplier and the penalty parameter are updated as follows:


refers to the maximum value of ‘’ and is an empirical constant (). The mathematical derivations to each sub-problems are provided in the supplementary material for reference. The closed form solution to the Eq:(17) is obtained by taking the derivative of Eq:(15) w.r.t variable ‘’ and equating it to zero ,


Similarly, rewriting the Eq:(15) treating as variable.


In contrast to the previous form, the solution to Eq:(20) is not straight forward. To obtain a closed form solution to this problem, lets define a soft-thresholding function . Also, let be the singular value decomposition of (), then the optimal solution to Eq:(20) is given by:


Here, is the weight assigned to the different singular values in the non-descending order based on its significance to the deformation data. For detail discussion on the initialization of weights kindly refer section §4.1.

4 Experiment and Discussion

To endorse our claim, we performed extensive experiments on both real and synthetic benchmark datasets [1, 14, 31]. We compared the performance of our algorithm against different state-of-the-art methods on these datasets [11, 22, 19]. Additionally, we unveil the substantial percentage boost in the reconstruction accuracy as high as 18% in comparison to the previous results reported for “simple prior-free" approach. For real-world applications to NRSfM, noisy data and missing feature tracks over frames are crucial, therefore, we also performed experiments to tackle such issues. To make the comparisons on noisy and missing data sequence, the experimental settings we used are same and consistent with Lee work [22]. Experimental results on dense datasets and more rigorous cases of missing trajectories are provided in the supplementary material. Before we provide details on our performance analysis, lets discuss the variable initialization.

4.1 Initialization

Our algorithm has few parameters and variables to initialize. For all our experiments on different datasets, we initialize , , , , and the ‘’ values are kept same as Dai method [9]. Practically, we considered the convergence of our optimization, if the gap or over iteration.

(a) Shark
(b) Pick-up
(c) Drink
(d) Face
(e) Dance
(f) Yoga
(g) Walking
(h) Stretch
Figure 5: Reconstruction results of our method on the NRSfM synthetic benchmark dataset [1, 2]. Ground-truth and reconstructed points are shown in filled(red) and non-filled circles respectively. Note: We used same ‘’ value as documented in [9] work for all the experiments.

1. Structure initialization: Using the result of Liu [24] on the uniqueness of minimizer for the rank minimization problem, we initialize the the 3D shape ‘’ as ‘’ = and . The pseudo-inverse solution to shape matrix provides a good enough initialization to our algorithm. Reader may refer to Dai [9] and Valmadre [32] work for detailed discussion on pseudo inverse solution to ‘’ in NRSfM.
2. Weight () initialization: It is well-known in NRSfM under factorization approach that the shape matrix lies in a low-rank space. Generally, the largest singular value of the shape matrix contains the most information about the non-rigid shape, therefore, while optimizing for the shape matrix, it’s illogical to treat each singular value equally. The singular values with major component must be penalized less and vice-versa. Using this inverse relation between singular values and its significance to the shape deformation modeling, we assign the weight () to be inversely propositional to the singular values of the shape matrix.


where, is a positive number and , a very small positive number to avoid division by zero as some singular values are likely to be zero (low rank). We initialized the weights by substituting the pseudo-inverse initialization of ‘’ using the relation = in the Eq:(22).

4.2 Performance Analysis

After a detailed discussion on the variable initialization and optimization, we present our experimental evaluation. We performed extensive experiments on both new and old benchmark datasets [1, 31, 14]. We report the quantitative result on the previous benchmark dataset using mean normalized 3D reconstruction error formulation


where, , are the estimated 3D shape and ground-truth 3D shape respectively. To keep our statistics consistent with the newly proposed NRSfM dataset, we used their error evaluation code to compute the robust root mean square error (RMSE) metric as proposed in Taylor work [29]

. For more details on NRSfM CVPR 2017 challenge dataset evaluation metric, kindly refer to Jensen work


Method PTA[1] CSF2[11] PND[22] BMM Ours
Drink 0.0287 0.0227 0.0037 0.0266 0.0119 (1.47%)
Pickup 0.1939 0.1791 0.0372 0.1731 0.0198 (15.3%)
Yoga 0.1243 0.1179 0.0140 0.1150 0.0129 (10.2%)
Stretch 0.1035 0.1136 0.0156 0.1034 0.0144 (8.90%)
Dance 0.2426 0.1877 0.1454 0.1864 0.1060 (8.04%)
Walking 0.3761 0.1938 0.0465 0.1298 0.0882 (4.16%)
Face 0.0489 0.0319 0.0165 0.0303 0.0179 (1.24%)
Shark 0.2933 0.1117 0.0135 0.2311 0.0551 (18.0%)
Table 1: Performance comparison in the shape recovery using our new approach with some of the state-of-the-art methods in single body NRSfM. The statistics clearly demonstrate our claim that we can achieve a significant improvement in the reconstruction accuracy without using complex mathematical formulation. The percentage value in the last column (red) show the improvements over the result documented by Dai original work (BMM) [9].

1. Benchmark datasets: Most of the methods proposed in non-rigid structure from motion often use it to evaluate the performance of the algorithm. Loosely speaking, this dataset is composed of eight standard sequences namely Drink, Pickup, Yoga, Stretch, Dance, Walking, Face and Shark. The number of frames () to number of points () set for these datasets are (1102, 41), (357, 41), (307, 41), (370, 41), (264, 75), (316, 40) and (240, 91) respectively. Table (1) show the statistical comparison of our approach in comparison to the other competing approaches for single body NRSfM. Our evaluation results clearly show a significant improvement in the reconstruction accuracy in comparison to the previously reported results for “prior-free” approach. Figure (5) show the qualitative reconstruction results w.r.t ground-truth on all of these sequences.
2. NRSfM challenge datasets: Jensen recently released this dataset as a part of NRSfM competition held at CVPR 2017 [14]. This is a high quality challenging dataset divided into five categories based on the deformation type, namely, Articulated, Balloon, Paper, Stretch and Tearing. Each of these categories is again shot using six different camera paths namely circle, flyby, line, semi-circle, tricky and zig-zag. This dataset is significantly larger and diverse to really test the performance of a NRSfM algorithm’s. However, the dataset provides only a single frame ground-truth 3D for each of the five categories to test the algorithm. To estimate the reliability of our approach, we compared our performance against the best performing algorithm on this dataset. Table (2) show the quantitative results of our method. The performance clearly demonstrates the significant improvement in the accuracy using “prior-free” idea under our modification. It also help infer that without using complex mathematical notions, we can reach performance accuracy close to the state-of-the-art. Figure (6) show some qualitative results using our method.

Method / Data Articulated Ballon Paper Stretch Tearing
Multibody [19] 10.15 10.64 15.78 9.96 14.17
BMM [9] 24.54 12.91 22.37 18.71 18.87
Ours 12.02 11.79 16.21 12.05 16.08
Table 2: Performance comparison of our method in comparison to the best performing algorithm (Multi-body) [19] on NRSfM challenge dataset [14]. The above statistics shows the average root-mean-square error in millimeters for the single test image on the orthogonal sequence available with the dataset. Our method shows a clear improvement over the originally proposed BMM approach and it’s accuracy got very close to the multi-body.
(a) Articulated
(b) Ballon
(c) Stretch
(d) 3D Articulated
(e) 3D Ballon
(f) 3D Stretch
Figure 6: Reconstruction results of our method on the NRSfM challenge dataset [14]. The results shown here are for the circular camera path. Ground-truth 3D and reconstructed 3D points are shown with filled and non-filled circles respectively.

3. Noisy data:

The feature tracks captured from a real-world motion capture system is noisy most of the time. Therefore, to test the reliability and robustness of our new approach, we performed experiments by re-synthesizing the trajectories added with Gaussian noise. We introduced the Gaussian noise with standard deviation set as

, where is varied from 0.05-0.25 [22]. Figure (7(a)) shows the variation in the normalized average 3D error for the stretch sequence using the performance of different algorithm recorded over 20 times. The plot clearly shows the robustness of our algorithm in comparison to other methods in the presence of large noise ratio’s.
4. Missing Data: In addition to the noisy data, the other problem with 3D reconstruction from a real video sequence is the missing trajectories over frames. We handle the missing trajectory quite robustly by incorporating a simple modification to the optimization proposed in Eq:(11). Let’s assume is the incomplete measurement matrix and is the mask matrix which indicates the presence or absence of the tracks over frames. Given , , we first find a complete matrix using the following optimization


The above optimization is a well studied optimization form. To keep things simple, we used Cabral work [6] to estimate . The motive is to first solve for complete ‘’ to estimate camera motion using our rectified approach §2, and then solve for shape using the following cost function:


Clearly, it’s just a minor adjustment to the proposed method based on the kind of data available in different situations. To evaluate our performance, we randomly set 30% of the data missing from the sequence same as Lee work [22] for comparison. Figure (7(b)) shows the performance of our algorithm with missing data. Supplementary video is provided to show results on real world missing data cases.

Discussion: Why not add the motion regularization in the final optimization and solve for both motion and shape? It’s definitely a valid argument. Nevertheless, we wanted to show the competence in a “prior free” way [9] which is to “solve for motion first using Intersection theorem and then solve for low-rank shape”, therefore, we avoided to add it in the final optimization. We showed that smooth solution [27] already exist within the solution to intersection theorem. Comprehensive analysis of our algorithm after adding motion regularization to solve the final optimization is left as an extension to the present idea.

Figure 7: (a) 3D reconstruction error comparison over noisy trajectories. (b) Comparison of our method performance with other competing methods with missing data in the measurement matrix. Note: BMM was not formulated for missing data case, therefore, its results are not present in the above figures.

5 Conclusion

With weighted nuclear norm minimization of the shape matrix and an analytic solution to the rotation matrix based on the smoothness of the camera motion [27], we witnessed that the prior-free idea performs almost as good as the best available algorithm’s. Without exploiting the “prior-free” idea [9] fully based on the well-known assumptions of smooth deformation of the non-rigid object and its low-rank shape, it may perform badly, which might be the reason that researchers have had poor results using it, even for the non-rigid objects that span a single linear subspace. Our work revealed the possibility of making “prior-free” [9] practically more accurate under the different conditions of measurement matrix with elementary modifications, and also conjecture some open problems. The accuracy of our algorithm on the benchmark datasets empirically validates that the “prior-free” theory is still a very powerful way to solve NRSfM and therefore, the proposition before the NRSfM researchers to consider is, it’s not the failure of the concept behind the prior-free idea for its inferior performance but, it’s possibly due to our inability to correctly cater, and cleverly exploit the arc of information and perspectives provided by it to solve NRSfM.


  • [1] I. Akhter, Y. Sheikh, S. Khan, and T. Kanade. Nonrigid structure from motion in trajectory space. In Advances in neural information processing systems, pages 41–48, 2009.
  • [2] I. Akhter, Y. Sheikh, S. Khan, and T. Kanade. Trajectory space: A dual representation for nonrigid structure from motion. IEEE Transactions on Pattern Analysis and Machine Intelligence, 33(7):1442–1456, 2011.
  • [3] S. Boyd, N. Parikh, E. Chu, B. Peleato, and J. Eckstein. Distributed optimization and statistical learning via the alternating direction method of multipliers.

    Foundations and Trends® in Machine Learning

    , 3(1):1–122, 2011.
  • [4] M. Brand. A direct method for 3d factorization of nonrigid motion observed in 2d. In

    2005 IEEE Computer Society Conference on Computer Vision and Pattern Recognition (CVPR’05)

    , volume 2, pages 122–128. IEEE, 2005.
  • [5] C. Bregler, A. Hertzmann, and H. 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.
  • [6] R. Cabral, F. De la Torre, J. P. Costeira, and A. Bernardino. Unifying nuclear norm and bilinear factorization approaches for low-rank matrix decomposition. In Proceedings of the IEEE International Conference on Computer Vision, pages 2488–2495, 2013.
  • [7] J.-F. Cai, E. J. Candès, and Z. Shen. A singular value thresholding algorithm for matrix completion. SIAM Journal on Optimization, 20(4):1956–1982, 2010.
  • [8] Y. Dai, H. Li, and M. He. A simple prior-free method for non-rigid structure-from-motion factorization. In Computer Vision and Pattern Recognition (CVPR), 2012 IEEE Conference on, pages 2018–2025. IEEE, 2012.
  • [9] Y. Dai, H. Li, and M. He. A simple prior-free method for non-rigid structure-from-motion factorization. International Journal of Computer Vision, 107(2):101–122, 2014.
  • [10] P. F. Gotardo and A. M. Martinez. Kernel non-rigid structure from motion. In IEEE International Conference on Computer Vision, pages 802–809. IEEE, 2011.
  • [11] P. F. Gotardo and A. 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.
  • [12] S. Gu, L. Zhang, W. Zuo, and X. Feng. Weighted nuclear norm minimization with application to image denoising. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pages 2862–2869, 2014.
  • [13] R. Hartley, J. Trumpf, Y. Dai, and H. Li. Rotation averaging. International journal of computer vision, 103(3):267–305, 2013.
  • [14] S. H. N. Jensen, A. Del Bue, M. E. B. Doest, and H. Aanæs. A benchmark and evaluation of non-rigid structure from motion. arXiv preprint arXiv:1801.08388, 2018.
  • [15] S. Kumar. Jumping manifolds: Geometry aware dense non-rigid structure from motion. arXiv preprint arXiv:1902.01077, 2019.
  • [16] S. Kumar, A. Cherian, Y. Dai, and H. Li. Scalable dense non-rigid structure-from-motion: A grassmannian perspective. arXiv preprint arXiv:1803.00233, 2018.
  • [17] S. Kumar, Y. Dai, and H. Li. Multi-body non-rigid structure-from-motion. In 3D Vision (3DV), 2016 Fourth International Conference on, pages 148–156. IEEE, 2016.
  • [18] S. Kumar, Y. Dai, and H. Li. Monocular dense 3d reconstruction of a complex dynamic scene from two perspective frames. In IEEE International Conference on Computer Vision, pages 4649–4657, Oct 2017.
  • [19] S. Kumar, Y. Dai, and H. Li. Spatio-temporal union of subspaces for multi-body non-rigid structure-from-motion. Pattern Recognition, 71:428–443, May 2017.
  • [20] S. Kumar, R. S. Ghorakavi, Y. Dai, and H. Li. Dense depth estimation of a complex dynamic scene without explicit 3d motion estimation. arXiv preprint arXiv:1902.03791, 2019.
  • [21] V. Larsson and C. Olsson. Compact matrix factorization with dependent subspaces. In Proceedings of the IEEE Conference on Computer Vision and Pattern Recognition, pages 280–289, 2017.
  • [22] M. Lee, J. Cho, C.-H. Choi, and S. Oh. Procrustean normal distribution for non-rigid structure from motion. In IEEE Conference on Computer Vision and Pattern Recognition, pages 1280–1287, 2013.
  • [23] M. Lee, J. Cho, and S. Oh. Consensus of non-rigid reconstructions. In IEEE Conference on Computer Vision and Pattern Recognition, pages 4670–4678, 2016.
  • [24] G. Liu, Z. Lin, S. Yan, J. Sun, Y. Yu, and Y. Ma. Robust recovery of subspace structures by low-rank representation. IEEE Transactions on Pattern Analysis and Machine Intelligence, 35(1):171–184, 2013.
  • [25] C. Lu, J. Tang, S. Yan, and Z. Lin. Nonconvex nonsmooth low-rank minimization via iteratively reweighted nuclear norm. arXiv preprint arXiv:1510.06895, 2015.
  • [26] S. Ma, D. Goldfarb, and L. Chen. Fixed point and bregman iterative methods for matrix rank minimization. Mathematical Programming, 128(1-2):321–353, 2011.
  • [27] V. Rabaud and S. Belongie. Re-thinking non-rigid structure from motion. In IEEE Conference on Computer Vision and Pattern Recognition, pages 1–8. IEEE, 2008.
  • [28] N. Srebro and T. Jaakkola. Weighted low-rank approximations. In Proceedings of the 20th International Conference on Machine Learning (ICML-03), pages 720–727, 2003.
  • [29] J. Taylor, A. D. Jepson, and K. N. Kutulakos. Non-rigid structure from locally-rigid motion. In IEEE Conference on Computer Vision and Pattern Recognition, pages 2761–2768. IEEE, 2010.
  • [30] C. Tomasi and T. Kanade. Shape and motion from image streams under orthography: a factorization method. International Journal of Computer Vision, 9(2):137–154, 1992.
  • [31] L. Torresani, A. Hertzmann, and C. Bregler. Nonrigid structure-from-motion: Estimating shape and motion with hierarchical priors. IEEE transactions on pattern analysis and machine intelligence, 30(5):878–892, 2008.
  • [32] J. Valmadre, S. Sridharan, S. Denman, C. Fookes, and S. Lucey. Closed-form solutions for low-rank non-rigid reconstruction. In Digital Image Computing: Techniques and Applications (DICTA), 2015 International Conference on, pages 1–6. IEEE, 2015.
  • [33] J. Xiao, J.-x. Chai, and T. Kanade. A closed-form solution to non-rigid shape and motion recovery. In European conference on computer vision, pages 573–587. Springer, 2004.
  • [34] Z. Zha, X. Zhang, Y. Wu, Q. Wang, Y. Bai, and L. Tang. Analyzing the weighted nuclear norm minimization and nuclear norm minimization based on group sparse representation. arXiv preprint arXiv:1702.04463, 2017.
  • [35] D. Zhang, Y. Hu, J. Ye, X. Li, and X. He. Matrix completion by truncated nuclear norm regularization. In Computer Vision and Pattern Recognition (CVPR), 2012 IEEE Conference on, pages 2192–2199. IEEE, 2012.
  • [36] Y. Zhu, D. Huang, F. De La Torre, and S. Lucey. Complex non-rigid motion 3d reconstruction by union of subspaces. In IEEE Conference on Computer Vision and Pattern Recognition, pages 1542–1549, 2014.