Optical Flow, KLT Feature Tracker.

Similar documents
Computer Vision I. Announcements

Motion estimation. Digital Visual Effects Yung-Yu Chuang. with slides by Michael Black and P. Anandan

Iterative Image Registration: Lucas & Kanade Revisited. Kentaro Toyama Vision Technology Group Microsoft Research

Optical flow. Subhransu Maji. CMPSCI 670: Computer Vision. October 20, 2016

Part I: Thin Converging Lens

Image Alignment and Mosaicing

Announcements. Tracking. Comptuer Vision I. The Motion Field. = ω. Pure Translation. Motion Field Equation. Rigid Motion: General Case

6.869 Advances in Computer Vision. Prof. Bill Freeman March 1, 2005

Lecture 12. Local Feature Detection. Matching with Invariant Features. Why extract features? Why extract features? Why extract features?

Motion Estimation (I)

Lecture Outline. Basics of Spatial Filtering Smoothing Spatial Filters. Sharpening Spatial Filters

Video and Motion Analysis Computer Vision Carnegie Mellon University (Kris Kitani)

Motion Estimation (I) Ce Liu Microsoft Research New England

8. THEOREM If the partial derivatives f x. and f y exist near (a, b) and are continuous at (a, b), then f is differentiable at (a, b).

Keypoint extraction: Corners Harris Corners Pkwy, Charlotte, NC

. This is the Basic Chain Rule. x dt y dt z dt Chain Rule in this context.

Lucas-Kanade Optical Flow. Computer Vision Carnegie Mellon University (Kris Kitani)

Primary dependent variable is fluid velocity vector V = V ( r ); where r is the position vector

Chapter 4 Imaging. Lecture 21. d (110) Chem 793, Fall 2011, L. Ma

Geometry review, part I

Department of Physics PHY 111 GENERAL PHYSICS I

y,z the subscript y, z indicating that the variables y and z are kept constant. The second partial differential with respect to x is written x 2 y,z

Shai Avidan Tel Aviv University

CS 335 Graphics and Multimedia. 2D Graphics Primitives and Transformation

AP Calculus Notes: Unit 1 Limits & Continuity. Syllabus Objective: 1.1 The student will calculate limits using the basic limit theorems.

Dynamics of the Atmosphere 11:670:324. Class Time: Tuesdays and Fridays 9:15-10:35

Chapter 8 Laminar Flows with Dependence on One Dimension

Feature extraction: Corners and blobs

Lesson 6: Apparent weight, Radial acceleration (sections 4:9-5.2)

A Unied Factorization Algorithm for Points, Line Segments and. Planes with Uncertainty Models. Daniel D. Morris and Takeo Kanade

MOTION IN 2-DIMENSION (Projectile & Circular motion And Vectors)

Chapter 11 Collision Theory

y2 = 0. Show that u = e2xsin(2y) satisfies Laplace's equation.

Tangent Line Approximations

Visual Object Recognition

Local Features (contd.)

Image Alignment and Mosaicing Feature Tracking and the Kalman Filter

DO PHYSICS ONLINE. WEB activity: Use the web to find out more about: Aristotle, Copernicus, Kepler, Galileo and Newton.

Invariant local features. Invariant Local Features. Classes of transformations. (Good) invariant local features. Case study: panorama stitching

Dynamics ( 동역학 ) Ch.2 Motion of Translating Bodies (2.1 & 2.2)

Mat 241 Homework Set 7key Due Professor David Schultz

SUBJECT:ENGINEERING MATHEMATICS-I SUBJECT CODE :SMT1101 UNIT III FUNCTIONS OF SEVERAL VARIABLES. Jacobians

Relating axial motion of optical elements to focal shift

Syllabus Objective: 2.9 The student will sketch the graph of a polynomial, radical, or rational function.

a by a factor of = 294 requires 1/T, so to increase 1.4 h 294 = h

An example of Lagrangian for a non-holonomic system

EDGES AND CONTOURS(1)

Chapter 3 Motion in a Plane

Differential Equations

III. Relative Velocity

Relating axial motion of optical elements to focal shift

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

Example: When describing where a function is increasing, decreasing or constant we use the x- axis values.

Lecture 13: Tracking mo3on features op3cal flow

Interest Operators. All lectures are from posted research papers. Harris Corner Detector: the first and most basic interest operator

New Functions from Old Functions

MATRIX TRANSFORMATIONS

Space Coordinates and Vectors in Space. Coordinates in Space

DYNAMICS. Kinematics of Particles Engineering Dynamics Lecture Note VECTOR MECHANICS FOR ENGINEERS: Eighth Edition CHAPTER

Lab on Taylor Polynomials. This Lab is accompanied by an Answer Sheet that you are to complete and turn in to your instructor.

CISE-301: Numerical Methods Topic 1:

1. Interference condition. 2. Dispersion A B. As shown in Figure 1, the path difference between interfering rays AB and A B is a(sin

Elaborazione delle Immagini Informazione multimediale - Immagini. Raffaella Lanzarotti

CAP 5415 Computer Vision

Quadratic Functions. The graph of the function shifts right 3. The graph of the function shifts left 3.

Math 20C. Lecture Examples.

Mathematical Preliminaries. Developed for the Members of Azera Global By: Joseph D. Fournier B.Sc.E.E., M.Sc.E.E.

6.1 The Linear Elastic Model

3. Several Random Variables

TFY4102 Exam Fall 2015

Instance-level l recognition. Cordelia Schmid INRIA

TRACKING and DETECTION in COMPUTER VISION

Introduction to motion correspondence

Section 6: PRISMATIC BEAMS. Beam Theory

CSE 252B: Computer Vision II

CEE598 - Visual Sensing for Civil Infrastructure Eng. & Mgmt.

Flip-Flop Functions KEY

Objectives. By the time the student is finished with this section of the workbook, he/she should be able

Interest Point Detection. Lecture-4

Chapter 2 Section 3. Partial Derivatives

Wavelets bases in higher dimensions

( x) f = where P and Q are polynomials.

Chapter 8: MULTIPLE CONTINUOUS RANDOM VARIABLES

Fluid Physics 8.292J/12.330J

Time Series Analysis for Quality Improvement: a Soft Computing Approach

Clicker questions. Clicker question 2. Clicker Question 1. Clicker question 2. Clicker question 1. the answers are in the lower right corner

EVALUATE: If the angle 40 is replaces by (cable B is vertical), then T = mg and

MATRIX REPRESENTATIONS OF THE COORDINATES OF THE ATOMS OF THE MOLECULES OF H 2 O AND NH 3 IN OPERATIONS OF SYMMETRY

2) If a=<2,-1> and b=<3,2>, what is a b and what is the angle between the vectors?

VECTORS IN THREE DIMENSIONS

Wave Motion A wave is a self-propagating disturbance in a medium. Waves carry energy, momentum, information, but not matter.

Instance-level recognition: Local invariant features. Cordelia Schmid INRIA, Grenoble

Mat 267 Engineering Calculus III Updated on 9/19/2010

z-axis SUBMITTED BY: Ms. Harjeet Kaur Associate Professor Department of Mathematics PGGCG 11, Chandigarh y-axis x-axis

Gaussian Plane Waves Plane waves have flat emag field in x,y Tend to get distorted by diffraction into spherical plane waves and Gaussian Spherical

Doppler shifts in astronomy

Polynomials, Linear Factors, and Zeros. Factor theorem, multiple zero, multiplicity, relative maximum, relative minimum

Optical Flow, Motion Segmentation, Feature Tracking

Chapter 6 Momentum Transfer in an External Laminar Boundary Layer

UNDERSTAND MOTION IN ONE AND TWO DIMENSIONS

Transcription:

Optical Flow, KL Feature racker E-mail: hogijung@hanang.ac.kr

Motion in Computer Vision Motion Structure rom motion Detection/segmentation with direction [1] E-mail: hogijung@hanang.ac.kr

Motion Field.s.. Optical Flow [2], [3] Motion Field: an ideal representation o 3D motion as it is projected onto a camera image. Optical Flow: the approimation (or estimate) o the motion ield which can be computed rom time-aring image sequences. Under the simpliing assumptions o 1) Lambertian surace, 2) pointwise light source at ininit, and 3) no photometric distortion. E-mail: hogijung@hanang.ac.kr

Motion Field [2] An ideal representation o 3D motion as it is projected onto a camera image. he time deriatie o the image position o all image points gien that the correspond to ied 3D points. ield : position ector he motion ield is deined as where P is a point in the scene where is the distance to that scene point. V is the relatie motion between the camera and the scene, is the translational component o the motion, and ω is the angular elocit o the motion. E-mail: hogijung@hanang.ac.kr

E-mail: hogijung@hanang.ac.kr Motion Field [2] Motion Field [2] P p (1) 3D point P (,,) and 2D point p (,), ocal length Motion ield can be obtained b taking the time deriatie o (1) 2 2 V V V V (2)

E-mail: hogijung@hanang.ac.kr Motion Field [2] Motion Field [2] P V (3) he motion o 3D point P, V is deined as low B substituting (3) into (2), the basic equations o the motion ield is acquired 2 2 (4) V V V

E-mail: hogijung@hanang.ac.kr Motion Field [2] Motion Field [2] he motion ield is the sum o two components, one o which depends on translation onl, the other on rotation onl. 2 2 2 2 ranslational components Rotational components (4) (5) (6)

E-mail: hogijung@hanang.ac.kr Motion Field: Pure ranslation [2] Motion Field: Pure ranslation [2] I there is no rotational motion, the resulting motion ield has a peculiar spatial structure. I (5) is regarded as a unction o 2D point position, (5) 0 0 I ( 0, 0 ) is deined as in (6) 0 0 (6) (7)

E-mail: hogijung@hanang.ac.kr Motion Field: Pure ranslation [2] Motion Field: Pure ranslation [2] Equation (7) sa that the motion ield o a pure translation is radial. In particular, i <0, the ectors point awa rom p 0 ( 0, 0 ), which is called the ocus o epansion (FOE). I >0, the motion ield ectors point towards p 0, which is called the ocus o contraction. I =0, rom (5), all the motion ield ectors are parallel. (5) (8)

Motion Field: Motion Paralla [6] Equation (8) sa that their lengths are inersel proportional to the depth o the corresponding 3D points. (8) http://upload.wikimedia.org /wikipedia/commons/a/ab/ Paralla.gi his animation is an eample o paralla. As the iewpoint moes side to side, the objects in the distance appear to moe more slowl than the objects close to the camera [6]. E-mail: hogijung@hanang.ac.kr

E-mail: hogijung@hanang.ac.kr Motion Field: Motion Paralla [2] Motion Field: Motion Paralla [2] I two 3D points are projected into one image point, that is coincident, rotational component will be the same. Notice that the motion ector V is about camera motion. (4) 2 2

Motion Field: Motion Paralla [2] he dierence o two points motion ield will be related with translation components. And, the will be radial w.r.t FOE or FOC. 1 1 1 1 0 1 2 1 2 1 1 1 1 0 1 2 1 2 FOC E-mail: hogijung@hanang.ac.kr

Motion Field: Motion Paralla [2] Motion Paralla he relatie motion ield o two instantaneousl coincident points: 1. Does not depend on the rotational component o motion 2. Points towards (awa rom) the point p0, the anishing point o the translation direction. E-mail: hogijung@hanang.ac.kr

Motion Field: Pure Rotation w.r.t -ais [7] I there is no translation motion and rotation w.r.t - and z- ais, rom (4) 2 E-mail: hogijung@hanang.ac.kr

Motion Field: Pure Rotation w.r.t -ais [7] ranslational Motion Distance to the point,, is constant. Rotational Motion 1 2 Distance to the point,, is changing. According to, is changing, too. E-mail: hogijung@hanang.ac.kr

Estimation o the Optical Flow [4] E-mail: hogijung@hanang.ac.kr

Estimation o the Optical Flow [4] he Image Brightness Constanc Equation [2] E-mail: hogijung@hanang.ac.kr

Estimation o the Optical Flow [4] I V I t V I t I Assumption he image brightness is continuous and dierentiable as man times as needed in both the spatial and temporal domain. he image brightness can be regarded as a plane in a small area. E-mail: hogijung@hanang.ac.kr

Optical Field: Aperture Problem [2], [4], [9] he component o the motion ield in the direction orthogonal to the spatial image gradient is not constrained b the image brightness constanc equation. Gien local inormation can determine component o optical low ector onl in direction o brightness gradient. E-mail: hogijung@hanang.ac.kr

Optical Field: Aperture Problem [2], [9] he aperture problem. he grating appears to be moing down and to the right, perpendicular to the orientation o the bars. But it could be moing in man other directions, such as onl down, or onl to the right. It is impossible to determine unless the ends o the bars become isible in the aperture. http://upload.wikimedia.org/wikipedia/commons//0/aperture_pr oblem_animated.gi E-mail: hogijung@hanang.ac.kr

Optical Field: Methods or Determining Optical Flow [4] E-mail: hogijung@hanang.ac.kr

Optical Field: Phase Correlation Method [10] E-mail: hogijung@hanang.ac.kr

Optical Field: Phase Correlation Method [10] E-mail: hogijung@hanang.ac.kr

Optical Field: Phase Correlation Method [10] E-mail: hogijung@hanang.ac.kr

Optical Field: Lucas-Kanade Method [8] Assuming that the optical low (V, ) is constant in a small window o size mm with m>1, which is center at (, ) and numbering the piels within as 1 n, n=m 2, a set o equations can be ound: I V I t E-mail: hogijung@hanang.ac.kr

Optical Field: Lucas-Kanade Method [8] c.) Harris corner detector E-mail: hogijung@hanang.ac.kr

KL Feature racker [11] E-mail: hogijung@hanang.ac.kr

KL Feature racker [11] F () 로 weighting E-mail: hogijung@hanang.ac.kr

KL Feature racker [11] I V I t KL: An Implementation o the Kanade-Lucas-omasi Feature racker http://www.ces.clemson.edu/~stb/klt/ E-mail: hogijung@hanang.ac.kr

Lucas-Kanade Method: A Uniing Framework [12] Since the Lucas-Kanade algorithm was proposed in 1981 image alignment has become one o the most widel used techniques in computer ision. Besides optical low, some o its other applications include - tracking (Black and Jepson, 1998; Hager and Belhumeur, 1998), - parametric and laered motion estimation (Bergen et al., 1992), - mosaic construction (Shum and Szeliski, 2000), - medical image registration (Christensen and Johnson, 2001), - ace coding (Baker and Matthews, 2001; Cootes et al., 1998). E-mail: hogijung@hanang.ac.kr

Lucas-Kanade Method: A Uniing Framework [12] E-mail: hogijung@hanang.ac.kr

Lucas-Kanade Method: A Uniing Framework [12] Minimizing the epression in Equation (1) is a non-linear optimization task een i W(; p) is linear in p because the piel alues I() are, in general, non-linear in. E-mail: hogijung@hanang.ac.kr

Lucas-Kanade Method: A Uniing Framework [12] E-mail: hogijung@hanang.ac.kr

Lucas-Kanade Method: A Uniing Framework [12] E-mail: hogijung@hanang.ac.kr

Lucas-Kanade Method: A Uniing Framework [12] Setting this epression to equal zero and soling gies the closed orm solution or the minimum o the epression in Equation (6) as: SD E-mail: hogijung@hanang.ac.kr

Lucas-Kanade Method: A Uniing Framework [12] E-mail: hogijung@hanang.ac.kr

Lucas-Kanade Method: A Uniing Framework [12] E-mail: hogijung@hanang.ac.kr

Reerences 1. Richard Szeliski, Dense motion estimation, Computer Vision: Algorithms and Applications, 19 June 2009 (drat), pp. 383-426. 2. Emanuele rucco, Alessandro Verri, 8. Motion, Introductor echniques or 3-D Computer Vision, Prentice Hall, New Jerse 1998, pp.177-218. 3. Wikipedia, Motion ield, aailable on www.wikipedia.org. 4. Wikipedia, Optical low, aailable on www.wikipedia.org. 5. Alessandro Verri, Emanuele rucco, Finding the Epipole rom Uncalibrated Optical Flow, BMVC 1997, aailable on http://www.bma.ac.uk/bmc/1997/papers/052/bmc.html. 6. Wikipedia, Paralle, aailable on www.wikipedia.org. 7. Jae Ku Suhr, Ho Gi Jung, Kwanghuk Bae, Jaihie Kim, Outlier rejection or cameras on intelligent ehicles, Pattern Recognition Letters 29 (2008) 828-840. 8. Wikipedia, Lucas-Kanade Optical Flow Method, aailable on www.wikipedia.org. 9. Wikipedia, Aperture Problem, aailable on www.wikipedia.org. 10. Wikipedia, Phase correlation, aailable on http://en.wikipedia.org/wiki/phase_correlation 11. Wikipedia, Kanade-Lucas-omasi eature tracker, aailable on http://en.wikipedia.org/wiki/kanade%e2%80%93lucas%e2%80%93omasi_eature_track er 12. Simon Baker, Iain Matthews, Lucas-Kanade 20 ears: A Uniing Framework, International Journal o Computer Vision, 56(3), 221-255, 2004. E-mail: hogijung@hanang.ac.kr