Camera: optical system

Similar documents
EPIPOLAR GEOMETRY WITH MANY DETAILS

Lines and points. Lines and points

Camera Projection Model

1. Projective geometry

A geometric interpretation of the homogeneous coordinates is given in the following Figure.

Camera Models and Affine Multiple Views Geometry

Single view metrology

Lecture 5. Epipolar Geometry. Professor Silvio Savarese Computational Vision and Geometry Lab. 21-Jan-15. Lecture 5 - Silvio Savarese

3D from Photographs: Camera Calibration. Dr Francesco Banterle

A Study of Kruppa s Equation for Camera Self-calibration

Multiple View Geometry in Computer Vision

Homogeneous Coordinates

EE Camera & Image Formation

Robotics - Homogeneous coordinates and transformations. Simone Ceriani

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 =

THE ENVELOPE OF LINES MEETING A FIXED LINE AND TANGENT TO TWO SPHERES

Math 302 Outcome Statements Winter 2013

Linear Algebra & Geometry why is linear algebra useful in computer vision?

Pose estimation from point and line correspondences

Trinocular Geometry Revisited

Conics and their duals

Geometric Modeling Summer Semester 2010 Mathematical Tools (1)

Camera calibration. Outline. Pinhole camera. Camera projection models. Nonlinear least square methods A camera calibration tool

Mathematics 309 Conic sections and their applicationsn. Chapter 2. Quadric figures. ai,j x i x j + b i x i + c =0. 1. Coordinate changes

A Practical Method for Decomposition of the Essential Matrix

Linear Algebra & Geometry why is linear algebra useful in computer vision?

8. Diagonalization.

Singular Value Decomposition

Optimisation on Manifolds

Extra Problems for Math 2050 Linear Algebra I

CSE 252B: Computer Vision II

Mobile Robotics 1. A Compact Course on Linear Algebra. Giorgio Grisetti

Exercises for Unit I I (Vector algebra and Euclidean geometry)

MATH PROBLEM SET 6

16. Analysis and Computation of the Intrinsic Camera Parameters

Outline. Linear Algebra for Computer Vision

APPLICATIONS The eigenvalues are λ = 5, 5. An orthonormal basis of eigenvectors consists of

Multiple View Geometry in Computer Vision

Hyperbolic Conformal Geometry with Clifford Algebra 1)

Augmented Reality VU Camera Registration. Prof. Vincent Lepetit

MAT1035 Analytic Geometry

TWO THEOREMS ON THE FOCUS-SHARING ELLIPSES: A THREE-DIMENSIONAL VIEW

6.801/866. Affine Structure from Motion. T. Darrell

Name of the Student:

Decomposition of Screw-motion Envelopes of Quadrics

Identifying second degree equations

QUESTION BANK ON STRAIGHT LINE AND CIRCLE

Final Exam Due on Sunday 05/06

PAijpam.eu EPIPOLAR GEOMETRY WITH A FUNDAMENTAL MATRIX IN CANONICAL FORM Georgi Hristov Georgiev 1, Vencislav Dakov Radulov 2

CLASS XII MATHEMATICS. Weightage (Marks) (i) Relations and Functions 10. Type of Questions Weightage of Number of Total Marks each question questions

Statistical Geometry Processing Winter Semester 2011/2012

Algorithms for Computing a Planar Homography from Conics in Correspondence

Lecture 2: Linear Algebra Review

MECH 576 Geometry in Mechanics September 16, 2009 Using Line Geometry

Kinematics. Chapter Multi-Body Systems

CSE 252B: Computer Vision II

ON THE INVERSE KINEMATICS OF A FRAGMENT OF PROTEIN BACKBONE

POINT. Preface. The concept of Point is very important for the study of coordinate

Basic Math for

CHAPTER 3. Gauss map. In this chapter we will study the Gauss map of surfaces in R 3.

Distance Formula in 3-D Given any two points P 1 (x 1, y 1, z 1 ) and P 2 (x 2, y 2, z 2 ) the distance between them is ( ) ( ) ( )

(D) (A) Q.3 To which of the following circles, the line y x + 3 = 0 is normal at the point ? 2 (A) 2

Two-View Segmentation of Dynamic Scenes from the Multibody Fundamental Matrix

1 Vectors and the Scalar Product

Algebra Workshops 10 and 11

The structure tensor in projective spaces

Image Transform. Panorama Image (Keller+Lind Hall)

Vectors Coordinate frames 2D implicit curves 2D parametric curves. Graphics 2008/2009, period 1. Lecture 2: vectors, curves, and surfaces

Multi-Frame Factorization Techniques

Linear algebra for computational statistics

arxiv: v6 [math.mg] 9 May 2014

DIRECTORATE OF EDUCATION GOVT. OF NCT OF DELHI

Methods in Computer Vision: Introduction to Matrix Lie Groups

11/29/2010. Propagation in Anisotropic Media 3. Introduction. Introduction. Gabriel Popescu

Part IB GEOMETRY (Lent 2016): Example Sheet 1

SIAM Conference on Applied Algebraic Geometry Daejeon, South Korea, Irina Kogan North Carolina State University. Supported in part by the

Chapter 3. Riemannian Manifolds - I. The subject of this thesis is to extend the combinatorial curve reconstruction approach to curves

Reading. 4. Affine transformations. Required: Watt, Section 1.1. Further reading:

MATH 1020 WORKSHEET 12.1 & 12.2 Vectors in the Plane

1 Differentiable manifolds and smooth maps. (Solutions)

Induced Planar Homologies in Epipolar Geometry

3. Interpret the graph of x = 1 in the contexts of (a) a number line (b) 2-space (c) 3-space

Linear Algebra Primer

Dynamic P n to P n Alignment

Matrix Algebra for Engineers Jeffrey R. Chasnov

Study guide for Exam 1. by William H. Meeks III October 26, 2012

Projection pencils of quadrics and Ivory s theorem

ROBOTICS 01PEEQW. Basilio Bona DAUIN Politecnico di Torino

Chapter 12. The cross ratio Klein s Erlanger program The projective line. Math 4520, Fall 2017

+ 2gx + 2fy + c = 0 if S

Introduction to pinhole cameras

Numerical Methods for Inverse Kinematics

Analytic Geometry MAT 1035

Chapter 8. Rigid transformations

Mathematical Fundamentals

APPENDIX A. Background Mathematics. A.1 Linear Algebra. Vector algebra. Let x denote the n-dimensional column vector with components x 1 x 2.

Q.2 A, B and C are points in the xy plane such that A(1, 2) ; B (5, 6) and AC = 3BC. Then. (C) 1 1 or

Definition: A vector is a directed line segment which represents a displacement from one point P to another point Q.

AOL Spring Wavefront Sensing. Figure 1: Principle of operation of the Shack-Hartmann wavefront sensor

Mathematics 1. Part II: Linear Algebra. Exercises and problems

Transcription:

Camera: optical system Y dαα α ρ ρ curvature radius Z thin lens small angles: α Y ρ α Y ρ

lens refraction inde: n Y incident light beam deviated beam θ '' α θ ' α sin( θ '' α ) sin( θ ' α ) n θ α θ α deviation angle? Δθ θ θ Δθ ( n ) Y ( + ) ρ ρ sin( θ α) sin( θ ' α ) ' n Z

a) Y Δθ b) f Δθ Y hin lens rules f ( n )( ρ + ) ρ parallel rays converge onto a focal plane beams through lens center: undeviated independent of y Y Δθ f

Hp: Z >> a r f the image of a point P belongs to the line (P,O) P image plane O p p image of P image plane line(o,p) interpretation line of p: line(o,p) locus of the scene points projecting onto image point p

Projective D geometry Notes based on di R.Hartley e A.Zisserman Multiple view geometry

Projective D Geometry Points, lines & conics ransformations & invariants D projective geometry and the Cross-ratio

Homogeneous coordinates Homogeneous representation of lines a + by + c ( a,b,c) ( ka ) + ( kb) y + kc, k ( a,b,c ) ~ k( a,b,c) equivalence class of vectors, any vector is representative Set of all equivalence classes in R 3 (,,) forms P Homogeneous representation of points, y on l a,b,c if and only if a + by + c ( ) ( ) (,y, )( a,b,c) (,y, ) l (, y,) ~ k(, y,), k he point lies on the line l if and only if ll ( ) Homogeneous coordinates,, Inhomogeneous coordinates (, y) 3 but only DOF

Points from lines and vice-versa Intersections of lines he intersection of two lines l and l' is l l' Line joining two points he line through two points and ' is Eample l Line joining two points: parametric equation A point on the line through two points and is y + θ ' y

Ideal points and the line at infinity Intersections of parallel lines l ( a, b, c) and l' ( a, b, c' ) l l' ( b, a,) Eample ( b, a) ( a,b) tangent vector normal direction Ideal points (,, ) Line at infinity l (,,) P R l Note that in P there is no distinction between ideal points and others

Duality l l l l l' l ' Duality principle: o any theorem of -dimensional projective geometry there corresponds a dual theorem, which may be derived by interchanging the role of points and lines in the original theorem

Conics Curve described by nd -degree equation in the plane a + by + cy + d + ey + f or homogenized a α, 3 y α + b + c + d3 + e3 + f3 or in matri form a b / d / C with C b / c e / d / e / f 3 5DOF: { a : b : c : d : e : f }

Five points define a conic For each point the conic passes through + + + + + f ey d cy y b a i i i i i i or ( ),,,,, c f y y y i i i i i i ( ) f e d c b a,,,,, c 5 5 5 5 5 5 4 4 4 4 4 4 3 3 3 3 3 3 c y y y y y y y y y y y y y y y stacking constraints yields

Pole-polar relationship he polar line lc of the point with respect to the conic C intersects the conic in two points. he two lines tangent to C at these points intersect at

Polarity: cross ratio Cross ratio of 4 colinear points y + θ (with i,..,4) ratio of ratios θ θ -θ -θ 3 4 θ θ i -θ -θ 3 4 i Harmonic 4-tuple of colinear points: such that CR-

angent lines to conics he line l tangent to C at point on C is given by lc l C

Dual conics A line tangent to the conic C satisfies l C * l C * In general (C full rank): in fact Dual conics line conics conic envelopes C Line l is the polar line of y : y C l, but since y Cy - - * - - l C C C l C C C -

Degenerate conics A conic is degenerate if matri C is not of full rank e.g. two lines (rank ) C lm + ml e.g. repeated line (rank ) m l C ll l Degenerate line conics: points (rank ), double point (rank) Note that for degenerate conics ( * C ) * C

Projective transformations Definition: A projectivity is an invertible mapping h from P to itself such that three points,, 3 lie on the same line if and only if h( ),h( ),h( 3 ) do. heorem: A mapping h:p P is a projectivity if and only if there eist a non-singular 33 matri H such that for any point in P represented by a vector it is true that h()h Definition: Projective transformation ' h ' h ' 3 h 3 h h h 3 h h h 3 3 33 3 or ' H 8DOF projectivitycollineationprojective transformationhomography

Mapping between planes central projection may be epressed by H (application of theorem)

Removing projective distortion select four points in a plane with known coordinates ' ' ' 3 h h 3 + h + h 3 y + h y + h 3 33 y' ' ' 3 h h 3 + h + h 3 y + h y + h 3 33 ( h3 + h3 y + h33 ) h + h y h3 ( h3 + h3 y + h33 ) h + h y h3 ' + y ' + (linear in h ij ) ( constraints/point, 8DOF 4 points needed) Remark: no calibration at all necessary, better ways to compute (see later)

More eamples

ransformation of lines and conics For a point transformation ' H ransformation for lines l' H - ransformation for conics l C ' H - CH - ransformation for dual conics * C' * HC H

A hierarchy of transformations Projective linear group Affine group (last row (,,)) Euclidean group (upper left orthogonal) Oriented Euclidean group (upper left det ) Alternative, characterize transformation in terms of elements or quantities that are preserved or invariant e.g. Euclidean transformations leave distances unchanged

Class I: Isometries (isosame, metricmeasure) ' y' ε cosθ ε sinθ sinθ cosθ t t y y ε ± orientation preserving: orientation reversing: ε ε ' H E R t 3DOF ( rotation, translation) R R special cases: pure rotation, pure translation I Invariants: length, angle, area

Class II: Similarities (isometry + scale) cos sin sin cos ' ' y t s s t s s y y θ θ θ θ ' t R H s S I R R also know as equi-form (shape preserving) metric structure structure up to similarity (in literature) 4DOF ( scale, rotation, translation) Invariants: ratios of length, angle, ratios of areas, parallel lines

Class III: Affine transformations ' ' y t a a t a a y y ' t A H A non-isotropic scaling! (DOF: scale ratio and orientation) 6DOF ( scale, rotation, translation) Invariants: parallel lines, ratios of parallel segment lengths, ratios of areas ( ) ( ) ( ) φ φ θ DR R R A λ λ D ) )( ( VDV UV UDV A where

Action of affinities and projectivities on line at infinity + v v v v A A t v A A t Line at infinity becomes finite, allows to observe vanishing points, horizon, Line at infinity stays at infinity, but points move along line

Class VI: Projective transformations ' H P A v t v v ( v v ), 8DOF ( scale, rotation, translation, line at infinity) Action: non-homogeneous over the plane Invariants: cross-ratio of four points on a line (ratio of ratios)

Projective geometry of D ( ), ' H 3DOF (-) he cross ratio Cross (,,, ) 3 4,, 3 3,, 4 4, i j i det i j j Invariant under projective transformations

Overview transformations Projective 8dof Affine 6dof Similarity 4dof Euclidean 3dof h h h3 a a sr sr r r h h h a a r r 3 sr sr h3 h 3 h 33 t t y t t y t t y Concurrency, collinearity, order of contact (intersection, tangency, inflection, etc.), cross ratio Parallellism, ratio of areas, ratio of lengths on parallel lines (e.g midpoints), linear combinations of vectors (centroids). he line at infinity l Ratios of lengths, angles. he circular points I,J lengths, areas.

Number of invariants? he number of functional invariants is equal to, or greater than, the number of degrees of freedom of the configuration less the number of degrees of freedom of the transformation e.g. configuration of 4 points in general position has 8 dof (/pt) and so 4 similarity, affinity and zero projective invariants

Recovering metric and affine properties from images Parallelism Parallel length ratios Angles Length ratios

he line at infinity l l l A t A H A he line at infinity l is a fied line under a projective transformation H if and only if H is an affinity Note: not fied pointwise

Affine properties from images projection (affine) rectification A P l l l H H 3 ' [ ], l' 3 3 l l l l in fact, any point on l is mapped to a point at the

Affine rectification v v l l l 3 l l 4 v l3 l4 v l l l v v

he circular points I i J i I cos sin sin cos I I i se i t s s t s s i y S θ θ θ θ θ H he circular points I, J are fied points under the projective transformation H iff H is a similarity

he circular points circular points l + + d3 + e3 + f3 + 3 I (, i,) Intersection points between any circle and l J, i, Algebraically, encodes orthogonal directions I (,, ) + i(,, ) ( )

Circular points invariance {I,J} l any circumference Similarity: circ circ Similarity: circ l circ l Similarity: {I,J} {I,J} circular points: invariant under similarity

* C Conic dual to the circular points * C IJ + JI : line conic set of lines through any of the circular points C * C * H C * S H S C he dual conic * is fied conic under the projective transformation H iff H is a similarity * C Note: has 4DOF l is the null vector

Angles ( )( ) cos m m l l m l m l + + + θ ( ) 3,, l l l l ( ) 3,, m m m m Euclidean: Projective: ( )( ) m m l l m l cos * * * C C C θ m l * C (orthogonal)

Metric properties from images ( ) ( ) ( ) ( ) ( ) ( ) v v v v ' * * * * K K KK H H C H H H H H C H H H H H H C H H H C A P A P A P S S A P S A P S A P U U C ' * Rectifying transformation from SVD H U

Normally: SVD (Singular Value Decomposition) U V C c b a ' * But is symmetric ' * C ' ' * * C VDU UDV C and SVD is unique V U with U and V orthogonal U U C ' * Why? Observation : HU orthogonal (33): not a P isometry

Once the image has been affinely rectified A H A C H C * * ' ' * t t KK K K C Metric from affine

Metric from affine ( ) 3 3 m m m l l l KK ( )( ),,,, + + k k k k k m l m l m l m l

Metric from projective ( ) ( ) ( ) ( ) c,,,,, 3 3 3 3 3 3 5..5.5 + + + m l m l m l m l m l m l m l m l m l ( ) v v v v 3 3 m m m l l l K K KK

Fied points and lines H e λ e (eigenvectors H fied points) (λ λ pointwise fied line) H l λ l (eigenvectors H - fied lines)

Projective 3D geometry

Singular Value Decomposition n n n m m m n m V Σ U A I U U n σ σ σ Λ I V V m n X V X V Σ X V UΣ n n U n V V U V U A σ σ σ + + + Λ Σ Λ Μ Μ Μ Λ Μ Ο Μ Μ Λ Λ σ n σ σ

Singular Value Decomposition Homogeneous least-squares min AX subject to X solution X Vn Span and null-space S L [ U U ]; N [ U ] L 3U4 V V ; R V3V 4 Closest rank r approimation A ~ UΣ ~ V [ ] [ ] S R N Pseudo inverse + A VΣ U A UΣ V ~ Σ diag Σ σ ( σ, σ, Λ, σ,,, ) r Λ + ( Σ + diag σ, σ, Λ, σ,,, ) σ r Λ

Projective 3D Geometry Points, lines, planes and quadrics ransformations П, ω and Ω

3D points 3D point ( X, Y, Z ) in R 3 X ( X, X, X X ) 3, 4 in P 3 X X X X 3,,, Z X 4 X 4 X 4 ( X, Y,,) ( X 4 ) projective transformation X' H X (44-5 dof)

3D plane Planes π X + πy + π3z + π4 π X + πx + π3x 3 + π4x 4 ransformation X' H X - π' H π π X Euclidean n.x ~ representation + d ( π ), π, π3 d n ( ) π 4 X ~ X, Y, Z X 4 d / n Dual: points planes, lines lines

Solve X X X Planes from points 3 π from 3 π [ X X X X ] det 3 X π, X π and X π X (solve π as right nullspace of X ) X3 Or implicitly from coplanarity condition X det X X X 3 4 ( X) ( X ) ( X ) 3 ( X) ( X ) ( X ) 3 ( X) ( X ) ( X 3) 3 3 ( X ) ( X ) ( X ) X D 34 X D 34 + X 3 D 4 X π ( D ) 4D3 34, D34, D4, D3 4 4 3 3 4

Points from planes X π π π 3 X M [ ] 3 X X X M π M M I p ( ) d c b a,,, π a d a c a b p,, X π X π, X π X 3 and from Solve (solve as right nullspace of ) X 3 π π π Representing a plane by its span

Lines W A B λa+ μb (4dof) W * P Q λp+ μq two points A and B two planes P and Q * * W W WW Eample: X-ais W W *

Points, lines and planes W M M π X X W W * M M X π * W π

Quadrics and dual quadrics X QX. 9 d.o.f. (Q : 44 symmetric matri). in general 9 points define quadric 3. det Q degenerate quadric 4. pole polar 5. (plane quadric)conic 6. transformation π * Q π π QX Q' H Q ο ο ο ο ο ο C M QM π : X M - QH - * - Q Q * * HQ H. relation to quadric (non-degenerate). transformation Q'

Quadric classification Rank Sign. Diagonal Equation Realization 4 4 (,,,) X + Y + Z + No real points (,,,-) X + Y + Z Sphere (,,-,-) X + Y Z + Hyperboloid (S) 3 3 (,,,) X + Y + Z Single point (,,-,) X + Y Z Cone (,,,) X + Y Single line (,-,,) X Y wo planes (,,,) X Single plane

Quadric classification Projectively equivalent to sphere: sphere Ruled quadrics: ellipsoid hyperboloid of two sheets paraboloid hyperboloids of one sheet Degenerate ruled quadrics: cone two planes

Hierarchy of transformations Projective 5dof A v t v Intersection and tangency Affine dof A t Parallellism of planes, Volume ratios, centroids, he plane at infinity π Similarity 7dof s R t he absolute conic Ω Euclidean 6dof R t Volume

Screw decomposition Any particular translation and rotation is equivalent to a rotation about a screw ais and a translation along the screw ais. screw ais // rotation ais t + t t //

he plane at infinity A π H A π π A t he plane at infinity π is a fied plane under a projective transformation H iff H is an affinity π (,,,). canical position. contains directions D ( X ), X, X 3, 3. two planes are parallel line of intersection in π 4. line // line (or plane) point of intersection in π

he absolute conic he absolute conic Ω is a (point) conic on π. In a metric frame: X or conic for directions: (with no real points) + X X 4 + X he absolute conic Ω is a fied conic under the projective transformation H iff H is a similarity 3 ( X, X, X ) I( X, X X ) 3, 3. Ω is only fied as a set. Circle intersect Ω in two points 3. Spheres intersect π in Ω

Absolute conic invariance Ω π any sphere Similarity: sphere sphere Similarity: sphere π sphere π Similarity: Ω Ω Ω : invariant under similarity

he absolute conic Euclidean: Projective: cosθ cosθ ( d ) d ( )( d ) d d d ( d ) Ω d ( d d )( Ω d Ω d ) d Ω d (orthogonalityconjugacy) normal plane

he absolute dual quadric * I Ω he absolute conic Ω * is a fied conic under the projective transformation H iff H is a similarity. 8 dof. plane at infinity π is the nullvector of Ω 3. Angles: * cosθ ( * )( * π Ω π π Ω π ) π Ω π

Y P X p c y f O Z y f f X Z Y Z perspective projection -nonlinear -not shape-preserving -not length-ratio preserving

Homogeneous coordinates In D: add a third coordinate, w Point [,y] epanded to [u,v,w] Any two sets of points [u,v,w ] and [u,v,w ] represent the same point if one is multiple of the other [u,v,w] [,y] with u/w, and yv/w [u,v,] is the point at the infinite along direction (u,v)

ransformations translation by vector [d,d y ] scaling (by different factors in and y) rotation by angle θ

Homogeneous coordinates In 3D: add a fourth coordinate, t Point [X,Y,Z] epanded to [,y,z,t] Any two sets of points [,y,z,t ] and [,y,z,t ] represent the same point if one is multiple of the other [,y,z,t] [X,Y,Z] with X/t, Yy/t, and Zz/t [,y,z,] is the point at the infinite along direction (,y,z)

ransformations translation scaling rotation Obs: rotation matri is an orthogonal matri i.e.: R - R

Pinhole camera model Z fy Z fx Z Y X ) /, / ( ),, ( α Z Y X f f Z fy fx Z Y X α

Scene->Image mapping: perspective transformation u v w f with f X Y Z y u w v w With ad hoc reference frames, for both image and scene

Let us recall them O image reference - centered on principal point - - and y-aes parallel to the sensor rows and columns - Euclidean reference c y f Y X scene reference - centered on lens center - Z-ais orthogonal to image plane - X- and Y-aes opposite to image - and y-aes Z

principal point c Actual references are generic y f O principal ais Y X image reference - centered on upper left corner - nonsquare piels (aspect ratio) noneuclidean reference Z scene reference - not attached to the camera

Principal point offset o o v Z fy u Z fx Z Y X ) /, / ( ),, ( + + α principal point o v o u ), ( + + Z Y X v f u f Z Zv fy Zu fx Z Y X o o o o α

CCD camera o y o v f u f K o o v f u f a K

Scene-image relationship wrt actual reference frames image u A u v w ( s a) a u v o o u v w normally, s X scene Y Z R t X

X t R X P u f f v a u o o t K R K t R I P o o v af u f K upper triangular: intrinsic camera parameters scene-camera tranformation etrinsic camera parameters orthogonal (3D rotation) matri P: - degrees of freedom ( if s)

z y o I M X o I M u ) ( o M o I M u i.e., defining [, y, z] o I M o M M m M t K R K P with R K M and m M o

Interpretation of o: u is image of if ( o) λm u i.e., if o + λm u he locus of the points whose image is u is a straight line through o having direction d M u o M m is independent of u o is the camera viewpoint (perspective projection center) line(o, d) Interpretation line of image point u

M K and R Intrinsic and etrinsic parameters from P M K R M R K R K RQ-decomposition of a matri: as the product between an orthogonal matri and an upper triangular matri M and m t o M m ( KR) Kt R K Kt R t t Ro

Camera center null-space camera projection matri P O X λa + ( λ)o PX λpa + ( λ)po For all A all points on AO project on image of A, therefore O is camera center Image of camera center is (,,), i.e. undefined o M m Finite cameras: O d Infinite cameras: o,md

Action of projective camera on point Forward projection PX [ M m] D Md PD Back-projection P O + X P P + P ( PP ) ( λ) P + λo X + - d M M X( λ) μ D - (pseudo-inverse) - M + O - m M - + PP I ( μ - m)

Camera matri decomposition Finding the camera center P O (use SVD to find null-space) X det( [ p,p3,m] ) Y det( [ p,p3,m] ) Z det( [ p,p,m] ) ([ p,p,m]) det Finding the camera orientation and internal parameters M KR (use RQ decomposition ~QR) (if only QR, invert) ( Q R ) - R - - Q

Radial distortion

Data normalization 3 X ~ i UX i u~ i u i (i) translate origin to gravity center (ii) (an)isotropic scaling

Eterior orientation Calibrated camera, position and orientation unkown Pose estimation 6 dof 3 points minimal (4 solutions in general)

Properties of perspective transformations ) vanishing points V image of the point at the along direction d u V M m d M d d M u V the interpretation line of V is parallel to d

d P O V he images of parallel lines are concurrent lines

Properties of perspective transformations ctd. ) cross ratio invariance Given four colinear points ( p, p, p p ) let (,, ) 3, 4 3, be their abscissae 4 CR ( p, p, p p ) 3, 4 3 4 3 4

Cross ratio invariance under perspective transformation a point on the line yz X [, y, z, t] [,,, t ] its image its coordinate u det u, u CR( u, u, u3, u4 ) det u, u det u [ u, v, w] 3 3 det u det u P X p u [ u, w] P4 p t 3, u, u 4 4 4 detp detp 4 4 4 det det,, 3 4 CR(,, 3, 4 ) det,, det det,, 3 belongs to a line detp detp 4 4 det det 3,, 4 4

Object localization : three colinear points geometric model of an object a perspective image of the object position and orientation of the object? calibrated camera: known P M known interpretation lines m C A B B A O C B A π C

a) orientation C A B C V O A B π Cross ratio invariance: a c solve CR( A', B', C', V ) CR( A, B, C, ) for V (image of ) b c V: vanishing point of the direction of (A,B,C) interpretation line of V parallel to (A,B,C) direction M uv

b) position (e.g., distance(o,a)) M u C γ C A B C V O α M u A A B M u V π interpretation lines OA AC sinγ sinα angles α and γ

Object localization : four coplanar points D C O E B A (i) orientation of (A,E,C) (ii) orientation of (B,E,D) (iii) distance (O,A)

Off-side Find vanishing point of the field-bottom line direction b b a a images of symmetric segments

a c b d a and b: c(a+b)/: d : abscissae of the endpoints of a segment abscissa of segment midpoint, point at the infinite along the segment direction CR a c a d b d a c b c ( a, b, c, d ) b c Harmonic 4-tuple (a,b,c,d) (a,b ) and (a, b ) are image of symmetric segments same image of the midpoint c, same vanishing point d

solve a' c' ( a', b', c', d' ) b' c' { CR for CR a' ', b'', c', d' a'' c' a' d' b' d' ( ) b'' c' a'' d' b'' d' c, d system of two linear equations in (c d ) and (c +d ) two degree equation, whose solutions are c and d among the two solutions, the one for d is the value eternal to the range [a,b ]

What can be told from a single image?

Action of projective camera on planes [ ] [ ] p p p p p p p PX 4 4 3 Y X Y X he most general transformation that can occur between a scene plane and an image plane under perspective imaging is a plane projective transformation

Action of projective camera on lines forward projection ( μ) P(A + μb) PA + μpb a μb X + back-projection l PX l PX with Π P l Π X Interpretation plane of line l

[ ] [ ] [ ] P Y X Y X 4 4 4 3 p p p p p p p p p p PX Image of a conic C C - P P therefore - C C' P P

Action of projective camera on conics back-projection of a conic C to cone Q co C Q co

back-projection of a conic C to cone Qco C PX C X P CPX with X Q co X Interpretation cone of a conic C Q co P CP eample: Q co K C [ K ] K CK

Images of smooth surfaces he contour generator Γ is the set of points X on S at which rays are tangent to the surface. he corresponding apparent contour γ is the set of points which are the image of X, i.e. γ is the image of Γ he contour generator Γ depends only on position of projection center, γ depends also on rest of P

Action of projective camera on quadrics apparent contour of a quadric Q dual quadric Q * Q is a plane quadric: the set of planes tangent to Q Π Q * Π Let us consider only those planes that are backprojection of image lines Π P l Π Q * Π l * PQ P l with C * * PQ P its dual is C C *

he plane containing the apparent contour Γ of a quadric Q from a camera center O follows from pole-polar relationship ΠQO he cone with verte V and tangent to the quadric Q is Q CO (V QV)Q - (QV)(QV) Q CO V back-projection to cone

What does calibration give? K[I ] d d K d d cos θ ( )( ) ( )( ) - - - - d d d d (K K ) (K K ) (K - K - ) An image line l defines a plane through the camera center with normal nk l measured in the camera s Euclidean frame. In fact the backprojection of l is P l nk l

he image of the absolute conic Ω d PX KR[I O] KRd mapping between π to an image is given by the planar homogaphy Hd, with HKR absolute conic (IAC), represented by I 3 within π (w) its image (IAC) ( ) - - KK K K ω ( C α H CH ) (i) IAC depends only on intrinsics (ii) angle between two rays cos θ (iii) DIACω * KK (iv) ω K (Cholesky factorization) (v) image of circular points belong to ω (image of absolute conic) ( ω )( ω ) ω

A simple calibration device (i) compute H i for each square (corners (,),(,),(,),(,)) (ii) compute the imaged circular points H i [,±i,] (iii) fit a conic ω to 6 imaged circular points (iv) compute K from ωk - K - through Cholesky factorization ( Zhang s calibration method)

Orthogonality relation cos θ ( )( ) v ωv v ωv v ωv v ωv * l ω l

Calibration from vanishing points and lines

Calibration from vanishing points and lines

wo-view geometry Epipolar geometry F-matri comp. 3D reconstruction Structure comp.

hree questions: (i) Correspondence geometry: Given an image point in the first view, how does this constrain the position of the corresponding point in the second image? (ii) Camera geometry (motion): Given a set of corresponding image points { i i }, i,,n, what are the cameras P and P for the two views? (iii) Scene geometry (structure): Given corresponding image points i i and cameras P, P, what is the position of (their pre-image) X in space?

he epipolar geometry C,C,, and X are coplanar

he epipolar geometry What if only C,C, are known?

he epipolar geometry All points on π project on l and l

he epipolar geometry Family of planes π and lines l and l Intersection in e and e

he epipolar geometry epipoles e,e intersection of baseline with image plane projection of projection center in other image vanishing point of camera motion direction an epipolar plane plane containing baseline (-D family) an epipolar line intersection of epipolar plane with image (always come in corresponding pairs)

Eample: converging cameras

Eample: motion parallel with image plane

Eample: forward motion e e

he fundamental matri F algebraic representation of epipolar geometry α l' we will see that mapping is (singular) correlation (i.e. projective mapping from points to lines) represented by the fundamental matri F

he fundamental matri F geometric derivation ' H π [ e' ] Hπ F l' e' ' mapping from -D to -D family (rank )

he fundamental matri F algebraic derivation ( λ) P λc X + + ( + P P I) l P'C P' P + P + X( λ) F [ ] + e' P' P (note: doesn t work for CC F)

he fundamental matri F correspondence condition he fundamental matri satisfies the condition that for any pair of corresponding points in the two images ' F ( ' l' )

he fundamental matri F F is the unique 33 rank matri that satisfies F for all (i) ranspose: if F is fundamental matri for (P,P ), then F is fundamental matri for (P,P) (ii) Epipolar lines: l F & lf (iii) Epipoles: on all epipolar lines, thus e F, e F, similarly Fe (iv) F has 7 d.o.f., i.e. 33-(homogeneous)-(rank) (v) F is a correlation, projective mapping from a point to a line l F (not a proper correlation, i.e. not invertible)

he epipolar line geometry l,l epipolar lines, k line not through e l F[k] l and symmetrically lf [k ] l l l k k e e' l Fk (pick ke, since e e ) F[ e] l l F [ e' ] l' l'

Fundamental matri for pure translation

Fundamental matri for pure translation

Fundamental matri for pure translation [ e' ] H [ ] F e' ( H K RK) eample: e' (,,) F - ' F y y' PX ' P' X K[I ]X K - K[I t] Z ( X,Y,Z) - K /Z ' + Kt/Z motion starts at and moves towards e, faster depending on Z pure translation: F only d.o.f., [e] auto-epipolar

General motion [ e' ] H [ e' ] ˆ ' ' ' K'RK - + K' t/z

Geometric representation of F F ( F F )/ ( F F F )/ S + A ( F + ) F S F A F ( F ) A F F s : Steiner conic, 5 d.o.f. F a [ a ] : pole of line ee w.r.t. F s, d.o.f. S

Geometric representation of F

Pure planar motion Steiner conic F s is degenerate (two lines)

Projective transformation and invariance Derivation based purely on projective concepts ˆ H, ˆ' H' ' Fˆ PX ( )( - PH H X) Pˆ Xˆ ( )( - P' H H X) Pˆ' Xˆ H' - FH - F invariant to transformations of projective 3-space ' P' X ( P, P' ) α F F α ( P,P' ) canonical form P' [I ] [M m] unique not unique P [ m] M F

Projective ambiguity of cameras given F previous slide: at least projective ambiguity this slide: not more! ~ ~ Show that if F is same for (P,P ) and (P,P ), there eists a projective transformation H so that ~ ~ PHP and P HP lemma: H P' H ~ P [I ] P' [A a] P [I ] P ~ ' F [ a] A [ ~ a] A ~ ~ A ~ a ka k A + av ( ) rank [A ~ ~ a] af a[] a A ~ af ~ a ka a A ~ a A ~ a ka ~ - A ka ~ - A av [ ] [ ] [ ] ( ) ( ) k I k v k [A a] k I k v k [ k ( A - av ) ka] P ~ ' (-57, ok)

Canonical cameras given F F matri corresponds to P,P iff P FP is skew-symmetric ( X P' FPX, X) F matri, S skew-symmetric matri P [I ] P' [SF e'] [SF e'] (fund.matrif) F[I ] F S F e' F F S F Possible choice: P [I ] P' [[e'] F e'] Canonical representation: P [I ] P' [[e'] F + e' v λe']

he essential matri ~fundamental matri for calibrated cameras (remove K) E ˆ ' [ t] R R[R t] Eˆ - - ( ˆ K ; ˆ' K ' ) E K' FK 5 d.o.f. (3 for R; for t up to scale) E is essential matri if and only if two singularvalues are equal (and third) E Udiag(,,)V

Four possible reconstructions from E (only one solution where points is in front of both cameras)

wo-view geometry Epipolar geometry F-matri comp. 3D reconstruction Structure comp.

hree questions: (i) Correspondence geometry: Given an image point in the first view, how does this constrain the position of the corresponding point in the second image? (ii) Camera geometry (motion): Given a set of corresponding image points { i i }, i,,n, what are the cameras P and P for the two views? (iii) Scene geometry (structure): Given corresponding image points i i and cameras P, P, what is the position of (their pre-image) X in space?

he epipolar geometry C,C,, and X are coplanar

he epipolar geometry What if only C,C, are known?

he epipolar geometry All points on π project on l and l

he epipolar geometry Family of planes π and lines l and l Intersection in e and e

he epipolar geometry epipoles e,e intersection of baseline with image plane projection of projection center in other image vanishing point of camera motion direction an epipolar plane plane containing baseline (-D family) an epipolar line intersection of epipolar plane with image (always come in corresponding pairs)

Eample: converging cameras

Eample: motion parallel with image plane

Eample: forward motion e e

he fundamental matri F algebraic representation of epipolar geometry α l' we will see that mapping is (singular) correlation (i.e. projective mapping from points to lines) represented by the fundamental matri F

he fundamental matri F geometric derivation ' H π [ e' ] Hπ F l' e' ' mapping from -D to -D family (rank )

l' he fundamental matri F algebraic derivation M + ( ) λc X λ - M [ ] [ ] [ ] M' m' C M' m' M' m' C M' M ' l' ' l' ' F P + X( λ) F [ ] e' M' M (note: doesn t work for CC F)

he fundamental matri F correspondence condition he fundamental matri satisfies the condition that for any pair of corresponding points in the two images ' F ( ' l' )

he fundamental matri F F is the unique 33 rank matri that satisfies F for all (i) ranspose: if F is fundamental matri for (P,P ), then F is fundamental matri for (P,P) (ii) Epipolar lines: l F & lf (iii) Epipoles: on all epipolar lines, thus e F, e F, similarly Fe (iv) F has 7 d.o.f., i.e. 33-(homogeneous)-(rank) (v) F is a correlation, projective mapping from a point to a line l F (not a proper correlation, i.e. not invertible)

he epipolar line geometry l,l epipolar lines, k line not through e l F[k] l and symmetrically lf [k ] l l l k k e e' l Fk (pick ke, since e e ) F[ e] l l F [ e' ] l' l'

Fundamental matri for pure translation

Fundamental matri for pure translation

Fundamental matri for pure translation [ e' ] H [ ] F e' ( H K RK) eample: e' (,,) F - ' F y y' PX ' P' X K[I ]X - K[I t] K Z ( X,Y,Z) - K /Z ' + Kt/Z motion starts at and moves towards e, faster depending on Z pure translation: F only d.o.f., [e] auto-epipolar

General motion [ e' ] H [ e' ] ˆ ' ' ' K'RK - + K' t/z

Geometric representation of F F ( F F )/ ( F F F )/ S + A ( F + ) F S F A F ( F ) A F F s : Steiner conic, 5 d.o.f. F a [ a ] : pole of line ee w.r.t. F s, d.o.f. S

Geometric representation of F

Pure planar motion Steiner conic F s is degenerate (two lines)

Projective transformation and invariance Derivation based purely on projective concepts ˆ H, ˆ' H' ' Fˆ PX ( )( - PH H X) Pˆ Xˆ ( )( - P' H H X) Pˆ' Xˆ H' - FH - F invariant to transformations of projective 3-space ' P' X ( P, P' ) α F F α ( P,P' ) canonical form P' [I ] [M m] unique not unique P [ m] M F

Projective ambiguity of cameras given F previous slide: at least projective ambiguity this slide: not more! ~ ~ Show that if F is same for (P,P ) and (P,P ), there eists a projective transformation H so that ~ ~ PHP and P HP lemma: H P' H ~ P [I ] P' [A a] P ~ [I ] P' F [ a] A [ ~ a] A ~ ~ A ~ a ka k A + av ( ) rank [A ~ ~ a] af a[] a A ~ af ~ a ka a A ~ a A ~ a ka ~ - A ka ~ - A av [ ] [ ] [ ] ( ) ( ) k I k v k [A a] k I k v k [ k ( A - av ) ka] P ~ ' (-57, ok)

Canonical cameras given F F matri corresponds to P,P iff P FP is skew-symmetric ( X P' FPX, X) F matri, S skew-symmetric matri P [I ] P' [SF e'] [SF e'] (fund.matrif) F[I ] F S F e' F F S F Possible choice: P [I ] P' [[e'] F e'] Canonical representation: P [I ] P' [[e'] F + e' v λe']

he essential matri ~fundamental matri for calibrated cameras (remove K) E ˆ ' [ t] R R[R t] Eˆ - - ( ˆ K ; ˆ' K ' ) E K' FK 5 d.o.f. (3 for R; for t up to scale) E is essential matri if and only if two singularvalues are equal (and third) SVD E Udiag(,,)V

Motion from E Given W Four solutions R' UW V and t' U 3 R' ' UWV t' ' U 3

Four possible reconstructions from E (only one solution where points is in front of both cameras)