Visual motion. Many slides adapted from S. Seitz, R. Szeliski, M. Pollefeys

Similar documents
Global Motion. Estimate motion using all pixels in the image. Parametric flow gives an equation, which describes optical flow for each pixel.

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

Lecture 8 Wrap-up Part1, Matlab

Matching patterns of line segments by eigenvector decomposition

1 Linear Least Squares

Multiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution

Matrices and Determinants

Chapter 16. Molecular Symmetry

Partial Derivatives. Limits. For a single variable function f (x), the limit lim

Reference. Vector Analysis Chapter 2

ARITHMETIC OPERATIONS. The real numbers have the following properties: a b c ab ac

g i fφdx dx = x i i=1 is a Hilbert space. We shall, henceforth, abuse notation and write g i f(x) = f

8.6 The Hyperbola. and F 2. is a constant. P F 2. P =k The two fixed points, F 1. , are called the foci of the hyperbola. The line segments F 1

NUMERICAL INTEGRATION. The inverse process to differentiation in calculus is integration. Mathematically, integration is represented by.

T 1 T 2 T 3 T 4 They may be illustrated by triangular patterns of numbers (hence their name) as shown:

Math 520 Final Exam Topic Outline Sections 1 3 (Xiao/Dumas/Liaw) Spring 2008

2A1A Vector Algebra and Calculus I

Satellite Retrieval Data Assimilation

Problems for HW X. C. Gwinn. November 30, 2009

2. VECTORS AND MATRICES IN 3 DIMENSIONS

Best Approximation in the 2-norm

SOLVING SYSTEMS OF EQUATIONS, ITERATIVE METHODS

Lecture 4 Single View Metrology

CS434a/541a: Pattern Recognition Prof. Olga Veksler. Lecture 2

Introduction. Definition of Hyperbola

P 1 (x 1, y 1 ) is given by,.

Thomas Whitham Sixth Form

Calculus - Activity 1 Rate of change of a function at a point.

Things to Memorize: A Partial List. January 27, 2017

Jack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah

Lecture Solution of a System of Linear Equation

4 The dynamical FRW universe

Section 14.3 Arc Length and Curvature

Purpose of the experiment

First midterm topics Second midterm topics End of quarter topics. Math 3B Review. Steve. 18 March 2009

Optical Flow, Motion Segmentation, Feature Tracking

A Signal-Level Fusion Model for Image-Based Change Detection in DARPA's Dynamic Database System

Some Methods in the Calculus of Variations

Non-Linear & Logistic Regression

Theoretical foundations of Gaussian quadrature

( dg. ) 2 dt. + dt. dt j + dh. + dt. r(t) dt. Comparing this equation with the one listed above for the length of see that

LINEAR ALGEBRA APPLIED

Chapter 3 Polynomials

A Planar Perspective Image Matching using Point Correspondences and Rectangle-to-Quadrilateral Mapping

Chapter 6 Polarization and Crystal Optics

Chapter 3 Solving Nonlinear Equations

Duality # Second iteration for HW problem. Recall our LP example problem we have been working on, in equality form, is given below.

19 Optimal behavior: Game theory

along the vector 5 a) Find the plane s coordinate after 1 hour. b) Find the plane s coordinate after 2 hours. c) Find the plane s coordinate

a * a (2,1) 1,1 0,1 1,1 2,1 hkl 1,0 1,0 2,0 O 2,1 0,1 1,1 0,2 1,2 2,2

Heat flux and total heat

Counting intersections of spirals on a torus

MATRICES AND VECTORS SPACE

Chapter 2. Vectors. 2.1 Vectors Scalars and Vectors

Is there an easy way to find examples of such triples? Why yes! Just look at an ordinary multiplication table to find them!

Operations with Polynomials

Multiple Integrals. Review of Single Integrals. Planar Area. Volume of Solid of Revolution

USA Mathematical Talent Search Round 1 Solutions Year 21 Academic Year

ADVANCEMENT OF THE CLOSELY COUPLED PROBES POTENTIAL DROP TECHNIQUE FOR NDE OF SURFACE CRACKS

2.4 Linear Inequalities and Interval Notation

In-Class Problems 2 and 3: Projectile Motion Solutions. In-Class Problem 2: Throwing a Stone Down a Hill

How do we solve these things, especially when they get complicated? How do we know when a system has a solution, and when is it unique?

KINEMATICS OF RIGID BODIES

Plate Theory. Section 11: PLATE BENDING ELEMENTS

Graduate Students do all problems. Undergraduate students choose three problems.

Data Assimilation. Alan O Neill Data Assimilation Research Centre University of Reading

Plate Theory. Section 13: PLATE BENDING ELEMENTS

u t = k 2 u x 2 (1) a n sin nπx sin 2 L e k(nπ/l) t f(x) = sin nπx f(x) sin nπx dx (6) 2 L f(x 0 ) sin nπx 0 2 L sin nπx 0 nπx

SECTION 9-4 Translation of Axes

APPLICATIONS OF THE DEFINITE INTEGRAL

AQA Further Pure 2. Hyperbolic Functions. Section 2: The inverse hyperbolic functions

03 Qudrtic Functions Completing the squre: Generl Form f ( x) x + x + c f ( x) ( x + p) + q where,, nd c re constnts nd 0. (i) (ii) (iii) (iv) *Note t

Fig. 1. Open-Loop and Closed-Loop Systems with Plant Variations

M344 - ADVANCED ENGINEERING MATHEMATICS

NOT TO SCALE. We can make use of the small angle approximations: if θ á 1 (and is expressed in RADIANS), then

Chapter 3. Vector Spaces

#6A&B Magnetic Field Mapping

Eigen Values and Eigen Vectors of a given matrix

fractions Let s Learn to

A Brief Review on Akkar, Sandikkaya and Bommer (ASB13) GMPE

Matrix Eigenvalues and Eigenvectors September 13, 2017

State space systems analysis (continued) Stability. A. Definitions A system is said to be Asymptotically Stable (AS) when it satisfies

The Properties of Stars

4.5 JACOBI ITERATION FOR FINDING EIGENVALUES OF A REAL SYMMETRIC MATRIX. be a real symmetric matrix. ; (where we choose θ π for.

The Islamic University of Gaza Faculty of Engineering Civil Engineering Department. Numerical Analysis ECIV Chapter 11

5.3 Nonlinear stability of Rayleigh-Bénard convection

Optimization Lecture 1 Review of Differential Calculus for Functions of Single Variable.

We will see what is meant by standard form very shortly

- 5 - TEST 2. This test is on the final sections of this session's syllabus and. should be attempted by all students.

Lecture 8. Newton s Laws. Applications of the Newton s Laws Problem-Solving Tactics. Physics 105; Fall Inertial Frames: T = mg

Name Solutions to Test 3 November 8, 2017

Scientific notation is a way of expressing really big numbers or really small numbers.

Numerical Linear Algebra Assignment 008

Jackson 2.26 Homework Problem Solution Dr. Christopher S. Baird University of Massachusetts Lowell

MATHEMATICS PART A. 1. ABC is a triangle, right angled at A. The resultant of the forces acting along AB, AC

CMDA 4604: Intermediate Topics in Mathematical Modeling Lecture 19: Interpolation and Quadrature

Energy Bands Energy Bands and Band Gap. Phys463.nb Phenomenon

Math 113 Exam 2 Practice

Physics 207 Lecture 5

We partition C into n small arcs by forming a partition of [a, b] by picking s i as follows: a = s 0 < s 1 < < s n = b.

Transcription:

Visul motion Mn slides dpted from S. Seitz, R. Szeliski, M. Pollefes

Outline Applictions of segmenttion to video Motion nd perceptul orgniztion Motion field Opticl flow Motion segmenttion with lers

Video A video is sequence of frmes cptured over time Now our imge dt is function of spce (, ) nd time (t)

Applictions of segmenttion to video Bckground subtrction A sttic cmer is observing scene Gol: seprte the sttic bckground from the moving foreground

Applictions of segmenttion to video Bckground subtrction Form n initil bckground estimte For ech frme: Updte estimte using moving verge Subtrct the bckground estimte from the frme Lbel s foreground ech piel where the mgnitude of the difference is greter thn some threshold Use medin filtering to clen up the results

Applictions of segmenttion to video Bckground subtrction Shot boundr detection Commercil video is usull composed of shots or sequences showing the sme objects or scene Gol: segment video into shots for summriztion nd browsing (ech shot cn be represented b single kefrme in user interfce) Difference from bckground subtrction: the cmer is not necessril sttionr

Applictions of segmenttion to video Bckground subtrction Shot boundr detection For ech frme Compute the distnce between the current frme nd the previous one» Piel-b-piel differences» Differences of color histogrms» Block comprison If the distnce is greter thn some threshold, clssif the frme s shot boundr

Applictions of segmenttion to video Bckground subtrction Shot boundr detection Motion segmenttion Segment the video into multiple coherentl moving objects

Motion nd perceptul orgniztion Sometimes, motion is the onl cue

Motion nd perceptul orgniztion Sometimes, motion is the onl cue

Motion nd perceptul orgniztion Even impoverished motion dt cn evoke strong percept

Motion nd perceptul orgniztion Even impoverished motion dt cn evoke strong percept

Motion nd perceptul orgniztion Even impoverished motion dt cn evoke strong percept

Uses of motion Estimting 3D structure Segmenting objects bsed on motion cues Lerning dnmicl models Recognizing events nd ctivities Improving video qulit (motion stbiliztion)

Motion estimtion techniques Direct methods Directl recover imge motion t ech piel from sptio-temporl imge brightness vritions Dense motion fields, but sensitive to ppernce vritions Suitble for video nd when imge motion is smll Feture-bsed methods Etrct visul fetures (corners, tetured res) nd trck them over multiple frmes Sprse motion fields, but more robust trcking Suitble when imge motion is lrge (10s of piels)

Motion field The motion field is the projection of the 3D scene motion into the imge

Motion field nd prll P(t) is moving 3D point Velocit of scene point: V = dp/dt p(t) = ((t),(t)) is the projection of P in the imge Apprent velocit v in the imge: given b components v = d/dt nd v = d/dt These components re known s the motion field of the imge P(t) V p(t) v P(tdt) p(tdt)

Motion field nd prll V = ( V, V, VZ ) p = f P Z P(t) V P(tdt) To find imge velocit v, differentite p with respect to t (using quotient rule): v = f Z V V 2 Z z P v p(tdt) v = f V V Z z v = f V V Z z p(t) Imge motion is function of both the 3D motion (V) nd the depth of the 3D point (Z)

Motion field nd prll Pure trnsltion: V is constnt everwhere Z V V f v z = Z V V f v z = ), ( 1 v 0 p v z V Z = ( ) V f V f, v 0 =

Motion field nd prll Pure trnsltion: V is constnt everwhere V z is nonzero: v = 1 Z v 0 = ( v 0 V p), Ever motion vector points towrd (or w from) v 0, the vnishing point of the trnsltion direction z ( f V, f ) V

Motion field nd prll Pure trnsltion: V is constnt everwhere V z is nonzero: Ever motion vector points towrd (or w from) v 0, the vnishing point of the trnsltion direction V z is zero: v = 1 Z v 0 = ( v 0 V p), Motion is prllel to the imge plne, ll the motion vectors re prllel The length of the motion vectors is inversel proportionl to the depth Z z ( f V, f ) V

Opticl flow Definition: opticl flow is the pprent motion of brightness ptterns in the imge Idell, opticl flow would be the sme s the motion field Hve to be creful: pprent motion cn be cused b lighting chnges without n ctul motion Think of uniform rotting sphere under fied lighting vs. sttionr sphere under moving illumintion

Estimting opticl flow I(,,t 1) I(,,t) Given two subsequent frmes, estimte the pprent motion field u(,) nd v(,) between them Ke ssumptions Brightness constnc: projection of the sme point looks the sme in ever frme Smll motion: points do not move ver fr Sptil coherence: points move like their neighbors

Brightness Constnc Eqution: ), ( 1),, ( ),, ( ), ( t v u I t I = ), ( ), ( ),, ( 1),, ( v I u I t I t I Linerizing the right side using Tlor epnsion: The brightness constnc constrint I(,,t 1) I(,,t) 0 t I v I u I Hence,

The brightness constnc constrint I u v How mn equtions nd unknowns per piel? One eqution, two unknowns I = 0 Intuitivel, wht does this constrint men? I The component of the flow perpendiculr to the grdient (i.e., prllel to the edge) is unknown I ( u, v) I = t t 0

The brightness constnc constrint I u v How mn equtions nd unknowns per piel? One eqution, two unknowns I = 0 Intuitivel, wht does this constrint men? I The component of the flow perpendiculr to the grdient (i.e., prllel to the edge) is unknown I ( u, v) I = t t 0 If (u, v) stisfies the eqution, so does (uu, vv ) if I ( u', v') = 0 grdient (u,v) (u,v ) (uu,vv ) edge

The perture problem Perceived motion

The perture problem Actul motion

The brber pole illusion http://en.wikipedi.org/wiki/brberpole_illusion

The brber pole illusion http://en.wikipedi.org/wiki/brberpole_illusion

The brber pole illusion http://en.wikipedi.org/wiki/brberpole_illusion

Solving the perture problem How to get more equtions for piel? Sptil coherence constrint: pretend the piel s neighbors hve the sme (u,v) If we use 55 window, tht gives us 25 equtions per piel B. Lucs nd T. Knde. An itertive imge registrtion technique with n ppliction to stereo vision. In Proceedings of the Interntionl Joint Conference on Artificil Intelligence, pp. 674 679, 1981.

Solving the perture problem Lest squres problem: When is this sstem solvble? Wht if the window contins just single stright edge? B. Lucs nd T. Knde. An itertive imge registrtion technique with n ppliction to stereo vision. In Proceedings of the Interntionl Joint Conference on Artificil Intelligence, pp. 674 679, 1981.

Conditions for solvbilit Bd cse: single stright edge

Conditions for solvbilit Good cse

Lucs-Knde flow Overconstrined liner sstem Lest squres solution for d given b The summtions re over ll piels in the K K window B. Lucs nd T. Knde. An itertive imge registrtion technique with n ppliction to stereo vision. In Proceedings of the Interntionl Joint Conference on Artificil Intelligence, pp. 674 679, 1981.

Conditions for solvbilit Optiml (u, v) stisfies Lucs-Knde eqution When is this solvble? A T A should be invertible A T A entries should not be too smll (noise) A T A should be well-conditioned

Eigenvectors of A T A Recll the Hrris corner detector: M = A T A is the second moment mtri The eigenvectors nd eigenvlues of M relte to edge direction nd mgnitude The eigenvector ssocited with the lrger eigenvlue points in the direction of fstest intensit chnge The other eigenvector is orthogonl to it

Interpreting the eigenvlues Clssifiction of imge points using eigenvlues of the second moment mtri: λ 2 Edge λ 2 >> λ 1 Corner λ 1 nd λ 2 re lrge, λ 1 ~ λ 2 λ 1 nd λ 2 re smll Flt region Edge λ 1 >> λ 2 λ 1

Edge grdients ver lrge or ver smll lrge λ 1, smll λ 2

Low-teture region grdients hve smll mgnitude smllλ 1, smll λ 2

High-teture region grdients re different, lrge mgnitudes lrge λ 1, lrge λ 2

Wht re good fetures to trck? Recll the Hrris corner detector Cn mesure qulit of fetures from just single imge

Motion models Trnsltion Affine Perspective 3D rottion 2 unknowns 6 unknowns 8 unknowns 3 unknowns

Substituting into the brightness constnc eqution: v u 6 5 4 3 2 1 ), ( ), ( = = 0 t I v I u I Affine motion

0 ) ( ) ( 6 5 4 3 2 1 t I I I Substituting into the brightness constnc eqution: v u 6 5 4 3 2 1 ), ( ), ( = = Ech piel provides 1 liner constrint in 6 unknowns [ ] 2 = t I I I Err ) ( ) ( ) ( 6 5 4 3 2 1 r Lest squres minimiztion: Affine motion

Errors in Lucs-Knde The motion is lrge (lrger thn piel) Itertive refinement, corse-to-fine estimtion A point does not move like its neighbors Motion segmenttion Brightness constnc does not hold Do ehustive neighborhood serch with normlized correltion

Itertive Refinement Estimte velocit t ech piel using one itertion of Lucs nd Knde estimtion Wrp one imge towrd the other using the estimted flow field Refine estimte b repeting the process

Deling with lrge motions

Reduce the resolution!

Corse-to-fine opticl flow estimtion u=1.25 piels u=2.5 piels u=5 piels imge H1 u=10 piels imge I imge 2 Gussin prmid of imge 1 Gussin prmid of imge 2

Corse-to-fine opticl flow estimtion run itertive L-K run itertive L-K wrp & upsmple... imge J1 imge I imge 2 Gussin prmid of imge 1 Gussin prmid of imge 2

Motion segmenttion How do we represent the motion in this scene? J. Wng nd E. Adelson. Lered Representtion for Motion Anlsis. CVPR 1993.

Lered motion Brek imge sequence into lers ech of which hs coherent motion J. Wng nd E. Adelson. Lered Representtion for Motion Anlsis. CVPR 1993.

Wht re lers? Ech ler is defined b n lph msk nd n ffine motion model J. Wng nd E. Adelson. Lered Representtion for Motion Anlsis. CVPR 1993.

v u 6 5 4 3 2 1 ), ( ), ( = = Locl flow estimtes Motion segmenttion with n ffine model J. Wng nd E. Adelson. Lered Representtion for Motion Anlsis. CVPR 1993.

Motion segmenttion with n ffine model v u 6 5 4 3 2 1 ), ( ), ( = = Eqution of plne (prmeters 1, 2, 3 cn be found b lest squres) J. Wng nd E. Adelson. Lered Representtion for Motion Anlsis. CVPR 1993.

Motion segmenttion with n ffine model u(, v(, ) ) = = 1 4 2 5 3 6 1D emple Eqution of plne (prmeters 1, 2, 3 cn be found b lest squres) u(,) True flow Locl flow estimte Foreground Bckground Segmented estimte Line fitting Occlusion J. Wng nd E. Adelson. Lered Representtion for Motion Anlsis. CVPR 1993.

How do we estimte the lers? Compute locl flow in corse-to-fine fshion Obtin set of initil ffine motion hpotheses Divide the imge into blocks nd estimte ffine motion prmeters in ech block b lest squres Eliminte hpotheses with high residul error Perform k-mens clustering on ffine motion prmeters Merge clusters tht re close nd retin the lrgest clusters to obtin smller set of hpotheses to describe ll the motions in the scene Iterte until convergence: Assign ech piel to best hpothesis Piels with high residul error remin unssigned Perform region filtering to enforce sptil constrints Re-estimte ffine motions in ech region J. Wng nd E. Adelson. Lered Representtion for Motion Anlsis. CVPR 1993.

Emple result J. Wng nd E. Adelson. Lered Representtion for Motion Anlsis. CVPR 1993.