Lecture 6 Stereo Systems Multi-view geometry

Size: px
Start display at page:

Download "Lecture 6 Stereo Systems Multi-view geometry"

Transcription

1 Lectre 6 Stereo Systems Mlti-view geometry Professor Silvio Savarese Comptational Vision and Geometry Lab Silvio Savarese Lectre 6-29-Jan-18

2 Lectre 6 Stereo Systems Mlti-view geometry Stereo systems Rectification Correspondence problem Mlti-view geometry The SFM problem Affine SFM Reading: [AZ] Chapter: 9 Epip. Geom. and the Fndam. Matrix Transf. [AZ] Chapter: 18 N view comptational methods [FP] Chapters: 7 Stereopsis [FP] Chapters: 8 Strctre from Motion Silvio Savarese Lectre 6-29-Jan-18

3 Epipolar geometry P p p e e O 1 O 2 Epipolar Plane Epipoles e, e Baseline Epipolar Lines = intersections of baseline with image planes = projections of the other camera center

4 Epipolar Constraint P p p O 1 O 2 p T E p = 0 E = Essential Matrix (Longet-Higgins, 1981) E = [ T ] R

5 Essential matrix [ ] R T E = " $ E = $ $ $ 0 T z T y T z 0 T x T y T x 0 % ' ' R ' '

6 Epipolar Constraint P p p O 1 O 2 p T F p = 0 F F = Fndamental Matrix (Fageras and Long, 1992) [ ] -1 = K - T T R K

7 Parallel image planes P v e p p e y z x O 1 O 2! p $! Epipolar lines are horizontal p = p Epipoles go to infinity v p' = v-coordinates are eqal " 1 % " p' p v 1 $ %

8 Parallel image planes P v e p p e y z x O 1 O 2 Hint : K E=? 1 =K 2 = known R = I T = (T, 0, 0) x parallel to O 1 O 2

9 Essential matrix for parallel images [ ] R T E = " $ E = $ $ $ 0 T z T y T z 0 T x T y T x 0 % " ' $ ' R = $ ' $ ' T 0 T 0 % ' ' ' T = [ T 0 0 ] R = I

10 Parallel image planes P v v e l p p l e y z O 1 What are the directions of epipolar lines? x l = E p' " $ = $ $ T 0 T 0 %" ' $ ' $ ' $ ' v' 1 % " ' $ ' = $ ' $ 0 T T v' % ' ' ' horizontal!

11 Parallel image planes P v v e p p e y z O 1 x How are p and p related? p T E p' = 0

12 Parallel image planes P v v e p p e y z O 1 x How are p and p related? é0 0 0 ùæ ö æ 0 ö ( v ê ) T úç ç = Þ v ( v ) T Tv Tv ê - ú = Þ - = Þ = ê 0 T 0 úç 1 ç Tv ë ûè ø è ø v = v" T pep

13 Parallel image planes P v v e p p e y z O 1 x Rectification: making two images parallel Why it is sefl? Epipolar constraint v = v New views can be synthesized by linear interpolation

14 Rectification: making two images parallel H Cortesy figre S. Lazebnik

15 Application: view morphing S. M. Seitz and C. R. Dyer, Proc. SIGGRAPH 96, 1996, 21-30

16 Rectification P H 1 I 1 I s î 1 v I 2 î S v C 1 v î 2 H 2 C s C 2

17

18

19

20 From its reflection!

21

22 Deep view morphing D. Ji, J. Kwon, M. McFarland, S. Savarese, CVPR 2017

23 Deep view morphing D. Ji, J. Kwon, M. McFarland, S. Savarese, CVPR 2017

24 Deep view morphing D. Ji, J. Kwon, M. McFarland, S. Savarese, CVPR 2017

25 Why are parallel images sefl? Makes trianglation easy Makes the correspondence problem easier

26 Point trianglation P v v e p p z e y z! p = " p p v 1 O $ % x! p' = " p' p v 1 $ % B O disparity = p p' B f z [Eq. 1] Disparity is inversely proportional to depth z!

27 Compting depth P (0, p v, 1) z f p p' O B O disparity = p p' B f [Eq. 1] z Disparity is inversely proportional to depth z!

28 Disparity maps p p' B f z p p [Eq. 1] Stereo pair Disparity map / depth map

29 Why are parallel images sefl? Makes trianglation easy Makes the correspondence problem easier

30 Correspondence problem P p belongs to l = F p p belongs to l = F T p p p O O Given a point in 3D, discover corresponding observations in left and right images [also called binoclar fsion problem]

31 Correspondence problem P p p O O When images are rectified, this problem is mch easier!

32 Correspondence problem A Cooperative Model (Marr and Poggio, 1976) Correlation Methods (1970--) Mlti-Scale Edge Matching (Marr, Poggio and Grimson, ) [FP] Chapters: 7

33 Correlation Methods (1970--) p Where is p' is? image 1 Image 2! p = " v 1 $ %! p' = " ' v 1 $ %

34 Correlation Methods (1970--) p Where is p' is? image 1 Image 2 100! p = " v 1 $ %! p' = " ' v 1 $ %

35 Correlation Methods (1970--) p Where is p' is? image 1 Image What s the problem with this?

36 Window-based correlation v W p image 1 image 2 Pick p a window W arond Bild vector w p = (, v )

37 Window-based correlation W v image 1 image 2 ' = + d = 1 Example: W is a 3x3 window in red w is a 9x1 vector w = [100, 100, 100, 90, 100, 20, 150, 150, 145] T Pick p a window W arond p = (, v ) Bild vector w Slide the window W along v = v in image 2 and compte w () for each Compte the dot prodct w T w () for each and retain the max vale

38 Window-based correlation W v image 1 image 2 ' = + d = 1 Example: W is a 3x3 window in red w is a 9x1 vector w = [100, 100, 100, 90, 100, 20, 150, 150, 145] T What s the problem with this?

39 Changes of brightness/exposre Changes in the mean and the variance of intensity vales in corresponding windows!

40 Normalized cross-correlation W v image 1 image 2 ' = + d = +1 Find that maximizes: (w w) T ( w "() w ") (w w) ( w "() w ") [Eq. 2] w = mean vale within W w'() = located at bar in image 1 mean vale within W located at in image 2

41 Example Image 1 Image 2 v scanline NCC Credit slide S. Lazebnik

42 Effect of the window s size Window size = 3 Window size = 20 Smaller window - More detail - More noise Larger window - Smoother disparity maps - Less prone to noise Credit slide S. Lazebnik

43 Isses Fore shortening effect O O Occlsions O O

44 Isses To redce the effect of foreshortening and occlsions, it is desirable to have small B / z ratio! However, when B/z is small, small errors in measrements imply large error in estimating depth Small B / z ratio Large B / z ratio O O O O e e

45 Isses Homogeneos regions mismatch

46 Repetitive patterns Isses

47 Correspondence problem is difficlt! - Occlsions - Fore shortening - Baseline trade-off - Homogeneos regions - Repetitive patterns Apply non-local constraints to help enforce the correspondences

48 Non-local constraints Uniqeness For any point in one image, there shold be at most one matching point in the other image Ordering Corresponding points shold be in the same order in both views Smoothness Disparity is typically a smooth fnction of x (except in occlding bondaries)

49 Lectre 6 Stereo Systems Mlti-view geometry Stereo systems Rectification Correspondence problem Mlti-view geometry The SFM problem Affine SFM Silvio Savarese Lectre 5-29-Jan-18

50 Strctre from motion problem Cortesy of Oxford Visal Geometry Grop

51 Strctre from motion problem X j x 1j M m M 1 x mj x 2j M 2 Given m images of n fixed 3D points x ij = M i X j, i = 1,, m, j = 1,, n

52 Strctre from motion problem X j x 1j M m M 1 x mj x 2j M 2 From the mxn observations x ij, estimate: m projection matrices M i n 3D points X j motion strctre

53 Affine strctre from motion (simpler problem) World point X j Image 1 x ij Image i From the mxn observations x ij, estimate: m projection matrices M i (affine cameras) n 3D points X j

54 Perspective x = M X! = " m 1 m 2 m 3 $! X = % " m 1 X m 2 X m 3 X $ %! M = A b " v 1 $ = % ém ê ê m êë m ù ú ú úû x E = ( m 1 X m 3 X, m 2 X m 3 X )T Affine x = M X! = " x E = (m 1 X, m 2 X) T magnification m 1 m 2 m 3 $ X % =! " A b! = " $ X =! " A b [Eq. 3] m 1 X m 2 X 1 $ %! % % $ % % " x y z 1! M = " m 1 m 2 m 3 = AX E + b $ $ %! = A 2 x3 b 2 x1 " 0 1x3 1! = " X E = m 1 m ! " x y z $ % $ % $ %

55 Affine cameras Image 1 x 1j X j x i j Image i For the affine case (in Eclidean space) x ij = A i X j + b i [Eq. 4] 2x1 2x3 3x1 2x1

56 The Affine Strctre-from-Motion Problem Given m images of n fixed points X j we can write x ij = A i X j + b i for i = 1,,m and j = 1,,n N. of cameras N. of points Problem: estimate m matrices A i, m matrices b i and the n positions X j from the m n observations x ij. How many eqations and how many nknown? 2m n eqations in 8m + 3n - 9 nknowns

57 The Affine Strctre-from-Motion Problem Two approaches: - Algebraic approach (affine epipolar geometry; estimate F; cameras; points) - Factorization method

58 Next lectre Mltiple view geometry: Affine and Perspective strctre from Motion

Lecture 5. Epipolar Geometry. Professor Silvio Savarese Computational Vision and Geometry Lab. 21-Jan-15. Lecture 5 - Silvio Savarese

Lecture 5. Epipolar Geometry. Professor Silvio Savarese Computational Vision and Geometry Lab. 21-Jan-15. Lecture 5 - Silvio Savarese Lecture 5 Epipolar Geometry Professor Silvio Savarese Computational Vision and Geometry Lab Silvio Savarese Lecture 5-21-Jan-15 Lecture 5 Epipolar Geometry Why is stereo useful? Epipolar constraints Essential

More information

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018

Geometric Image Manipulation. Lecture #4 Wednesday, January 24, 2018 Geometric Image Maniplation Lectre 4 Wednesda, Janar 4, 08 Programming Assignment Image Maniplation: Contet To start with the obvios, an image is a D arra of piels Piel locations represent points on the

More information

CS4495/6495 Introduction to Computer Vision. 3D-L3 Fundamental matrix

CS4495/6495 Introduction to Computer Vision. 3D-L3 Fundamental matrix CS4495/6495 Introduction to Computer Vision 3D-L3 Fundamental matrix Weak calibration Main idea: Estimate epipolar geometry from a (redundant) set of point correspondences between two uncalibrated cameras

More information

CSE 252B: Computer Vision II

CSE 252B: Computer Vision II CSE 252B: Computer Vision II Lecturer: Serge Belongie Scribe: Hamed Masnadi Shirazi, Solmaz Alipour LECTURE 5 Relationships between the Homography and the Essential Matrix 5.1. Introduction In practice,

More information

F O R SOCI AL WORK RESE ARCH

F O R SOCI AL WORK RESE ARCH 7 TH EUROPE AN CONFERENCE F O R SOCI AL WORK RESE ARCH C h a l l e n g e s i n s o c i a l w o r k r e s e a r c h c o n f l i c t s, b a r r i e r s a n d p o s s i b i l i t i e s i n r e l a t i o n

More information

Simplified Identification Scheme for Structures on a Flexible Base

Simplified Identification Scheme for Structures on a Flexible Base Simplified Identification Scheme for Strctres on a Flexible Base L.M. Star California State University, Long Beach G. Mylonais University of Patras, Greece J.P. Stewart University of California, Los Angeles

More information

Modelling by Differential Equations from Properties of Phenomenon to its Investigation

Modelling by Differential Equations from Properties of Phenomenon to its Investigation Modelling by Differential Eqations from Properties of Phenomenon to its Investigation V. Kleiza and O. Prvinis Kanas University of Technology, Lithania Abstract The Panevezys camps of Kanas University

More information

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises

1. State-Space Linear Systems 2. Block Diagrams 3. Exercises LECTURE 1 State-Space Linear Sstems This lectre introdces state-space linear sstems, which are the main focs of this book. Contents 1. State-Space Linear Sstems 2. Block Diagrams 3. Exercises 1.1 State-Space

More information

M 408 K Fall 2005 Inverse Trig Functions Important Decimal Approximations and Useful Trig Identities Decimal Approximations: p

M 408 K Fall 2005 Inverse Trig Functions Important Decimal Approximations and Useful Trig Identities Decimal Approximations: p M 408 K Fall 005 Inverse Trig Fnctions Imortant Decimal Aroimations an Usefl Trig Ientities Decimal Aroimations: 0 0000 0 0 0000 054 0500 6 0577 ( æ ö ç ø è 4 0785 0707 ; 0866 047 4 000 57 6 55 094 8 44

More information

A fundamental inverse problem in geosciences

A fundamental inverse problem in geosciences A fndamental inverse problem in geosciences Predict the vales of a spatial random field (SRF) sing a set of observed vales of the same and/or other SRFs. y i L i ( ) + v, i,..., n i ( P)? L i () : linear

More information

probability of k samples out of J fall in R.

probability of k samples out of J fall in R. Nonparametric Techniques for Density Estimation (DHS Ch. 4) n Introduction n Estimation Procedure n Parzen Window Estimation n Parzen Window Example n K n -Nearest Neighbor Estimation Introduction Suppose

More information

Structure from Motion. rot

Structure from Motion. rot Motion and Flow II Strctre from Motion tr rot r t, w F tr ˆ z t r ( zˆ ( r ([ ] r) )) rot? Passive Navigation and Strctre from Motion Thesystem moves with a rigid motion with translational velocity

More information

Formal Methods for Deriving Element Equations

Formal Methods for Deriving Element Equations Formal Methods for Deriving Element Eqations And the importance of Shape Fnctions Formal Methods In previos lectres we obtained a bar element s stiffness eqations sing the Direct Method to obtain eact

More information

Homogeneous Liner Systems with Constant Coefficients

Homogeneous Liner Systems with Constant Coefficients Homogeneos Liner Systems with Constant Coefficients Jly, 06 The object of stdy in this section is where A is a d d constant matrix whose entries are real nmbers. As before, we will look to the exponential

More information

A Tutorial on Data Reduction. Principal Component Analysis Theoretical Discussion. By Shireen Elhabian and Aly Farag

A Tutorial on Data Reduction. Principal Component Analysis Theoretical Discussion. By Shireen Elhabian and Aly Farag A Tutorial on Data Reduction Principal Component Analysis Theoretical Discussion By Shireen Elhabian and Aly Farag University of Louisville, CVIP Lab November 2008 PCA PCA is A backbone of modern data

More information

! " # $! % & '! , ) ( + - (. ) ( ) * + / 0 1 2 3 0 / 4 5 / 6 0 ; 8 7 < = 7 > 8 7 8 9 : Œ Š ž P P h ˆ Š ˆ Œ ˆ Š ˆ Ž Ž Ý Ü Ý Ü Ý Ž Ý ê ç è ± ¹ ¼ ¹ ä ± ¹ w ç ¹ è ¼ è Œ ¹ ± ¹ è ¹ è ä ç w ¹ ã ¼ ¹ ä ¹ ¼ ¹ ±

More information

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski

Advanced topics in Finite Element Method 3D truss structures. Jerzy Podgórski Advanced topics in Finite Element Method 3D trss strctres Jerzy Podgórski Introdction Althogh 3D trss strctres have been arond for a long time, they have been sed very rarely ntil now. They are difficlt

More information

EPIPOLAR GEOMETRY WITH MANY DETAILS

EPIPOLAR GEOMETRY WITH MANY DETAILS EPIPOLAR GEOMERY WIH MANY DEAILS hank ou for the slides. he come mostl from the following source. Marc Pollefes U. of North Carolina hree questions: (i) Correspondence geometr: Given an image point in

More information

Study of the diffusion operator by the SPH method

Study of the diffusion operator by the SPH method IOSR Jornal of Mechanical and Civil Engineering (IOSR-JMCE) e-issn: 2278-684,p-ISSN: 2320-334X, Volme, Isse 5 Ver. I (Sep- Oct. 204), PP 96-0 Stdy of the diffsion operator by the SPH method Abdelabbar.Nait

More information

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation:

Math 263 Assignment #3 Solutions. 1. A function z = f(x, y) is called harmonic if it satisfies Laplace s equation: Math 263 Assignment #3 Soltions 1. A fnction z f(x, ) is called harmonic if it satisfies Laplace s eqation: 2 + 2 z 2 0 Determine whether or not the following are harmonic. (a) z x 2 + 2. We se the one-variable

More information

Upper Bounds on the Spanning Ratio of Constrained Theta-Graphs

Upper Bounds on the Spanning Ratio of Constrained Theta-Graphs Upper Bonds on the Spanning Ratio of Constrained Theta-Graphs Prosenjit Bose and André van Renssen School of Compter Science, Carleton University, Ottaa, Canada. jit@scs.carleton.ca, andre@cg.scs.carleton.ca

More information

Chapter 6. Inverse Circular Functions and Trigonometric Equations. Section 6.1 Inverse Circular Functions y = 0

Chapter 6. Inverse Circular Functions and Trigonometric Equations. Section 6.1 Inverse Circular Functions y = 0 Chapter Inverse Circlar Fnctions and Trigonometric Eqations Section. Inverse Circlar Fnctions. onetoone. range. cos... = tan.. Sketch the reflection of the graph of f across the line =. 7. (a) [, ] é ù

More information

Differentiating Functions & Expressions - Edexcel Past Exam Questions

Differentiating Functions & Expressions - Edexcel Past Exam Questions - Edecel Past Eam Questions. (a) Differentiate with respect to (i) sin + sec, (ii) { + ln ()}. 5-0 + 9 Given that y =, ¹, ( -) 8 (b) show that = ( -). (6) June 05 Q. f() = e ln, > 0. (a) Differentiate

More information

Multiple View Geometry in Computer Vision

Multiple View Geometry in Computer Vision Multiple View Geometry in Computer Vision Prasanna Sahoo Department of Mathematics University of Louisville 1 Scene Planes & Homographies Lecture 19 March 24, 2005 2 In our last lecture, we examined various

More information

FRTN10 Exercise 12. Synthesis by Convex Optimization

FRTN10 Exercise 12. Synthesis by Convex Optimization FRTN Exercise 2. 2. We want to design a controller C for the stable SISO process P as shown in Figre 2. sing the Yola parametrization and convex optimization. To do this, the control loop mst first be

More information

Vision 3D articielle Session 2: Essential and fundamental matrices, their computation, RANSAC algorithm

Vision 3D articielle Session 2: Essential and fundamental matrices, their computation, RANSAC algorithm Vision 3D articielle Session 2: Essential and fundamental matrices, their computation, RANSAC algorithm Pascal Monasse monasse@imagine.enpc.fr IMAGINE, École des Ponts ParisTech Contents Some useful rules

More information

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University

Lecture Notes: Finite Element Analysis, J.E. Akin, Rice University 9. TRUSS ANALYSIS... 1 9.1 PLANAR TRUSS... 1 9. SPACE TRUSS... 11 9.3 SUMMARY... 1 9.4 EXERCISES... 15 9. Trss analysis 9.1 Planar trss: The differential eqation for the eqilibrim of an elastic bar (above)

More information

PAijpam.eu EPIPOLAR GEOMETRY WITH A FUNDAMENTAL MATRIX IN CANONICAL FORM Georgi Hristov Georgiev 1, Vencislav Dakov Radulov 2

PAijpam.eu EPIPOLAR GEOMETRY WITH A FUNDAMENTAL MATRIX IN CANONICAL FORM Georgi Hristov Georgiev 1, Vencislav Dakov Radulov 2 International Journal of Pure and Applied Mathematics Volume 105 No. 4 2015, 669-683 ISSN: 1311-8080 (printed version); ISSN: 1314-3395 (on-line version) url: http://www.ijpam.eu doi: http://dx.doi.org/10.12732/ijpam.v105i4.8

More information

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction

PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS. 1. Introduction PHASE PLANE DIAGRAMS OF DIFFERENCE EQUATIONS TANYA DEWLAND, JEROME WESTON, AND RACHEL WEYRENS Abstract. We will be determining qalitatie featres of a discrete dynamical system of homogeneos difference

More information

MATH 174: Numerical Analysis I. Math Division, IMSP, UPLB 1 st Sem AY

MATH 174: Numerical Analysis I. Math Division, IMSP, UPLB 1 st Sem AY MATH 74: Numerical Analysis I Math Division, IMSP, UPLB st Sem AY 0809 Eample : Prepare a table or the unction e or in [0,]. The dierence between adjacent abscissas is h step size. What should be the step

More information

CHAPTER 4 PRINCIPAL COMPONENT ANALYSIS-BASED FUSION

CHAPTER 4 PRINCIPAL COMPONENT ANALYSIS-BASED FUSION 59 CHAPTER 4 PRINCIPAL COMPONENT ANALYSIS-BASED FUSION 4. INTRODUCTION Weighted average-based fusion algorithms are one of the widely used fusion methods for multi-sensor data integration. These methods

More information

Prandl established a universal velocity profile for flow parallel to the bed given by

Prandl established a universal velocity profile for flow parallel to the bed given by EM 0--00 (Part VI) (g) The nderlayers shold be at least three thicknesses of the W 50 stone, bt never less than 0.3 m (Ahrens 98b). The thickness can be calclated sing Eqation VI-5-9 with a coefficient

More information

1 JAXA Special Pblication JAXA-SP-1-E Small-scale trblence affects flow fields arond a blff body and therefore it governs characteristics of cross-sec

1 JAXA Special Pblication JAXA-SP-1-E Small-scale trblence affects flow fields arond a blff body and therefore it governs characteristics of cross-sec First International Symposim on Fltter and its Application, 1 11 IEXPERIMENTAL STUDY ON TURBULENCE PARTIAL SIMULATION FOR BLUFF BODY Hiroshi Katschi +1 and Hitoshi Yamada + +1 Yokohama National University,

More information

Edexcel GCE A Level Maths Further Maths 3 Matrices.

Edexcel GCE A Level Maths Further Maths 3 Matrices. Edexcel GCE A Level Maths Further Maths 3 Matrices. Edited by: K V Kumaran kumarmathsweebly.com kumarmathsweebly.com 2 kumarmathsweebly.com 3 kumarmathsweebly.com 4 kumarmathsweebly.com 5 kumarmathsweebly.com

More information

IE 400 Principles of Engineering Management. Graphical Solution of 2-variable LP Problems

IE 400 Principles of Engineering Management. Graphical Solution of 2-variable LP Problems IE 400 Principles of Engineering Management Graphical Solution of 2-variable LP Problems Graphical Solution of 2-variable LP Problems Ex 1.a) max x 1 + 3 x 2 s.t. x 1 + x 2 6 - x 1 + 2x 2 8 x 1, x 2 0,

More information

Hybrid modelling and model reduction for control & optimisation

Hybrid modelling and model reduction for control & optimisation Hybrid modelling and model redction for control & optimisation based on research done by RWTH-Aachen and TU Delft presented by Johan Grievink Models for control and optimiation market and environmental

More information

UNIQUE FJORDS AND THE ROYAL CAPITALS UNIQUE FJORDS & THE NORTH CAPE & UNIQUE NORTHERN CAPITALS

UNIQUE FJORDS AND THE ROYAL CAPITALS UNIQUE FJORDS & THE NORTH CAPE & UNIQUE NORTHERN CAPITALS Q J j,. Y j, q.. Q J & j,. & x x. Q x q. ø. 2019 :. q - j Q J & 11 Y j,.. j,, q j q. : 10 x. 3 x - 1..,,. 1-10 ( ). / 2-10. : 02-06.19-12.06.19 23.06.19-03.07.19 30.06.19-10.07.19 07.07.19-17.07.19 14.07.19-24.07.19

More information

Induced Planar Homologies in Epipolar Geometry

Induced Planar Homologies in Epipolar Geometry Global Journal of Pure and Applied Mathematics. ISSN 0973-1768 Volume 12, Number 4 (2016), pp. 3759 3773 Research India Publications http://www.ripublication.com/gjpam.htm Induced Planar Homologies in

More information

Mechanisms and topology determination of complex chemical and biological network systems from first-passage theoretical approach

Mechanisms and topology determination of complex chemical and biological network systems from first-passage theoretical approach Mechanisms and topology determination of complex chemical and biological network systems from first-passage theoretical approach Xin Li and Anatoly B. Kolomeisky Citation: J. Chem. Phys. 39, 4406 (203);

More information

The Dual of the Maximum Likelihood Method

The Dual of the Maximum Likelihood Method Department of Agricltral and Resorce Economics University of California, Davis The Dal of the Maximm Likelihood Method by Qirino Paris Working Paper No. 12-002 2012 Copyright @ 2012 by Qirino Paris All

More information

Lecture: Corporate Income Tax

Lecture: Corporate Income Tax Lectre: Corporate Income Tax Ltz Krschwitz & Andreas Löffler Disconted Cash Flow, Section 2.1, Otline 2.1 Unlevered firms Similar companies Notation 2.1.1 Valation eqation 2.1.2 Weak atoregressive cash

More information

Bayes and Naïve Bayes Classifiers CS434

Bayes and Naïve Bayes Classifiers CS434 Bayes and Naïve Bayes Classifiers CS434 In this lectre 1. Review some basic probability concepts 2. Introdce a sefl probabilistic rle - Bayes rle 3. Introdce the learning algorithm based on Bayes rle (ths

More information

Elements of Coordinate System Transformations

Elements of Coordinate System Transformations B Elements of Coordinate System Transformations Coordinate system transformation is a powerfl tool for solving many geometrical and kinematic problems that pertain to the design of gear ctting tools and

More information

Chapter 1: Differential Form of Basic Equations

Chapter 1: Differential Form of Basic Equations MEG 74 Energ and Variational Methods in Mechanics I Brendan J. O Toole, Ph.D. Associate Professor of Mechanical Engineering Howard R. Hghes College of Engineering Universit of Nevada Las Vegas TBE B- (7)

More information

), σ is a parameter, is the Euclidean norm in R d.

), σ is a parameter, is the Euclidean norm in R d. WEIGHTED GRAPH LAPLACIAN AND IMAGE INPAINTING ZUOQIANG SHI, STANLEY OSHER, AND WEI ZHU Abstract. Inspired by the graph Laplacian and the point integral method, we introdce a novel weighted graph Laplacian

More information

Research Article Uniqueness of Solutions to a Nonlinear Elliptic Hessian Equation

Research Article Uniqueness of Solutions to a Nonlinear Elliptic Hessian Equation Applied Mathematics Volme 2016, Article ID 4649150, 5 pages http://dx.doi.org/10.1155/2016/4649150 Research Article Uniqeness of Soltions to a Nonlinear Elliptic Hessian Eqation Siyan Li School of Mathematics

More information

arxiv:quant-ph/ v4 14 May 2003

arxiv:quant-ph/ v4 14 May 2003 Phase-transition-like Behavior of Qantm Games arxiv:qant-ph/0111138v4 14 May 2003 Jiangfeng D Department of Modern Physics, University of Science and Technology of China, Hefei, 230027, People s Repblic

More information

Model Predictive Control Lecture VIa: Impulse Response Models

Model Predictive Control Lecture VIa: Impulse Response Models Moel Preictive Control Lectre VIa: Implse Response Moels Niet S. Kaisare Department of Chemical Engineering Inian Institte of Technolog Maras Ingreients of Moel Preictive Control Dnamic Moel Ftre preictions

More information

Self-Testing Polynomial Functions Efficiently and over Rational Domains

Self-Testing Polynomial Functions Efficiently and over Rational Domains Chapter 1 Self-Testing Polynomial Functions Efficiently and over Rational Domains Ronitt Rubinfeld Madhu Sudan Ý Abstract In this paper we give the first self-testers and checkers for polynomials over

More information

Affine Invariant Total Variation Models

Affine Invariant Total Variation Models Affine Invariant Total Variation Models Helen Balinsky, Alexander Balinsky Media Technologies aboratory HP aboratories Bristol HP-7-94 Jne 6, 7* Total Variation, affine restoration, Sobolev ineqality,

More information

FUZZY BOUNDARY ELEMENT METHODS: A NEW MULTI-SCALE PERTURBATION APPROACH FOR SYSTEMS WITH FUZZY PARAMETERS

FUZZY BOUNDARY ELEMENT METHODS: A NEW MULTI-SCALE PERTURBATION APPROACH FOR SYSTEMS WITH FUZZY PARAMETERS MODELOWANIE INŻYNIERSKIE ISNN 896-77X 3, s. 433-438, Gliwice 6 FUZZY BOUNDARY ELEMENT METHODS: A NEW MULTI-SCALE PERTURBATION APPROACH FOR SYSTEMS WITH FUZZY PARAMETERS JERZY SKRZYPCZYK HALINA WITEK Zakład

More information

Feature extraction: Corners and blobs

Feature extraction: Corners and blobs Featre etraction: Corners and blobs Wh etract featres? Motiation: panorama stitching We hae two images how do we combine them? Wh etract featres? Motiation: panorama stitching We hae two images how do

More information

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by

Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introduction The transmission line equations are given by Reflections on a mismatched transmission line Reflections.doc (4/1/00) Introdction The transmission line eqations are given by, I z, t V z t l z t I z, t V z, t c z t (1) (2) Where, c is the per-nit-length

More information

Final Exam Due on Sunday 05/06

Final Exam Due on Sunday 05/06 Final Exam Due on Sunday 05/06 The exam should be completed individually without collaboration. However, you are permitted to consult with the textbooks, notes, slides and even internet resources. If you

More information

Lecture: Corporate Income Tax - Unlevered firms

Lecture: Corporate Income Tax - Unlevered firms Lectre: Corporate Income Tax - Unlevered firms Ltz Krschwitz & Andreas Löffler Disconted Cash Flow, Section 2.1, Otline 2.1 Unlevered firms Similar companies Notation 2.1.1 Valation eqation 2.1.2 Weak

More information

Outline. Linear Algebra for Computer Vision

Outline. Linear Algebra for Computer Vision Outline Linear Algebra for Computer Vision Introduction CMSC 88 D Notation and Basics Motivation Linear systems of equations Gauss Elimination, LU decomposition Linear Spaces and Operators Addition, scalar

More information

Gradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How

Gradient Projection Anti-windup Scheme on Constrained Planar LTI Systems. Justin Teo and Jonathan P. How 1 Gradient Projection Anti-windp Scheme on Constrained Planar LTI Systems Jstin Teo and Jonathan P. How Technical Report ACL1 1 Aerospace Controls Laboratory Department of Aeronatics and Astronatics Massachsetts

More information

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator

Nonparametric Identification and Robust H Controller Synthesis for a Rotational/Translational Actuator Proceedings of the 6 IEEE International Conference on Control Applications Mnich, Germany, October 4-6, 6 WeB16 Nonparametric Identification and Robst H Controller Synthesis for a Rotational/Translational

More information

Discontinuous Fluctuation Distribution for Time-Dependent Problems

Discontinuous Fluctuation Distribution for Time-Dependent Problems Discontinos Flctation Distribtion for Time-Dependent Problems Matthew Hbbard School of Compting, University of Leeds, Leeds, LS2 9JT, UK meh@comp.leeds.ac.k Introdction For some years now, the flctation

More information

Convex Hypersurfaces of Constant Curvature in Hyperbolic Space

Convex Hypersurfaces of Constant Curvature in Hyperbolic Space Srvey in Geometric Analysis and Relativity ALM 0, pp. 41 57 c Higher Edcation Press and International Press Beijing-Boston Convex Hypersrfaces of Constant Crvatre in Hyperbolic Space Bo Gan and Joel Sprck

More information

Thermal balance of a wall with PCM-enhanced thermal insulation

Thermal balance of a wall with PCM-enhanced thermal insulation Thermal balance of a wall with PCM-enhanced thermal inslation E. Kossecka Institte of Fndamental Technological esearch of the Polish Academy of Sciences, Warsaw, Poland J. Kośny Oak idge National aboratory;

More information

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018

Lecture 3. (2) Last time: 3D space. The dot product. Dan Nichols January 30, 2018 Lectre 3 The dot prodct Dan Nichols nichols@math.mass.ed MATH 33, Spring 018 Uniersity of Massachsetts Janary 30, 018 () Last time: 3D space Right-hand rle, the three coordinate planes 3D coordinate system:

More information

M3: Multiple View Geometry

M3: Multiple View Geometry M3: Multiple View Geometry L18: Projective Structure from Motion: Iterative Algorithm based on Factorization Based on Sections 13.4 C. V. Jawahar jawahar-at-iiit.net Mar 2005: 1 Review: Reconstruction

More information

Trinocular Geometry Revisited

Trinocular Geometry Revisited Trinocular Geometry Revisited Jean Pounce and Martin Hebert 报告人 : 王浩人 2014-06-24 Contents 1. Introduction 2. Converging Triplets of Lines 3. Converging Triplets of Visual Rays 4. Discussion 1. Introduction

More information

MATH140 Exam 2 - Sample Test 1 Detailed Solutions

MATH140 Exam 2 - Sample Test 1 Detailed Solutions www.liontutors.com 1. D. reate a first derivative number line MATH140 Eam - Sample Test 1 Detailed Solutions cos -1 0 cos -1 cos 1 cos 1/ p + æp ö p æp ö ç è 4 ø ç è ø.. reate a second derivative number

More information

CSE 252B: Computer Vision II

CSE 252B: Computer Vision II CSE 252B: Computer Vision II Lecturer: Serge Belongie Scribe: Tasha Vanesian LECTURE 3 Calibrated 3D Reconstruction 3.1. Geometric View of Epipolar Constraint We are trying to solve the following problem:

More information

CHAPTER 3 Describing Relationships

CHAPTER 3 Describing Relationships CHAPTER 3 Describing Relationships 3.1 Scatterplots and Correlation The Practice of Statistics, 5th Edition Starnes, Tabor, Yates, Moore Bedford Freeman Worth Publishers Scatterplots and Correlation Learning

More information

TrustSVD: Collaborative Filtering with Both the Explicit and Implicit Influence of User Trust and of Item Ratings

TrustSVD: Collaborative Filtering with Both the Explicit and Implicit Influence of User Trust and of Item Ratings Proceedings of the Twenty-Ninth AAAI Conference on Artificial Intelligence TrstSVD: Collaborative Filtering with Both the Explicit and Implicit Inflence of User Trst and of Item Ratings Gibing Go Jie Zhang

More information

MODELING FREQUENCY SELECTIVE SURFACES WITH MULTIPLE METALLIZATION AND CONTINUOUSLY LAYERED DIELECTRIC SUBSTRATES

MODELING FREQUENCY SELECTIVE SURFACES WITH MULTIPLE METALLIZATION AND CONTINUOUSLY LAYERED DIELECTRIC SUBSTRATES ' # MODELING FREQUENCY SELECTIVE SURFACES WITH MULTIPLE METALLIZATION AND CONTINUOUSLY LAYERED DIELECTRIC SUBSTRATES B.J. Morsink, S.W. Rienstra and A.G. Tijhis (1) Thales Naal Nederland, Bsiness Unit

More information

CDS 110b: Lecture 1-2 Introduction to Optimal Control

CDS 110b: Lecture 1-2 Introduction to Optimal Control CDS 110b: Lectre 1-2 Introdction to Optimal Control Richard M. Mrray 4 Janary 2006 Goals: Introdce the problem of optimal control as method of trajectory generation State the maimm principle and give eamples

More information

CMSC427 Parametric curves: Hermite, Catmull-Rom, Bezier

CMSC427 Parametric curves: Hermite, Catmull-Rom, Bezier CMSC427 Parametric curves: Hermite, Catmull-Rom, Bezier Modeling Creating 3D objects How to construct complicated surfaces? Goal Specify objects with few control points Resulting object should be visually

More information

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty

Technical Note. ODiSI-B Sensor Strain Gage Factor Uncertainty Technical Note EN-FY160 Revision November 30, 016 ODiSI-B Sensor Strain Gage Factor Uncertainty Abstract Lna has pdated or strain sensor calibration tool to spport NIST-traceable measrements, to compte

More information

PCA FACE RECOGNITION

PCA FACE RECOGNITION PCA FACE RECOGNITION The slides are from several sources through James Hays (Brown); Srinivasa Narasimhan (CMU); Silvio Savarese (U. of Michigan); Shree Nayar (Columbia) including their own slides. Goal

More information

Decision making is the process of selecting

Decision making is the process of selecting Jornal of Advances in Compter Engineering and Technology, (4) 06 A New Mlti-Criteria Decision Making Based on Fzzy- Topsis Theory Leila Yahyaie Dizaji, Sohrab khanmohammadi Received (05-09-) Accepted (06--)

More information

An Investigation into Estimating Type B Degrees of Freedom

An Investigation into Estimating Type B Degrees of Freedom An Investigation into Estimating Type B Degrees of H. Castrp President, Integrated Sciences Grop Jne, 00 Backgrond The degrees of freedom associated with an ncertainty estimate qantifies the amont of information

More information

LINEAR COMBINATIONS AND SUBSPACES

LINEAR COMBINATIONS AND SUBSPACES CS131 Part II, Linear Algebra and Matrices CS131 Mathematics for Compter Scientists II Note 5 LINEAR COMBINATIONS AND SUBSPACES Linear combinations. In R 2 the vector (5, 3) can be written in the form

More information

The Multibody Trifocal Tensor: Motion Segmentation from 3 Perspective Views

The Multibody Trifocal Tensor: Motion Segmentation from 3 Perspective Views The Multibody Trifocal Tensor: Motion Segmentation from 3 Perspective Views Richard Hartley 1,2 and RenéVidal 2,3 1 Dept. of Systems Engineering 3 Center for Imaging Science Australian National University

More information

Curves - Foundation of Free-form Surfaces

Curves - Foundation of Free-form Surfaces Crves - Fondation of Free-form Srfaces Why Not Simply Use a Point Matrix to Represent a Crve? Storage isse and limited resoltion Comptation and transformation Difficlties in calclating the intersections

More information

Affine Structure from Motion

Affine Structure from Motion Affie Structure from Motio EECS 598-8 Fall 24! Foudatios of Computer Visio!! Istructor: Jaso Corso (jjcorso)! web.eecs.umich.edu/~jjcorso/t/598f4!! Readigs: FP 8.2! Date: /5/4!! Materials o these slides

More information

DEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS (H/U < 20) AND CONSEQUENCES ON CRITICALITY SAFETY

DEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS (H/U < 20) AND CONSEQUENCES ON CRITICALITY SAFETY DEFINITION OF A NEW UO 2 F 2 DENSITY LAW FOR LOW- MODERATED SOLUTIONS ( < 20) AND CONSEQUENCES ON CRITICALITY SAFETY N. Leclaire, S. Evo, I.R.S.N., France Introdction In criticality stdies, the blk density

More information

Approximate Solution for the System of Non-linear Volterra Integral Equations of the Second Kind by using Block-by-block Method

Approximate Solution for the System of Non-linear Volterra Integral Equations of the Second Kind by using Block-by-block Method Astralian Jornal of Basic and Applied Sciences, (1): 114-14, 008 ISSN 1991-8178 Approximate Soltion for the System of Non-linear Volterra Integral Eqations of the Second Kind by sing Block-by-block Method

More information

Reducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation

Reducing Conservatism in Flutterometer Predictions Using Volterra Modeling with Modal Parameter Estimation JOURNAL OF AIRCRAFT Vol. 42, No. 4, Jly Agst 2005 Redcing Conservatism in Fltterometer Predictions Using Volterra Modeling with Modal Parameter Estimation Rick Lind and Joao Pedro Mortaga University of

More information

Applications of Light-Front Dynamics in Hadron Physics

Applications of Light-Front Dynamics in Hadron Physics Applications of Light-Front Dynamics in Hadron Physics Chueng-Ryong Ji North Carolina State University 1. What is light-front dynamics (LFD)? 2. Why is LFD useful in hadron physics? 3. Any first principle

More information

Citation Osaka Journal of Mathematics. 43(2)

Citation Osaka Journal of Mathematics. 43(2) TitleIrreducible representations of the Author(s) Kosuda, Masashi Citation Osaka Journal of Mathematics. 43(2) Issue 2006-06 Date Text Version publisher URL http://hdl.handle.net/094/0396 DOI Rights Osaka

More information

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum

The Real Stabilizability Radius of the Multi-Link Inverted Pendulum Proceedings of the 26 American Control Conference Minneapolis, Minnesota, USA, Jne 14-16, 26 WeC123 The Real Stabilizability Radis of the Mlti-Link Inerted Pendlm Simon Lam and Edward J Daison Abstract

More information

Chapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of

Chapter 2. Linear Algebra. rather simple and learning them will eventually allow us to explain the strange results of Chapter 2 Linear Algebra In this chapter, we study the formal structure that provides the background for quantum mechanics. The basic ideas of the mathematical machinery, linear algebra, are rather simple

More information

Control Systems Design

Control Systems Design ELEC4410 Control Systems Design Lectre 16: Controllability and Observability Canonical Decompositions Jlio H. Braslavsky jlio@ee.newcastle.ed.a School of Electrical Engineering and Compter Science Lectre

More information

Matrix-Vector Products and the Matrix Equation Ax = b

Matrix-Vector Products and the Matrix Equation Ax = b Matrix-Vector Products and the Matrix Equation Ax = b A. Havens Department of Mathematics University of Massachusetts, Amherst January 31, 2018 Outline 1 Matrices Acting on Vectors Linear Combinations

More information

Examination paper for TFY4240 Electromagnetic theory

Examination paper for TFY4240 Electromagnetic theory Department of Physics Examination paper for TFY4240 Electromagnetic theory Academic contact during examination: Associate Professor John Ove Fjærestad Phone: 97 94 00 36 Examination date: 16 December 2015

More information

FINITE ELEMENT MODELING OF EDDY CURRENT PROBES FOR EDGE EFFECT

FINITE ELEMENT MODELING OF EDDY CURRENT PROBES FOR EDGE EFFECT FIITE ELEMET MODELIG OF EDDY CURRET PROBES FOR EDGE EFFECT REDUCTIO Sarit Sharma, Ibrahim Elshafiey, Lalita Udpa, and Satish Udpa Department of Electrical and Compter Engineering Iowa State University

More information

Section 7.4: Integration of Rational Functions by Partial Fractions

Section 7.4: Integration of Rational Functions by Partial Fractions Section 7.4: Integration of Rational Fnctions by Partial Fractions This is abot as complicated as it gets. The Method of Partial Fractions Ecept for a few very special cases, crrently we have no way to

More information

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers

Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers Theoretical and Experimental Implementation of DC Motor Nonlinear Controllers D.R. Espinoza-Trejo and D.U. Campos-Delgado Facltad de Ingeniería, CIEP, UASLP, espinoza trejo dr@aslp.mx Facltad de Ciencias,

More information

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation

A Survey of the Implementation of Numerical Schemes for Linear Advection Equation Advances in Pre Mathematics, 4, 4, 467-479 Pblished Online Agst 4 in SciRes. http://www.scirp.org/jornal/apm http://dx.doi.org/.436/apm.4.485 A Srvey of the Implementation of Nmerical Schemes for Linear

More information

A Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations

A Computational Study with Finite Element Method and Finite Difference Method for 2D Elliptic Partial Differential Equations Applied Mathematics, 05, 6, 04-4 Pblished Online November 05 in SciRes. http://www.scirp.org/jornal/am http://d.doi.org/0.46/am.05.685 A Comptational Stdy with Finite Element Method and Finite Difference

More information

When Closed Graph Manifolds are Finitely Covered by Surface Bundles Over S 1

When Closed Graph Manifolds are Finitely Covered by Surface Bundles Over S 1 Acta Mathematica Sinica, English Series 1999, Jan, Vol15, No1, p 11 20 When Closed Graph Manifolds are Finitely Covered by Srface Bndles Over S 1 Yan Wang Fengchn Y Department of Mathematics, Peking University,

More information

Move Blocking Strategies in Receding Horizon Control

Move Blocking Strategies in Receding Horizon Control Move Blocking Strategies in Receding Horizon Control Raphael Cagienard, Pascal Grieder, Eric C. Kerrigan and Manfred Morari Abstract In order to deal with the comptational brden of optimal control, it

More information

Experimental Study of an Impinging Round Jet

Experimental Study of an Impinging Round Jet Marie Crie ay Final Report : Experimental dy of an Impinging Rond Jet BOURDETTE Vincent Ph.D stdent at the Rovira i Virgili University (URV), Mechanical Engineering Department. Work carried ot dring a

More information

Segmentation of Subspace Arrangements III Robust GPCA

Segmentation of Subspace Arrangements III Robust GPCA Segmentation of Subspace Arrangements III Robust GPCA Berkeley CS 294-6, Lecture 25 Dec. 3, 2006 Generalized Principal Component Analysis (GPCA): (an overview) x V 1 V 2 (x 3 = 0)or(x 1 = x 2 = 0) {x 1x

More information

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration

A New Method for Calculating of Electric Fields Around or Inside Any Arbitrary Shape Electrode Configuration Proceedings of the 5th WSEAS Int. Conf. on Power Systems and Electromagnetic Compatibility, Corf, Greece, Agst 3-5, 005 (pp43-48) A New Method for Calclating of Electric Fields Arond or Inside Any Arbitrary

More information

Linear and Nonlinear Model Predictive Control of Quadruple Tank Process

Linear and Nonlinear Model Predictive Control of Quadruple Tank Process Linear and Nonlinear Model Predictive Control of Qadrple Tank Process P.Srinivasarao Research scholar Dr.M.G.R.University Chennai, India P.Sbbaiah, PhD. Prof of Dhanalaxmi college of Engineering Thambaram

More information