Estimation of complex motions in video signals. Erhardt Barth. Institute for Neuro- and Bioinformatics University of Lübeck, Germany

Similar documents
Theory of! Partial Differential Equations!

Theory of! Partial Differential Equations-I!

Lecture #6: Continuous-Time Signals

Math 221: Mathematical Notation

Rotation, Scale and Translation Resilient Public Watermarking for Images

Tracking. Announcements

Solutions for homework 12

The motions of the celt on a horizontal plane with viscous friction

3, so θ = arccos

Chapter 1 Fundamental Concepts

Some Basic Information about M-S-D Systems

EELE Lecture 8 Example of Fourier Series for a Triangle from the Fourier Transform. Homework password is: 14445

Scientific Research of the Institute of Mathematics and Computer Science DIFFERENT VARIANTS OF THE BOUNDARY ELEMENT METHOD FOR PARABOLIC EQUATIONS

CS376 Computer Vision Lecture 6: Optical Flow

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

Mean Square Projection Error Gradient-based Variable Forgetting Factor FAPI

UNIVERSITY OF TRENTO MEASUREMENTS OF TRANSIENT PHENOMENA WITH DIGITAL OSCILLOSCOPES. Antonio Moschitta, Fabrizio Stefani, Dario Petri.

Announcements: Warm-up Exercise:

dt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.

Predator - Prey Model Trajectories and the nonlinear conservation law

Kinematics and kinematic functions

Traveling Waves. Chapter Introduction

Lecture 9: Advanced DFT concepts: The Exchange-correlation functional and time-dependent DFT

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...

Wall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.

A Limit Symmetry of Modified KdV Equation and Its Applications

Tutorial Sheet #2 discrete vs. continuous functions, periodicity, sampling

SOUND AND VIBRATION SIGNAL ANALYSIS USING IMPROVED SHORT-TIME FOURIER REPRESENTATION. June-Yule Lee

Rapid Termination Evaluation for Recursive Subdivision of Bezier Curves

On the Separation Theorem of Stochastic Systems in the Case Of Continuous Observation Channels with Memory

Fractional Method of Characteristics for Fractional Partial Differential Equations

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

Optical Flow I. Guido Gerig CS 6320, Spring 2015

di Bernardo, M. (1995). A purely adaptive controller to synchronize and control chaotic systems.

Let ( α, β be the eigenvector associated with the eigenvalue λ i

SUFFICIENT CONDITIONS FOR EXISTENCE SOLUTION OF LINEAR TWO-POINT BOUNDARY PROBLEM IN MINIMIZATION OF QUADRATIC FUNCTIONAL

t + t sin t t cos t sin t. t cos t sin t dt t 2 = exp 2 log t log(t cos t sin t) = Multiplying by this factor and then integrating, we conclude that

Chapter 6. Systems of First Order Linear Differential Equations

From Particles to Rigid Bodies

Extracting a Fluid Dynamic Texture and the Background from Video

Institute for Mathematical Methods in Economics. University of Technology Vienna. Singapore, May Manfred Deistler

Dirac Particles Emission from Reissner-Nordstrom-Vaidya Black Hole

4.5 Constant Acceleration

Single and Double Pendulum Models

ME 391 Mechanical Engineering Analysis

Chapter 3 Boundary Value Problem

Brock University Physics 1P21/1P91 Fall 2013 Dr. D Agostino. Solutions for Tutorial 3: Chapter 2, Motion in One Dimension

Improved Video Coding Using Long-term Global Motion Compensation

Augmented Reality II - Kalman Filters - Gudrun Klinker May 25, 2004

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

Trajectory planning in Cartesian space

Plasma Astrophysics Chapter 3: Kinetic Theory. Yosuke Mizuno Institute of Astronomy National Tsing-Hua University

Fitting Parameters in a Differential Equation With a Non-Analytical Solution

Lecture 1 Overview. course mechanics. outline & topics. what is a linear dynamical system? why study linear systems? some examples

Basilio Bona ROBOTICA 03CFIOR 1

Simulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010

Stability and Bifurcation in a Neural Network Model with Two Delays

Welcome Back to Physics 215!

Optimal trajectory generation for wheeled mobile robot

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

k The function Ψ(x) is called wavelet function and shows band-pass behavior. The wavelet coefficients d a,b

1. Kinematics I: Position and Velocity

Bifurcation Analysis of a Stage-Structured Prey-Predator System with Discrete and Continuous Delays

Section 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients

EECS 2602 Winter Laboratory 3 Fourier series, Fourier transform and Bode Plots in MATLAB

A Note on Fractional Electrodynamics. Abstract

Using the Kalman filter Extended Kalman filter

The physics most important to 20th Century. The most useful interactions - Electricity & PHYSICS. » Physics on the atomic scale.

6.003: Signal Processing

Ordinary Differential Equations

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

Physics 20 Lesson 5 Graphical Analysis Acceleration

Bianchi Type-II, VIII & IX Universe Filled with Wet Dark Fluid in f ( R,

Human-Robot Cooperative Manipulation with Motion Estimation

11/16/99 (T.F. Weiss) Lecture #18: Continuous time periodic signals. Outline: Fourier series of periodic functions. Motivation:

System of Linear Differential Equations

Mean-square Stability Control for Networked Systems with Stochastic Time Delay

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

I. OBJECTIVE OF THE EXPERIMENT.

Ch.1. Group Work Units. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

STAD57 Time Series Analysis. Lecture 5

Two Coupled Oscillators / Normal Modes

ENGINEERING FOR RURAL DEVELOPMENT Jelgava, MECHANISM MOTION STUDIES WITH COLLISIONS AT SEVERAL POINTS

2.1: What is physics? Ch02: Motion along a straight line. 2.2: Motion. 2.3: Position, Displacement, Distance

Production Analysis Reserves Estimation in Unconventional Reservoirs New Rate-Time Relations

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

Physics for Scientists and Engineers I

Signals and Systems Prof. Brian L. Evans Dept. of Electrical and Computer Engineering The University of Texas at Austin

On Safe Navigation Support System using GPS Compass

Reduced-order Observer Design for Structure and Motion Estimation

Accurate RMS Calculations for Periodic Signals by. Trapezoidal Rule with the Least Data Amount

ln 2 1 ln y x c y C x

EMS SCM joint meeting. On stochastic partial differential equations of parabolic type

Lecture 10: Wave equation, solution by spherical means

CH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC

IMPROVED HYPERBOLIC FUNCTION METHOD AND EXACT SOLUTIONS FOR VARIABLE COEFFICIENT BENJAMIN-BONA-MAHONY-BURGERS EQUATION

Anti-Disturbance Control for Multiple Disturbances

Lecture 10 Estimating Nonlinear Regression Models

Multi-scale 2D acoustic full waveform inversion with high frequency impulsive source

Transcription:

Esimaion o complex moions in video signals Erhard Barh Insiue or Neuro- and Bioinormaics Universiy o Lübec, Germany Nonlinear analysis o mulidimensional signals: LOCOMOTOR LOcal adapive esimaion o COmplex MOTion and ORienaion paerns LOCOMOTOR Franur Heidelberg Im DFG Schwerpun 4: Mahemaical mehods or ime series analysis and digial image processing hp://www.mah.uni-bremen.de/zeem/dfg-schwerpun/ Projec Ba-76/7 Lübec

Inrinsic dimension in 3D i0d: consan in all direcions: id: consan in direcions: id: consan in one direcion: i3d: no consan direcion: (x, y,) = cons. (x, y,) = g(ξ) (x, y,) = g(ξ,ζ) ( x,y,) = g(ξ,ζ,τ) Inrinsic dimension id and i3d regions are he leas requen bu mos signiican. i0d id id i3d

The classical moion model (x,y,) image-sequence inensiy v = (v,v ) T moion vecor α (v) = v + v + = V = 0 x y r = ( v,v, ) Why complex moions? Basic problems wih naural image sequences: muliple moions, occlusions, noisy daa, eaures a dieren scales,... 3

Transparen moions Movie (op) and x plo (below). Moion model or ransparen moions Opical low or one moion (x,y,) image sequence, v = (v,v ) T moion vecor α (v) = v + v + = V = 0 x y wih he derivaive operaor: α (v) = v + v + x y Opical low or N moions: = g (x-v )+...+g N (x-v N ) α(v ) α(v ) α (v ) = 0 N 4

Mixed-moion parameers Example or wo moions u, v : (x,)=g (x-v )+g (x-u ) α(v) α(u) = u v + ( u v xx + u v + u v ) + ( u + v ) x yy xy + ( u + v ) We deine he mixed-moion parameers as: y + c = vu c = vu c = uv +uv xx yy xy c = u+v c = u+v c = x y Solving or he mixed-moion parameers Wih he mixed-moion parameers we obain: α(v) α (u) = c + c + c xx xx yy yy xy xy + c + c + c x x y y = V L= 0 For one moion we had V = 0 The above consrain can be used in a number o ways o derive he mixed moion parameers in V rom, e.g. by deining he generalized srucure ensor e.g. JV= N 0 and solving he sysem m i im i(m ) i J N = T h LL V (M, M,...,( ) M ) M ij, i =,...,m are he minors o J N 5

Separaion o he mixed-moion parameers We inerpre he moion vecors u and v as complex numbers (v = v +i v, u = u +i u ) and observe ha u v = A0 = cxx cyy + icxy u+ v= A = c + ic x y A 0 and A are homogeneous symmerical uncions o he coordinaes u and v and hereore by Viea s rule he coeicien o he complex polynomial Q(z) = (z u)(z v) = z Az + A0 ha has he complex roos u and v. In case o N moions Q(z) = z A z + + ( ) A N N N N 0 Fourier analysis o muliple ransparen moions Coninuous ime ( x, ) = g α( u) α( v) ( x u) + g = 0 F = G δ( u ω+ ω ) + G ( x v) δ( v ω+ ω ) Discree ime ( x) = g F ( x u) + g ( ω) = φ G ( ω) + φ ( x v) G ( ω) α (u ) = u + φ -ju ω = e 6

Bloc-maching or ransparen moions F ( ω) = φ G( ω) + φg ( ω) Layer eliminaion F a ( ) ( ) 0( ) ω a F ω + a F ω = 0 = φ +φ a = φ φ Bloc-maching consrain Bac o space domain x) ( x u) ( x v) + ( x u v) = 0 ( 0 Synheic resuls or bloc-maching SNR : 35 db Bloc size : 5 x 5 pixels 7

Fourier analysis o muliple occluded moions ( x, ) = χ( x u) g( x u) + χ( x u) g ( x v) g ( x) g ( x) occluding objec occluded objec χ characerisic uncion o he occluding objec χ = χ How does he Fourier ransorm o (x,) deviae rom he wo moion planes? α ( u ) α ( v ) ( x, )? Occluded moion α( u) α( v) ( x, ) = q( x, u, v, ) δ( B( x u)) q( x, u, v, ) = ( u v) N( x u) ( u v) g ( x v) Foreground moion Bacground moion Normal N o he border B(x)=0 The equaion is he same as or ransparen moions excep a he occluding boundary. q( x, u, v, ) is maximal i u - v is orhogonal o he boundary is null i u - v is angen o he boundary 8

Occluded moion in he Fourier domain F( ω, ω ) = A( ω) δ( u ω+ω ) + B( ω) δ( v ω+ω ) + C( ω, ω ) C( ω, ω ) d /( v ω+ω ) along he plane u ω+ω = c The disorion C has hyperbolic decay Log(+F) Log(+F) The hyperbolic decay o C has been irs recognized by Beauchemin a al. or he paricular case o sraigh line boundary Reerences C. Moa, I. Sue, T. Aach, and E. Barh: Spaio-emporal moion esimaion or ransparency and occlusions. Signal Processing: Image Communicaion. Elsevier Science, 005. C. Moa, I. Sue, T. Aach, and E. Barh: Divide-and-Conquer Sraegies or Esimaing Muliple Transparen Moions. In: Proceedings o he s Inernaional Worshop on Complex Moion, Schloss Reisensburg, Germany. Lecure Noes on Compuer Science, LNCS Vol. 347, 005. C. Moa, I. Sue, T. Aach, and E. Barh: Esimaion o muliple orienaions a corners and juncions. In: 6h Paern Recogniion Symposium (DAGM'04), Tübingen, 63-70, 004. C. Moa, T. Aach, I. Sue, and E. Barh: Esimaion o muliple orienaions in muli-dimensional signals. In: IEEE Inernaional Conerence on Image Processing, 665-668, 004. I. Sue, T. Aach, E. Barh, and C. Moa: Muliple-moion-esimaion by bloc-maching using MRF. Inernaional Journal o Compuer & Inormaion Science, Vol. 5, No., pp. 4-5, 004. I. Sue, T. Aach, E. Barh, and C. Moa: Muliple-moion esimaion by bloc-maching using Marov random ields. In: Inernaional Journal o Compuer & Inormaion Science, Vol. 5, No., pp. 4-5, June 004. E. Barh, I. Sue, T. Aach, and C. Moa: Spaio-Temporal Moion Esimaion or Transparency and Occlusions.In: Proceedings o IEEE Inernaional Conerence on Image Processing (ICIP 03), Barcelona, Spain, Sepember 4-7, 003, Vol. III, pp. 69-7. E. Barh, I. Sue, and C. Moa: Analysis o moion and curvaure in image sequences (invied paper). In: Proc. o he 5h IEEE Souhwes Symposium on Image Analysis and Inerpreaion, 06-0, 00. C. Moa, I. Sue, and E. Barh: Analyic Soluions or Muliple Moions. In: Proceedings o he IEEE Inernaional Conerence on Image Processing (ICIP 0), Thessalonii, Greece, Ocober 7-0, 00, Vol. II, pp. 97-90. 9