Lecture 3. Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab. 13- Jan- 15.

Similar documents
Geometric Camera Calibration

Lecture 3 Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab

Lecture 8: Camera Calibration

Lecture 8: Camera Calibra0on

CS4495/6495 Introduction to Computer Vision. 3C-L3 Calibrating cameras

AUTO-CALIBRATION. FACTORIZATION. STRUCTURE FROM MOTION.

[ ] T Journal of Software. Vol.13, No /2002/13(05)

THE PRACTICE OF AUTOMATIC SATELLITE IMAGE REGISTRATION

Structure from Motion. Forsyth&Ponce: Chap. 12 and 13 Szeliski: Chap. 7

Trifocal Tensor for Heterogeneous Cameras

Solutions for Homework #9

y new = M x old Feature Selection: Linear Transformations Constraint Optimization (insertion)

Least Squares Fitting of Data

Physics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum. Physics 3A: Linear Momentum

Linear Momentum. Center of Mass.

Gradient Descent Learning and Backpropagation

Least Squares Fitting of Data

total If no external forces act, the total linear momentum of the system is conserved. This occurs in collisions and explosions.

Lecture 21: Numerical methods for pricing American type derivatives

Perceptual Organization (IV)

LECTURE :FACTOR ANALYSIS

Excess Error, Approximation Error, and Estimation Error

Physics 231. Topic 8: Rotational Motion. Alex Brown October MSU Physics 231 Fall

1.3 Hence, calculate a formula for the force required to break the bond (i.e. the maximum value of F)

Slobodan Lakić. Communicated by R. Van Keer

Multiple view geometry

Final Exam Solutions, 1998

Demonstrate Teaching FEM, Particle Dynamics, and Simple MD Codes Using Matlab & Mathematica

1 Matrix representations of canonical matrices

Multipoint Analysis for Sibling Pairs. Biostatistics 666 Lecture 18

Description of the Force Method Procedure. Indeterminate Analysis Force Method 1. Force Method con t. Force Method con t

Quantum Mechanics I - Session 4

BAYESIAN CURVE FITTING USING PIECEWISE POLYNOMIALS. Dariusz Biskup

XII.3 The EM (Expectation-Maximization) Algorithm

CHAPTER 10 ROTATIONAL MOTION

Our focus will be on linear systems. A system is linear if it obeys the principle of superposition and homogenity, i.e.

System in Weibull Distribution

COMP th April, 2007 Clement Pang

Recap: the SVM problem

Lecture Topics VMSC Prof. Dr.-Ing. habil. Hermann Lödding Prof. Dr.-Ing. Wolfgang Hintze. PD Dr.-Ing. habil.

On the Eigenspectrum of the Gram Matrix and the Generalisation Error of Kernel PCA (Shawe-Taylor, et al. 2005) Ameet Talwalkar 02/13/07

p p +... = p j + p Conservation Laws in Physics q Physical states, process, and state quantities: Physics 201, Lecture 14 Today s Topics

10/2/2003 PHY Lecture 9 1

Important Instructions to the Examiners:

Physics 231. Topic 8: Rotational Motion. Alex Brown October MSU Physics 231 Fall

On Pfaff s solution of the Pfaff problem

ENVIRONMENT OBSERVATION BY STRUCTURE FROM MOTION WITH AN OMNI-DIRECTIONAL CAMERA. Tomoaki Harada, Atsushi Yamashita and Toru Kaneko

Linear Momentum. Center of Mass.

Chemical Engineering 160/260 Polymer Science and Engineering. Lecture 10 - Phase Equilibria and Polymer Blends February 7, 2001

Einstein Summation Convention

Rectilinear motion. Lecture 2: Kinematics of Particles. External motion is known, find force. External forces are known, find motion

Applied Mathematics Letters

COS 511: Theoretical Machine Learning

Report on Image warping

Minimizing Algebraic Error in Geometric Estimation Problems

p 1 c 2 + p 2 c 2 + p 3 c p m c 2

An Optimal Bound for Sum of Square Roots of Special Type of Integers

Computational and Statistical Learning theory Assignment 4

Denote the function derivatives f(x) in given points. x a b. Using relationships (1.2), polynomials (1.1) are written in the form

Error Bars in both X and Y

Xiangwen Li. March 8th and March 13th, 2001

11. Dynamics in Rotating Frames of Reference

Geometric Registration for Deformable Shapes. 2.1 ICP + Tangent Space optimization for Rigid Motions

Estimating the Fundamental Matrix by Transforming Image Points in Projective Space 1

Moments of Inertia. and reminds us of the analogous equation for linear momentum p= mv, which is of the form. The kinetic energy of the body is.

Collaborative Filtering Recommendation Algorithm

3. Be able to derive the chemical equilibrium constants from statistical mechanics.

ITERATIVE ESTIMATION PROCEDURE FOR GEOSTATISTICAL REGRESSION AND GEOSTATISTICAL KRIGING

PHYS 1443 Section 002 Lecture #20

Chapter 12 Lyes KADEM [Thermodynamics II] 2007

Page 1. SPH4U: Lecture 7. New Topic: Friction. Today s Agenda. Surface Friction... Surface Friction...

Efficient, General Point Cloud Registration with Kernel Feature Maps

Newton s Method for One - Dimensional Optimization - Theory

CS 468 Lecture 16: Isometry Invariance and Spectral Techniques

For all questions, answer choice E) NOTA" means none of the above answers is correct.

IV. Performance Optimization

1 Definition of Rademacher Complexity

Application of Kalman Filter in Track Prediction of Shuttlecock

From Biot-Savart Law to Divergence of B (1)

PHYS 705: Classical Mechanics. Calculus of Variations II

The Impact of the Earth s Movement through the Space on Measuring the Velocity of Light

1. Statement of the problem

Minimization of l 2 -Norm of the KSOR Operator

A Sub-Pixel Image Registration Technique with Applications to Defect Detection

The Geometry of Logit and Probit

First Law: A body at rest remains at rest, a body in motion continues to move at constant velocity, unless acted upon by an external force.

Spring 2002 Lecture #13

The Feynman path integral

Chapter 8. Momentum Impulse and Collisions. Analysis of motion: 2 key ideas. Newton s laws of motion. Conservation of Energy

Estimating the Odometry Error of a Mobile Robot during Navigation

MTH 263 Practice Test #1 Spring 1999

Support Vector Machines. Vibhav Gogate The University of Texas at dallas

Ch 35 Images. Eunil Won Department of Physics Korea University. Fundamentals of Physics by Eunil Won, Korea University 1

Preference and Demand Examples

Transfer Functions. Convenient representation of a linear, dynamic model. A transfer function (TF) relates one input and one output: ( ) system

Efficiency Improvement of Measurement Pose Selection Techniques in Robot Calibration

where λ = Z/f. where a 3 4 projection matrix represents a map from 3D to 2D. Part I: Single and Two View Geometry Internal camera parameters

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

Implicit scaling of linear least squares problems

Introductory Optomechanical Engineering. 2) First order optics

Transcription:

Lecture Caera Models Caera Calbraton rofessor Slvo Savarese Coputatonal Vson and Geoetry Lab Slvo Savarese Lecture - - Jan- 5

Lecture Caera Models Caera Calbraton Recap of caera odels Caera calbraton proble Caera calbraton wth radal dstorton Exaple Readng: [F] Chapter Geoetrc Caera Calbraton [HZ] Chapter 7 Coputaton of Caera Matrx Soe sldes n ths lecture are courtesy to rofs. J. once, F-F L Slvo Savarese Lecture - 4- Jan- 5

rojectve caera nhole perspectve projecton f f focal length

rojectve caera nhole perspectve projecton f y y c x c f focal length u o, v o offset (note a dfferent conventon w.r.t. lecture ) C [u o, v o ] x

rojectve caera f Unts: k,l [pxel/] α, β f [] Non-square pxels [pxel] f focal length u o, v o offset α, β non-square pxels

rojectve caera f y y c θ C [u o, v o ] x c x f focal length u o, v o offset α, β non-square pxels θ skew angle

rojectve caera f ( ( ( ( '( α α cotθ u o β snθ v o K has 5 degrees of freedo ) + ( + ( + ( + ( * + ' x y z ) + + + + * f focal length u o, v o offset α, β non-square pxels θ skew angle

rojectve caera j w R, f O w k w ' R 4 4 w f focal length u o, v o offset α, β non-square pxels θ skew angle R, rotaton, translaton

rojectve caera j w R, f O w k w M w K[ R ] w Internal paraeters External paraeters f focal length u o, v o offset α, β non-square pxels θ skew angle R, rotaton, translaton

rojectve caera j w R, f O w k w ' M K R w W W W W w4 4 M E ' E ( w, w ) w w

heore (Faugeras, 99) [ R ] [ K R K ] [ A b] M K [Eq.; lect. ] A a a a

Exercse R, j w f k w O w w à M K R ' E ( w w, w w ) K I ( f x w z w, f y w z w ) f f w x w y w z w

rojectve caera p q r f O Q R

Weak perspectve projecton When the relatve scene depth s sall copared to ts dstance fro the caera p f R R_ q r Q Q_ O zo _ π

Weak perspectve projecton f zo p R_ q r Q Q_ O R x' f ' z x y' f ' z y x' f ' z x y' f ' z y _ π Magnfcaton

Weak perspectve projecton f zo p R_ q r Q Q_ O π R _ rojectve (perspectve) Weak perspectve M K[ R ] A b v à M A b

W W W ), ( w w E M w b A M agnfcaton M w v b A M W W W W ), ( w w w w E erspectve Weak perspectve

Orthographc (affne) projecton Dstance fro center of projecton to age plane s nfnte x' f ' z x y' f ' z y x' x y' y

ros and Cons of hese Models Weak perspectve uch spler ath. Accurate when object s sall and dstant. Most useful for recognton. nhole perspectve uch ore accurate for scenes. Used n structure fro oton or SLAM.

Weak perspectve projecton he Kangx Eperor's Southern Inspecton our (69-698) By Wang Hu

Weak perspectve projecton he Kangx Eperor's Southern Inspecton our (69-698) By Wang Hu

Lecture Caera Calbraton Recap of caera odels Caera calbraton proble Caera calbraton wth radal dstorton Exaple Readng: [F] Chapter Geoetrc Caera Calbraton [HZ] Chapter 7 Coputaton of Caera Matrx Soe sldes n ths lecture are courtesy to rofs. J. once, F-F L Slvo Savarese Lecture - 4- Jan- 5

rojectve caera v u cot K o o snθ β θ α α R r r r z y x t t t 4 M w [ ] w K R Internal paraeters External paraeters

Goal of calbraton M w [ ] w K R Estate ntrnsc and extrnsc paraeters fro or ultple ages Change notaton: w p v u cot K o o snθ β θ α α R r r r z y x t t t 4

Calbraton roble Calbraton rg j C n wth known postons n [O w, w,j w,k w ]

Calbraton roble Calbraton rg age p j C n wth known postons n [O w, w,j w,k w ] p, p n known postons n the age Goal: copute ntrnsc and extrnsc paraeters

Calbraton roble Calbraton rg age p j C How any correspondences do we need? M has unknown We need equatons 6 correspondences would do t

Calbraton roble Calbraton rg age p j C In practce, usng ore than 6 correspondences enables ore robust results

Calbraton roble j C M v u p M n pxels [Eq. ]

Calbraton roble u ) ( v ) ( u v v u ) ( v u ( ) [Eq. ] [Eqs. ]

Calbraton roble ) ( v ) ( u ) ( v ) ( u ) ( n n n v ) ( n n n u [Eqs. ]

Block Matrx Multplcaton B B B B B A A A A A What s AB? + + + + B A B A B A B A B A B A B A B A AB

Calbraton roble u ( ) + v ( ) + known unknown [Eq. 4] u v n n ( n ) + n ( n ) + n Hoogenous lnear syste u v n un n n un n ' def x4 n x def 4x x

Hoogeneous M x N Lnear Systes Mnuber of equatons n Nnuber of unknown M N Rectangular syste (M>N) s always a soluton o fnd non-zero soluton Mnze under the constrant

Calbraton roble How do we solve ths hoogenous lnear syste? Va SVD decoposton

Calbraton roble SVD decoposton of U D V n Last colun of V gves Why? See pag 59 of HZ def M

Extractng caera paraeters M ρ

Extractng caera paraeters A A a a a [ ] K R a ± ρ b b b b Estated values ( ) a a ρ u o ) ( v o a a ρ ( ) ( ) cos a a a a a a a a θ Intrnsc b v u cot K o o snθ β θ α α ρ M Box

heore (Faugeras, 99)

Extractng caera paraeters A a a a [ ] K R b b b b Estated values Intrnsc θ ρ α sn a a θ ρ β sn a a ρ A b

Extractng caera paraeters A A a a a b [ ] K R b b b b Estated values Extrnsc ( ) a a a a r a a r ± r r r b K ρ ρ

Degenerate cases s cannot le on the sae plane onts cannot le on the ntersecton curve of two quadrc surfaces

Lecture Caera Calbraton Recap of projectve caeras Caera calbraton proble Caera calbraton wth radal dstorton Exaple Readng: [F] Chapter Geoetrc Caera Calbraton [HZ] Chapter 7 Coputaton of Caera Matrx Soe sldes n ths lecture are courtesy to rofs. J. once, F-F L Slvo Savarese Lecture - 4- Jan- 5

Radal Dstorton Iage agnfcaton (n)decreases wth dstance fro the optcal axs Caused by perfect lenses Devatons are ost notceable for rays that pass through the edge of the lens No dstorton n cushon Barrel

Radal Dstorton p v u M λ λ d v v u c v b u a d + + u o odel radal behavor ± p p κ p d λ olynoal functon Dstorton coeffcent [Eq. 5] [Eq. 6] S λ Iage agnfcaton decreases wth dstance fro the optcal center

v u p Q q q q q q q q p v u M λ λ Q v u q q q q Is ths a lnear syste of equatons? Radal Dstorton No why? [Eqs.7]

General Calbraton roble u v q q q q n X f (Q) [Eq.8] paraeters easureents f( ) s the nonlnear appng -Newton Method -Levenberg-Marquardt Algorth Iteratve, starts fro ntal soluton May be slow f ntal soluton far fro real soluton Estated soluton ay be functon of the ntal soluton Newton requres the coputaton of J, H Levenberg-Marquardt doesn t requre the coputaton of H

General Calbraton roble u v q q q q n X f (Q) [Eq.8] paraeters easureents f( ) s the nonlnear appng A possble algorth. Solve lnear part of the syste to fnd approxated soluton. Use ths soluton as ntal condton for the full syste. Solve full syste usng Newton or L.M.

ypcal assuptons: - zero-skew, square pxel - u o, v o known center of the age General Calbraton roble X f (Q) easureent paraeter f( ) s nonlnear q q q q v u

Can we estate and and gnore the radal dstorton? Radal Dstorton v u p λ d v u Hnt: slope v u q q q q

Estatng and Radal Dstorton v u p λ ) ( ) ( u v ) ( ) ( u v ) ( ) ( n n n n u v L n n sa [87] v u ) ( ) ( ) ( ) ( Get and by SVD [Eq.9] [Eq.] [Eq.]

Once that and are estated Radal Dstorton v u p λ s non lnear functon of λ here are soe degenerate confguratons for whch and cannot be coputed,,

Lecture Caera Calbraton Recap of projectve caeras Caera calbraton proble Caera calbraton wth radal dstorton Exaple Readng: [F] Chapter Geoetrc Caera Calbraton [HZ] Chapter 7 Coputaton of Caera Matrx Soe sldes n ths lecture are courtesy to rofs. J. once, F-F L Slvo Savarese Lecture - 4- Jan- 5

Calbraton rocedure Caera Calbraton oolbox for Matlab J. Bouguet [998-] http://www.vson.caltech.edu/bouguetj/calb_doc/ndex.htlexaples

Calbraton rocedure

Calbraton rocedure

Calbraton rocedure

Calbraton rocedure

Calbraton rocedure

Calbraton rocedure

Calbraton rocedure

Next lecture Sngle vew reconstructon