CS376 Computer Vision Lecture 6: Optical Flow

Similar documents
Optical Flow I. Guido Gerig CS 6320, Spring 2015

Tracking. Announcements

Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle

Online Convex Optimization Example And Follow-The-Leader

Two Coupled Oscillators / Normal Modes

WEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x

Ground Rules. PC1221 Fundamentals of Physics I. Kinematics. Position. Lectures 3 and 4 Motion in One Dimension. A/Prof Tay Seng Chuan

Lecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.

CS 4495 Computer Vision Tracking 1- Kalman,Gaussian

Using the Kalman filter Extended Kalman filter

translational component of a rigid motion appear to originate.

Single and Double Pendulum Models

t is a basis for the solution space to this system, then the matrix having these solutions as columns, t x 1 t, x 2 t,... x n t x 2 t...

3.1.3 INTRODUCTION TO DYNAMIC OPTIMIZATION: DISCRETE TIME PROBLEMS. A. The Hamiltonian and First-Order Conditions in a Finite Time Horizon

An introduction to the theory of SDDP algorithm

Computer Vision. Motion Extraction

Chapter 3 Boundary Value Problem

Non-uniform circular motion *

L07. KALMAN FILTERING FOR NON-LINEAR SYSTEMS. NA568 Mobile Robotics: Methods & Algorithms

Solutions for homework 12

Numerical Dispersion

Today: Falling. v, a

NEWTON S SECOND LAW OF MOTION

Week 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)

IB Physics Kinematics Worksheet

Technical Report Doc ID: TR March-2013 (Last revision: 23-February-2016) On formulating quadratic functions in optimization models.

Supplement for Stochastic Convex Optimization: Faster Local Growth Implies Faster Global Convergence

a. Show that these lines intersect by finding the point of intersection. b. Find an equation for the plane containing these lines.

Y 0.4Y 0.45Y Y to a proper ARMA specification.

KEY. Math 334 Midterm III Winter 2008 section 002 Instructor: Scott Glasgow

Then. 1 The eigenvalues of A are inside R = n i=1 R i. 2 Union of any k circles not intersecting the other (n k)

Math Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.

1. VELOCITY AND ACCELERATION

Lecture 9: September 25

Some Ramsey results for the n-cube

PENALIZED LEAST SQUARES AND PENALIZED LIKELIHOOD

Hamilton- J acobi Equation: Explicit Formulas In this lecture we try to apply the method of characteristics to the Hamilton-Jacobi equation: u t

Let us start with a two dimensional case. We consider a vector ( x,

In this chapter the model of free motion under gravity is extended to objects projected at an angle. When you have completed it, you should

State-Space Models. Initialization, Estimation and Smoothing of the Kalman Filter

Math Wednesday March 3, , 4.3: First order systems of Differential Equations Why you should expect existence and uniqueness for the IVP

Best test practice: Take the past test on the class website

PHYSICS 149: Lecture 9

Notes for Lecture 17-18

3, so θ = arccos

Kinematics Vocabulary. Kinematics and One Dimensional Motion. Position. Coordinate System in One Dimension. Kinema means movement 8.

Exercises: Similarity Transformation

Sequential Importance Resampling (SIR) Particle Filter

Hamilton- J acobi Equation: Weak S olution We continue the study of the Hamilton-Jacobi equation:

Linear Dynamic Models

Suggested Practice Problems (set #2) for the Physics Placement Test

Econ107 Applied Econometrics Topic 7: Multicollinearity (Studenmund, Chapter 8)

Lecture 10 Estimating Nonlinear Regression Models

Linear Response Theory: The connection between QFT and experiments

Linear Algebra Primer

Two Popular Bayesian Estimators: Particle and Kalman Filters. McGill COMP 765 Sept 14 th, 2017

CSE 3802 / ECE Numerical Methods in Scientific Computation. Jinbo Bi. Department of Computer Science & Engineering

ODEs II, Lecture 1: Homogeneous Linear Systems - I. Mike Raugh 1. March 8, 2004

Elements of Computer Graphics

Course Notes for EE227C (Spring 2018): Convex Optimization and Approximation

Some Basic Information about M-S-D Systems

Outline of Topics. Analysis of ODE models with MATLAB. What will we learn from this lecture. Aim of analysis: Why such analysis matters?

HW6: MRI Imaging Pulse Sequences (7 Problems for 100 pts)

Today: Graphing. Note: I hope this joke will be funnier (or at least make you roll your eyes and say ugh ) after class. v (miles per hour ) Time

Position, Velocity, and Acceleration

Lecture 4 Kinetics of a particle Part 3: Impulse and Momentum

Zürich. ETH Master Course: L Autonomous Mobile Robots Localization II

Lecture 2 October ε-approximation of 2-player zero-sum games

KINEMATICS IN ONE DIMENSION

Concourse Math Spring 2012 Worked Examples: Matrix Methods for Solving Systems of 1st Order Linear Differential Equations

Chapter 4. Truncation Errors

2001 November 15 Exam III Physics 191

STAD57 Time Series Analysis. Lecture 14

Seminar 4: Hotelling 2

Module 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II

Testing What You Know Now

Ordinary Differential Equations

GMM - Generalized Method of Moments

Non-parametric techniques. Instance Based Learning. NN Decision Boundaries. Nearest Neighbor Algorithm. Distance metric important

Retrieval Models. Boolean and Vector Space Retrieval Models. Common Preprocessing Steps. Boolean Model. Boolean Retrieval Model

!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)

Solutions from Chapter 9.1 and 9.2

Math 23 Spring Differential Equations. Final Exam Due Date: Tuesday, June 6, 5pm

Presentation Overview

Introduction to Physical Oceanography Homework 5 - Solutions

Probabilistic Robotics

Open loop vs Closed Loop. Example: Open Loop. Example: Feedforward Control. Advanced Control I

ME 391 Mechanical Engineering Analysis

KEY. Math 334 Midterm III Fall 2008 sections 001 and 003 Instructor: Scott Glasgow

15. Which Rule for Monetary Policy?

Appendix to Online l 1 -Dictionary Learning with Application to Novel Document Detection

T L. t=1. Proof of Lemma 1. Using the marginal cost accounting in Equation(4) and standard arguments. t )+Π RB. t )+K 1(Q RB

= ( ) ) or a system of differential equations with continuous parametrization (T = R

Announcements: Warm-up Exercise:

On-line Adaptive Optimal Timing Control of Switched Systems

Robotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.

Roller-Coaster Coordinate System

Lecture 33: November 29

Introduction D P. r = constant discount rate, g = Gordon Model (1962): constant dividend growth rate.

Transcription:

CS376 Compuer Vision Lecure 6: Opical Flow Qiing Huang Feb. 11 h 2019 Slides Credi: Krisen Grauman and Sebasian Thrun, Michael Black, Marc Pollefeys

Opical Flow mage racking 3D compuaion mage sequence (single camera) Tracked sequence 3D srucure + 3D rajecory 2

Wha is Opical Flow? p 2 p 3 v 2 v 3 p 1 Opical Flow v 1 p 4 v 4 ( i ),{ p } ( +1) Velociy vecors { v i } Opical flow is he 2D projecion of he physical movemen of poins relaive o he observer A common assumpion is brighness consancy: ( p, ) = ( p + v, + 1) i i i 3

When does Brighness Assumpion Break down? TV is based on illusory moion he se is saionary ye hings seem o move A uniform roaing sphere nohing seems o move, ye i is roaing Changing direcions or inensiies of lighing can make hings seem o move for eample, if he specular highligh on a roaing sphere moves Muscle movemen can make some spos on a cheeah move opposie direcion of moion

Opical Flow Assumpions: Brighness Consancy 5 * Slide from Michael Black, CS143 2003

Opical Flow Assumpions: Neighboring piels end o have similar moions When does his break down? 6 * Slide from Michael Black, CS143 2003

Opical Flow Assumpions: The image moion of a surface pah changes gradually over ime 7 * Slide from Michael Black, CS143 2003

1D Opical Flow

9 Opical Flow: 1D Case Brighness Consancy Assumpion: ) ), ( ( ) ), ( ( ) ( d d f + + = 0 ) ( = + v v = { 0 ) ( = f Because no change in brighness wih ime

2D Opical Flow

11 From 1D o 2D racking 0 ) ( = + 1D: 0 ) ( = + + y y 2D: 0 ) ( = + + v y u One equaion bu wo velociy (u,v) unknowns

How does his show up visually? Known as he Aperure Problem 12

Aperure Problem in Real Life 13

From 1D o 2D racking y (, y, + 1) (, ) (, + 1) v v 1 v 2 v 3 v 4 (, y, ) The Mah is very similar: v Aperure problem v G = b = Window size here ~ 1111 G 1 window around p window around p b 2 y y y 2 y 14

More Deail: Solving he aperure problem How o ge more equaions for a piel? -- impose addiional consrains mos common is o assume ha he flow field is smooh locally one mehod: preend he piel s neighbors have he same (u,v) Suppose a 55 window 15 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Lukas-Kanade flow Prob: we have more equaions han unknowns Soluion: solve leas squares problem minimum leas squares soluion given by soluion (in d) of: The summaions are over all piels in he K K window This echnique was firs proposed by Lukas & Kanade (1981) described in Trucco & Verri reading 16 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Condiions for solvabiliy Opimal (u, v) saisfies Lucas-Kanade equaion When is This Solvable? A T A should be inverible A T A should no be oo small due o noise eigenvalues l 1 and l 2 of A T A should no be oo small A T A should be well-condiioned l 1 / l 2 should no be oo large (l 1 = larger eigenvalue) 17 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Eigenvecors of A T A Suppose (,y) is on an edge. Wha is A T A? gradiens along edge all poin he same direcion gradiens away from edge have small magniude is an eigenvecor wih eigenvalue Wha s he oher eigenvecor of A T A? le N be perpendicular o N is he second eigenvecor wih eigenvalue 0 The eigenvecors of A T A relae o edge direcion and magniude 18 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Edge large gradiens, all he same large l 1, small l 2 19 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Low eure region gradiens have small magniude small l 1, small l 2 20 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

High eured region gradiens are differen, large magniudes large l 1, large l 2 21 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Observaion This is a wo image problem BUT Can measure sensiiviy by jus looking a one of he images! This ells us which piels are easy o rack, which are hard very useful laer on when we do feaure racking... 22 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Errors in Lukas-Kanade Wha are he poenial causes of errors in his procedure? Suppose A T A is easily inverible Suppose here is no much noise in he image When our assumpions are violaed Brighness consancy is no saisfied The moion is no small A poin does no move like is neighbors window size is oo large wha is he ideal window size? 23 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

mproving accuracy Recall our small moion assumpion -1 (,y) -1 (,y) This is no eac To do beer, we need o add higher order erms back in: -1 (,y) This is a polynomial roo finding problem Can solve using Newon s mehod Also known as Newon-Raphson mehod Lukas-Kanade mehod does one ieraion of Newon s mehod Beer resuls are obained via more ieraions 24 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

eraive Refinemen eraive Lukas-Kanade Algorihm 1. Esimae velociy a each piel by solving Lucas- Kanade equaions 2. Warp (-1) owards () using he esimaed flow field - use image warping echniques 3. Repea unil convergence 25 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Revisiing he small moion assumpion s his moion small enough? Probably no i s much larger han one piel (2 nd order erms dominae) How migh we solve his problem? 26 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Reduce he resoluion! 27 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Opical Flow Resuls 28 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

Opical Flow Resuls 29 * From Khurram Hassan-Shafique CAP5415 Compuer Vision 2003

30

31 * From Marc Pollefeys COMP 256 2003

Affine Flow 32 * Slide from Michael Black, CS143 2003

Horn & Schunck algorihm Addiional smoohness consrain : 2 2 2 es (( u + uy ) + ( v + = v 2 y )) ddy, besides Op. Flow consrain equaion erm ec ( u + yv + = ) 2 ddy, minimize es+aec 33 * From Marc Pollefeys COMP 256 2003

Horn & Schunck algorihm n simpler erms: f we wan dense flow, we need o regularize wha happens in ill condiioned (rank deficien) areas of he image. We ake he old cos funcion: d = arg min d ( (, ) ( + d, + 1) ) N And add a regularizaion erm o he cos: 2 d = arg min d N 2 ( (, ) ( + d, + 1) ) + a d Conve Program! We will see a lo of such formulaions in in robus regression! 34

Discussion: Wha are he oher mehods o improve opical flows?