Augmented Reality VU Camera Registration. Prof. Vincent Lepetit
|
|
- Antony White
- 6 years ago
- Views:
Transcription
1 Augmented Reality VU Camera Registration Prof. Vincent Lepetit
2 Different Approaches to Vision-based 3D Tracking [From D. Wagner] [From Drummond PAMI02] [From Davison ICCV01] Consider natural features Consider natural features Make use of visual markers No a priori 3D knowledge Use some 3D knowledge Use some 3D knowledge 13
3 Camera Model From World to Images 14
4 From the World to the Image World coordinate system M m What is the relation between the 3D coordinates of a point M and its correspondent m in the image captured by the camera? 15
5 Perspective Projection World coordinate system M m C Camera center The image formation is modeled as a perspective projection, which is realistic for standard cameras: The rays passing through a 3D point M and its correspondent m in the image all intersect at a single point C, the camera center. 16
6 Expressing M in the Camera Coordinates System World coordinate system M M cam m Z C Y X Camera coordinate system Step 1: Express the coordinates of M in the camera coordinates system as M cam. This transformation corresponds to a Euclidean displacement (a rotation plus a translation): M cam = RM + T where: R is a 3x3 rotation matrix, and T is a 3- vector. 18
7 Homogeneous Coordinates World coordinate system X X M = Y M Y = Z Z 1 M cam C Y Z Lets replace M by the 4- homogeneous vector : Just add a 1 as the fourth coordinate. m X Camera coordinate system Now, the Euclidean displacement can be expressed as an linear transformation instead of an affine one: M cam = RM + T X cam Y cam Z cam X ( ( = R ( Y ( + T Z M X cam Y cam Z cam X ( ( ( = ( R T Y ) ( Z( ( 1 ( ) M cam = R T M 19 (R T ) is a 3x4 matrix.
8 The matrix ( R T) = The External Calibration Matrix R 11 R 13 R 13 T 1 R 21 R 22 R 23 T 2 R 31 R 32 R 33 T 3 is called the External calibration matrix, or the extrinsic calibration matrix, or the matrix of external parameters. It can be parameterized by 6 values: 3 for the rotation, 3 for the translation. The parameterization of the translation is trivial, while parametrizing rotations is not. We will detail the possible parameterization of rotations in the 3D space. 20
9 21 Homogeneous Coordinates A mathematical trick to express non-linear transformations such as translation or projection as linear transformations. A 3- vector expressed in homogeneous coordinates becomes a 4-vector: To go back to the 3-vector, one just has to divide the 3 first coordinates by the last one. à A homogeneous vector is defined up to a scale factor: X Y Z kx ky kz k All the homogeneous vectors of the form kx ky kz k correspond to the same 3 vector X Y Z. The homogeneous vector u v w t corresponds to the 3 vector u t v t w t.
10 22 Homogeneous Coordinates In homogeneous coordinates, an Euclidean displacement can be written using a 4x4 matrix: The Computer Vision community uses a 3x4 matrix and considers that the 4th coordinate of is 1: k X k Y k Z k ( ( ( ( = R T ( ( ( ( kx ky kz k ( ( ( ( X cam Y cam Z cam = R T ( ) X Y Z 1 M
11 Projection M cam C Y Z X f Camera coordinate system m Z m X M cam Computation of the coordinates of m in the image plane, from M cam (expressed in the camera coordinates system): Simply use Thales theorem: f m X f = X Z m X = f X Z C m X X 23
12 From Projection to Image Coordinates of m in pixels? Image coordinate u m f m X = f X Z, m Y = f Y Z system m u =?, m v =? v 0 v u 0 1 k u 1 1 pixel k v C Camera coordinate system 24
13 - (u 0, v 0 ): Projection of the camera center C in the image coordinate system. - k u : Inverse of the size of one pixel along X axis. - k v : Same for Y axis. Image coordinate system v 0 u v f u 0 m m X = f X Z, m = f Y Y Z m u = u 0 + k u m X, m v = v 0 + k v m Y C Camera coordinate system 25
14 Putting the perspective projection and the transformation into pixel coordinates together: m X = f X Z, m Y = f Y Z m u = u 0 + k u m X, m v = v 0 + k v m Y In matrix form: u k u f 0 u 0 X v = 0 k v f v 0 Y w Z u v defines m in homogeneous coordinates w m u = u w = u + k f X 0 u Z m v = v w = v + k f Y 0 v Z 26
15 The matrix The Internal Calibration Matrix k u f 0 u 0 0 k v f v is called the Internal calibration matrix, or the intrinsic calibration matrix, or the calibration matrix, or the matrix of internal parameters. It can be parameterized by 4 values: α u, α v, u 0, v 0 : α u 0 u 0 ( 0 α v v 0 ( The number of parameters can be reduced under some assumptions: (u 0, v 0 ) is sometimes taken at the center of the image; Taking α u = α v assumes that the pixels are square. 27
16 28 The Full Transformation The two transformations are chained to form the full transformation from a 3D point in the world coordinate system to its projection in the image: The product of the internal calibration matrix and the external calibration matrix is a 3x4 matrix called the projection matrix. The projection matrix is defined up to a scale factor. u v w = k u f 0 u 0 0 k v f v R 11 R 13 R 13 T 1 R 21 R 22 R 23 T 2 R 31 R 32 R 33 T 3 X Y Z 1 = P 11 P 12 P 13 P 14 P 21 P 22 P 23 P 24 P 31 P 32 P 33 P 34 X Y Z 1 projection matrix
17 Notations R, T World coordinate system M C Camera center m m = A[ R T] M = PM A: Internal calibration matrix (3x3); R: Rotation matrix (3x3); T: Translation vector (3-); P: Projection matrix (3x4); M : world 3D point in homogenous coordinates (4-); m: image 2D point in homogenous coordinates (3-), projection of M in the image; 29
18 World coordinate system M R, t m Camera center C Camera coordinate system
19 Internal Parameters (A) A : Often assumed to be known in the following. Can be estimated beforehand using a calibration grid. m i = A[R T]M i 31
20 Camera Distortion The previous camera model does not take into account distortions due to fish eye lenses or bad quality lenses. Fish eye lenses are interesting because images have a larger field of view. Distortion can be removed before applying perspective projection model. 32
21 Distortion Estimation (u,v) (u c,v c ) ( u, v ) Radial + tangential distortions. Tangential distortion component has much less influenceà neglected. (u c,v c ) : Center of radial distortion. The correction is written: ( ) ( ) u = u c + L(r) u u c v = v c + L(r) v v c with ( ) 2 + ( v v c ) 2 and r 2 = u u c L(r) =1+ κ 1 r + κ 2 r 2 + κ 3 r Once the distortion parameters are estimated, the distortion can be compensated in two ways: Undistort the images, using a look-up table for efficiency; Distort the model projection when computing a reprojection error. 33
22 Camera Pose (Rotation + translation) Estimation Algorithms From 3D-2D points correspondences: DLT algorithm + Extraction of A, R, t; P3P algorithm; POSIT and others; Non-linear minimization. From a planar structure. 34
23 Linear Estimation of P from 3D/2D point correspondences (DLT algorithm) Each correspondence gives us 2 equations: m i = P P u i = 11 X i + P 12 Y i + P 13 Z i + P 14 P M 31 X i + P 32 Y i + P 33 Z i + P 34 i v i = P X + P Y + P Z + P 21 i 22 i 23 i 24 P 31 X i + P 32 Y i + P 33 Z i + P 34 P 11 X i + P 12 Y i + P 13 Z i + P 14 P 31 X i u i P 32 Y i u i P 33 Z i u i P 34 u i = 0 P 21 X i + P 22 Y i + P 23 Z i + P 24 P 31 X i v i P 32 Y i v i P 33 Z i v i P 34 v i = 0 X i Y i Z i X i u i Y i u i Z i u i u i ( X i Y i Z i 1 X i v i Y i v i Z i v i v i P P 34 ( ( = 0 Stacking all the equations into matrix B yields the linear system: B P P 34 = 0 35
24 DLT algorithm P P 34 ( ( B Question: Why is the null vector not a valid solution? P P 34 = 0 is the eigen vector of B corresponding to the smallest eigenvalue of B. 11 coefficients to estimate; each 3D-2D correspondence gives 2 equations: 6 3D-2D correspondences must be known. In practice, much more are needed to get a good estimation. 36
25 Extraction of A, R, and T from P 1. A is known. 2. A is not known. Trick: Consider the matrix made of the first 3 columns of P: P 3 = AR P = ( P 3 c 4 ) And compute P 3 P 3 T P 3 P 3 T = (AR)(AR) T = ARR T A T = AA T à since A is a upper-triangular matrix, it can be found using Choleski factorization of P 3 P 3 T ( AA T ) 1 is called the image of the absolute conic. 37
26 Extraction of R, and T from P Once A is known: [R t] A 1 P No guarantee that R is a rotation matrix: R still must be ortho-normalized. 38
27 Ortho-Normalization of R Lets Q = A -1 P 3 a first approximation of R. No guarantee that Q is a rotation matrix. How to find the best rotation matrix that approximates Q? 1. Singular Value Decomposition of Q: Q = USV T, where S = diag(s 1, s 2, s 3 ) U and V have orthogonal columns, so that U T U = V T V = I 2. R = UV T [From Zhang TR98] 39
28 Ortho-Normalization of R: Proof Lets Q = A -1 P 3 a first approximation of R. No guarantee that Q is a rotation matrix. How to find the best rotation matrix that approximates Q? - Take best in the sense of the smallest Frobenius norm of the difference R Q. R Q F 2 min R = trace R Q R Q F Minimizing R - Q 2 is equivalent to maximize trace(r T Q)) Let the SVD of Q be USV T, where S = diag(s 1, s 2, s 3 ): with Z =V T R T U. Maximum reached with Z = I R = UV T 2 subject to RR T = I (( ) T ( R Q) ) (definition of Frobenius norm) = trace(r T R Q T R R T Q + Q T Q) = 3 + trace Q T Q trace(r T Q) ( ) 2trace( R T Q) = trace(r T USV) = trace(v T RUS) = trace(zs) = z ii s i s i i i 40 [From Zhang TR98]
29 DLT : Summary Not a good idea when the internal parameters are known. 41
30 P-n-P Problem PnP: Perspective-n-Point Estimation of the Position and Orientation (R and T) from 2D/3D point correspondences when the internal parameters (A) are known. m i M i R, T 42
31 P-n-P Problem PnP: Perspective-n-Point Estimation of the Position and Orientation (R and T) from 2D/3D point correspondences when the internal parameters (A) are known. 3 correspondences are sufficient (6 equations for 6 parameters), BUT yield up to 4 solutions, generally only two. A fourth correspondence is needed to remove the ambiguity. M i R, T 43
32 The P3P Problem M 3 M 1 M2 m 1 m 3 m 2 C P3P: Pose estimation from 3 correspondences. Yield up to 4 solutions, generally only two. A fourth correspondence will be needed to remove the ambiguity. 44
33 M 3 M 1 R, t M2 m 1 m 3 m 2 C
34 The P3P Problem M i d ij Mj m i θ ij m j x i C Each pair of correspondences M i m i and M j m j gives a constraint on the (unknown) camera-point distances x i = M i - C and x j = M j - C : d ij 2 = x i 2 + x j 2 2 x i x j cosθ ij, where: d ij = M i - M j is the (known) distance between M i and M j ; θ ij is the (also known) angle sustended at the camera center by M i and M j. 46
35 Overview of the P3P M i d ij Mj m i θ ij m j x i C 1. Solve for the distances x i ; 2. The positions M i C of the points M i in the camera coordinates system can be computed; 3. R and T are computed as the Euclidean displacement from the M i to the M i C.. 47
36 Computation of cosθ ij m i θ ij m j C cosθ ij must be computed. It depends only on the (known) 2D positions m i and m j. ( cosθ ij = Cm j ) T ( Cm i ) ( Cm j ) T ( Cm i ) = Cm j Cm i Cm i = A 1 m i cosθ ij = m j T ωm i ( m T j ωm j ) 1 2 ( m T i ωm i ) ( Cm j ) T ( Cm j ) ( Cm i ) T ( Cm i ) 1 2 with ω = AAT ( ) -1 the image of the absolute conic. 49
37 Quadratic System d ij 2 = x i 2 + x j 2 2 x i x j cosθ ij f ij (x i, x j ) = x i 2 + x j 2 2 x i x j cosθ ij - d ij 2 = 0 f 12 (x 1, x 2 ) = 0 f 13 (x 1, x 3 ) = 0 Resultant g(x 1,x 2 ) [degree 4] f 23 (x 2,x 3 ) = 0 Resultant h(x 1 ) [degree 8] h(x 1 ): Polynomial of degree 8 in x 1 with (fortunately) only even terms Degree polynomial of degree 4 in x=x 12 : At most 4 solutions for x; Can be solved in closed form. x 1 is positive and x 1 = x. x 2 and x 3 are uniquely determined from x 1. To obtain a unique solution, add one more point: Solve the P3P problem for each subset of 3 points, and keep the common solution. 51
38 Resultant of Two Polynomials Sylvester resultant (for two polynomials of a single unknown): p 1 (x) = a m x m + a m-1 x m a 1 x + a 0 p 2 (x) = b n x n + b n-1 x n b 1 x + b 0 Syl(p 1, p 2 ) = a m a a m a 0 a m a 0 b n b b n b 0 b n b 0 p 1 and p 2 have a common root iff determinant(syl(p 1, p 2 )) = 0. determinant(syl(p 1, p 2 )) is called the Sylvester resultant of p 1 and p 2. 52
39 Example: System of two polynomials of two unknowns p 1 (x, y) = 6x 2 + 3xy - xy 2 + y + 1 p 2 (x, y) = x 2 y + 5x + 4y - 1 p 1 (x, y) = 0 p 2 (x, y) = 0 Lets considers y as a constant, and write these two polynomials as polynomials in x: p 1 (x, y) = 6x 2 + (3y - y 2 )x + (y + 1) p 2 (x, y) = y x 2 + 5x + (4y - 1) 6 0 y (3y y 2 ) 6 5 (y +1) (3y y 2 ) (4y 1) 0 (y +1) 0 0 y 5 (4y 1) determinant(syl(p 1, p 2 )) = 4y 6-20y y 4-310y y y + 36 is a polynomial in y only! = 0 First, solve 4y 6-20y y 4-310y y y + 36 = 0 Once y is known, x can be found from p 1 (x, y) or p 2 (x, y). 54
40 Estimating R and T The coordinates of the 3D points M i in the camera coordinates system can now be computed: M i M i C = x i A 1 m i m i x i C 55
41 Estimating R and T The positions M i C of the points M i in the camera coordinates system can now be computed: M i W R,T? M i C n M W = 1 W M i N W i = M W i M W M C = 1 n n i=1 v t.( Ru) = trace R t uv t C ( ) ( N i ) t W.( RN i ) i max R rotation matrix i C ( N i ) t W. RN i = trace R t C N i i ( ) ( ) t.n i W R = argmax trace(r t L) R rotation matrix R = UV t where L = USV t T = M C RM W n i=1 M i C ) = trace R t L ( N i C = M i C M C C ( ) where L = ( N i ) t W.N i i 56
42 What are the possible solutions? Special case: the center of projection lies on the plane perpendicular to the triangle: C M 2 M 1 M 3 [From The Perspective View of Three Points, Wolfe et al. PAMI91] A second solution can be obtained by rotating the triangle about the side M 1 M 2 57
43 What are the possible solutions? General case: M 2 C M 1 58
44 What are the possible solutions? General case: M 3 M 3 M 2 C M 1 [From The Perspective View of Three Points, Wolfe et al. PAMI91] The set of circles generates a surface. Each intersection between the line of sight [CM 3 ) and this surface corresponds to one solution to the P3P problem. 59
45 P3P: Summary Well adapted to the problem Closed-form solution. Only 3 points + one constraint needed to compute the pose. Inconvenient: Pose computed with only 3 points, can be inaccurate. à Pose refined by a non-linear estimation with all the points available. 60
46 PnP Lu-Hager-Mjolsness: One on the most efficient algorithms in the general case. Iterative algorithm. Requires an initialization. Pick a starting R 0 Repeat Compute T(R k ) [can be computed in closed-form] Compute N i (R k ) Compute L(R k ) Compute L(R k ) = U D V t Set R k+1 = V U t Until convergence Can converge into a local minimum. 61
47 PnP algorithm: Application [Moreno-Noguer et al. ICCV07] 62
48 64 Pose Estimation from a Plane Let assume that the homography between the plane and the plane image is known. H m= Hm ku kv k = a b c d e f g h i u v 1 m = PM with M a point on to the plane = A R 1 R 2 R 3 T [ ] X Y 0 1 = A R 1 R 2 T [ ] X Y 1 = H X Y 1
49 Pose Estimation from a Plane We have: H A[ R 1 R 2 T] R 1, R 2 and T can be retrieved from the product G = A -1 H G is defined up to a scale factor and R 1 and R 2 are not perfectly orthonormal. Normalization : l = G 1 G 2 R 1 = G 1 l, R 2 = G 2 l,t = G 3 l R 2 c c = R 1 + R 2, p = R 1 R 2, d = c p R 1 = 1 2 c c + d d (, R 2 = 1 2 c c d d ( p R 1 R 3 = R 1 R 2 + Refinement with a non-linear optimization. d 65
50 Pose Estimation from a Plane: Application [Benhimane et al. CVPR06] 66
51 Non-Linear Optimization Non-linear least-squares minimization: min 2 A[R T]M i m i R,T i Minimization of a physical, meaningful error (reprojection error, in pixels) Handle an arbitrary number of points. Gauss-Newton, Levenberg Marquardt... minimization (next lecture). 67
Pose estimation from point and line correspondences
Pose estimation from point and line correspondences Giorgio Panin October 17, 008 1 Problem formulation Estimate (in a LSE sense) the pose of an object from N correspondences between known object points
More informationCamera Calibration The purpose of camera calibration is to determine the intrinsic camera parameters (c 0,r 0 ), f, s x, s y, skew parameter (s =
Camera Calibration The purpose of camera calibration is to determine the intrinsic camera parameters (c 0,r 0 ), f, s x, s y, skew parameter (s = cotα), and the lens distortion (radial distortion coefficient
More informationAugmented Reality VU numerical optimization algorithms. Prof. Vincent Lepetit
Augmented Reality VU numerical optimization algorithms Prof. Vincent Lepetit P3P: can exploit only 3 correspondences; DLT: difficult to exploit the knowledge of the internal parameters correctly, and the
More information3D from Photographs: Camera Calibration. Dr Francesco Banterle
3D from Photographs: Camera Calibration Dr Francesco Banterle francesco.banterle@isti.cnr.it 3D from Photographs Automatic Matching of Images Camera Calibration Photographs Surface Reconstruction Dense
More informationCamera calibration Triangulation
Triangulation Perspective projection in homogenous coordinates ~x img I 0 apple R t 0 T 1 ~x w ~x img R t ~x w Matrix transformations in 2D ~x img K R t ~x w K = 2 3 1 0 t u 40 1 t v 5 0 0 1 Translation
More informationTheory of Bouguet s MatLab Camera Calibration Toolbox
Theory of Bouguet s MatLab Camera Calibration Toolbox Yuji Oyamada 1 HVRL, University 2 Chair for Computer Aided Medical Procedure (CAMP) Technische Universität München June 26, 2012 MatLab Camera Calibration
More informationPose Tracking II! Gordon Wetzstein! Stanford University! EE 267 Virtual Reality! Lecture 12! stanford.edu/class/ee267/!
Pose Tracking II! Gordon Wetzstein! Stanford University! EE 267 Virtual Reality! Lecture 12! stanford.edu/class/ee267/!! WARNING! this class will be dense! will learn how to use nonlinear optimization
More informationOutline. Linear Algebra for Computer Vision
Outline Linear Algebra for Computer Vision Introduction CMSC 88 D Notation and Basics Motivation Linear systems of equations Gauss Elimination, LU decomposition Linear Spaces and Operators Addition, scalar
More informationAlgorithms for Computing a Planar Homography from Conics in Correspondence
Algorithms for Computing a Planar Homography from Conics in Correspondence Juho Kannala, Mikko Salo and Janne Heikkilä Machine Vision Group University of Oulu, Finland {jkannala, msa, jth@ee.oulu.fi} Abstract
More informationLecture 2: Linear Algebra Review
EE 227A: Convex Optimization and Applications January 19 Lecture 2: Linear Algebra Review Lecturer: Mert Pilanci Reading assignment: Appendix C of BV. Sections 2-6 of the web textbook 1 2.1 Vectors 2.1.1
More informationIntroduction to pinhole cameras
Introduction to pinhole cameras Lesson given by Sébastien Piérard in the course Introduction aux techniques audio et vidéo (ULg, Pr JJ Embrechts) INTELSIG, Montefiore Institute, University of Liège, Belgium
More informationLecture 4.3 Estimating homographies from feature correspondences. Thomas Opsahl
Lecture 4.3 Estimating homographies from feature correspondences Thomas Opsahl Homographies induced by central projection 1 H 2 1 H 2 u uu 2 3 1 Homography Hu = u H = h 1 h 2 h 3 h 4 h 5 h 6 h 7 h 8 h
More informationBasic Math for
Basic Math for 16-720 August 23, 2002 1 Linear Algebra 1.1 Vectors and Matrices First, a reminder of a few basic notations, definitions, and terminology: Unless indicated otherwise, vectors are always
More informationMulti-Frame Factorization Techniques
Multi-Frame Factorization Techniques Suppose { x j,n } J,N j=1,n=1 is a set of corresponding image coordinates, where the index n = 1,...,N refers to the n th scene point and j = 1,..., J refers to the
More informationCamera calibration by Zhang
Camera calibration by Zhang Siniša Kolarić September 2006 Abstract In this presentation, I present a way to calibrate a camera using the method by Zhang. NOTE. This
More informationHomographies and Estimating Extrinsics
Homographies and Estimating Extrinsics Instructor - Simon Lucey 16-423 - Designing Computer Vision Apps Adapted from: Computer vision: models, learning and inference. Simon J.D. Prince Review: Motivation
More informationLecture 5. Epipolar Geometry. Professor Silvio Savarese Computational Vision and Geometry Lab. 21-Jan-15. Lecture 5 - Silvio Savarese
Lecture 5 Epipolar Geometry Professor Silvio Savarese Computational Vision and Geometry Lab Silvio Savarese Lecture 5-21-Jan-15 Lecture 5 Epipolar Geometry Why is stereo useful? Epipolar constraints Essential
More informationUniversity of Houston, Department of Mathematics Numerical Analysis, Fall 2005
3 Numerical Solution of Nonlinear Equations and Systems 3.1 Fixed point iteration Reamrk 3.1 Problem Given a function F : lr n lr n, compute x lr n such that ( ) F(x ) = 0. In this chapter, we consider
More informationLines and points. Lines and points
omogeneous coordinates in the plane Homogeneous coordinates in the plane A line in the plane a + by + c is represented as (a, b, c). A line is a subset of points in the plane. All vectors (ka, kb, kc)
More informationSingular Value Decomposition
Chapter 6 Singular Value Decomposition In Chapter 5, we derived a number of algorithms for computing the eigenvalues and eigenvectors of matrices A R n n. Having developed this machinery, we complete our
More informationLinear Subspace Models
Linear Subspace Models Goal: Explore linear models of a data set. Motivation: A central question in vision concerns how we represent a collection of data vectors. The data vectors may be rasterized images,
More informationThe Multibody Trifocal Tensor: Motion Segmentation from 3 Perspective Views
The Multibody Trifocal Tensor: Motion Segmentation from 3 Perspective Views Richard Hartley 1,2 and RenéVidal 2,3 1 Dept. of Systems Engineering 3 Center for Imaging Science Australian National University
More informationLinear Algebra & Geometry why is linear algebra useful in computer vision?
Linear Algebra & Geometry why is linear algebra useful in computer vision? References: -Any book on linear algebra! -[HZ] chapters 2, 4 Some of the slides in this lecture are courtesy to Prof. Octavia
More information3D Computer Vision - WT 2004
3D Computer Vision - WT 2004 Singular Value Decomposition Darko Zikic CAMP - Chair for Computer Aided Medical Procedures November 4, 2004 1 2 3 4 5 Properties For any given matrix A R m n there exists
More informationMachine Learning. B. Unsupervised Learning B.2 Dimensionality Reduction. Lars Schmidt-Thieme, Nicolas Schilling
Machine Learning B. Unsupervised Learning B.2 Dimensionality Reduction Lars Schmidt-Thieme, Nicolas Schilling Information Systems and Machine Learning Lab (ISMLL) Institute for Computer Science University
More informationMultiple View Geometry in Computer Vision
Multiple View Geometry in Computer Vision Prasanna Sahoo Department of Mathematics University of Louisville 1 Scene Planes & Homographies Lecture 19 March 24, 2005 2 In our last lecture, we examined various
More informationLinear Algebra Massoud Malek
CSUEB Linear Algebra Massoud Malek Inner Product and Normed Space In all that follows, the n n identity matrix is denoted by I n, the n n zero matrix by Z n, and the zero vector by θ n An inner product
More informationA Study of Kruppa s Equation for Camera Self-calibration
Proceedings of the International Conference of Machine Vision and Machine Learning Prague, Czech Republic, August 14-15, 2014 Paper No. 57 A Study of Kruppa s Equation for Camera Self-calibration Luh Prapitasari,
More informationMobile Robotics 1. A Compact Course on Linear Algebra. Giorgio Grisetti
Mobile Robotics 1 A Compact Course on Linear Algebra Giorgio Grisetti SA-1 Vectors Arrays of numbers They represent a point in a n dimensional space 2 Vectors: Scalar Product Scalar-Vector Product Changes
More informationCS 231A Computer Vision (Fall 2011) Problem Set 2
CS 231A Computer Vision (Fall 2011) Problem Set 2 Solution Set Due: Oct. 28 th, 2011 (5pm) 1 Some Projective Geometry Problems (15 points) Suppose there are two parallel lines that extend to infinity in
More informationLinear Algebra for Machine Learning. Sargur N. Srihari
Linear Algebra for Machine Learning Sargur N. srihari@cedar.buffalo.edu 1 Overview Linear Algebra is based on continuous math rather than discrete math Computer scientists have little experience with it
More informationAffine invariant Fourier descriptors
Affine invariant Fourier descriptors Sought: a generalization of the previously introduced similarityinvariant Fourier descriptors H. Burkhardt, Institut für Informatik, Universität Freiburg ME-II, Kap.
More informationPreliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012
Instructions Preliminary/Qualifying Exam in Numerical Analysis (Math 502a) Spring 2012 The exam consists of four problems, each having multiple parts. You should attempt to solve all four problems. 1.
More information7. Dimension and Structure.
7. Dimension and Structure 7.1. Basis and Dimension Bases for Subspaces Example 2 The standard unit vectors e 1, e 2,, e n are linearly independent, for if we write (2) in component form, then we obtain
More information13. Nonlinear least squares
L. Vandenberghe ECE133A (Fall 2018) 13. Nonlinear least squares definition and examples derivatives and optimality condition Gauss Newton method Levenberg Marquardt method 13.1 Nonlinear least squares
More informationNumerical computation II. Reprojection error Bundle adjustment Family of Newtonʼs methods Statistical background Maximum likelihood estimation
Numerical computation II Reprojection error Bundle adjustment Family of Newtonʼs methods Statistical background Maximum likelihood estimation Reprojection error Reprojection error = Distance between the
More informationLinear Algebra in Computer Vision. Lecture2: Basic Linear Algebra & Probability. Vector. Vector Operations
Linear Algebra in Computer Vision CSED441:Introduction to Computer Vision (2017F Lecture2: Basic Linear Algebra & Probability Bohyung Han CSE, POSTECH bhhan@postech.ac.kr Mathematics in vector space Linear
More informationCSE 554 Lecture 7: Alignment
CSE 554 Lecture 7: Alignment Fall 2012 CSE554 Alignment Slide 1 Review Fairing (smoothing) Relocating vertices to achieve a smoother appearance Method: centroid averaging Simplification Reducing vertex
More informationLinear Algebra Review. Vectors
Linear Algebra Review 9/4/7 Linear Algebra Review By Tim K. Marks UCSD Borrows heavily from: Jana Kosecka http://cs.gmu.edu/~kosecka/cs682.html Virginia de Sa (UCSD) Cogsci 8F Linear Algebra review Vectors
More informationOPPA European Social Fund Prague & EU: We invest in your future.
OPPA European Social Fund Prague & EU: We invest in your future. 1D Projective Coordinates The 1-D projective coordinate of a point P : [P ] = [P P 0 P I P ] = [p p 0 p I p] = p pi p 0 p I p 0 p p p p
More informationMath Review: parameter estimation. Emma
Math Review: parameter estimation Emma McNally@flickr Fitting lines to dots: We will cover how Slides provided by HyunSoo Park 1809, Carl Friedrich Gauss What about fitting line on a curved surface? Least
More informationSparse Levenberg-Marquardt algorithm.
Sparse Levenberg-Marquardt algorithm. R. I. Hartley and A. Zisserman: Multiple View Geometry in Computer Vision. Cambridge University Press, second edition, 2004. Appendix 6 was used in part. The Levenberg-Marquardt
More informationLinear Algebra (Review) Volker Tresp 2017
Linear Algebra (Review) Volker Tresp 2017 1 Vectors k is a scalar (a number) c is a column vector. Thus in two dimensions, c = ( c1 c 2 ) (Advanced: More precisely, a vector is defined in a vector space.
More informationMathematical Properties of Stiffness Matrices
Mathematical Properties of Stiffness Matrices CEE 4L. Matrix Structural Analysis Department of Civil and Environmental Engineering Duke University Henri P. Gavin Fall, 0 These notes describe some of the
More informationExercise Sheet 1.
Exercise Sheet 1 You can download my lecture and exercise sheets at the address http://sami.hust.edu.vn/giang-vien/?name=huynt 1) Let A, B be sets. What does the statement "A is not a subset of B " mean?
More informationLINEAR ALGEBRA W W L CHEN
LINEAR ALGEBRA W W L CHEN c W W L Chen, 1997, 2008. This chapter is available free to all individuals, on the understanding that it is not to be used for financial gain, and may be downloaded and/or photocopied,
More informationM3: Multiple View Geometry
M3: Multiple View Geometry L18: Projective Structure from Motion: Iterative Algorithm based on Factorization Based on Sections 13.4 C. V. Jawahar jawahar-at-iiit.net Mar 2005: 1 Review: Reconstruction
More informationCSE 252B: Computer Vision II
CSE 252B: Computer Vision II Lecturer: Serge Belongie Scribe: Tasha Vanesian LECTURE 3 Calibrated 3D Reconstruction 3.1. Geometric View of Epipolar Constraint We are trying to solve the following problem:
More informationLinear Algebra & Geometry why is linear algebra useful in computer vision?
Linear Algebra & Geometry why is linear algebra useful in computer vision? References: -Any book on linear algebra! -[HZ] chapters 2, 4 Some of the slides in this lecture are courtesy to Prof. Octavia
More informationEE/ACM Applications of Convex Optimization in Signal Processing and Communications Lecture 17
EE/ACM 150 - Applications of Convex Optimization in Signal Processing and Communications Lecture 17 Andre Tkacenko Signal Processing Research Group Jet Propulsion Laboratory May 29, 2012 Andre Tkacenko
More informationA Theory of Multi-Layer Flat Refractive Geometry Supplementary Materials
A Theory of Multi-Layer Flat Refractive Geometry Supplementary Materials Amit Agrawal, Srikumar Ramalingam, Yuichi Taguchi Mitsubishi Electric Research Labs MERL) [agrawal,ramalingam,taguchi] at merl.com
More informationCamera Projection Model
amera Projection Model 3D Point Projection (Pixel Space) O ( u, v ) ccd ccd f Projection plane (, Y, Z ) ( u, v ) img Z img (0,0) w img ( u, v ) img img uimg f Z p Y v img f Z p x y Zu f p Z img img x
More informationMysteries of Parameterizing Camera Motion - Part 2
Mysteries of Parameterizing Camera Motion - Part 2 Instructor - Simon Lucey 16-623 - Advanced Computer Vision Apps Today Gauss-Newton and SO(3) Alternative Representations of SO(3) Exponential Maps SL(3)
More informationThere are six more problems on the next two pages
Math 435 bg & bu: Topics in linear algebra Summer 25 Final exam Wed., 8/3/5. Justify all your work to receive full credit. Name:. Let A 3 2 5 Find a permutation matrix P, a lower triangular matrix L with
More informationLecture: Linear algebra. 4. Solutions of linear equation systems The fundamental theorem of linear algebra
Lecture: Linear algebra. 1. Subspaces. 2. Orthogonal complement. 3. The four fundamental subspaces 4. Solutions of linear equation systems The fundamental theorem of linear algebra 5. Determining the fundamental
More informationReview of Linear Algebra
Review of Linear Algebra Definitions An m n (read "m by n") matrix, is a rectangular array of entries, where m is the number of rows and n the number of columns. 2 Definitions (Con t) A is square if m=
More informationAssignment 1 Math 5341 Linear Algebra Review. Give complete answers to each of the following questions. Show all of your work.
Assignment 1 Math 5341 Linear Algebra Review Give complete answers to each of the following questions Show all of your work Note: You might struggle with some of these questions, either because it has
More informationLecture 20: 6.1 Inner Products
Lecture 0: 6.1 Inner Products Wei-Ta Chu 011/1/5 Definition An inner product on a real vector space V is a function that associates a real number u, v with each pair of vectors u and v in V in such a way
More informationAPPLICATIONS The eigenvalues are λ = 5, 5. An orthonormal basis of eigenvectors consists of
CHAPTER III APPLICATIONS The eigenvalues are λ =, An orthonormal basis of eigenvectors consists of, The eigenvalues are λ =, A basis of eigenvectors consists of, 4 which are not perpendicular However,
More informationMultivariate Statistical Analysis
Multivariate Statistical Analysis Fall 2011 C. L. Williams, Ph.D. Lecture 4 for Applied Multivariate Analysis Outline 1 Eigen values and eigen vectors Characteristic equation Some properties of eigendecompositions
More informationUncertainty Models in Quasiconvex Optimization for Geometric Reconstruction
Uncertainty Models in Quasiconvex Optimization for Geometric Reconstruction Qifa Ke and Takeo Kanade Department of Computer Science, Carnegie Mellon University Email: ke@cmu.edu, tk@cs.cmu.edu Abstract
More information6.801/866. Affine Structure from Motion. T. Darrell
6.801/866 Affine Structure from Motion T. Darrell [Read F&P Ch. 12.0, 12.2, 12.3, 12.4] Affine geometry is, roughly speaking, what is left after all ability to measure lengths, areas, angles, etc. has
More informationOptimisation on Manifolds
Optimisation on Manifolds K. Hüper MPI Tübingen & Univ. Würzburg K. Hüper (MPI Tübingen & Univ. Würzburg) Applications in Computer Vision Grenoble 18/9/08 1 / 29 Contents 2 Examples Essential matrix estimation
More informationElementary linear algebra
Chapter 1 Elementary linear algebra 1.1 Vector spaces Vector spaces owe their importance to the fact that so many models arising in the solutions of specific problems turn out to be vector spaces. The
More informationTwo-View Segmentation of Dynamic Scenes from the Multibody Fundamental Matrix
Two-View Segmentation of Dynamic Scenes from the Multibody Fundamental Matrix René Vidal Stefano Soatto Shankar Sastry Department of EECS, UC Berkeley Department of Computer Sciences, UCLA 30 Cory Hall,
More informationMath 302 Outcome Statements Winter 2013
Math 302 Outcome Statements Winter 2013 1 Rectangular Space Coordinates; Vectors in the Three-Dimensional Space (a) Cartesian coordinates of a point (b) sphere (c) symmetry about a point, a line, and a
More informationMath Camp Lecture 4: Linear Algebra. Xiao Yu Wang. Aug 2010 MIT. Xiao Yu Wang (MIT) Math Camp /10 1 / 88
Math Camp 2010 Lecture 4: Linear Algebra Xiao Yu Wang MIT Aug 2010 Xiao Yu Wang (MIT) Math Camp 2010 08/10 1 / 88 Linear Algebra Game Plan Vector Spaces Linear Transformations and Matrices Determinant
More informationEXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. 1. Determinants
EXERCISES ON DETERMINANTS, EIGENVALUES AND EIGENVECTORS. Determinants Ex... Let A = 0 4 4 2 0 and B = 0 3 0. (a) Compute 0 0 0 0 A. (b) Compute det(2a 2 B), det(4a + B), det(2(a 3 B 2 )). 0 t Ex..2. For
More informationMathematical foundations - linear algebra
Mathematical foundations - linear algebra Andrea Passerini passerini@disi.unitn.it Machine Learning Vector space Definition (over reals) A set X is called a vector space over IR if addition and scalar
More informationThroughout these notes we assume V, W are finite dimensional inner product spaces over C.
Math 342 - Linear Algebra II Notes Throughout these notes we assume V, W are finite dimensional inner product spaces over C 1 Upper Triangular Representation Proposition: Let T L(V ) There exists an orthonormal
More informationIntro Vectors 2D implicit curves 2D parametric curves. Graphics 2011/2012, 4th quarter. Lecture 2: vectors, curves, and surfaces
Lecture 2, curves, and surfaces Organizational remarks Tutorials: Tutorial 1 will be online later today TA sessions for questions start next week Practicals: Exams: Make sure to find a team partner very
More informationMachine Learning. Principal Components Analysis. Le Song. CSE6740/CS7641/ISYE6740, Fall 2012
Machine Learning CSE6740/CS7641/ISYE6740, Fall 2012 Principal Components Analysis Le Song Lecture 22, Nov 13, 2012 Based on slides from Eric Xing, CMU Reading: Chap 12.1, CB book 1 2 Factor or Component
More informationCOMP 558 lecture 18 Nov. 15, 2010
Least squares We have seen several least squares problems thus far, and we will see more in the upcoming lectures. For this reason it is good to have a more general picture of these problems and how to
More informationCamera Models and Affine Multiple Views Geometry
Camera Models and Affine Multiple Views Geometry Subhashis Banerjee Dept. Computer Science and Engineering IIT Delhi email: suban@cse.iitd.ac.in May 29, 2001 1 1 Camera Models A Camera transforms a 3D
More informationProbabilistic Latent Semantic Analysis
Probabilistic Latent Semantic Analysis Seungjin Choi Department of Computer Science and Engineering Pohang University of Science and Technology 77 Cheongam-ro, Nam-gu, Pohang 37673, Korea seungjin@postech.ac.kr
More informationA geometric interpretation of the homogeneous coordinates is given in the following Figure.
Introduction Homogeneous coordinates are an augmented representation of points and lines in R n spaces, embedding them in R n+1, hence using n + 1 parameters. This representation is useful in dealing with
More informationExtra Problems for Math 2050 Linear Algebra I
Extra Problems for Math 5 Linear Algebra I Find the vector AB and illustrate with a picture if A = (,) and B = (,4) Find B, given A = (,4) and [ AB = A = (,4) and [ AB = 8 If possible, express x = 7 as
More informationEXERCISE SET 5.1. = (kx + kx + k, ky + ky + k ) = (kx + kx + 1, ky + ky + 1) = ((k + )x + 1, (k + )y + 1)
EXERCISE SET 5. 6. The pair (, 2) is in the set but the pair ( )(, 2) = (, 2) is not because the first component is negative; hence Axiom 6 fails. Axiom 5 also fails. 8. Axioms, 2, 3, 6, 9, and are easily
More informationLinear Algebra. Session 12
Linear Algebra. Session 12 Dr. Marco A Roque Sol 08/01/2017 Example 12.1 Find the constant function that is the least squares fit to the following data x 0 1 2 3 f(x) 1 0 1 2 Solution c = 1 c = 0 f (x)
More informationIntroduction to Mobile Robotics Compact Course on Linear Algebra. Wolfram Burgard, Cyrill Stachniss, Kai Arras, Maren Bennewitz
Introduction to Mobile Robotics Compact Course on Linear Algebra Wolfram Burgard, Cyrill Stachniss, Kai Arras, Maren Bennewitz Vectors Arrays of numbers Vectors represent a point in a n dimensional space
More informationPCA & ICA. CE-717: Machine Learning Sharif University of Technology Spring Soleymani
PCA & ICA CE-717: Machine Learning Sharif University of Technology Spring 2015 Soleymani Dimensionality Reduction: Feature Selection vs. Feature Extraction Feature selection Select a subset of a given
More informationOptimization Methods. Lecture 23: Semidenite Optimization
5.93 Optimization Methods Lecture 23: Semidenite Optimization Outline. Preliminaries Slide 2. SDO 3. Duality 4. SDO Modeling Power 5. Barrier lgorithm for SDO 2 Preliminaries symmetric matrix is positive
More informationLecture 7 Spectral methods
CSE 291: Unsupervised learning Spring 2008 Lecture 7 Spectral methods 7.1 Linear algebra review 7.1.1 Eigenvalues and eigenvectors Definition 1. A d d matrix M has eigenvalue λ if there is a d-dimensional
More informationMath Introduction to Numerical Analysis - Class Notes. Fernando Guevara Vasquez. Version Date: January 17, 2012.
Math 5620 - Introduction to Numerical Analysis - Class Notes Fernando Guevara Vasquez Version 1990. Date: January 17, 2012. 3 Contents 1. Disclaimer 4 Chapter 1. Iterative methods for solving linear systems
More informationLinear Algebra and Eigenproblems
Appendix A A Linear Algebra and Eigenproblems A working knowledge of linear algebra is key to understanding many of the issues raised in this work. In particular, many of the discussions of the details
More information15 Singular Value Decomposition
15 Singular Value Decomposition For any high-dimensional data analysis, one s first thought should often be: can I use an SVD? The singular value decomposition is an invaluable analysis tool for dealing
More informationAPPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.
APPENDIX A Background Mathematics A. Linear Algebra A.. Vector algebra Let x denote the n-dimensional column vector with components 0 x x 2 B C @. A x n Definition 6 (scalar product). The scalar product
More informationLinear Algebra in Actuarial Science: Slides to the lecture
Linear Algebra in Actuarial Science: Slides to the lecture Fall Semester 2010/2011 Linear Algebra is a Tool-Box Linear Equation Systems Discretization of differential equations: solving linear equations
More informationLecture II: Linear Algebra Revisited
Lecture II: Linear Algebra Revisited Overview Vector spaces, Hilbert & Banach Spaces, etrics & Norms atrices, Eigenvalues, Orthogonal Transformations, Singular Values Operators, Operator Norms, Function
More information1 9/5 Matrices, vectors, and their applications
1 9/5 Matrices, vectors, and their applications Algebra: study of objects and operations on them. Linear algebra: object: matrices and vectors. operations: addition, multiplication etc. Algorithms/Geometric
More informationHomogeneous Coordinates
Homogeneous Coordinates Basilio Bona DAUIN-Politecnico di Torino October 2013 Basilio Bona (DAUIN-Politecnico di Torino) Homogeneous Coordinates October 2013 1 / 32 Introduction Homogeneous coordinates
More informationMultiple View Geometry in Computer Vision
Multiple View Geometry in Computer Vision Prasanna Sahoo Department of Mathematics University of Louisville 1 Basic Information Instructor: Professor Ron Sahoo Office: NS 218 Tel: (502) 852-2731 Fax: (502)
More informationDuke University, Department of Electrical and Computer Engineering Optimization for Scientists and Engineers c Alex Bronstein, 2014
Duke University, Department of Electrical and Computer Engineering Optimization for Scientists and Engineers c Alex Bronstein, 2014 Linear Algebra A Brief Reminder Purpose. The purpose of this document
More informationAdvances in Computer Vision. Prof. Bill Freeman. Image and shape descriptors. Readings: Mikolajczyk and Schmid; Belongie et al.
6.869 Advances in Computer Vision Prof. Bill Freeman March 3, 2005 Image and shape descriptors Affine invariant features Comparison of feature descriptors Shape context Readings: Mikolajczyk and Schmid;
More information14 Singular Value Decomposition
14 Singular Value Decomposition For any high-dimensional data analysis, one s first thought should often be: can I use an SVD? The singular value decomposition is an invaluable analysis tool for dealing
More informationAlgebra C Numerical Linear Algebra Sample Exam Problems
Algebra C Numerical Linear Algebra Sample Exam Problems Notation. Denote by V a finite-dimensional Hilbert space with inner product (, ) and corresponding norm. The abbreviation SPD is used for symmetric
More information0, otherwise. Find each of the following limits, or explain that the limit does not exist.
Midterm Solutions 1, y x 4 1. Let f(x, y) = 1, y 0 0, otherwise. Find each of the following limits, or explain that the limit does not exist. (a) (b) (c) lim f(x, y) (x,y) (0,1) lim f(x, y) (x,y) (2,3)
More informationCorners, Blobs & Descriptors. With slides from S. Lazebnik & S. Seitz, D. Lowe, A. Efros
Corners, Blobs & Descriptors With slides from S. Lazebnik & S. Seitz, D. Lowe, A. Efros Motivation: Build a Panorama M. Brown and D. G. Lowe. Recognising Panoramas. ICCV 2003 How do we build panorama?
More informationVision par ordinateur
Vision par ordinateur Géométrie épipolaire Frédéric Devernay Avec des transparents de Marc Pollefeys Epipolar geometry π Underlying structure in set of matches for rigid scenes C1 m1 l1 M L2 L1 l T 1 l
More informationCamera calibration. Outline. Pinhole camera. Camera projection models. Nonlinear least square methods A camera calibration tool
Outline Camera calibration Camera projection models Camera calibration i Nonlinear least square methods A camera calibration tool Applications Digital Visual Effects Yung-Yu Chuang with slides b Richard
More information