Two Efficient Algorithms for Outlier Removal in Multi-view Geometry Using L Norm
|
|
- Christal Crawford
- 5 years ago
- Views:
Transcription
1 Two Efficient Algorithms for Outlier Removal in Multi-view Geometry Using L Norm Yuchao Dai 1, Mingyi He 1 Hongdong Li 3, 1 School of Electronics and Information, Northwestern Polytechnical University Shaanxi Key Laboratory of Information Acquisition and Processing Xi an, 7119, China Research School of Information Sciences and Engineering, Australian National University Canberra, ACT 1, Australia 3 Canberra Research Lab, NICTA Abstract L norm has been recently introduced to multi-view geometry computation to achieve globally optimal computation. It however suffers from a serious sensitivity to outliers. A few remedies have been proposed but with high computational complexity. This paper presents two efficient algorithms to overcome these problems. Our first algorithm is based on a cheap and effective local descent method (as opposed to the conventional but expensive SOCP(Second Order Cone Programming)). The second algorithm further improves the first one by using a Depth-first search heuristics. Both algorithms retain the nice property of global optimality of the L scheme, while at cost only a small fraction of the original computation. Experiments on both synthetic data and real images have validated the proposed algorithms. 1. Introduction A recent trend in Multi-view geometry computation in computer vision is to seek truly globally optimal (hence unique) solution. The L -norm based methods is such a scheme with the desired property. Under the L norm, an optimal solution with geometrically meaningful merit is guaranteed. The L scheme has therefore been recently applied to a variety of geometric vision problems such as triangulation, camera resection(uncalibrated case), homography, structure and motion recovery with known camera orientation [5], camera motion recovery [1], outlier detection [7], camera pose(calibrated case), two-view relative orientation [], multi-camera system [] and etc. Outlier Removing. A major criticism of using L minimization is, however, its sensitiveness to outliers. In order to apply the L optimization framework, the data should be free of outliers or with as few outliers as possible. To deal with this issue, researchers have tried various approaches. RANSAC [] is the most popular method. But, as it relies on random sampling to find a set of measurements containing only inliers, it may easily miss a few outliers out. This is already fatal for the L scheme because in which case even a single remaining outlier may destroy the estimation. Sim and Hartley [13] proposed a method by throwing away bad points with the maximal residuals. However, as their experimental results showed, this remedy is generally not very reliable for discriminating true inliers from outliers as it will throw outliers and inliers at the same time. Li [7] proposed a theoretic algorithm for L triangulation with outliers. The algorithm aimed to remove small number of outliers after RANSAC polish. However, the algorithm is generally expensive, particularly when the ratio of outliers is large. Another reason is that it uses SOCP as its solver, which we argue that SOCP is relatively expensive. Olsson et al. [11] proposed a framework that first finding local maxima subset through RANSAC then verifying whether the model is global maxima and they claimed that the algorithm had a polynomial-time bound. Martinec et al. [9] proposed to use L to delete the outliers in relative motion measurements or epipolar geometry(eg) rather than image points. In this paper we present two algorithms to solve the outlier removal problem more efficiently. The first one is a unified framework for outlier removal based on the so-called LP-type problem theory along with fast localdescent primitive solver. The second is a heuristic algorithm using depth-first priority search. Both algorithms are
2 able to speed up the computation significantly even when the faction of outliers is large. The algorithms proposed in this paper can be employed after the RANSAC in order to further remove any remaining outliers.. L with bisection method For completeness, this section briefly reviews the standard bisection SOCP method for solving L problems. Many problems in geometric vision may be written in the following form: A i x + b i min max x i c T i x + d i subject to c T i x + d i (1) It is easily seen that the problem (1) can be transformed, by introducing an additional variable γ, into an equivalent problem of the following form: min γ subject to A ix+b i γ(c T i x+d i ) and γ > γ,x () For any fixed value of γ, the set of points satisfying the constraint A i x + b i γ(c T i x + d i) is a convex set. Thus () is an SOCP problem which we may easily solve via a bisection scheme as follows. There has been recent work for speeding up the computations using fractional programming theory [1]. Require: Initial interval [γ l, γ u ] known to contain the optimal value of γ and tolerance ɛ >. 1: repeat : γ := (γ l + γ u )/. 3: Solve the SOCP feasibility problem to get x. : if feasible then 5: x = x, γ u := max i A ix+b i c T i x+di. : else 7: γ l := γ : end if 9: until γ u γ l ɛ. Algorithm 1: Improved Bisection Algorithm 3. Solving L optimization as LP-type Problem It is easily to observe that L minimization is determined by only part of the constraints thus implies that we could solve the problem efficiently. In this section, we describe the work on solving the L optimization problem treated as an LP-type problem LP-type Problem Paper [7] introduced a very original perspective to revisit the L problems in multi-view geometry,which is based on the so-called LP-type problems theory originally due to [1, 3]. Consider optimization problems specified by pairs (H, w), where H is a finite set, and w : H W is a function with values in a linearly ordered set (W, ); we assume that W has a minimum value. The elements of H are called constraints, and for G H, w(g) is called the value of G. Definition An abstract optimization problem is called an LP-type problem if the monotonicity and locality axioms below are satisfied. 1. Monotonicity axiom: For any F, G with F G H, we have w(f ) w(g).. Locality axiom: For any F G H with w(f ) = w(g) and for any h H, w(g) < w(g {h}) implies w(f ) < w(f {h}). The goal of LP-type computation is to compute a minimal subset B H H with the same value as H that is ω(b H ) = ω(h) (from which, in general, the value of H is easy to determine). Here we refer some basic terms to describe the LP-type problem. A set B of constraints is called a basis if w(b ) < w(b) for every proper subset of B. Given a set G of constraints, a subset B G is called a basis of set of constraints if it is a basis and w(b) = w(g)(i.e., an inclusion-minimal subset of G with equal w-value). Combinatorial dimension: The maximum cardinality of any basis in LP-type problem (H, w), denoted by δ = δ(h, w). Violation test: Decides whether or not w(b) < w(b {h}), for a basis B and a constraint h. Basis computation: Delivers a basis of B {h}, for a basis B and a constraint h. Support set: Given a minmax problem as with a solution Θ opt, the support set I supp is the set I supp = {i I f i (Θ opt ) = δ opt } From the definition of basis and support set, we can easily draw the conclusion that the basis for a set of constraints is the subset of the support set. We say that a constraint h H violates a subset G H if we have ω(g h) > ω(g). The violation set V (G) is all the constraints that violate G. The cardinality k of the violation set V (G) of G is called the level of the set G, k = V (G). For more details about the application of the LP-type theory to L computation, the reader is referred to [7] and references therein.. Algorithm-1: an exact approach In this section, we present our first algorithm for outlier removal under L norm. The key motivation is to find the optimal solution in the presence of outliers.
3 We first formulate the outlier removal with the example of triangulation under L norm. In this case we are given camera projection matrices P i, i = 1,, n and correspondences x i across multiple views. The aim is to infer the 3D point X in 3D space. Under L norm, we seek to attain the smallest maximal residual that is δ opt = min X max i f i (X). The definition of outlier can be stated that some measurements (the inliers) fit a parameterized model within a given maximum residual δ in while the other measurements are outliers. Denote the index set I of all the constraints made up of two subsets I in and I out as the inlier and outlier sets, respectively where the inlier set satisfies min X max i Iin f i (X) < δ in. According to this definition, the inlier set may be taken to be the largest subset, or any maximal (or even non-maximal) subset of I satisfying the inlier condition. On the other side, when we seeking model for the constraints, we aim to find the largest possible subset with a preset tolerance δ in. Thus we have the following mathematical model as: min max f i (X) X i I in max I in (3) Subject to f i (X) < δ in This is a multiple objective optimization problem. To solve the problem, a direct way is to change one of the objective function to constraint thus obtain the following model as: max I in Subject to f i (X) < δ in () which aims at a model with maximum inliers, whose residuals are all under given threshold. One of the key innovations that we use to speed up the computation is to accelerate the LP-type basis computation via a fast local method applied to primitive (sub)problems (to be defined later). First, we define the operation of Basis Change. Given the basis B k for the constraint set H, for each constraint in the basis b i, compute the new basis for the set H\{b i }. Thank to the locality, starting from the basis B k \{b i }, we are easily arrive at a new basis in fewer iteration. Second, for the tree structure contains all the basis at different levels, along all the paths, the function value will decease. So when we need to find the basis with the smallest function value, we only need to access the bases at the leaves of the tree. Algorithm-1 1. Initial basis finding: Finding the basis B satisfying all the constraints H = H, B = B(H ) and k =, n = 1;. Generate new bases: Generate new bases by performing the following procedures. Choose the current basis B n in the queue with level k; For each b B n, generate a basis B i = B(H n \b), H n = H n \b; (c) Add the new basis to the front of the queue, n = n + 1; 3. Selection Choose from the bases at level k with least L residual..1. Some details of the algorithm Primitive problem reduction. In [] [11], a main reduction theorem is proven which also gives an upper bound of the size of LP-type bases. Theorem.1 (main reduction theorem) Consider L problem min x max i f i (x), where x R n, i = 1,..., N, N > (n + 1), each f i (x) is pseudoconvex. Denote f I (x) = max i I f i (x), where I is a subset of N = 1,..., N. Also denote fi as the minimum of f I (x). Then there must exists a proper subset B which is of size B n + 1 such that fb = f N and x B = x N are their(equal) optimizers. This theorem predicts that the solution of the big problem is dominated by a small basis set of at most (n + 1) constraints. For the multi-view triangulation problem, the combinational dimension is, so only distinct primitive problems exist: 1-view which is trivial, -view, 3-view and -view. The left three primitive problems can be solved using local gradient method applying the KKT conditions. After obtaining the KKT conditions, the problem reduces to the following very small primitive problems and can be solved effectively by any local gradient descent method or Newton method, as did in []. -view case { Ci (x) = A i x + b i γ(c T i x + d i) =, i = 1, λ C 1 (x) + (1 λ) C (x) =, λ > 3-view case { Ci (x) = A i x + b i γ(c T i x + d i) =, i = 1,, 3 i λ i C i (x) =, i λ i = 1, λ i > -view case C i (x) = A i x + b i γ(c T i x + d i ) =, i = 1,...,.. Basis Change and Finding Now, we detail the steps for basis change using as an example the triangulation problem. The complete algorithm is given as follows.
4 1. Initialization: Randomly choose a pair of constraints as the initial basis B, set the number of elements in the basis as k =.. Violation Test: Solve the basis B t to obtain the 3D point X t, compute and sort the L residual for each constraint. If all the constraints are met, return the basis as B t ; else, select the constraint b which violates basis B t with the largest residual. 3. Compute new basis for subset B t b If the current basis consists of constraints denoted as constraints 1 and, while the candidate constraint to add is 3, test the bases {1, 3},{, 3} and {1,, 3} to meet all the constraints 1, and 3, set the corresponding basis size, the new basis is B t+1 ; If the current basis consists of 3 constraints denoted as constraints 1, and 3, while the candidate constraint to add is, test the bases {1, },{, }, {3, }, {1,, },{1, 3, }, {, 3, } and {1,, 3, } to meet all the constraints 1,,3 and, set the corresponding basis size, the new basis is B t+1 ; (c) If the current basis consists of constraints denoted as constraints 1,,3 and, while the candidate constraint to add is 5, test the bases {1, 5},{, 5}, {3, 5}, {, 5}, {1,, 5},{1, 3, 5}, {, 3, 5}, {,, 5}, {1,, 3, 5}, {1,,, 5},{1, 3,, 5} and {, 3,, 5} to meet all the constraints 1,,3, and 5, set the corresponding basis size, the new basis is B t+1 ;. Iteration: Do the above step until a basis satisfies all the constraints is found. Convergence Analysis. During each iteration, we add the constraint violate the basis with the largest residual thus the size of inlier set for the current basis B t will increase and the size of violation set will decrease. So the algorithm will terminate in finite steps and output a basis for the whole set I. The number of iteration is upper bounded by N the size of the whole set I however in practice the iteration will generally be far smaller than this bound. 5. Algorithm-: a heuristic depth-first search algorithm The above algorithm can obtain all the k-level basis thus output the optimal result. The upper bound for basis number is k = k Thus the complexity is bound by O( k+1 ) in the worst case. When k is small, we could obtain the optimal solution efficiently. However with the increase of k, the complexity become incredible to implement. Thus an efficient algorithm is in need. Let us consider the following model: min max f i (X) < δ in X i I in I out k where N is number of measurements and k is the upper bound of outliers in the measurements. To speed up the implementation, one choice is to relax the constraint on the size of inlier set. In Thus the model change to: min max f i (X) < δ in X i I in I out k + ɛ Compared with model 5, the number of outlier has been increased by ɛ that means we could remove a little more inliers constraints as outliers but it does not matter as at the cost of a few inliers removed we can speed up the implementation thus a compromise between efficiency and quality. Our analysis is based on some heuristic observation. As the basic property of outlier shown, with the outlier added, the L residual will change dramatically. So with the removing of the outliers, the L residual will decrease dramatically, when all outliers have been removed, the residual will not change more by removing more constraints. Based on this observation, we could apply a local gradient like algorithm. Given a basis and the constraints, we know that there is at least one outlier in the basis, removing each constraint in the basis will decrease the residual to different extent. To determine which constraint is the outlier with the largest possibility, we choose the one which cause the residual decrease with the largest amplitude. Repeated this operation, we could obtain a sequence of bases B, B 1,, B k. If we are lucky enough, we will obtain the optimal solution this way. If not we could remove a few more constraint to make the residual under the prefix threshold. When we do not obtain the feasible solution after the above iterations, or with more constraints removed we still can not obtain the feasible solution with residual under threshold, then we could apply backtrace operation to revisit the tree. The complete algorithm description is given below: Algorithm-: Depth-first search 1. Obtain the basis B at level ; (5) (). For each constraint in the basis, remove b B k, generate the new basis originally from the basis B k. Find the new basis with the smallest L residual. 3. Do the above step until the following cases are met:
5 The residual is lower than a given threshold δ; The depth is the larger than the upper bound of outlier number; The implementation process is illustrated in Figure 1 where B is the basis for all the constraints also the basis at level. The path that finds the basis for k = 3 is B B 1 B B 3. We can conclude that by solving bases number have linear relationship with level we can obtain the solution k= k=5 k= k= k= k=5 k= k=7 Figure. Synthetic experiment result for view with views corrupted by outliers. Log of residual for the view triangulation problem after removing,,5,,7 outliers respectively; Sorted residual for removing,,5,,7 outliers respectively. Log of k= k= k= k= Figure 1. Illustration of the operation on tree of bases Complexity Analysis For the problem of outlier removing, when applying the method of exhaustive search, the total number of feasible test is k i= Ci N. In [7] the total number of feasible test has been decreased to O( k ) which is independent of size of all the constraints N thus efficient reduce the number. In our new algorithm, the number of feasible tests is (k + ɛ).. Experiments To evaluate the performance of the proposed algorithms, we conduct both synthetic and real image experiments..1. Synthetic Data First to illustrate the success of the algorithm, we conduct two groups of experiments on synthetic data set. In the first group, we generate view triangulation problem where 3D points is randomly generated and projected on the image plane with Gaussian noise added and to simulate the effect of outlier, random image is corrupted by Gaussian noise with variance times the image noise. Figure illustrate the result which is obtained on 1 experiments. From this figure, we see that when we have removed k = constraint, almost all of the L residual is under.1, with further removing of constraints, the residual does decrease but with small amplitude thus show almost all the outliers have been removed successfully. The second group of synthetic experiments are to illustrate the performance for large number of outliers. In this experiment, outliers occurs in the 1 views. The result is shown as Figure 3. The figure shows the algorithm works well on large number of views. Figure shows the number of primitive problems for outlier dealing with the 5th, th, 7th outlier for the view 1 Figure 3. Synthetic experiment result for 1 view with views corrupted by outliers. Log of residual for the 1 view triangulation problem after removing,,1,,3 outliers respectively; Sorted residual for removing,,1,,3 outliers respectively. outlier problem. From the histogram, we see the number of primitive problems almost constant for different outlier Remving the 5th outlier Removing the th outlier Removing the 7th outlier Figure. Histogram of number of primitive problems for outlier. Removing the 5th outlier; Removing the th outlier; (c) Removing the 7th outlier. Experiments are conducted times for view triangulation with views is set as outlier... Real Image Experiments In order to test the proposed algorithm, we have used two publicly available sequences. 1 The first sequence consists of 11 images(of size 51 51) in a corridor. There are 1 points correspondence visible in all images. The other image set is a turnable sequence of a dinosaur, containing 3 images(of size 57 7) and, in total, 3 images correspondences with lots of occlusions. For this real image multiple view geometry triangulation problem, views of image points are perturbed by 5 pixel Gaussian noise as a simulation of the outliers. The result is shown as Figure and Figure 7. From the figures, we see the residual change dramatically from outlier to 1 outlier 1 Available at vgg/data.html (c)
6 resection, structure and motion known rotation and etc. Acknowledgment This work is supported by National Natural Science Foundation of China under key project number 737. The first author would like to thank the China Scholarship Council. NICTA is funded by the Australian Government as represented by the Department of Broadband, Communications and the Digital Economy and the Australian Research Council through the ICT Centre of Excellence program. (c) (d) Figure 5. Real Image: Corridor and Dinosaur until outlier, then the residual does not change more and under 1 pixel which means we have successfully removed almost all of the outliers as there is outliers in fact. 1 1 k= k= 1 3 Figure. Experimental Results on Corridor. Log of residual for the 1 view triangulation problem; for the 1 view triangulation problem. 1 1 k= k= Figure 7. Experimental Results on Dinosaur. Log of residual for the 1 view triangulation problem; for the 1 view triangulation problem. 7. Conclusion In this paper, two efficient algorithms are proposed to remove outliers using L norm. The first algorithm is based on a cheap and effective local descent method (as opposed to the conventional but expensive SOCP. The second algorithm further improves the first one by using a Depth-first search heuristics. Experiments on both synthetic and real data have proven the performance of the algorithms. For future work, rigorous analysis should be discussed and the same strategy can be applied to other outlier removing problem under L norm such as homograph, camera References [1] S. Agarwal, N. Snavely, and S. Seitz. Fast algorithms for L problems in multiview geometry. In CVPR, pages 1, June. [] M. A. Fischler and R. C. Bolles. Random sample consensus: a paradigm for model fitting with applications to image analysis and automated cartography. Commun. ACM, ():31 395, [3] J. E. Goodman and J. O Rourke, editors. Handbook of discrete and computational geometry. CRC Press, Inc., Boca Raton, FL, USA, [] R. I. Hartley and F. Kahl. Global optimization through rotation space search. IJCV, (1): 79, 9. 1 [5] F. Kahl and R. Hartley. Multiple-view geometry under the L -norm. PAMI, 3(9):13 117,. 1 [] J.-H. Kim, H. Li, and R. Hartley. Motion estimation for multi-camera systems using global optimization. In CVPR, pages 1,. 1 [7] H. Li. A practical algorithm for L triangulation with outliers. In CVPR, pages 1, 7. 1,, 5 [] H. Li. Efficient reduction for solving L-infinity problems in multiview geometry. In CVPR, pages 1, June 9. 3 [9] D. Martinec and T. Pajdla. Robust rotation and translation estimation in multiview reconstruction. In CVPR, volume 1, pages 1, June 7. 1 [1] J. Matousek, M. Sharir, and E. Welzl. A subexponential bound for linear programming. Algorithmica, 1(/5):9 51, 199. [11] C. Olsson, O. Enqvist, and F. Kahl. A polynomial-time bound for matching and registration with outliers. In CVPR, pages 1, June. 1, 3 [1] K. Sim and R. Hartley. Recovering camera motion using L minimization. In CVPR, pages ,. 1 [13] K. Sim and R. Hartley. Removing outliers using the L norm. In CVPR, pages 5 9,. 1
A practical algorithm for L triangulation with outliers
A practical algorithm for L triangulation with outliers Hongdong Li RSISE, Australian National University VISTA, National ICT Australia Abstract This paper addresses the problem of robust optimal multi-view
More informationEfficient Reduction of L-infinity Geometry Problems
Efficient Reduction of L-infinity Geometry Problems Hongdong Li Research School of Information Sciences and Engineering Australian National University and NICTA hongdong.li@anu.edu.au Abstract This paper
More informationOutlier Removal by Convex Optimization for L-Infinity Approaches
Outlier Removal by Convex Optimization for L-Infinity Approaches Yongduek Seo, Hyunjung Lee, and Sang Wook Lee Department of Media Technology, Sogang University, Korea {yndk,whitetobi,slee}@sogang.ac.kr
More informationOptimisation in Multiple View Geometry: The L-infinity Way (Part 2)
Optimisation in Multiple View Geometry: The L-infinity Way (Part 2) Tat-Jun Chin The University of Adelaide CVPR 2018 Tutorial Tat-Jun Chin (The University of Adelaide) Optimisation in Multiple View Geometry
More informationMotion Estimation for Multi-Camera Systems using Global Optimization
Motion Estimation for Multi-Camera Systems using Global Optimization Jae-Hak Kim, Hongdong Li, Richard Hartley The Australian National University and NICTA {Jae-Hak.Kim, Hongdong.Li, Richard.Hartley}@anu.edu.au
More informationVisual Odometry for Non-Overlapping Views Using Second-Order Cone Programming
Visual Odometry for Non-Overlapping Views Using Second-Order Cone Programming Jae-Hak Kim 1, Richard Hartley 1, Jan-Michael Frahm 2 and Marc Pollefeys 2 1 Research School of Information Sciences and Engineering
More informationGeneralized Convexity in Multiple View Geometry
Noname manuscript No. (will be inserted by the editor) Generalized Convexity in Multiple View Geometry Carl Olsson Fredrik Kahl Received: date / Accepted: date Abstract Recent work on geometric vision
More informationRotation Averaging with Application to Camera-Rig Calibration
Rotation Averaging with Application to Camera-Rig Calibration Yuchao Dai 1,2, Jochen Trumpf 2, Hongdong Li 3,2, Nick Barnes 3,2, and Richard Hartley 2,3 1 School of Electronics and Information,Northwestern
More informationFast Outlier Rejection by Using Parallax-Based Rigidity Constraint for Epipolar Geometry Estimation
Fast Outlier Rejection by Using Parallax-Based Rigidity Constraint for Epipolar Geometry Estimation Engin Tola 1 and A. Aydın Alatan 2 1 Computer Vision Laboratory, Ecóle Polytechnique Fédéral de Lausanne
More informationEfficient Optimization for L -problems using Pseudoconvexity
Efficient Optimization for L -problems using Pseudoconvexity Carl Olsson calle@maths.lth.se Anders P. Eriksson anderspe@maths.lth.se Centre for Mathematical Sciences Lund University, Sweden Fredrik Kahl
More informationA Summary of Projective Geometry
A Summary of Projective Geometry Copyright 22 Acuity Technologies Inc. In the last years a unified approach to creating D models from multiple images has been developed by Beardsley[],Hartley[4,5,9],Torr[,6]
More informationStereo and Epipolar geometry
Previously Image Primitives (feature points, lines, contours) Today: Stereo and Epipolar geometry How to match primitives between two (multiple) views) Goals: 3D reconstruction, recognition Jana Kosecka
More informationNonparametric estimation of multiple structures with outliers
Nonparametric estimation of multiple structures with outliers Wei Zhang and Jana Kosecka Department of Computer Science, George Mason University, 44 University Dr. Fairfax, VA 223 USA {wzhang2,kosecka}@cs.gmu.edu
More informationImage correspondences and structure from motion
Image correspondences and structure from motion http://graphics.cs.cmu.edu/courses/15-463 15-463, 15-663, 15-862 Computational Photography Fall 2017, Lecture 20 Course announcements Homework 5 posted.
More informationAutomatic Image Alignment (feature-based)
Automatic Image Alignment (feature-based) Mike Nese with a lot of slides stolen from Steve Seitz and Rick Szeliski 15-463: Computational Photography Alexei Efros, CMU, Fall 2006 Today s lecture Feature
More informationFast and Reliable Two-View Translation Estimation
Fast and Reliable Two-View Translation Estimation Johan Fredriksson 1 1 Centre for Mathematical Sciences Lund University, Sweden johanf@maths.lth.se Olof Enqvist 2 Fredrik Kahl 1,2 2 Department of Signals
More informationOptimal Relative Pose with Unknown Correspondences
Optimal Relative Pose with Unknown Correspondences Johan Fredriksson 1,Viktor Larsson 1, Carl Olsson 1 and Fredrik Kahl 1,2 1 Lund University, Sweden 2 Chalmers University of Technology, Sweden {johanf,
More informationarxiv: v1 [cs.cv] 28 Sep 2018
Camera Pose Estimation from Sequence of Calibrated Images arxiv:1809.11066v1 [cs.cv] 28 Sep 2018 Jacek Komorowski 1 and Przemyslaw Rokita 2 1 Maria Curie-Sklodowska University, Institute of Computer Science,
More informationNonparametric estimation of multiple structures with outliers
Nonparametric estimation of multiple structures with outliers Wei Zhang and Jana Kosecka George Mason University, 44 University Dr. Fairfax, VA 223 USA Abstract. Common problem encountered in the analysis
More informationHomographies and RANSAC
Homographies and RANSAC Computer vision 6.869 Bill Freeman and Antonio Torralba March 30, 2011 Homographies and RANSAC Homographies RANSAC Building panoramas Phototourism 2 Depth-based ambiguity of position
More informationVision par ordinateur
Epipolar geometry π Vision par ordinateur Underlying structure in set of matches for rigid scenes l T 1 l 2 C1 m1 l1 e1 M L2 L1 e2 Géométrie épipolaire Fundamental matrix (x rank 2 matrix) m2 C2 l2 Frédéric
More informationLeow Wee Kheng CS4243 Computer Vision and Pattern Recognition. Robust Methods. CS4243 Robust Methods 1
Leow Wee Kheng CS4243 Computer Vision and Pattern Recognition Robust Methods CS4243 Robust Methods 1 Consider this data set Fitting a line to blue points give blue line. Outliers cause fitting error (red
More informationRANSAC and some HOUGH transform
RANSAC and some HOUGH transform Thank you for the slides. They come mostly from the following source Dan Huttenlocher Cornell U Matching and Fitting Recognition and matching are closely related to fitting
More informationQuasiconvex Optimization for Robust Geometric Reconstruction
Quasiconvex Optimization for Robust Geometric Reconstruction Qifa Ke and Takeo Kanade, Computer Science Department, Carnegie Mellon University {Qifa.Ke,tk}@cs.cmu.edu Abstract Geometric reconstruction
More informationMultiple View Geometry. Frank Dellaert
Multiple View Geometry Frank Dellaert Outline Intro Camera Review Stereo triangulation Geometry of 2 views Essential Matrix Fundamental Matrix Estimating E/F from point-matches Why Consider Multiple Views?
More informationROBUST ESTIMATION TECHNIQUES IN COMPUTER VISION
ROBUST ESTIMATION TECHNIQUES IN COMPUTER VISION Half-day tutorial at ECCV 14, September 7 Olof Enqvist! Fredrik Kahl! Richard Hartley! Robust Optimization Techniques in Computer Vision! Olof Enqvist! Chalmers
More informationarxiv: v1 [cs.cv] 28 Sep 2018
Extrinsic camera calibration method and its performance evaluation Jacek Komorowski 1 and Przemyslaw Rokita 2 arxiv:1809.11073v1 [cs.cv] 28 Sep 2018 1 Maria Curie Sklodowska University Lublin, Poland jacek.komorowski@gmail.com
More informationCS 231A Computer Vision (Winter 2018) Problem Set 3
CS 231A Computer Vision (Winter 2018) Problem Set 3 Due: Feb 28, 2018 (11:59pm) 1 Space Carving (25 points) Dense 3D reconstruction is a difficult problem, as tackling it from the Structure from Motion
More informationROBUST LINE-BASED CALIBRATION OF LENS DISTORTION FROM A SINGLE VIEW
ROBUST LINE-BASED CALIBRATION OF LENS DISTORTION FROM A SINGLE VIEW Thorsten Thormählen, Hellward Broszio, Ingolf Wassermann thormae@tnt.uni-hannover.de University of Hannover, Information Technology Laboratory,
More informationInstance-level recognition part 2
Visual Recognition and Machine Learning Summer School Paris 2011 Instance-level recognition part 2 Josef Sivic http://www.di.ens.fr/~josef INRIA, WILLOW, ENS/INRIA/CNRS UMR 8548 Laboratoire d Informatique,
More informationFeature Based Registration - Image Alignment
Feature Based Registration - Image Alignment Image Registration Image registration is the process of estimating an optimal transformation between two or more images. Many slides from Alexei Efros http://graphics.cs.cmu.edu/courses/15-463/2007_fall/463.html
More informationAn Improved Evolutionary Algorithm for Fundamental Matrix Estimation
03 0th IEEE International Conference on Advanced Video and Signal Based Surveillance An Improved Evolutionary Algorithm for Fundamental Matrix Estimation Yi Li, Senem Velipasalar and M. Cenk Gursoy Department
More informationModel Fitting, RANSAC. Jana Kosecka
Model Fitting, RANSAC Jana Kosecka Fitting: Overview If we know which points belong to the line, how do we find the optimal line parameters? Least squares What if there are outliers? Robust fitting, RANSAC
More informationTwo View Geometry Estimation with Outliers
ENQVIST, KAHL: TWO VIEW GEOMETRY ESTIMATION WITH OUTLIERS 1 Two View Geometry Estimation with Outliers Olof Enqvist www.maths.lth.se/vision/people/olofe Fredrik Kahl www.maths.lth.se/vision/people/fredrik
More informationFitting. Fitting. Slides S. Lazebnik Harris Corners Pkwy, Charlotte, NC
Fitting We ve learned how to detect edges, corners, blobs. Now what? We would like to form a higher-level, more compact representation of the features in the image by grouping multiple features according
More information/ Approximation Algorithms Lecturer: Michael Dinitz Topic: Linear Programming Date: 2/24/15 Scribe: Runze Tang
600.469 / 600.669 Approximation Algorithms Lecturer: Michael Dinitz Topic: Linear Programming Date: 2/24/15 Scribe: Runze Tang 9.1 Linear Programming Suppose we are trying to approximate a minimization
More informationParameter estimation. Christiano Gava Gabriele Bleser
Parameter estimation Christiano Gava Christiano.Gava@dfki.de Gabriele Bleser gabriele.bleser@dfki.de Introduction Previous lectures: P-matrix 2D projective transformations Estimation (direct linear transform)
More informationPractical Global Optimization for Multiview Geometry
Practical Global Optimization for Multiview Geometry Sameer Agarwal 1, Manmohan Krishna Chandraker 1, Fredrik Kahl 2, David Kriegman 1, and Serge Belongie 1 1 University of California, San Diego, CA 92093,
More informationFitting (LMedS, RANSAC)
Fitting (LMedS, RANSAC) Thursday, 23/03/2017 Antonis Argyros e-mail: argyros@csd.uoc.gr LMedS and RANSAC What if we have very many outliers? 2 1 Least Median of Squares ri : Residuals Least Squares n 2
More informationB553 Lecture 12: Global Optimization
B553 Lecture 12: Global Optimization Kris Hauser February 20, 2012 Most of the techniques we have examined in prior lectures only deal with local optimization, so that we can only guarantee convergence
More informationA Systems View of Large- Scale 3D Reconstruction
Lecture 23: A Systems View of Large- Scale 3D Reconstruction Visual Computing Systems Goals and motivation Construct a detailed 3D model of the world from unstructured photographs (e.g., Flickr, Facebook)
More informationPractical Global Optimization for Multiview Geometry
Practical Global Optimization for Multiview Geometry Fredrik Kahl 1, Sameer Agarwal 2, Manmohan Krishna Chandraker 2, David Kriegman 2 and Serge Belongie 2 1 Centre for Mathematical Sciences Lund University,
More informationEE795: Computer Vision and Intelligent Systems
EE795: Computer Vision and Intelligent Systems Spring 2012 TTh 17:30-18:45 FDH 204 Lecture 14 130307 http://www.ee.unlv.edu/~b1morris/ecg795/ 2 Outline Review Stereo Dense Motion Estimation Translational
More informationModule 1 Lecture Notes 2. Optimization Problem and Model Formulation
Optimization Methods: Introduction and Basic concepts 1 Module 1 Lecture Notes 2 Optimization Problem and Model Formulation Introduction In the previous lecture we studied the evolution of optimization
More informationStable Structure from Motion for Unordered Image Collections
Stable Structure from Motion for Unordered Image Collections Carl Olsson and Olof Enqvist Centre for Mathematical Sciences, Lund University, Sweden Abstract. We present a non-incremental approach to structure
More informationExercise session using MATLAB: Quasiconvex Optimixation
Optimization in Computer Vision, May 2008 Exercise session using MATLAB: Quasiconvex Optimixation Overview In this laboratory session you are going to use matlab to study structure and motion estimation
More informationSrikumar Ramalingam. Review. 3D Reconstruction. Pose Estimation Revisited. School of Computing University of Utah
School of Computing University of Utah Presentation Outline 1 2 3 Forward Projection (Reminder) u v 1 KR ( I t ) X m Y m Z m 1 Backward Projection (Reminder) Q K 1 q Q K 1 u v 1 What is pose estimation?
More informationFeature Detectors and Descriptors: Corners, Lines, etc.
Feature Detectors and Descriptors: Corners, Lines, etc. Edges vs. Corners Edges = maxima in intensity gradient Edges vs. Corners Corners = lots of variation in direction of gradient in a small neighborhood
More informationCombining Appearance and Topology for Wide
Combining Appearance and Topology for Wide Baseline Matching Dennis Tell and Stefan Carlsson Presented by: Josh Wills Image Point Correspondences Critical foundation for many vision applications 3-D reconstruction,
More informationRANSAC: RANdom Sampling And Consensus
CS231-M RANSAC: RANdom Sampling And Consensus Roland Angst rangst@stanford.edu www.stanford.edu/~rangst CS231-M 2014-04-30 1 The Need for RANSAC Why do I need RANSAC? I know robust statistics! Robust Statistics
More informationWeek 2: Two-View Geometry. Padua Summer 08 Frank Dellaert
Week 2: Two-View Geometry Padua Summer 08 Frank Dellaert Mosaicking Outline 2D Transformation Hierarchy RANSAC Triangulation of 3D Points Cameras Triangulation via SVD Automatic Correspondence Essential
More informationPlanar Scene Modeling from Quasiconvex Subproblems
Planar Scene Modeling from Quasiconvex Subproblems by Visesh Chari, Anil Kumar Nelakantji, Chetan Jakkoju, C V Jawahar in ACCV 29 (Asian Conference on Computer Vision) Report No: IIIT/TR/29/228 Centre
More information1 Projective Geometry
CIS8, Machine Perception Review Problem - SPRING 26 Instructions. All coordinate systems are right handed. Projective Geometry Figure : Facade rectification. I took an image of a rectangular object, and
More informationClustering in Registration of 3D Point Clouds
Sergey Arkhangelskiy 1 Ilya Muchnik 2 1 Google Moscow 2 Rutgers University, New Jersey, USA International Workshop Clusters, orders, trees: Methods and applications in honor of Professor Boris Mirkin December
More informationAccurate Image Registration from Local Phase Information
Accurate Image Registration from Local Phase Information Himanshu Arora, Anoop M. Namboodiri, and C.V. Jawahar Center for Visual Information Technology, IIIT, Hyderabad, India { himanshu@research., anoop@,
More informationObservations. Basic iteration Line estimated from 2 inliers
Line estimated from 2 inliers 3 Observations We need (in this case!) a minimum of 2 points to determine a line Given such a line l, we can determine how well any other point y fits the line l For example:
More informationSrikumar Ramalingam. Review. 3D Reconstruction. Pose Estimation Revisited. School of Computing University of Utah
School of Computing University of Utah Presentation Outline 1 2 3 Forward Projection (Reminder) u v 1 KR ( I t ) X m Y m Z m 1 Backward Projection (Reminder) Q K 1 q Presentation Outline 1 2 3 Sample Problem
More informationA novel point matching method for stereovision measurement using RANSAC affine transformation
A novel point matching method for stereovision measurement using RANSAC affine transformation Naiguang Lu, Peng Sun, Wenyi Deng, Lianqing Zhu, Xiaoping Lou School of Optoelectronic Information & Telecommunication
More informationMultiple View Geometry in Computer Vision Second Edition
Multiple View Geometry in Computer Vision Second Edition Richard Hartley Australian National University, Canberra, Australia Andrew Zisserman University of Oxford, UK CAMBRIDGE UNIVERSITY PRESS Contents
More informationOn-line Convex Optimization based Solution for Mapping in VSLAM
On-line Convex Optimization based Solution for Mapping in VSLAM by Abdul Hafeez, Shivudu Bhuvanagiri, Madhava Krishna, C.V.Jawahar in IROS-2008 (Intelligent Robots and Systems) Report No: IIIT/TR/2008/178
More informationL minimization in geometric reconstruction problems. Abstract. 1 The triangulation problem. Richard Hartley 1,2 and Frederik Schaffalitzky 2,3
L minimization in geometric reconstruction problems Richard Hartley 1,2 and Frederik Schaffalitzky 2,3 1 National ICT Australia 2 The Australian National University 3 Oxford University Abstract We investigate
More informationCS 664 Image Matching and Robust Fitting. Daniel Huttenlocher
CS 664 Image Matching and Robust Fitting Daniel Huttenlocher Matching and Fitting Recognition and matching are closely related to fitting problems Parametric fitting can serve as more restricted domain
More informationStep-by-Step Model Buidling
Step-by-Step Model Buidling Review Feature selection Feature selection Feature correspondence Camera Calibration Euclidean Reconstruction Landing Augmented Reality Vision Based Control Sparse Structure
More informationCS664 Lecture #16: Image registration, robust statistics, motion
CS664 Lecture #16: Image registration, robust statistics, motion Some material taken from: Alyosha Efros, CMU http://www.cs.cmu.edu/~efros Xenios Papademetris http://noodle.med.yale.edu/~papad/various/papademetris_image_registration.p
More informationOptimal Algorithms in Multiview Geometry
Optimal Algorithms in Multiview Geometry Richard Hartley 1 and Fredrik Kahl 2 1 Research School of Information Sciences and Engineering, The Australian National University National ICT Australia (NICTA)
More informationVisual Odometry. Features, Tracking, Essential Matrix, and RANSAC. Stephan Weiss Computer Vision Group NASA-JPL / CalTech
Visual Odometry Features, Tracking, Essential Matrix, and RANSAC Stephan Weiss Computer Vision Group NASA-JPL / CalTech Stephan.Weiss@ieee.org (c) 2013. Government sponsorship acknowledged. Outline The
More informationA GLOBAL APPROACH FOR IMAGE ORIENTATION USING LIE ALGEBRAIC ROTATION AVERAGING AND CONVEX L MINIMISATION
The International Archives of the Photogrammetry, Remote Sensing and Spatial Information Sciences, Volume XL-3, 214 ISPRS Technical Commission III Symposium, 7 September 214, Zurich, Switzerland A GLOBAL
More informationNonrigid Surface Modelling. and Fast Recovery. Department of Computer Science and Engineering. Committee: Prof. Leo J. Jia and Prof. K. H.
Nonrigid Surface Modelling and Fast Recovery Zhu Jianke Supervisor: Prof. Michael R. Lyu Committee: Prof. Leo J. Jia and Prof. K. H. Wong Department of Computer Science and Engineering May 11, 2007 1 2
More informationOptimal Essential Matrix Estimation via Inlier-Set Maximization
Optimal Essential Matrix Estimation via Inlier-Set Maximization Jiaolong Yang 1,2, Hongdong Li 2, and Yunde Jia 1 1 Beijing Laboratory of Intelligent Information Technology, Beijing Institute of Technology,
More informationInteger Programming Theory
Integer Programming Theory Laura Galli October 24, 2016 In the following we assume all functions are linear, hence we often drop the term linear. In discrete optimization, we seek to find a solution x
More informationMathematical and Algorithmic Foundations Linear Programming and Matchings
Adavnced Algorithms Lectures Mathematical and Algorithmic Foundations Linear Programming and Matchings Paul G. Spirakis Department of Computer Science University of Patras and Liverpool Paul G. Spirakis
More informationCSE 527: Introduction to Computer Vision
CSE 527: Introduction to Computer Vision Week 5 - Class 1: Matching, Stitching, Registration September 26th, 2017 ??? Recap Today Feature Matching Image Alignment Panoramas HW2! Feature Matches Feature
More informationInstance-level recognition II.
Reconnaissance d objets et vision artificielle 2010 Instance-level recognition II. Josef Sivic http://www.di.ens.fr/~josef INRIA, WILLOW, ENS/INRIA/CNRS UMR 8548 Laboratoire d Informatique, Ecole Normale
More informationCS675: Convex and Combinatorial Optimization Spring 2018 The Simplex Algorithm. Instructor: Shaddin Dughmi
CS675: Convex and Combinatorial Optimization Spring 2018 The Simplex Algorithm Instructor: Shaddin Dughmi Algorithms for Convex Optimization We will look at 2 algorithms in detail: Simplex and Ellipsoid.
More informationRANSAC RANdom SAmple Consensus
Talk Outline importance for computer vision principle line fitting epipolar geometry estimation RANSAC RANdom SAmple Consensus Tomáš Svoboda, svoboda@cmp.felk.cvut.cz courtesy of Ondřej Chum, Jiří Matas
More informationReal-Time Scene Reconstruction. Remington Gong Benjamin Harris Iuri Prilepov
Real-Time Scene Reconstruction Remington Gong Benjamin Harris Iuri Prilepov June 10, 2010 Abstract This report discusses the implementation of a real-time system for scene reconstruction. Algorithms for
More informationAccurate Localization and Pose Estimation for Large 3D Models
Accurate Localization and Pose Estimation for Large 3D Models Linus Svärm 1 Olof Enqvist 2 1 Centre for Mathematical Sciences Lund University, Sweden {linus,magnuso}@maths.lth.se Magnus Oskarsson 1 Fredrik
More informationLarge Scale 3D Reconstruction by Structure from Motion
Large Scale 3D Reconstruction by Structure from Motion Devin Guillory Ziang Xie CS 331B 7 October 2013 Overview Rome wasn t built in a day Overview of SfM Building Rome in a Day Building Rome on a Cloudless
More informationFeature Matching and RANSAC
Feature Matching and RANSAC Recognising Panoramas. [M. Brown and D. Lowe,ICCV 2003] [Brown, Szeliski, Winder, CVPR 2005] with a lot of slides stolen from Steve Seitz, Rick Szeliski, A. Efros Introduction
More informationResearch Interests Optimization:
Mitchell: Research interests 1 Research Interests Optimization: looking for the best solution from among a number of candidates. Prototypical optimization problem: min f(x) subject to g(x) 0 x X IR n Here,
More informationDETECTION AND ROBUST ESTIMATION OF CYLINDER FEATURES IN POINT CLOUDS INTRODUCTION
DETECTION AND ROBUST ESTIMATION OF CYLINDER FEATURES IN POINT CLOUDS Yun-Ting Su James Bethel Geomatics Engineering School of Civil Engineering Purdue University 550 Stadium Mall Drive, West Lafayette,
More informationImage stitching. Digital Visual Effects Yung-Yu Chuang. with slides by Richard Szeliski, Steve Seitz, Matthew Brown and Vaclav Hlavac
Image stitching Digital Visual Effects Yung-Yu Chuang with slides by Richard Szeliski, Steve Seitz, Matthew Brown and Vaclav Hlavac Image stitching Stitching = alignment + blending geometrical registration
More informationCOMPARATIVE STUDY OF DIFFERENT APPROACHES FOR EFFICIENT RECTIFICATION UNDER GENERAL MOTION
COMPARATIVE STUDY OF DIFFERENT APPROACHES FOR EFFICIENT RECTIFICATION UNDER GENERAL MOTION Mr.V.SRINIVASA RAO 1 Prof.A.SATYA KALYAN 2 DEPARTMENT OF COMPUTER SCIENCE AND ENGINEERING PRASAD V POTLURI SIDDHARTHA
More informationGuided Image Super-Resolution: A New Technique for Photogeometric Super-Resolution in Hybrid 3-D Range Imaging
Guided Image Super-Resolution: A New Technique for Photogeometric Super-Resolution in Hybrid 3-D Range Imaging Florin C. Ghesu 1, Thomas Köhler 1,2, Sven Haase 1, Joachim Hornegger 1,2 04.09.2014 1 Pattern
More informationFast and stable algebraic solution to L 2 three-view triangulation
Fast and stable algebraic solution to L 2 three-view triangulation Zuzana Kukelova, Tomas Pajdla Czech Technical University, Faculty of Electrical Engineering, Karlovo namesti 13, Prague, Czech Republic
More informationImage Features: Local Descriptors. Sanja Fidler CSC420: Intro to Image Understanding 1/ 58
Image Features: Local Descriptors Sanja Fidler CSC420: Intro to Image Understanding 1/ 58 [Source: K. Grauman] Sanja Fidler CSC420: Intro to Image Understanding 2/ 58 Local Features Detection: Identify
More informationCS 395T Lecture 12: Feature Matching and Bundle Adjustment. Qixing Huang October 10 st 2018
CS 395T Lecture 12: Feature Matching and Bundle Adjustment Qixing Huang October 10 st 2018 Lecture Overview Dense Feature Correspondences Bundle Adjustment in Structure-from-Motion Image Matching Algorithm
More informationCS 231A Computer Vision (Winter 2014) Problem Set 3
CS 231A Computer Vision (Winter 2014) Problem Set 3 Due: Feb. 18 th, 2015 (11:59pm) 1 Single Object Recognition Via SIFT (45 points) In his 2004 SIFT paper, David Lowe demonstrates impressive object recognition
More informationUsing Subspace Constraints to Improve Feature Tracking Presented by Bryan Poling. Based on work by Bryan Poling, Gilad Lerman, and Arthur Szlam
Presented by Based on work by, Gilad Lerman, and Arthur Szlam What is Tracking? Broad Definition Tracking, or Object tracking, is a general term for following some thing through multiple frames of a video
More informationBayesian Methods in Vision: MAP Estimation, MRFs, Optimization
Bayesian Methods in Vision: MAP Estimation, MRFs, Optimization CS 650: Computer Vision Bryan S. Morse Optimization Approaches to Vision / Image Processing Recurring theme: Cast vision problem as an optimization
More informationMATH3016: OPTIMIZATION
MATH3016: OPTIMIZATION Lecturer: Dr Huifu Xu School of Mathematics University of Southampton Highfield SO17 1BJ Southampton Email: h.xu@soton.ac.uk 1 Introduction What is optimization? Optimization is
More informationMulti-view stereo. Many slides adapted from S. Seitz
Multi-view stereo Many slides adapted from S. Seitz Beyond two-view stereo The third eye can be used for verification Multiple-baseline stereo Pick a reference image, and slide the corresponding window
More informationAn LMI Approach for Reliable PTZ Camera Self-Calibration
An LMI Approach for Reliable PTZ Camera Self-Calibration Hongdong Li 1, Chunhua Shen 2 RSISE, The Australian National University 1 ViSTA, National ICT Australia, Canberra Labs 2. Abstract PTZ (Pan-Tilt-Zoom)
More informationAUTOMATIC RECTIFICATION OF LONG IMAGE SEQUENCES. Kenji Okuma, James J. Little, David G. Lowe
AUTOMATIC RECTIFICATION OF LONG IMAGE SEQUENCES Kenji Okuma, James J. Little, David G. Lowe The Laboratory of Computational Intelligence The University of British Columbia Vancouver, British Columbia,
More informationGeometry for Computer Vision
Geometry for Computer Vision Lecture 5b Calibrated Multi View Geometry Per-Erik Forssén 1 Overview The 5-point Algorithm Structure from Motion Bundle Adjustment 2 Planar degeneracy In the uncalibrated
More informationEdge and corner detection
Edge and corner detection Prof. Stricker Doz. G. Bleser Computer Vision: Object and People Tracking Goals Where is the information in an image? How is an object characterized? How can I find measurements
More informationHough Transform and RANSAC
CS4501: Introduction to Computer Vision Hough Transform and RANSAC Various slides from previous courses by: D.A. Forsyth (Berkeley / UIUC), I. Kokkinos (Ecole Centrale / UCL). S. Lazebnik (UNC / UIUC),
More informationMultiple-View Structure and Motion From Line Correspondences
ICCV 03 IN PROCEEDINGS OF THE IEEE INTERNATIONAL CONFERENCE ON COMPUTER VISION, NICE, FRANCE, OCTOBER 003. Multiple-View Structure and Motion From Line Correspondences Adrien Bartoli Peter Sturm INRIA
More informationHartley - Zisserman reading club. Part I: Hartley and Zisserman Appendix 6: Part II: Zhengyou Zhang: Presented by Daniel Fontijne
Hartley - Zisserman reading club Part I: Hartley and Zisserman Appendix 6: Iterative estimation methods Part II: Zhengyou Zhang: A Flexible New Technique for Camera Calibration Presented by Daniel Fontijne
More informationComputer Vision 6 Segmentation by Fitting
Computer Vision 6 Segmentation by Fitting MAP-I Doctoral Programme Miguel Tavares Coimbra Outline The Hough Transform Fitting Lines Fitting Curves Fitting as a Probabilistic Inference Problem Acknowledgements:
More information