Single View Geometry Camera Model & Orientation + Position Estimation

Total Page:16

File Type:pdf, Size:1020Kb

Single View Geometry Camera Model & Orientation + Position Estimation Single View Geometry Camera model & Orientation + Position estimation What am I? Vanishing point Mapping from 3D to 2D Point & Line Goal: Homogeneous coordinates Point – represent coordinates in 2 dimensions with a 3-vector &x# &x# homogeneous coords $y! $ ! ''''''→$ ! %y" %$1"! The projective plane • Why do we need homogeneous coordinates? – represent points at infinity, homographies, perspective projection, multi-view relationships • What is the geometric intuition? – a point in the image is a ray in projective space y (sx,sy,s) (x,y,1) (0,0,0) z x image plane • Each point (x,y) on the plane is represented by a ray (sx,sy,s) – all points on the ray are equivalent: (x, y, 1) ≅ (sx, sy, s) Projective Lines Projective lines • What does a line in the image correspond to in projective space? • A line is a plane of rays through origin – all rays (x,y,z) satisfying: ax + by + cz = 0 ⎡x⎤ in vector notation : 0 a b c ⎢y⎥ = [ ]⎢ ⎥ ⎣⎢z⎦⎥ l p • A line is also represented as a homogeneous 3-vector l Line Representation • a line is • is the distance from the origin to the line • is the norm direction of the line • It can also be written as Example of Line Example of Line (2) 0.42*pi Homogeneous representation Line in Is represented by a point in : But correspondence of line to point is not unique We define set of equivalence class of vectors in R^3 - (0,0,0) As projective space Projective lines from two points Line passing through two points Two points: x Define a line l is the line passing two points Proof: Line passing through two points • More specifically, Matlab codes • function l = get_line_by_two_points(x, y) • x1 = [x(1), y(1), 1]'; • x2 = [x(2), y(2), 1]'; • l = cross(x1, x2); • l = l / sqrt(l(1)*l(1)+l(2)*l(2)); Example of Line Verify p2 lies on the line P2 >> 218*0.6+13*0.8 ans = 141.2000 (close) P1 Example of Line P2 Test line: >> p1 = [42,142,1]; >> p2=[218,13,1]; >> l = cross(p1,p2) l = 129 176 -30410 P1 >> l = l/sqrt(l(1)^2+l(2)^2) l = 0.5912 0.8066 -139.3591 >> x = [221;23]; >> y = [218;268]; >> l = get_line_by_two_points(x,y); >> l l = -0.2448 -0.9696 265.4744 Projective lines from two points When does the line has the form (a, b, 0)? -all lines passing through the optical center When does the line has the form (0, 0, 1)? -line at infinity Points from two lines Intersection of lines Given two lines: l , Define a point X is the intersection of the two lines Intersection of lines • More specifically, Matlab codes • function x0 = get_point_by_two_line(l, l1) • x0 = cross(l, l1); • x0 = [x0(1)/x0(3); x0(2)/x0(3)]; Example of Lines Intersection Point and line duality – A line l is a homogeneous 3-vector – It is ⊥ to every point (ray) p on the line: l p=0 p p2 l1 p l 1 l2 What is the line l spanned by rays p1 and p2 ? • l is ⊥ to p1 and p2 ⇒ l = p1 × p2 • l is the plane normal What is the intersection of two lines l1 and l2 ? • p is ⊥ to l1 and l2 ⇒ p = l1 × l2 Points and lines are dual in projective space • given any formula, can switch the meanings of points and lines to get another formula Points from two lines When P has the form (x,y,0)? Point at infinity Example: Consider two parallel horizontal lines: Intersection = Point at infinity in the direction of y Point at infinity, Ideal points Intersection: Any point (x1,x2,0) is intersection of lines at infinity Points at infinity • Under projective transformation, – All parallel lines intersects at the point at infinity – One point at infinity ó one parallel line direction • Where are the points at infinity in the image plane? Line at infinity • A line passing all points at infinity: • Because : Ideal points and lines (a,b,0) y (sx,sy,0) y z image plane z image plane x x • Ideal point (“point at infinity”) – p ≅ (x, y, 0) – parallel to image plane – It has infinite image coordinates Ideal line • l ≅ (a, b, 0) – parallel to image plane • Corresponds to a line in the image (finite coordinates) Comparing heights Vanishing Point Measuring height 5.4 5 Camera height 4 3.3 3 2.8 2 1 Measuring height vz r vanishing line (horizon) t0 t vx v vy H R H b0 b Focal Length • A Lighting Reproduction Approach to Live-Action Compositing, P. Debevec 2002 FOV depends of Focal Length f f Smaller fOV = larger focal Length Field of View (Zoom) Field of View (Zoom) Large Focal Length compresses depth 400 mm 200 mm 100 mm 50 mm 28 mm 17 mm © 1995-2005 Michael Reichmann http://www.creativepro.com/img/story/102604_fg3b.jpg Nikon’s 1200-1700mm f/5.6-8P lens. Nikon’s 1200-1700mm f/5.6-8P lens. " /012'+$3+/4$&'()5$2 35'0:+$3 +>5'?+!6 ;<=6 $%&'()*+,+-+. !! ;<=66 36+!+3$(10+0'27)8 # 366+!+3$(10+0'27)8 9.+!+:5*)12('+)$+$%&'()+. 9,+!+:5*)12('+)$+$%&'()+, Canon D30 takes image with 2160x1440 pixels. What is the pixel size? What f to make tower appear half the size of image? Fix f=50mm, D2 = ? in order to make the tower 1/3 size of image? D1=? to make human the same height as tower? What f’ to make tower appear twice as big as before (2/3 the height)? Setting f’=100mm, D1=? in order to keep the original size of human? (1/3 the height) Lens distortion with focal length Lens distortion with focal length Camera projection model The overall goal is to compute 3D geometry of the scene from just 2D images. We will first study 3D->2D camera projection using linear algebra. The Pinhole Camera – Light enters a darkened chamber through a pinhole opening and forms an image on the further surface The Pinhole Camera Model Ideal case: yc y Projection equation: y’ x’ = f X / Z p C y’ = f Y / Z z zc f xc Zx’ = f X Zy’ = f Y Z = Z Step 1: Camera projection matrix Zx’ = f X Zy’ = f Y y Z = Z c y y’ p C f 0 0 0 Xc x’ z Z zc f 0 f 0 0 Yc = y’ c Z 0 0 1 0 c 1 x 1 c P0 X = x Step 2: Intrinsic camera parameters: map camera coordinate to pixel coordinate yc αx s px x’ x 0 p y’ y αy y = y 0 0 1 z’ z y’ p C z zc f K Optical world Pixel world (3x3 submatrix) xc αx, αy is pixel scaling factor px, py is the principle point (where optical axis hits image plane) s is the slant factor, when the image plane is not normal to the optical axis Combine the Intrinsic camera parameters X αx s px f 0 0 0 c x Y 0 αy py 0 f 0 0 c = y 0 0 1 0 0 1 0 Zc 1 1 x αxf s px 0 c x y 0 αyf py 0 c = y 0 0 1 0 zc z 1 K (Calibration matrix) P0 X = x P = [K , 0] = K [I, 0] Step 3: External parameters: rotation and translation map the world to camera coordinates 3x3 3x1 (R,t) Camera coordinate yc World coordinate y y’ p C z zc f xc Combining Internal and External parameters (R,t) Camera coordinate 1) yc World coordinate y y’ p C z zc f 2) xc 1) Translate the world coordinate into the camera coordinate 2) Translate the camera coordinate into the pixel coordinate Combining Internal and External parameters (R,t) Camera coordinate yc World coordinate y y’ p C z After simplication: zc f xc x = K [R, t] X pixel world Special case, planar world, homograph (R,t) X= (xw,yw,zw=0,1) yc World coordinate y y’ p C x = K [R, t] X z zc f Expand: xc xw Camera coordinate yw x = K r1 r2 r3 t zw=0 1 Special case, planar world, homograph (R,t) X= (xw,yw,zw=0,1) yc World coordinate y After simplication: y’ p C z zc f w x xc x = K r1 r2 t yw Camera coordinate 1 We have the homographic mapping: X = K H3x3 X Special case: rotating camera (R,t) x = K [R, t] X yc World coordinate with t = 0 y y’ p C Expand: z x = K [R] X zc f xc Camera coordinate xw w This is also a homography with X = y H = R zw Panoramas x = K [R] X Visual Odometry: Overview 1) The transformation (R,t) tells us (R,t) how to rotate + move the world coordinate to the camera yc World coordinate. coordinate y 2) In the following slides, we will see y’ p C z 1) how to get the world reference’s zc f frame orientation+position from (R,t), xc 2) how to get the pan/tilt/yaw angles from R 3) If we want to know the camera orientation+position, we need to look at its inverse transform, which we should see is , 1) The computation of the camera pan/tilt/yaw angle is the same as in the previous case Recover Camera orientation Case 1: single vanishing point (z- direction) Pan/tilt/yaw angles Pan : rotation Yaw :rotation Tilt : rotation around y-axis around z-axis around x-axis z y y x x z Rotation: Pan/Tilt/Yaw Given the 3x3 rotation matrix R, we can recover, pan/tilt/yaw angle, and vise.
Recommended publications
  • Projective Geometry: a Short Introduction
    Projective Geometry: A Short Introduction Lecture Notes Edmond Boyer Master MOSIG Introduction to Projective Geometry Contents 1 Introduction 2 1.1 Objective . .2 1.2 Historical Background . .3 1.3 Bibliography . .4 2 Projective Spaces 5 2.1 Definitions . .5 2.2 Properties . .8 2.3 The hyperplane at infinity . 12 3 The projective line 13 3.1 Introduction . 13 3.2 Projective transformation of P1 ................... 14 3.3 The cross-ratio . 14 4 The projective plane 17 4.1 Points and lines . 17 4.2 Line at infinity . 18 4.3 Homographies . 19 4.4 Conics . 20 4.5 Affine transformations . 22 4.6 Euclidean transformations . 22 4.7 Particular transformations . 24 4.8 Transformation hierarchy . 25 Grenoble Universities 1 Master MOSIG Introduction to Projective Geometry Chapter 1 Introduction 1.1 Objective The objective of this course is to give basic notions and intuitions on projective geometry. The interest of projective geometry arises in several visual comput- ing domains, in particular computer vision modelling and computer graphics. It provides a mathematical formalism to describe the geometry of cameras and the associated transformations, hence enabling the design of computational ap- proaches that manipulates 2D projections of 3D objects. In that respect, a fundamental aspect is the fact that objects at infinity can be represented and manipulated with projective geometry and this in contrast to the Euclidean geometry. This allows perspective deformations to be represented as projective transformations. Figure 1.1: Example of perspective deformation or 2D projective transforma- tion. Another argument is that Euclidean geometry is sometimes difficult to use in algorithms, with particular cases arising from non-generic situations (e.g.
    [Show full text]
  • Algebraic Geometry Michael Stoll
    Introductory Geometry Course No. 100 351 Fall 2005 Second Part: Algebraic Geometry Michael Stoll Contents 1. What Is Algebraic Geometry? 2 2. Affine Spaces and Algebraic Sets 3 3. Projective Spaces and Algebraic Sets 6 4. Projective Closure and Affine Patches 9 5. Morphisms and Rational Maps 11 6. Curves — Local Properties 14 7. B´ezout’sTheorem 18 2 1. What Is Algebraic Geometry? Linear Algebra can be seen (in parts at least) as the study of systems of linear equations. In geometric terms, this can be interpreted as the study of linear (or affine) subspaces of Cn (say). Algebraic Geometry generalizes this in a natural way be looking at systems of polynomial equations. Their geometric realizations (their solution sets in Cn, say) are called algebraic varieties. Many questions one can study in various parts of mathematics lead in a natural way to (systems of) polynomial equations, to which the methods of Algebraic Geometry can be applied. Algebraic Geometry provides a translation between algebra (solutions of equations) and geometry (points on algebraic varieties). The methods are mostly algebraic, but the geometry provides the intuition. Compared to Differential Geometry, in Algebraic Geometry we consider a rather restricted class of “manifolds” — those given by polynomial equations (we can allow “singularities”, however). For example, y = cos x defines a perfectly nice differentiable curve in the plane, but not an algebraic curve. In return, we can get stronger results, for example a criterion for the existence of solutions (in the complex numbers), or statements on the number of solutions (for example when intersecting two curves), or classification results.
    [Show full text]
  • Projective Geometry Lecture Notes
    Projective Geometry Lecture Notes Thomas Baird March 26, 2014 Contents 1 Introduction 2 2 Vector Spaces and Projective Spaces 4 2.1 Vector spaces and their duals . 4 2.1.1 Fields . 4 2.1.2 Vector spaces and subspaces . 5 2.1.3 Matrices . 7 2.1.4 Dual vector spaces . 7 2.2 Projective spaces and homogeneous coordinates . 8 2.2.1 Visualizing projective space . 8 2.2.2 Homogeneous coordinates . 13 2.3 Linear subspaces . 13 2.3.1 Two points determine a line . 14 2.3.2 Two planar lines intersect at a point . 14 2.4 Projective transformations and the Erlangen Program . 15 2.4.1 Erlangen Program . 16 2.4.2 Projective versus linear . 17 2.4.3 Examples of projective transformations . 18 2.4.4 Direct sums . 19 2.4.5 General position . 20 2.4.6 The Cross-Ratio . 22 2.5 Classical Theorems . 23 2.5.1 Desargues' Theorem . 23 2.5.2 Pappus' Theorem . 24 2.6 Duality . 26 3 Quadrics and Conics 28 3.1 Affine algebraic sets . 28 3.2 Projective algebraic sets . 30 3.3 Bilinear and quadratic forms . 31 3.3.1 Quadratic forms . 33 3.3.2 Change of basis . 33 1 3.3.3 Digression on the Hessian . 36 3.4 Quadrics and conics . 37 3.5 Parametrization of the conic . 40 3.5.1 Rational parametrization of the circle . 42 3.6 Polars . 44 3.7 Linear subspaces of quadrics and ruled surfaces . 46 3.8 Pencils of quadrics and degeneration . 47 4 Exterior Algebras 52 4.1 Multilinear algebra .
    [Show full text]
  • Chapter 12 the Cross Ratio
    Chapter 12 The cross ratio Math 4520, Fall 2017 We have studied the collineations of a projective plane, the automorphisms of the underlying field, the linear functions of Affine geometry, etc. We have been led to these ideas by various problems at hand, but let us step back and take a look at one important point of view of the big picture. 12.1 Klein's Erlanger program In 1872, Felix Klein, one of the leading mathematicians and geometers of his day, in the city of Erlanger, took the following point of view as to what the role of geometry was in mathematics. This is from his \Recent trends in geometric research." Let there be given a manifold and in it a group of transforma- tions; it is our task to investigate those properties of a figure belonging to the manifold that are not changed by the transfor- mation of the group. So our purpose is clear. Choose the group of transformations that you are interested in, and then hunt for the \invariants" that are relevant. This search for invariants has proved very fruitful and useful since the time of Klein for many areas of mathematics, not just classical geometry. In some case the invariants have turned out to be simple polynomials or rational functions, such as the case with the cross ratio. In other cases the invariants were groups themselves, such as homology groups in the case of algebraic topology. 12.2 The projective line In Chapter 11 we saw that the collineations of a projective plane come in two \species," projectivities and field automorphisms.
    [Show full text]
  • Where Parallel Lines Meet
    Where Parallel Lines WhereMeet parallel lines meet ...a geometric love story Renzo Cavalieri Renzo Cavalieri Colorado State University MATH DAY 2009 University of Northern Colorado Oct 22, 2014 Renzo Cavalieri Where parallel lines meet... Renzo Cavalieri Where Parallel Lines Meet Geometry γ"!µετρια \Geos = earth" + \Metron = measure" Was developed by King Sesostris as a way to counter tax fraud! Renzo Cavalieri Where Parallel Lines Meet Abstraction Abstraction In order to study in broader generalityIn order to the study properties in broader of spaces, __________ mathematiciansgenerality the properties created concepts of spaces, thatmathematicians abstract and created generalize concepts our physicalthat abstract experience. and generalize our Forphysical example, experience. a line or a plane do not reallyFor example, exist in our a line world!or a plane do not really exist in our world! Renzo Cavalieri Where Parallel Lines Meet Renzo Cavalieri Where Parallel Lines Meet Euclid's Posulates 1 A straight line segment can be drawn joining any two points. 2 Any straight line segment can be extended indefinitely in a straight line. 3 Given any straight line segment, a circle can be drawn having the segment as radius and one endpoint as center. 4 All right angles are congruent. 5 Given a line m and a point P not belonging to m, there is a unique line through P that never intersects m. Renzo Cavalieri Where Parallel Lines Meet The parallel postulate The parallel postulate was first shaken by the following puzzle: A hunter gets off her pickup and walks south for one mile. Then she makes a 90◦ left turn and walks east for another mile.
    [Show full text]
  • 6. the Riemann Sphere It Is Sometimes Convenient to Add a Point at Infinity ∞ to the Usual Complex Plane to Get the Extended Complex Plane
    6. The Riemann sphere It is sometimes convenient to add a point at infinity 1 to the usual complex plane to get the extended complex plane. Definition 6.1. The extended complex plane, denoted P1, is simply the union of C and the point at infinity. It is somewhat curious that when we add points at infinity to the reals we add two points ±∞ but only only one point for the complex numbers. It is rare in geometry that things get easier as you increase the dimension. One very good way to understand the extended complex plane is to realise that P1 is naturally in bijection with the unit sphere: Definition 6.2. The Riemann sphere is the unit sphere in R3: 2 3 2 2 2 S = f (x; y; z) 2 R j x + y + z = 1 g: To make a correspondence with the sphere and the plane is simply to make a map (a real map not a function). We define a function 2 3 F : S − f(0; 0; 1)g −! C ⊂ R as follows. Pick a point q = (x; y; z) 2 S2, a point of the unit sphere, other than the north pole p = N = (0; 0; 1). Connect the point p to the point q by a line. This line will meet the plane z = 0, 2 f (x; y; 0) j (x; y) 2 R g in a unique point r. We then identify r with a point F (q) = x + iy 2 C in the usual way. F is called stereographic projection.
    [Show full text]
  • Euclidean Versus Projective Geometry
    Projective Geometry Projective Geometry Euclidean versus Projective Geometry n Euclidean geometry describes shapes “as they are” – Properties of objects that are unchanged by rigid motions » Lengths » Angles » Parallelism n Projective geometry describes objects “as they appear” – Lengths, angles, parallelism become “distorted” when we look at objects – Mathematical model for how images of the 3D world are formed. Projective Geometry Overview n Tools of algebraic geometry n Informal description of projective geometry in a plane n Descriptions of lines and points n Points at infinity and line at infinity n Projective transformations, projectivity matrix n Example of application n Special projectivities: affine transforms, similarities, Euclidean transforms n Cross-ratio invariance for points, lines, planes Projective Geometry Tools of Algebraic Geometry 1 n Plane passing through origin and perpendicular to vector n = (a,b,c) is locus of points x = ( x 1 , x 2 , x 3 ) such that n · x = 0 => a x1 + b x2 + c x3 = 0 n Plane through origin is completely defined by (a,b,c) x3 x = (x1, x2 , x3 ) x2 O x1 n = (a,b,c) Projective Geometry Tools of Algebraic Geometry 2 n A vector parallel to intersection of 2 planes ( a , b , c ) and (a',b',c') is obtained by cross-product (a'',b'',c'') = (a,b,c)´(a',b',c') (a'',b'',c'') O (a,b,c) (a',b',c') Projective Geometry Tools of Algebraic Geometry 3 n Plane passing through two points x and x’ is defined by (a,b,c) = x´ x' x = (x1, x2 , x3 ) x'= (x1 ', x2 ', x3 ') O (a,b,c) Projective Geometry Projective Geometry
    [Show full text]
  • A Gallery of Conics by Five Elements
    Forum Geometricorum Volume 14 (2014) 295–348. FORUM GEOM ISSN 1534-1178 A Gallery of Conics by Five Elements Paris Pamfilos Abstract. This paper is a review on conics defined by five elements, i.e., either lines to which the conic is tangent or points through which the conic passes. The gallery contains all cases combining a number (n) of points and a number (5 − n) of lines and also allowing coincidences of some points with some lines. Points and/or lines at infinity are handled as particular cases. 1. Introduction In the following we review the construction of conics by five elements: points and lines, briefly denoted by (αP βT ), with α + β =5. In these it is required to construct a conic passing through α given points and tangent to β given lines. The six constructions, resulting by giving α, β the values 0 to 5 and considering the data to be in general position, are considered the most important ([18, p. 387]), and are to be found almost on every book about conics. It seems that constructions for which some coincidences are allowed have attracted less attention, though they are related to many interesting theorems of the geometry of conics. Adding to the six main cases those with the projectively different possible coincidences we land to 12 main constructions figuring on the first column of the classifying table below. The six added cases can be considered as limiting cases of the others, in which a point tends to coincide with another or with a line. The twelve main cases are the projectively inequivalent to which every other case can be reduced by means of a projectivity of the plane.
    [Show full text]
  • 1 Affine and Projective Coordinate Notation 2 Transformations
    CS348a: Computer Graphics Handout #9 Geometric Modeling Original Handout #9 Stanford University Tuesday, 3 November 1992 Original Lecture #2: 6 October 1992 Topics: Coordinates and Transformations Scribe: Steve Farris 1 Affine and Projective Coordinate Notation Recall that we have chosen to denote the point with Cartesian coordinates (X,Y) in affine coordinates as (1;X,Y). Also, we denote the line with implicit equation a + bX + cY = 0as the triple of coefficients [a;b,c]. For example, the line Y = mX + r is denoted by [r;m,−1],or equivalently [−r;−m,1]. In three dimensions, points are denoted by (1;X,Y,Z) and planes are denoted by [a;b,c,d]. This transfer of dimension is natural for all dimensions. The geometry of the line — that is, of one dimension — is a little strange, in that hyper- planes are points. For example, the hyperplane with coefficients [3;7], which represents all solutions of the equation 3 + 7X = 0, is the same as the point X = −3/7, which we write in ( −3) coordinates as 1; 7 . 2 Transformations 2.1 Affine Transformations of a line Suppose F(X) := a + bX and G(X) := c + dX, and that we want to compose these functions. One way to write the two possible compositions is: G(F(X)) = c + d(a + bX)=(c + da)+bdX and F(G(X)) = a + b(c + d)X =(a + bc)+bdX. Note that it makes a difference which function we apply first, F or G. When writing the compositions above, we used prefix notation, that is, we wrote the func- tion on the left and its argument on the right.
    [Show full text]
  • Dirichlet Problem at Infinity for P-Harmonic Functions on Negatively
    CORE Metadata, citation and similar papers at core.ac.uk Provided by Helsingin yliopiston digitaalinen arkisto Dirichlet problem at infinity for p-harmonic functions on negatively curved spaces Aleksi Vähäkangas Department of Mathematics and Statistics Faculty of Science University of Helsinki 2008 Academic dissertation To be presented for public examination with the permission of the Faculty of Science of the University of Helsinki in Auditorium XIII of the University Main Building on 13 December 2008 at 10 a.m. ISBN 978-952-92-4874-2 (pbk) Helsinki University Print ISBN 978-952-10-5170-8 (PDF) http://ethesis.helsinki.fi/ Helsinki 2008 Acknowledgements I am in greatly indebted to my advisor Ilkka Holopainen for his continual support and encouragement. His input into all aspects of writing this thesis has been invaluable to me. I would also like to thank the staff and students at Kumpula for many interesting discussions in mathematics. The financial support of the Academy of Finland and the Graduate School in Mathe- matical Analysis and Its Applications is gratefully acknowledged. Finally, I would like to thank my brother Antti for advice and help with this work and my parents Simo and Anna-Kaisa for support. Helsinki, September 2008 Aleksi Vähäkangas List of included articles This thesis consists of an introductory part and the following three articles: [A] A. Vähäkangas: Dirichlet problem at infinity for A-harmonic functions. Potential Anal. 27 (2007), 27–44. [B] I. Holopainen and A. Vähäkangas: Asymptotic Dirichlet problem on negatively curved spaces. Proceedings of the International Conference on Geometric Function Theory, Special Functions and Applications (ICGFT) (R.
    [Show full text]
  • Projective Geometry Yingkun Li 11/20/2011
    UCLA Math Circle: Projective Geometry Yingkun Li 11/20/2011 This handout follows the x9 of [2] and x1 of [1]. For further information, you can check out those two wonderful books. 1 Homogeneous Coordinate On the real number line, a point is represented by a real number x if x is finite. Infinity in this setting could be represented by the symbol 1. However, it is more difficult to calculate with infinity. One way to resolve this issue is through a different coordinate system call homogeneous coordinate. For the real line, this looks a lot like fractions. Instead of using one number x 2 R, we will use two real numbers [x : y], with at least one nonzero, to represent a point on the number line. Furthermore, we define the following equivalent relationship ∼ 0 0 0 0 [x : y] ∼ [x : y ] , there exists 2 Rnf0g s.t. x = x and y = y : Now two points [x : y] and [x0 : y0] represent the same point if and only if they are equivalent. The set of coordinates f[x : y] j x; y 2 R not both zerog with the equivalence relation ∼ represents all the points on the real number line including 1, i.e. [1 : 0]. Furthermore, the operation on the real number line, such as addition, multiplication, inverse, can still be done in this new coordinate system. The line described by this coordinate system is called the projective line, denoted by the 1 symbol RP . Now it is not too difficult to extend this coordinate system to the plane.
    [Show full text]
  • Announcements the Course Equation of Perspective Projection
    Announcements • Assignment 0: “Getting Started with Image Formation, Matlab” is due today Cameras (cont.) • Read Chapters 1 & 2 of Forsyth & Ponce Computer Vision I CSE 252A Lecture 4 CS252A, Fall 2010 Computer Vision I CS252A, Fall 2010 Computer Vision I The course Pinhole Camera: Perspective projection • Part 1: The physics of imaging • Abstract camera model - box • Part 2: Early vision with a small hole in it • Part 3: Reconstruction • Part 4: Recognition CS252A, Fall 2010 Computer Vision I CS252A, Fall 2010 Forsyth&Ponce Computer Vision I Equation of Perspective Projection Geometric properties of projection • 3-D points map to points • 3-D lines map to lines • Planes map to whole image or half-plane • Polygons map to polygons • Important point to note: Angles & distances not preserved, nor are inequalities of angles & distances. Cartesian coordinates: • We have, by similar triangles, that (x, y, z) -> (f’ x/z, f’ y/z, f’) • Degenerate cases: • Establishing an image plane coordinate system at C’ aligned with i – line through focal point project to point and j, we get – plane through focal point projects to a line CS252A, Fall 2010 Computer Vision I CS252A, Fall 2010 Computer Vision I 1 Homogenous coordinates • Our usual coordinate system is called a Euclidean or affine coordinate system A Digression • Rotations, translations and projection in Homogenous coordinates can be expressed linearly as matrix multiplies Projective Geometry Projection Convert Convert and Homogenous Coordinates Euclidean Homogenous Homogenous Euclidean World World Image World 3D 3D 2D 2D CS252A, Fall 2010 Computer Vision I CS252A, Fall 2010 Computer Vision I What is the intersection of Do two lines in the plane always two lines in a plane? intersect at a point? No, Parallel lines don’t A Point meet at a point.
    [Show full text]