WORCESTER POLYTECHNIC INSTITUTE

Size: px
Start display at page:

Download "WORCESTER POLYTECHNIC INSTITUTE"

Transcription

1 WORCESTER POLYTECHNIC INSTITUTE MECHANICAL ENGINEERING DEPARTMENT Optical Metrology and NDT ME-593L, C 2018 Introduction: Fringe Skeletonization February 2018

2 Quantitative analysis Fringe skeletonization Recording o intererograms Identiication o boundary area (AOI) Preprocessing, e.g., iltering Identiication o ringe centers Numbering o intererence ringes: ringe ordering Reconstruction o phase using interpolation Review reerence papers (list is on next page)

3 Reerence papers: Quantitative analysis Fringe skeletonization J. Novak, Techniques or automatic identiication and numbering o intererence ringes using MATLAB. T. Merz, D. Paulus, and H. Niemann, Line segmentation or intererograms o continuously deorming objects. L. Z. Cai, Q. Liu, and X. L. Yang, A simple method o contrast enhancement and extremum extraction or intererence ringes. J. Budzinski, SNOP: a method or skeletonization o a ringe pattern along the ringe direction.

4 Recording o intererograms Classic intererometry Michelson Mach-Zehnder Sagnac Holographic intererometry Digital holographic intererometry Speckle pattern intererometry White-light intererometry

5 Identiication o boundary area Identiication o area o interest (AOI) Area o operation Selection o convolution Kernel Minimization o power leakage

6 Pre-processing, e.g., iltering: convolution Application o digital spatial convolution to the pixel located at the (m,n) image plane position Image convolution: m Q D ( m, n) H Q' ( m, n) kl D m 1 m m+1 n 1 n * (0,0) = n n+1 ( m, n) Q D H kl ( m, n) Q ' D Original image Kernel Convolved image

7 Pre-processing, e.g., iltering: convolution Digital spatial convolution by scanning the convolution kernel line by line over the entire image plane m n Kernel (3 x 3, in this case)

8 Pre-processing, e.g., iltering: convolution Low-pass iltering Image convolution: ' 1 Q D ( m, n) QD ( m k, n l) ( 2r 1) ( 2r 1) r k r r l r Weighting actors: H kl 1 ( 2r 1) ( 2r 1)

9 Pre-processing, e.g., iltering: convolution Median iltering Digital spatial convolution to median iltering the pixel value located at the (m,n) image plane position Original image Convolved image m m n n (a) (b) Kernel (3 x 3, in this case) Sorted list

10 Fringe skeletonization Pre-processing, e.g., iltering: convolution Fourier iltering Fourier transormation: )}, ( { ), ( n m Q Q D v u D Convolution in the requency domain: ), ( ), ( ), ( v u v u D v u T D W Q Q ), ( v W u is the weighting unction Inverse Fourier transormation: )}, ( { ), ( 1 v u D ' D Q m n Q T

11 Pre-processing, e.g., iltering: convolution Fourier iltering Fourier transormation: Q D (, ) { Q ( m, n)} u v D Original intererogram Original intererogram: requency domain

12 Weighting actor Fringe skeletonization Pre-processing, e.g., iltering: convolution Fourier iltering W (, v) u, weighting unction: Gaussian unction 2D representation 3D representation

13 Pre-processing, e.g., iltering: convolution Fourier iltering Convolution in the requency domain: QT D( u, v ) QD ( u, v ) W ( u, v ) ' 1 T Inverse Fourier transormation: Q D ( m, n) { QD( u, v)} Original intererogram Filtered intererogram

14 Pre-processing: thinning Multiple algorithms exist, e.g., contrast enhancement, edge detection, etc. Original intererogram Thinned ringes: contrast enhancement

15 Identiication o ringe centers Fringe centers Thinned ringes: contrast enhancement

16 Semi-quantitative analysis 1 ringe is 2/: ringe ordering and counting Intererogram o a turbine blade: contouring Intererogram o a turbine blade: vibrations Fringe-locus unction (ringe localization): ( x, y, z) K L ( x, y, z) 2 n n = is the ringe order A ringe represents a contour o constant phase n is the ringe order or number o waves

17 Sample: identiication o ringe centers Closed ringe pattern and continuous deormations:

18 Fringe ordering: n-th order Identiy zero-order: understanding o the physical phenomena Use sequential ordering Fringe orders: n Zero order

19 Fringe ordering Opened ringe pattern and continuous deormations: A A Intensity line: A-A Fringe ordering:

20 Phase reconstruction Fringe interpolation, e.g., Lagrange, Zernike polynomials, etc. Lagrange polynomials o degree n-1 x o

21 Phase reconstruction Fringe interpolation, e.g., Lagrange, Zernike polynomials, etc. Lagrange polynomials o degree n-1 x o. x 2

22 Phase reconstruction: 1-ringe = /2 Reconstruction along speciic lines: Interpolation along speciic lines:

Two-step self-tuning phase-shifting interferometry

Two-step self-tuning phase-shifting interferometry Two-step sel-tuning phase-shiting intererometry J. Vargas, 1,* J. Antonio Quiroga, T. Belenguer, 1 M. Servín, 3 J. C. Estrada 3 1 Laboratorio de Instrumentación Espacial, Instituto Nacional de Técnica

More information

WORCESTER POLYTECHNIC INSTITUTE

WORCESTER POLYTECHNIC INSTITUTE WORCESTER POLYTECHNIC INSTITUTE MECHANICAL ENGINEERING DEPARTMENT Optical Metrology and NDT ME-593L, C 018 Introduction: Wave Optics January 018 Wave optics: light waves Wave equation An optical wave --

More information

Lecture 13: Applications of Fourier transforms (Recipes, Chapter 13)

Lecture 13: Applications of Fourier transforms (Recipes, Chapter 13) Lecture 13: Applications o Fourier transorms (Recipes, Chapter 13 There are many applications o FT, some o which involve the convolution theorem (Recipes 13.1: The convolution o h(t and r(t is deined by

More information

Lab 3: The FFT and Digital Filtering. Slides prepared by: Chun-Te (Randy) Chu

Lab 3: The FFT and Digital Filtering. Slides prepared by: Chun-Te (Randy) Chu Lab 3: The FFT and Digital Filtering Slides prepared by: Chun-Te (Randy) Chu Lab 3: The FFT and Digital Filtering Assignment 1 Assignment 2 Assignment 3 Assignment 4 Assignment 5 What you will learn in

More information

Notes on Wavelets- Sandra Chapman (MPAGS: Time series analysis) # $ ( ) = G f. y t

Notes on Wavelets- Sandra Chapman (MPAGS: Time series analysis) # $ ( ) = G f. y t Wavelets Recall: we can choose! t ) as basis on which we expand, ie: ) = y t ) = G! t ) y t! may be orthogonal chosen or appropriate properties. This is equivalent to the transorm: ) = G y t )!,t )d 2

More information

Filtering in Frequency Domain

Filtering in Frequency Domain Dr. Praveen Sankaran Department of ECE NIT Calicut February 4, 2013 Outline 1 2D DFT - Review 2 2D Sampling 2D DFT - Review 2D Impulse Train s [t, z] = m= n= δ [t m T, z n Z] (1) f (t, z) s [t, z] sampled

More information

Digital Image Processing. Lecture 6 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 2009

Digital Image Processing. Lecture 6 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 2009 Digital Image Processing Lecture 6 (Enhancement) Bu-Ali Sina University Computer Engineering Dep. Fall 009 Outline Image Enhancement in Spatial Domain Spatial Filtering Smoothing Filters Median Filter

More information

Midterm Summary Fall 08. Yao Wang Polytechnic University, Brooklyn, NY 11201

Midterm Summary Fall 08. Yao Wang Polytechnic University, Brooklyn, NY 11201 Midterm Summary Fall 8 Yao Wang Polytechnic University, Brooklyn, NY 2 Components in Digital Image Processing Output are images Input Image Color Color image image processing Image Image restoration Image

More information

CAP 5415 Computer Vision Fall 2011

CAP 5415 Computer Vision Fall 2011 CAP 545 Computer Vision Fall 2 Dr. Mubarak Sa Univ. o Central Florida www.cs.uc.edu/~vision/courses/cap545/all22 Oice 247-F HEC Filtering Lecture-2 General Binary Gray Scale Color Binary Images Y Row X

More information

CAP 5415 Computer Vision

CAP 5415 Computer Vision CAP 545 Computer Vision Dr. Mubarak Sa Univ. o Central Florida Filtering Lecture-2 Contents Filtering/Smooting/Removing Noise Convolution/Correlation Image Derivatives Histogram Some Matlab Functions General

More information

Reading. 3. Image processing. Pixel movement. Image processing Y R I G Q

Reading. 3. Image processing. Pixel movement. Image processing Y R I G Q Reading Jain, Kasturi, Schunck, Machine Vision. McGraw-Hill, 1995. Sections 4.-4.4, 4.5(intro), 4.5.5, 4.5.6, 5.1-5.4. 3. Image processing 1 Image processing An image processing operation typically defines

More information

Least-Squares Spectral Analysis Theory Summary

Least-Squares Spectral Analysis Theory Summary Least-Squares Spectral Analysis Theory Summary Reerence: Mtamakaya, J. D. (2012). Assessment o Atmospheric Pressure Loading on the International GNSS REPRO1 Solutions Periodic Signatures. Ph.D. dissertation,

More information

Additional Pointers. Introduction to Computer Vision. Convolution. Area operations: Linear filtering

Additional Pointers. Introduction to Computer Vision. Convolution. Area operations: Linear filtering Additional Pointers Introduction to Computer Vision CS / ECE 181B andout #4 : Available this afternoon Midterm: May 6, 2004 W #2 due tomorrow Ack: Prof. Matthew Turk for the lecture slides. See my ECE

More information

Design of Coprime DFT Arrays and Filter Banks

Design of Coprime DFT Arrays and Filter Banks Design o Coprime DFT Arrays and Filter Banks Chun-Lin Liu clliu@caltechedu P P Vaidyanathan ppvnath@systemscaltechedu Digital Signal Processing Group Electrical Engineering Caliornia Institute o Technology

More information

Edge Detection. CS 650: Computer Vision

Edge Detection. CS 650: Computer Vision CS 650: Computer Vision Edges and Gradients Edge: local indication of an object transition Edge detection: local operators that find edges (usually involves convolution) Local intensity transitions are

More information

TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall Problem 1.

TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall TLT-5200/5206 COMMUNICATION THEORY, Exercise 3, Fall Problem 1. TLT-5/56 COMMUNICATION THEORY, Exercise 3, Fall Problem. The "random walk" was modelled as a random sequence [ n] where W[i] are binary i.i.d. random variables with P[W[i] = s] = p (orward step with probability

More information

Presenters. Time Domain and Statistical Model Development, Simulation and Correlation Methods for High Speed SerDes

Presenters. Time Domain and Statistical Model Development, Simulation and Correlation Methods for High Speed SerDes JANUARY 28-31, 2013 SANTA CLARA CONVENTION CENTER Time Domain and Statistical Model Development, Simulation and Correlation Methods or High Speed SerDes Presenters Xingdong Dai Fangyi Rao Shiva Prasad

More information

Spatial Enhancement Region operations: k'(x,y) = F( k(x-m, y-n), k(x,y), k(x+m,y+n) ]

Spatial Enhancement Region operations: k'(x,y) = F( k(x-m, y-n), k(x,y), k(x+m,y+n) ] CEE 615: Digital Image Processing Spatial Enhancements 1 Spatial Enhancement Region operations: k'(x,y) = F( k(x-m, y-n), k(x,y), k(x+m,y+n) ] Template (Windowing) Operations Template (window, box, kernel)

More information

Edge Detection. Image Processing - Computer Vision

Edge Detection. Image Processing - Computer Vision Image Processing - Lesson 10 Edge Detection Image Processing - Computer Vision Low Level Edge detection masks Gradient Detectors Compass Detectors Second Derivative - Laplace detectors Edge Linking Image

More information

The Fourier Transform

The Fourier Transform The Fourier Transorm Fourier Series Fourier Transorm The Basic Theorems and Applications Sampling Bracewell, R. The Fourier Transorm and Its Applications, 3rd ed. New York: McGraw-Hill, 2. Eric W. Weisstein.

More information

A Comparative Study of Non-separable Wavelet and Tensor-product. Wavelet; Image Compression

A Comparative Study of Non-separable Wavelet and Tensor-product. Wavelet; Image Compression Copyright c 007 Tech Science Press CMES, vol., no., pp.91-96, 007 A Comparative Study o Non-separable Wavelet and Tensor-product Wavelet in Image Compression Jun Zhang 1 Abstract: The most commonly used

More information

Application of Wavelet Transform Modulus Maxima in Raman Distributed Temperature Sensors

Application of Wavelet Transform Modulus Maxima in Raman Distributed Temperature Sensors PHOTONIC SENSORS / Vol. 4, No. 2, 2014: 142 146 Application o Wavelet Transorm Modulus Maxima in Raman Distributed Temperature Sensors Zongliang WANG, Jun CHANG *, Sasa ZHANG, Sha LUO, Cuanwu JIA, Boning

More information

Image Recognition Using Modified Zernike Moments

Image Recognition Using Modified Zernike Moments Sensors & Transducers 204 by IFSA Publishing S. L. http://www.sensorsportal.com Image Recognition Using Modified ernike Moments Min HUANG Yaqiong MA and Qiuping GONG School of Computer and Communication

More information

Random Error Analysis of Inertial Sensors output Based on Allan Variance Shaochen Li1, a, Xiaojing Du2,b and Junyi Zhai3,c

Random Error Analysis of Inertial Sensors output Based on Allan Variance Shaochen Li1, a, Xiaojing Du2,b and Junyi Zhai3,c International Conerence on Civil, Transportation and Environment (ICCTE 06) Random Error Analysis o Inertial Sensors output Based on Allan Variance Shaochen Li, a, Xiaojing Du, and Junyi Zhai3,c School

More information

Low order continuous-time filters for approximation of the ISO human vibration sensitivity weightings

Low order continuous-time filters for approximation of the ISO human vibration sensitivity weightings Journal o Sound and Vibration 265 (23) 459 465 JOURNAL OF SOUND AND VIBRATION www.elsevier.com/locate/jsvi Letter to the Editor Low order continuous-time ilters or approximation o the ISO 263- human vibration

More information

Image Acquisition and Sampling Theory

Image Acquisition and Sampling Theory Image Acquisition and Sampling Theory Electromagnetic Spectrum The wavelength required to see an object must be the same size of smaller than the object 2 Image Sensors 3 Sensor Strips 4 Digital Image

More information

Introduction to Sparsity in Signal Processing

Introduction to Sparsity in Signal Processing 1 Introduction to Sparsity in Signal Processing Ivan Selesnick Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 212 2 Under-determined linear equations Consider a system

More information

ITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms

ITK Filters. Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms ITK Filters Thresholding Edge Detection Gradients Second Order Derivatives Neighborhood Filters Smoothing Filters Distance Map Image Transforms ITCS 6010:Biomedical Imaging and Visualization 1 ITK Filters:

More information

Introduction to Computer Vision. 2D Linear Systems

Introduction to Computer Vision. 2D Linear Systems Introduction to Computer Vision D Linear Systems Review: Linear Systems We define a system as a unit that converts an input function into an output function Independent variable System operator or Transfer

More information

Feature Extraction and Image Processing

Feature Extraction and Image Processing Feature Extraction and Image Processing Second edition Mark S. Nixon Alberto S. Aguado :*авш JBK IIP AMSTERDAM BOSTON HEIDELBERG LONDON NEW YORK OXFORD PARIS SAN DIEGO SAN FRANCISCO SINGAPORE SYDNEY TOKYO

More information

Chapter 4 Imaging. Lecture 21. d (110) Chem 793, Fall 2011, L. Ma

Chapter 4 Imaging. Lecture 21. d (110) Chem 793, Fall 2011, L. Ma Chapter 4 Imaging Lecture 21 d (110) Imaging Imaging in the TEM Diraction Contrast in TEM Image HRTEM (High Resolution Transmission Electron Microscopy) Imaging or phase contrast imaging STEM imaging a

More information

CSE 421 Algorithms. T(n) = at(n/b) + n c. Closest Pair Problem. Divide and Conquer Algorithms. What you really need to know about recurrences

CSE 421 Algorithms. T(n) = at(n/b) + n c. Closest Pair Problem. Divide and Conquer Algorithms. What you really need to know about recurrences CSE 421 Algorithms Richard Anderson Lecture 13 Divide and Conquer What you really need to know about recurrences Work per level changes geometrically with the level Geometrically increasing (x > 1) The

More information

Chirp images in 2-D fractional Fourier transform domain. Presenter: Ming-Feng Lu Beijing Institute of Technology July 20, 2016

Chirp images in 2-D fractional Fourier transform domain. Presenter: Ming-Feng Lu Beijing Institute of Technology July 20, 2016 Chirp images in -D fractional Fourier transform domain Presenter: Ming-Feng Lu lumingfeng@bit.edu.cn Beijing Institute of Technology July 0, 016 目录 Introduction of chirp images Chirp images in FRFT domain

More information

Systems & Signals 315

Systems & Signals 315 1 / 15 Systems & Signals 315 Lecture 13: Signals and Linear Systems Dr. Herman A. Engelbrecht Stellenbosch University Dept. E & E Engineering 2 Maart 2009 Outline 2 / 15 1 Signal Transmission through a

More information

Introduction to Sparsity in Signal Processing

Introduction to Sparsity in Signal Processing 1 Introduction to Sparsity in Signal Processing Ivan Selesnick Polytechnic Institute of New York University Brooklyn, New York selesi@poly.edu 212 2 Under-determined linear equations Consider a system

More information

In many diverse fields physical data is collected or analysed as Fourier components.

In many diverse fields physical data is collected or analysed as Fourier components. 1. Fourier Methods In many diverse ields physical data is collected or analysed as Fourier components. In this section we briely discuss the mathematics o Fourier series and Fourier transorms. 1. Fourier

More information

Time Frequency Aggregation Performance Optimization of Power Quality Disturbances Based on Generalized S Transform

Time Frequency Aggregation Performance Optimization of Power Quality Disturbances Based on Generalized S Transform Time Frequency Aggregation Perormance Otimization o Power Quality Disturbances Based on Generalized S Transorm Mengda Li Shanghai Dianji University, Shanghai 01306, China limd @ sdju.edu.cn Abstract In

More information

Exam 2. Average: 85.6 Median: 87.0 Maximum: Minimum: 55.0 Standard Deviation: Numerical Methods Fall 2011 Lecture 20

Exam 2. Average: 85.6 Median: 87.0 Maximum: Minimum: 55.0 Standard Deviation: Numerical Methods Fall 2011 Lecture 20 Exam 2 Average: 85.6 Median: 87.0 Maximum: 100.0 Minimum: 55.0 Standard Deviation: 10.42 Fall 2011 1 Today s class Multiple Variable Linear Regression Polynomial Interpolation Lagrange Interpolation Newton

More information

Lecture Outline. Basics of Spatial Filtering Smoothing Spatial Filters. Sharpening Spatial Filters

Lecture Outline. Basics of Spatial Filtering Smoothing Spatial Filters. Sharpening Spatial Filters 1 Lecture Outline Basics o Spatial Filtering Smoothing Spatial Filters Averaging ilters Order-Statistics ilters Sharpening Spatial Filters Laplacian ilters High-boost ilters Gradient Masks Combining Spatial

More information

Edge Detection. Introduction to Computer Vision. Useful Mathematics Funcs. The bad news

Edge Detection. Introduction to Computer Vision. Useful Mathematics Funcs. The bad news Edge Detection Introduction to Computer Vision CS / ECE 8B Thursday, April, 004 Edge detection (HO #5) Edge detection is a local area operator that seeks to find significant, meaningful changes in image

More information

Therefore the new Fourier coefficients are. Module 2 : Signals in Frequency Domain Problem Set 2. Problem 1

Therefore the new Fourier coefficients are. Module 2 : Signals in Frequency Domain Problem Set 2. Problem 1 Module 2 : Signals in Frequency Domain Problem Set 2 Problem 1 Let be a periodic signal with fundamental period T and Fourier series coefficients. Derive the Fourier series coefficients of each of the

More information

Automatic detection and of dipoles in large area SQUID magnetometry

Automatic detection and of dipoles in large area SQUID magnetometry Automatic detection and of dipoles in large area SQUID magnetometry Lisa Qian December 4, Introduction. Scanning SQUID magnetometry Scanning SQUID magnetometry is a powerful tool for metrology of individual

More information

Outline. Convolution. Filtering

Outline. Convolution. Filtering Filtering Outline Convolution Filtering Logistics HW1 HW2 - out tomorrow Recall: what is a digital (grayscale) image? Matrix of integer values Images as height fields Let s think of image as zero-padded

More information

Enhancement Using Local Histogram

Enhancement Using Local Histogram Enhancement Using Local Histogram Used to enhance details over small portions o the image. Deine a square or rectangular neighborhood hose center moves rom piel to piel. Compute local histogram based on

More information

Strain Analysis by Digital Shearography on an Aluminium Plate with a Geometric Defect under Thermal Loading

Strain Analysis by Digital Shearography on an Aluminium Plate with a Geometric Defect under Thermal Loading Journal of Stress Analysis Vol. 1, No. 1, Spring Summer 2016 Strain Analysis by Digital Shearography on an Aluminium Plate with a Geometric Defect under Thermal Loading R. Moharrami, M. Taghiloo, A. Darud

More information

Optical Sciences Center, Rm 704 University of Arizona Tucson, AZ Office Hours: Call for appointment or see after class

Optical Sciences Center, Rm 704 University of Arizona Tucson, AZ Office Hours: Call for appointment or see after class Term: Spring 2000 Course #: OPTI 505 Course Title: Diffraction and Interferometry Instructor: James C. Wyant Optical Sciences Center, Rm 704 University of Arizona Tucson, AZ 85721 Phone: 520-621-2448 E-Mail:

More information

Laplacian Filters. Sobel Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters

Laplacian Filters. Sobel Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters. Laplacian Filters Sobel Filters Note that smoothing the image before applying a Sobel filter typically gives better results. Even thresholding the Sobel filtered image cannot usually create precise, i.e., -pixel wide, edges.

More information

Metrology and Sensing

Metrology and Sensing Metrology and Sensing Lecture 5: Interferometry I 06--09 Herbert Gross Winter term 06 www.iap.uni-jena.de Preliminary Schedule No Date Subject Detailed Content 8.0. Introduction Introduction, optical measurements,

More information

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals.

Z - Transform. It offers the techniques for digital filter design and frequency analysis of digital signals. Z - Transform The z-transform is a very important tool in describing and analyzing digital systems. It offers the techniques for digital filter design and frequency analysis of digital signals. Definition

More information

Image preprocessing in spatial domain

Image preprocessing in spatial domain Image preprocessing in spatial domain Sampling theorem, aliasing, interpolation, geometrical transformations Revision: 1.4, dated: May 25, 2006 Tomáš Svoboda Czech Technical University, Faculty of Electrical

More information

Lecture 04 Image Filtering

Lecture 04 Image Filtering Institute of Informatics Institute of Neuroinformatics Lecture 04 Image Filtering Davide Scaramuzza 1 Lab Exercise 2 - Today afternoon Room ETH HG E 1.1 from 13:15 to 15:00 Work description: your first

More information

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is:

Computed Tomography Notes, Part 1. The equation that governs the image intensity in projection imaging is: Noll 3 CT Notes : Page Compute Tomograph Notes Part Challenges with Projection X-ra Sstems The equation that governs the image intensit in projection imaging is: z I I ep µ z Projection -ra sstems are

More information

TRACKING and DETECTION in COMPUTER VISION Filtering and edge detection

TRACKING and DETECTION in COMPUTER VISION Filtering and edge detection Technischen Universität München Winter Semester 0/0 TRACKING and DETECTION in COMPUTER VISION Filtering and edge detection Slobodan Ilić Overview Image formation Convolution Non-liner filtering: Median

More information

Image preprocessing in spatial domain

Image preprocessing in spatial domain Image preprocessing in spatial domain Sampling theorem, aliasing, interpolation, geometrical transformations Revision: 1.3, dated: December 7, 2005 Tomáš Svoboda Czech Technical University, Faculty of

More information

Probabilistic Model of Error in Fixed-Point Arithmetic Gaussian Pyramid

Probabilistic Model of Error in Fixed-Point Arithmetic Gaussian Pyramid Probabilistic Model o Error in Fixed-Point Arithmetic Gaussian Pyramid Antoine Méler John A. Ruiz-Hernandez James L. Crowley INRIA Grenoble - Rhône-Alpes 655 avenue de l Europe 38 334 Saint Ismier Cedex

More information

SIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( )

SIO 211B, Rudnick. We start with a definition of the Fourier transform! ĝ f of a time series! ( ) SIO B, Rudnick! XVIII.Wavelets The goal o a wavelet transorm is a description o a time series that is both requency and time selective. The wavelet transorm can be contrasted with the well-known and very

More information

Fourier transform spectroscopy: an introduction. David Naylor University of Lethbridge

Fourier transform spectroscopy: an introduction. David Naylor University of Lethbridge Fourier transform spectroscopy: an introduction David Naylor University of Lethbridge Outline History Ideal vs real FTS Pros/cons Extension to ifts Examples: Sitelle, SPIRE, Safari Michelson s original

More information

COMPUTER GENERATED HOLOGRAMS Optical Sciences W.J. Dallas PART II: CHAPTER ONE HOLOGRAPHY IN A NUTSHELL

COMPUTER GENERATED HOLOGRAMS Optical Sciences W.J. Dallas PART II: CHAPTER ONE HOLOGRAPHY IN A NUTSHELL What is a Hologram? Holography in a Nutshell: Page 1 o 1 C:\_Dallas\_Courses\3_OpSci_67 8\ MsWord\_TheCgh\1_MSWord\_1 Holography.doc Version: Wednesday, September 4, 8, 8: AM COMPUTER GENERATED HOLOGRAMS

More information

Fourier Analysis and Imaging Ronald Bracewell L.M. Terman Professor of Electrical Engineering Emeritus Stanford University Stanford, California

Fourier Analysis and Imaging Ronald Bracewell L.M. Terman Professor of Electrical Engineering Emeritus Stanford University Stanford, California Fourier Analysis and Imaging Ronald Bracewell L.M. Terman Professor of Electrical Engineering Emeritus Stanford University Stanford, California 4u Springer Contents fa 1 PREFACE INTRODUCTION xiii 1 11

More information

Case Studies of Logical Computation on Stochastic Bit Streams

Case Studies of Logical Computation on Stochastic Bit Streams Case Studies of Logical Computation on Stochastic Bit Streams Peng Li 1, Weikang Qian 2, David J. Lilja 1, Kia Bazargan 1, and Marc D. Riedel 1 1 Electrical and Computer Engineering, University of Minnesota,

More information

Role of Assembling Invariant Moments and SVM in Fingerprint Recognition

Role of Assembling Invariant Moments and SVM in Fingerprint Recognition 56 Role of Assembling Invariant Moments SVM in Fingerprint Recognition 1 Supriya Wable, 2 Chaitali Laulkar 1, 2 Department of Computer Engineering, University of Pune Sinhgad College of Engineering, Pune-411

More information

Optical Metrology and Sensing

Optical Metrology and Sensing Optical Metrology and Sensing Lecture 1: Introduction 2016-10-18 Herbert Gross Winter term 2016 www.iap.uni-jena.de 2 Preliminary Schedule No Date Subject Detailed Content 1 18.10. Introduction Introduction,

More information

5.6 Convolution and FFT

5.6 Convolution and FFT 5.6 Convolution and FFT Fast Fourier Transform: Applications Applications. Optics, acoustics, quantum physics, telecommunications, control systems, signal processing, speech recognition, data compression,

More information

Various Shape Descriptors in Image Processing A Review

Various Shape Descriptors in Image Processing A Review Various Shape Descriptors in Image Processing A Review Preetika D Silva 1, P. Bhuvaneswari 2 1 PG Student, Department of Electronics and Communication, Rajarajeswari College of Engineering, Bangalore,

More information

Double-slit interference of biphotons generated in spontaneous parametric downconversion from a thick crystal

Double-slit interference of biphotons generated in spontaneous parametric downconversion from a thick crystal INSTITUTE OF PHYSICS PUBLISHING JOURNAL OF OPTICS B: QUANTUM AND SEMICLASSICAL OPTICS J. Opt. B: Quantum Semiclass. Opt. 3 (2001 S50 S54 www.iop.org/journals/ob PII: S1464-4266(0115159-1 Double-slit intererence

More information

APPLICATION OF A SCANNING VIBROMETER FOR THE PERIODIC CALIBRATION OF FORCE TRANSDUCERS

APPLICATION OF A SCANNING VIBROMETER FOR THE PERIODIC CALIBRATION OF FORCE TRANSDUCERS XX IMEKO World Congress Metrology or Green Growth September 9 14, 01, Busan, Republic o Korea APPLICATION OF A SCANNING VIBROMETER FOR THE PERIODIC CALIBRATION OF FORCE TRANSDUCERS Christian Schlegel,

More information

Nyquist sampling a bandlimited function. Sampling: The connection from CT to DT. Impulse Train sampling. Interpolation: Signal Reconstruction

Nyquist sampling a bandlimited function. Sampling: The connection from CT to DT. Impulse Train sampling. Interpolation: Signal Reconstruction Sampling: he connection rom C to D Nyquist sampling a bandlimited unction ( Consider as an example sampling the bandlimited signal sinc 2 (5t F 1 5 Λ 5 An ininite number o waveorms could exactly pass through

More information

3. Lecture. Fourier Transformation Sampling

3. Lecture. Fourier Transformation Sampling 3. Lecture Fourier Transformation Sampling Some slides taken from Digital Image Processing: An Algorithmic Introduction using Java, Wilhelm Burger and Mark James Burge Separability ² The 2D DFT can be

More information

Image Filtering. Slides, adapted from. Steve Seitz and Rick Szeliski, U.Washington

Image Filtering. Slides, adapted from. Steve Seitz and Rick Szeliski, U.Washington Image Filtering Slides, adapted from Steve Seitz and Rick Szeliski, U.Washington The power of blur All is Vanity by Charles Allen Gillbert (1873-1929) Harmon LD & JuleszB (1973) The recognition of faces.

More information

Interference. Reminder: Exam 2 and Review quiz, more details on the course website

Interference. Reminder: Exam 2 and Review quiz, more details on the course website Chapter 9 Interference Phys 322 Lecture 25 Reminder: Exam 2 and Review quiz, more details on the course website Interferometers Wavefront-splitting interferometers Amplitude-splitting interferometers ed

More information

Open Access Measuring Method for Diameter of Bearings Based on the Edge Detection Using Zernike Moment

Open Access Measuring Method for Diameter of Bearings Based on the Edge Detection Using Zernike Moment Send Orders for Reprints to reprints@benthamscience.ae 114 The Open Automation and Control Systems Journal, 2015, 7, 114-119 Open Access Measuring Method for Diameter of Bearings Based on the Edge Detection

More information

Comparative Measurement in Speckle Interferometry using Fourier Transform

Comparative Measurement in Speckle Interferometry using Fourier Transform Egypt. J. Solids, Vol. (30), No. (2), (2007) 253 Comparative Measurement in Speckle Interferometry using Fourier Transform Nasser A. Moustafa Physics Department, Faculty of Science, Helwan University,

More information

Department of Electrical Engineering, Polytechnic University, Brooklyn Fall 05 EL DIGITAL IMAGE PROCESSING (I) Final Exam 1/5/06, 1PM-4PM

Department of Electrical Engineering, Polytechnic University, Brooklyn Fall 05 EL DIGITAL IMAGE PROCESSING (I) Final Exam 1/5/06, 1PM-4PM Department of Electrical Engineering, Polytechnic University, Brooklyn Fall 05 EL512 --- DIGITAL IMAGE PROCESSING (I) Y. Wang Final Exam 1/5/06, 1PM-4PM Your Name: ID Number: Closed book. One sheet of

More information

Computer Vision & Digital Image Processing

Computer Vision & Digital Image Processing Computer Vision & Digital Image Processing Image Restoration and Reconstruction I Dr. D. J. Jackson Lecture 11-1 Image restoration Restoration is an objective process that attempts to recover an image

More information

Efficient Inference in Fully Connected CRFs with Gaussian Edge Potentials

Efficient Inference in Fully Connected CRFs with Gaussian Edge Potentials Efficient Inference in Fully Connected CRFs with Gaussian Edge Potentials Philipp Krähenbühl and Vladlen Koltun Stanford University Presenter: Yuan-Ting Hu 1 Conditional Random Field (CRF) E x I = φ u

More information

Engine fault feature extraction based on order tracking and VMD in transient conditions

Engine fault feature extraction based on order tracking and VMD in transient conditions Engine fault feature extraction based on order tracking and VMD in transient conditions Gang Ren 1, Jide Jia 2, Jian Mei 3, Xiangyu Jia 4, Jiajia Han 5 1, 4, 5 Fifth Cadet Brigade, Army Transportation

More information

Maximal Entropy for Reconstruction of Back Projection Images

Maximal Entropy for Reconstruction of Back Projection Images Maximal Entropy for Reconstruction of Back Projection Images Tryphon Georgiou Department of Electrical and Computer Engineering University of Minnesota Minneapolis, MN 5545 Peter J Olver Department of

More information

MB-JASS D Reconstruction. Hannes Hofmann March 2006

MB-JASS D Reconstruction. Hannes Hofmann March 2006 MB-JASS 2006 2-D Reconstruction Hannes Hofmann 19 29 March 2006 Outline Projections Radon Transform Fourier-Slice-Theorem Filtered Backprojection Ramp Filter 19 29 Mar 2006 Hannes Hofmann 2 Parallel Projections

More information

Kalman filtering based probabilistic nowcasting of object oriented tracked convective storms

Kalman filtering based probabilistic nowcasting of object oriented tracked convective storms Kalman iltering based probabilistic nowcasting o object oriented traced convective storms Pea J. Rossi,3, V. Chandrasear,2, Vesa Hasu 3 Finnish Meteorological Institute, Finland, Eri Palménin Auio, pea.rossi@mi.i

More information

Metrology and Sensing

Metrology and Sensing Metrology and Sensing Lecture 5: Interferometry I 08--6 Herbert Gross Winter term 08 www.iap.uni-jena.de Schedule Optical Metrology and Sensing 08 No Date Subject Detailed Content 6.0. Introduction Introduction,

More information

The Poisson summation formula, the sampling theorem, and Dirac combs

The Poisson summation formula, the sampling theorem, and Dirac combs The Poisson summation ormula, the sampling theorem, and Dirac combs Jordan Bell jordanbell@gmailcom Department o Mathematics, University o Toronto April 3, 24 Poisson summation ormula et S be the set o

More information

Metrology and Sensing

Metrology and Sensing Metrology and Sensing Lecture 5: Interferometry I 017-11-16 Herbert Gross Winter term 017 www.iap.uni-jena.de Preliminary Schedule No Date Subject Detailed Content 1 19.10. Introduction Introduction, optical

More information

Effective Fresnel-number concept for evaluating the relative focal shift in focused beams

Effective Fresnel-number concept for evaluating the relative focal shift in focused beams Martíne-Corral et al. Vol. 15, No. / February 1998/ J. Opt. Soc. Am. A 449 Eective Fresnel-number concept or evaluating the relative ocal shit in ocused beams Manuel Martíne-Corral, Carlos J. Zapata-Rodrígue,

More information

Products and Convolutions of Gaussian Probability Density Functions

Products and Convolutions of Gaussian Probability Density Functions Tina Memo No. 003-003 Internal Report Products and Convolutions o Gaussian Probability Density Functions P.A. Bromiley Last updated / 9 / 03 Imaging Science and Biomedical Engineering Division, Medical

More information

! Downsampling/Upsampling. ! Practical Interpolation. ! Non-integer Resampling. ! Multi-Rate Processing. " Interchanging Operations

! Downsampling/Upsampling. ! Practical Interpolation. ! Non-integer Resampling. ! Multi-Rate Processing.  Interchanging Operations Lecture Outline ESE 531: Digital Signal Processing Lec 10: February 14th, 2017 Practical and Non-integer Sampling, Multirate Sampling! Downsampling/! Practical Interpolation! Non-integer Resampling! Multi-Rate

More information

Optical determination of field angular correlation for transmission through three-dimensional turbid media

Optical determination of field angular correlation for transmission through three-dimensional turbid media 1040 J. Opt. Soc. Am. A/Vol. 16, No. 5/May 1999 Brian G. Hoover Optical determination o ield angular correlation or transmission through three-dimensional turbid media Brian G. Hoover Department o Electrical

More information

OBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL

OBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL OBSERVER/KALMAN AND SUBSPACE IDENTIFICATION OF THE UBC BENCHMARK STRUCTURAL MODEL Dionisio Bernal, Burcu Gunes Associate Proessor, Graduate Student Department o Civil and Environmental Engineering, 7 Snell

More information

Using a Membrane DM to Generate Zernike Modes

Using a Membrane DM to Generate Zernike Modes Using a Membrane DM to Generate Zernike Modes Author: Justin D. Mansell, Ph.D. Active Optical Systems, LLC Revision: 12/23/08 Membrane DMs have been used quite extensively to impose a known phase onto

More information

Fourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012)

Fourier Transform in Image Processing. CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) Fourier Transform in Image Processing CS/BIOEN 6640 U of Utah Guido Gerig (slides modified from Marcel Prastawa 2012) Basis Decomposition Write a function as a weighted sum of basis functions f ( x) wibi(

More information

Advanced Signal Processing Techniques for Feature Extraction in Data Mining

Advanced Signal Processing Techniques for Feature Extraction in Data Mining Volume 9 o.9, April 0 Advanced Signal Processing Techniques or Feature Extraction in Data Mining Maya ayak Proessor Orissa Engineering College BPUT, Bhubaneswar Bhawani Sankar Panigrahi Asst. Proessor

More information

Introduction to time-frequency analysis Centre for Doctoral Training in Healthcare Innovation

Introduction to time-frequency analysis Centre for Doctoral Training in Healthcare Innovation Introduction to time-frequency analysis Centre for Doctoral Training in Healthcare Innovation Dr. Gari D. Clifford, University Lecturer & Director, Centre for Doctoral Training in Healthcare Innovation,

More information

Fourier Sampling. Fourier Sampling. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2005 Linear Systems Lecture 3.

Fourier Sampling. Fourier Sampling. Bioengineering 280A Principles of Biomedical Imaging. Fall Quarter 2005 Linear Systems Lecture 3. Bioengineering 280A Principles of Biomedical Imaging Fall Quarter 2005 Linear Sstems Lecture 3 Fourier Sampling F Instead of sampling the signal, we sample its Fourier Transform Sample??? F -1 Fourier

More information

An Alternative Poincaré Section for Steady-State Responses and Bifurcations of a Duffing-Van der Pol Oscillator

An Alternative Poincaré Section for Steady-State Responses and Bifurcations of a Duffing-Van der Pol Oscillator An Alternative Poincaré Section or Steady-State Responses and Biurcations o a Duing-Van der Pol Oscillator Jang-Der Jeng, Yuan Kang *, Yeon-Pun Chang Department o Mechanical Engineering, National United

More information

DESIGN OF MULTI-DIMENSIONAL DERIVATIVE FILTERS. Eero P. Simoncelli

DESIGN OF MULTI-DIMENSIONAL DERIVATIVE FILTERS. Eero P. Simoncelli Published in: First IEEE Int l Conf on Image Processing, Austin Texas, vol I, pages 790--793, November 1994. DESIGN OF MULTI-DIMENSIONAL DERIVATIVE FILTERS Eero P. Simoncelli GRASP Laboratory, Room 335C

More information

References Ideal Nyquist Channel and Raised Cosine Spectrum Chapter 4.5, 4.11, S. Haykin, Communication Systems, Wiley.

References Ideal Nyquist Channel and Raised Cosine Spectrum Chapter 4.5, 4.11, S. Haykin, Communication Systems, Wiley. Baseand Data Transmission III Reerences Ideal yquist Channel and Raised Cosine Spectrum Chapter 4.5, 4., S. Haykin, Communication Systems, iley. Equalization Chapter 9., F. G. Stremler, Communication Systems,

More information

A STUDY OF SHOСK WAVE PROCESSES IN THE COMBUSTION CHAMBER AND THE ESTIMATION METHOD OF THE KNOCK INTENSITY BASED ON ION CURRENT SIGNAL ANALYSIS

A STUDY OF SHOСK WAVE PROCESSES IN THE COMBUSTION CHAMBER AND THE ESTIMATION METHOD OF THE KNOCK INTENSITY BASED ON ION CURRENT SIGNAL ANALYSIS A STUDY OF SHOСK WAVE PROCESSES IN THE COMBUSTION CHAMBER AND THE ESTIMATION METHOD OF THE KNOCK INTENSITY BASED ON ION CURRENT SIGNAL ANALYSIS Artem Yurevich Budko, Anatolij Sergeevich Nazarkin and Mikhail

More information

Computer-Generated Holographic Gratings in Soft Matter

Computer-Generated Holographic Gratings in Soft Matter Computer-Generated Holographic Gratings in Sot Matter Gianluigi Zito 1,2, Antigone Marino 1,3, Bruno Piccirillo 1,2, Volodymyr Tkachenko 3, Enrico Santamato 1,2, Giancarlo Abbate 1,3,4 1 University Federico

More information

Spectroscopic Instruments

Spectroscopic Instruments Spectroscopic Instruments 95 Spectroscopic Instruments by division of amplitude Mach-Zehnder (division of amplitude) Michelson Fringe localisation LIGO Fabry-Perot (FPI) Multi-layer coatings 96 Mach-Zehnder

More information

Topic 2: Imaging Process

Topic 2: Imaging Process A/TP T V Topic 2: maging Process Aim This lecture covers the parts of optical image formation required for this course. Contents: 1. maging as a Linear Process 2. P of ideal system 3. T of ideal system

More information

Algorithm Design and Analysis

Algorithm Design and Analysis Algorithm Design and Analysis LECTURE 14 Divide and Conquer Fast Fourier Transform Sofya Raskhodnikova 10/7/2016 S. Raskhodnikova; based on slides by K. Wayne. 5.6 Convolution and FFT Fast Fourier Transform:

More information