ECE 472/572 - Digital Image Processing. Roadmap. Questions. Lecture 6 Geometric and Radiometric Transformation 09/27/11

Similar documents
Report on Image warping

Lecture Topics VMSC Prof. Dr.-Ing. habil. Hermann Lödding Prof. Dr.-Ing. Wolfgang Hintze. PD Dr.-Ing. habil.

Fourier Transform. Additive noise. Fourier Tansform. I = S + N. Noise doesn t depend on signal. We ll consider:

Tutorial 2. COMP4134 Biometrics Authentication. February 9, Jun Xu, Teaching Asistant

The Fourier Transform

Bruce A. Draper & J. Ross Beveridge, January 25, Geometric Image Manipulation. Lecture #1 January 25, 2013

1 Matrix representations of canonical matrices

CSci 6974 and ECSE 6966 Math. Tech. for Vision, Graphics and Robotics Lecture 21, April 17, 2006 Estimating A Plane Homography

Autonomous On-orbit Calibration Approaches for Star Trackers. D. Todd Griffith

Change Detection: Current State of the Art and Future Directions

Transform Coding. Transform Coding Principle

CS4495/6495 Introduction to Computer Vision. 3C-L3 Calibrating cameras

Invariant deformation parameters from GPS permanent networks using stochastic interpolation

Structure from Motion. Forsyth&Ponce: Chap. 12 and 13 Szeliski: Chap. 7

ADAPTIVE IMAGE FILTERING

Formulas for the Determinant

THE ASTER IMAGES FOR THE ENVIRONMENTAL MONITORING

Review: Fit a line to N data points

CS 468 Lecture 16: Isometry Invariance and Spectral Techniques

Lectures - Week 4 Matrix norms, Conditioning, Vector Spaces, Linear Independence, Spanning sets and Basis, Null space and Range of a Matrix

Image Processing for Bubble Detection in Microfluidics

Video Data Analysis. Video Data Analysis, B-IT. Lecture plan:

Integrals and Invariants of Euler-Lagrange Equations

AERODYNAMICS I LECTURE 6 AERODYNAMICS OF A WING FUNDAMENTALS OF THE LIFTING-LINE THEORY

Convexity preserving interpolation by splines of arbitrary degree

Module 3 LOSSY IMAGE COMPRESSION SYSTEMS. Version 2 ECE IIT, Kharagpur

Originated from experimental optimization where measurements are very noisy Approximation can be actually more accurate than

ρ some λ THE INVERSE POWER METHOD (or INVERSE ITERATION) , for , or (more usually) to

A NOTE ON CES FUNCTIONS Drago Bergholt, BI Norwegian Business School 2011

COMPUTATIONALLY EFFICIENT WAVELET AFFINE INVARIANT FUNCTIONS FOR SHAPE RECOGNITION. Erdem Bala, Dept. of Electrical and Computer Engineering,

Instituto Tecnológico de Aeronáutica FINITE ELEMENTS I. Class notes AE-245

Error Bars in both X and Y

The Geometry of Logit and Probit

ACTM State Calculus Competition Saturday April 30, 2011

Problem Do any of the following determine homomorphisms from GL n (C) to GL n (C)?

ANALYSIS OF ELECTROMAGNETIC FIELD USING THE CONSTRAINED INTERPOLATION PROFILE METHOD PHÂN TÍCH TRƯỜNG ĐIỆN TỪ SỬ DỤNG PHƯƠNG PHÁP CIP

IMGS-261 Solutions to Homework #9

Numerical Heat and Mass Transfer

Problem adapted reduced models based on Reaction-Diffusion Manifolds (REDIMs)

Rigid body simulation

ALGORITHM FOR THE CALCULATION OF THE TWO VARIABLES CUBIC SPLINE FUNCTION

Tensor Analysis. For orthogonal curvilinear coordinates, ˆ ˆ (98) Expanding the derivative, we have, ˆ. h q. . h q h q

Modelli Clamfim Equazione del Calore Lezione ottobre 2014

GEMINI GEneric Multimedia INdexIng

Modeling curves. Graphs: y = ax+b, y = sin(x) Implicit ax + by + c = 0, x 2 +y 2 =r 2 Parametric:

Quantum Mechanics I - Session 4

8.1 Arc Length. What is the length of a curve? How can we approximate it? We could do it following the pattern we ve used before

9 Characteristic classes

A NEW DISCRETE WAVELET TRANSFORM

MEM 255 Introduction to Control Systems Review: Basics of Linear Algebra

JAB Chain. Long-tail claims development. ASTIN - September 2005 B.Verdier A. Klinger

APPROXIMATE PRICES OF BASKET AND ASIAN OPTIONS DUPONT OLIVIER. Premia 14

PHYS 705: Classical Mechanics. Calculus of Variations II

The Fundamental Theorem of Algebra. Objective To use the Fundamental Theorem of Algebra to solve polynomial equations with complex solutions

Affine transformations and convexity

5.04, Principles of Inorganic Chemistry II MIT Department of Chemistry Lecture 32: Vibrational Spectroscopy and the IR

Problems & Techniques

Centre of Rotation Determination Using Projection Data in X-ray Micro Computed Tomography

5.62 Physical Chemistry II Spring 2008

Important Instructions to the Examiners:

Numerical Methods. ME Mechanical Lab I. Mechanical Engineering ME Lab I

Some Comments on Accelerating Convergence of Iterative Sequences Using Direct Inversion of the Iterative Subspace (DIIS)

Impulse Noise Removal Technique Based on Fuzzy Logic

= = = (a) Use the MATLAB command rref to solve the system. (b) Let A be the coefficient matrix and B be the right-hand side of the system.

However, since P is a symmetric idempotent matrix, of P are either 0 or 1 [Eigen-values

THE CHINESE REMAINDER THEOREM. We should thank the Chinese for their wonderful remainder theorem. Glenn Stevens

Turbulence. Lecture 21. Non-linear Dynamics. 30 s & 40 s Taylor s work on homogeneous turbulence Kolmogorov.

Affine and Riemannian Connections

COMPLEX NUMBERS AND QUADRATIC EQUATIONS

12 MATH 101A: ALGEBRA I, PART C: MULTILINEAR ALGEBRA. 4. Tensor product

THE VIBRATIONS OF MOLECULES II THE CARBON DIOXIDE MOLECULE Student Instructions

MAE140 - Linear Circuits - Winter 16 Final, March 16, 2016

RECONSTRUCTION OF NON- CARTESIAN DATA USING BURS/RBURS ALGORITHM

Causal Diamonds. M. Aghili, L. Bombelli, B. Pilgrim

A New Design Approach for Recursive Diamond-Shaped Filters

NON-LINEAR CONVOLUTION: A NEW APPROACH FOR THE AURALIZATION OF DISTORTING SYSTEMS

Lecture 6/7 (February 10/12, 2014) DIRAC EQUATION. The non-relativistic Schrödinger equation was obtained by noting that the Hamiltonian 2

Feb 14: Spatial analysis of data fields

CS 523: Computer Graphics, Spring Shape Modeling. PCA Applications + SVD. Andrew Nealen, Rutgers, /15/2011 1

Geometrical Optics Mirrors and Prisms

Linear Approximation with Regularization and Moving Least Squares

Spatial Statistics and Analysis Methods (for GEOG 104 class).

New Method for Solving Poisson Equation. on Irregular Domains

Signal space Review on vector space Linear independence Metric space and norm Inner product

COMPOSITE BEAM WITH WEAK SHEAR CONNECTION SUBJECTED TO THERMAL LOAD

Lecture 13 APPROXIMATION OF SECOMD ORDER DERIVATIVES

OPTIMISATION. Introduction Single Variable Unconstrained Optimisation Multivariable Unconstrained Optimisation Linear Programming

x yi In chapter 14, we want to perform inference (i.e. calculate confidence intervals and perform tests of significance) in this setting.

Computational Astrophysics

Canonical transformations

Solutions to Homework 7, Mathematics 1. 1 x. (arccos x) (arccos x) 1

Gaussian Conditional Random Field Network for Semantic Segmentation - Supplementary Material

Modal Identification of the Elastic Properties in Composite Sandwich Structures

Development of a Semi-Automated Approach for Regional Corrector Surface Modeling in GPS-Levelling

Workshop: Approximating energies and wave functions Quantum aspects of physical chemistry

NMT EE 589 & UNM ME 482/582 ROBOT ENGINEERING. Dr. Stephen Bruder NMT EE 589 & UNM ME 482/582

C/CS/Phy191 Problem Set 3 Solutions Out: Oct 1, 2008., where ( 00. ), so the overall state of the system is ) ( ( ( ( 00 ± 11 ), Φ ± = 1

Assessing inter-annual and seasonal variability Least square fitting with Matlab: Application to SSTs in the vicinity of Cape Town

Lecture 3. Camera Models 2 & Camera Calibration. Professor Silvio Savarese Computational Vision and Geometry Lab. 13- Jan- 15.

Inner Product. Euclidean Space. Orthonormal Basis. Orthogonal

Transcription:

ECE 472/572 - Dgtal Image Processng Lecture 6 Geometrc and Radometrc Transformaton 09/27/ Roadmap Introducton Image format vector vs. btmap IP vs. CV vs. CG HLIP vs. LLIP Image acquston Percepton Structure of human ee Brghtness adaptaton and Dscrmnaton Image resoluton Image enhancement Enhancement vs. restoraton Spatal doman methods Pont-based methods Log trans. vs. Power-law Contrast stretchng vs. HE Gra-level vs. Bt plane slcng Image averagng prncple Mask-based methods - spatal flter Smoothng vs. Sharpenng flter Lnear vs. Non-lnear flter Smoothng average vs. Gaussan vs. medan Sharpenng UM vs. st vs. 2nd dervatves Frequenc doman methods Understandng Fourer transform Implementaton n the frequenc doman Low-pass flters vs. hgh-pass flters vs. homomorphc flter Geometrc correcton Affne vs. Perspectve Homogeneous coordnates Inverse vs. forward transform Composte General Model dstorton wth polnomal Least square soluton 2 Questons Affne vs. Perspectve Forward vs. Inverse Composte vs. Sequental Homogeneous coordnate General geometrc s 3

Usage Image correcton Color nterpolaton Forensc analss Entertanment effect http://www.mp-sb.mpg.de/resources/fam/demos.html 4 http://w3.mpa.br/~morph/ Affne s " u " a Preserve lnes and parallel lnes a 2 a 3 " x Homogeneous coordnates v a 2 a 22 a 23 0 0 General form a a2 a3! a2 a22 a23! 0 0!" Specal matrces R: rotaton, S: scalng, T: translaton, H: shear cosθ snθ 0 sx 0 0 0 tx hx 0 R snθ cosθ 0,S 0 s 0,T 0 t,h h 0 0 0 0 0 0 0 0 0 5 6 2

3 7 Composte vs. Sequental Orgnal Image fx, R H T S Transformed Image gu, v u v " S T H R x " R H T S C 8 Forward vs. Inverse transforms u v " C x " Forward transform C u v x Inverse transform 9 Examples - Shear h 0.2 hx 0.2 hx h 0.2

Examples Translaton + Rotaton theta PI/4 theta PI/4 tx -40, t 60 0 Perspectve Preserve parallel lnes onl when the are parallel to the projecton plane. Otherwse, lnes converge to a vanshng pont General form " a a 2 a 3 a 2 a 22 a 23 a 3 0, a 32 0 a 3 a 32 u " a a 2 a 3 x v " a 2 a 22 a 23,u u " w ",v v " w " w " a 3 a 32 Determne the coeffcents u " a a 2 a 3 x v " a 2 a 22 a 23,u u " w ",v v " w " w " a 3 a 32 8 unknowns, 4-pont least squares 0,0 0,255 0,0 0,255 255,0 255,255 255,0 5,5 2 4

Example - PT 3 General approaches Fnd teponts Spatal 4 Example CCD buttng msalgnment greater than 50 mcron x-ra senstve scntllator fber optcs CCD arra 242 x 52 9/27/ 5 5

Sources of dstortons defects n the producton of fber-optc tapers mperfect compresson and cuttng dfferent lght transfer effcenc across the whole surface 9/27/ 6 Geometrc correcton map close approxmaton P u, v x, control pont nterpolaton T u, v x, x, u, v map exactl 9/27/ 7 Spatal Blnear equaton x ˆ r u,v a u + a 2 v + a 3 uv + a 4 ˆ s u,v b u + b 2 v + b 3 uv + b 4 n-th degree polnomal " x ˆ " ˆ P xu,v d " a krs u r s v P u,v b krs u r s k 0 r+s k v Use nformaton from teponts to solve coeffcents Exact soluton Least square soluton ε mn m 2 2 a, b [ x xˆ + ˆ ] 0 8 6

How s t appled? Step : Choose a set of te ponts x, : coordnates of te ponts n the orgnal or dstorted mage u,v : coordnates of te ponts n the corrected mage Step 2: Decde on whch degree of polnomal to use to model the nverse of the dstorton, e.g., x ˆ r u,v a u + a 2 v + a 3 uv + a 4 ˆ s u,v b u + b 2 v + b 3 uv + b 4 Step 3: Solve the coeffcents of the polnomal usng least-squares approach Step 4: Use the derved polnomal model to correct the entre orgnal mage X [ x 0 x x 24 ] T Y [ 0 24 ] T " u 0 v 0 u 0 v 0 u W v u v u 24 v 24 u 24 v 24 A W X,B W Y W W T W W T A [ a a 2 a 3 a 4 ] T B [ b b 2 b 3 b 4 ] T For each u,v n the corrected mage, fnd the correspondng x, n the orgnal mage and use ts ntenst as the ntenst at u,v. 9 Example Geometrc correcton Geometrc correcton of mages from butted CCD arras 20 Example - Color correcton Teponts are colors R, G, B, nstead of spatal coordnates 2 7

Example - Image warpng 22 Image warpng Two-pass mesh warpng b Douglas Smthe Reference: G. Wolberg, Dgtal Image Warpng, 990 23 Example From Joe Howell and Cor McKa, ECE472, Fall 2000 24 8

Example 2 From Adam Mller, Truman Bonds, Randal Waldrop, ECE472, Fall 2000 25 9