Today s lecture. Local neighbourhood processing. The convolution. Removing uncorrelated noise from an image The Fourier transform
|
|
- Derick Brooks
- 5 years ago
- Views:
Transcription
1 Cris Luengo TD396 fall 4 cris@cbuuse Today s lecture Local neighbourhood processing smoothing an image sharpening an image The convolution What is it? What is it useful for? How can I compute it? Removing uncorrelated noise from an image The Fourier transform What is it? What is it useful for?
2 Cris Luengo TD396 fall 4 cris@cbuuse Neighbourhoods
3 Cris Luengo TD396 fall 4 cris@cbuuse Local neighbourhood operation For each pixel, examine its neighbourhood and compute an output value (mean)
4 Cris Luengo TD396 fall 4 cris@cbuuse Local neighbourhood operation Possible operations to do for each neighbourhood: Neighbourhood size and shape is very important average (mean, median, etc) weighted average other statistics (variance, maximum, etc) difference (to compute derivative) round neighbourhood gives rotation invariance Adaptive filtering: changing size, shape and/or operation depending on local image properties
5 Cris Luengo TD396 fall 4 cris@cbuuse Smoothing an image Input image mean filter median filter
6 Cris Luengo TD396 fall 4 cris@cbuuse Smoothing an image Input image mean filter weighted mean filter
7 Cris Luengo TD396 fall 4 cris@cbuuse How to define Gaussian weights σ determines the amount of smoothing the neighbourhood size should be large enough to contain the whole Gaussian bell! rule of thumb: ceil(3σ) + sum of all weights normalised to x + y exp πσ σ ( ) ceil(3σ) +
8 Cris Luengo TD396 fall 4 cris@cbuuse Weighted mean filter For each pixel, multiply the values in its neighbourhood with the corresponding weights, then sum /9 /9 /9 /9 /9 / /9 /9 /9
9 Cris Luengo TD396 fall 4 cris@cbuuse Applications? Write down as many applications of a smoothing filter as you can come up with
10 Cris Luengo TD396 fall 4 cris@cbuuse Application: noise reduction input image Normally distributed noise Salt & pepper noise 3x3 mean filter 3x3 median filter
11 Cris Luengo TD396 fall 4 cris@cbuuse Application: abstraction Sometimes you just don t want all those details
12 Cris Luengo TD396 fall 4 cris@cbuuse Application: shading correction Gaussian smoothing, σ = pixels
13 Cris Luengo TD396 fall 4 cris@cbuuse Sharpening an image Unsharp masking original smoothed (3x3) sharpened (α = 9) sharpened = (+α) original α smoothed
14 Cris Luengo TD396 fall 4 cris@cbuuse Sharpening an image sharpened = (+α) original α smoothed sharpened = original + α ( original smoothed ) original /9 /9 /9 smoothed /9 /9 /9-8 - /9 /9 / diff
15 Cris Luengo TD396 fall 4 cris@cbuuse Laplace filter Laplace operator: Δ= = + x y -4-8 sharpened = original + 9 ( original smoothed ) sharpened = original - Laplace
16 Cris Luengo TD396 fall 4 cris@cbuuse Sobel filter Approximates the first derivatives: - Sx, x y Sy
17 Cris Luengo TD396 fall 4 cris@cbuuse Detecting edges Approximates the gradient magnitude: ( + x y sqrt ( Sx^ + Sy^ ) ) ( )
18 Cris Luengo TD396 fall 4 cris@cbuuse Adaptive filtering Many non-linear filters are meant to reduce noise without blurring the edges One common technique is to adapt the kernel so that it does not extend across any edges The bilateral filter is the most common one input image median filter bilateral filter
19 Cris Luengo TD396 fall 4 cris@cbuuse Bilateral filter A new kernel is designed for each output pixel Kernel weights are reduced if the corresponding pixel in the input image has a large difference in intensity with the central pixel h x ( x ) = Gσ ( x x ) Gσ ( f ( x ) f ( x )) x f
20 Cris Luengo TD396 fall 4 cris@cbuuse What happens at the image edge?
21 Cris Luengo TD396 fall 4 cris@cbuuse What happens at the image edge? Write down as many different ways of extending the edge as you can think of
22 Cris Luengo TD396 fall 4 cris@cbuuse What happens at the image edge? Mean padding f[end+x] = mean(f) Zero order hold f[end+x] = f[end]
23 Cris Luengo TD396 fall 4 cris@cbuuse What happens at the image edge? Periodic boundary condition f[end+x] = f[x] Symmetric boundary condition f[end+x] = f[end-x]
24 Cris Luengo TD396 fall 4 cris@cbuuse Linear neighbourhood operation For each pixel, multiply the values in its neighbourhood with the corresponding weights, then sum /9 /9 /9 /9 /9 / /9 /9 /9
25 Cris Luengo TD396 fall 4 cris@cbuuse Linear neighbourhood operation For each pixel, multiply the values in its neighbourhood with the corresponding weights, then sum (-,-) (,-) (,-) f(x,y) h(i,j) g(x,y) (-,) (,) (,) (-,) (,) (,) (x,y) (x,y)
26 Cris Luengo TD396 fall 4 cris@cbuuse Convolution h is: impulse response function point-spread function convolution kernel g (t ) = f (t ) h(t ) g (t ) = f (t τ) h( τ) d τ b g [n] = f [n k ] h[k ] k =a [a,b] is the interval where h is defined, eg [-,]
27 Cris Luengo TD396 fall 4 cris@cbuuse Convolution properties Linear: Scaling invariant: C f h = C f h Distributive: f g h = f h g h Time Invariant: shift f h = shift f h Commutative: f h = h f Associative: f h h = f h h (= shift invariant)
28 Cris Luengo TD396 fall 4 cris@cbuuse Associativity of convolution f (h h ) = (f h ) h if h = h h then f h = (f h ) h thus: you can decompose h to speed up the operation! Eg the Gaussian can be decomposed into two one-dimensional filters: G( x, y ) = e π σ x + y σ = e π σ x σ e π σ y σ
29 Cris Luengo TD396 fall 4 cris@cbuuse Kernel decomposition G = G x G y original convolved with Gx Gx and Gy are both a kernel with 3x values G is a kernel with 3x3 values convolved with Gy 3+3 = 6 MADs 3x3 = 96 MADs
30 Cris Luengo TD396 fall 4 cris@cbuuse Sequence of filters f (h h h3 ) = (((f h ) h ) h3 ) 3+3 ops 4(3+3) ops = 4 ops 9+9 ops = 8 ops
31 Cris Luengo TD396 fall 4 cris@cbuuse Sequence of filters 3+3 ops 4(3+3) ops = 4 ops 9+9 ops = 8 ops
32 Cris Luengo TD396 fall 4 cris@cbuuse Sequence of filters a b e f (example from Section 37) b = Laplace(a) c=a-b d = Sobel(a) e = smooth(d) f=c*e g=a+f h = gamma(g) c d g h
33 Cris Luengo TD396 fall 4 cris@cbuuse Jean Baptiste Joseph Fourier Born March 768, Auxerre (Bourgogne region) Died 6 May 83, Paris Same age as Napoleon Bonaparte Permanent Secretary of the French Academy of Sciences (8-83) Foreign member of the Royal Swedish Academy of Sciences (83)
34 Cris Luengo TD396 fall 4 cris@cbuuse The Fourier transform = + + +
35 Cris Luengo TD396 fall 4 cris@cbuuse The Fourier transform Can you think of a function that cannot be decomposed into Fourier basis functions?
36 Cris Luengo TD396 fall 4 cris@cbuuse Fourier analysis
37 Cris Luengo TD396 fall 4 cris@cbuuse Fourier example
38 Cris Luengo TD396 fall 4 cris@cbuuse Fourier example
39 Cris Luengo TD396 fall 4 cris@cbuuse Fourier example
40 Cris Luengo TD396 fall 4 cris@cbuuse Fourier example
41 Cris Luengo TD396 fall 4 cris@cbuuse Complex numbers i= i i= x =a+ i b x * =a i b (complex conjugate) * (Euler s formula) x x =a + b = x a= x cos( x ) b x=arctan( ) a b= x sin( x ) e i φ=cos φ+ i sin φ x = x cos( x ) + i x sin( x ) = x e i x
42 Cris Luengo TD396 fall 4 cris@cbuuse Fourier basis function e iωx = cos(ω x )+ i sin(ω x ) ω= π f = T π T
43 Cris Luengo TD396 fall 4 cris@cbuuse Fourier transform F (ω) = f ( x ) e i ω x d x f (x) = = F (ω )e + i ω x F (ω )e + i ω x F (ω3 )e + i ω3 x + + +
44 Cris Luengo TD396 fall 4 cris@cbuuse Fourier transform F (ω) = complex value f ( x ) e i ω x d x f (x) = complex function F (ω )e + i ω x F (ω )e + i ω x F (ω3 )e + i ω3 x complex function???
45 Cris Luengo TD396 fall 4 cris@cbuuse Fourier basis function A ei ω x + A* e i ω x is a real-valued function Thus: we need negative frequencies! For real-valued images: At frequency ω we have weight A At frequency -ω we have weight A* F ( ω) = F * (ω)
46 Cris Luengo TD396 fall 4 cris@cbuuse Inverse Fourier transform f (x) = π F (ω) e iωx dω normalization no minus sign Compare with the forward transform: F (ω) = i ω x f (x) e dx
47 Cris Luengo TD396 fall 4 cris@cbuuse Fourier transform pairs Spatial cosine impulses sine impulses box sinc sinc box Gaussian white noise Notice the symmetry! Gaussian white noise Frequency impulse
48 Cris Luengo TD396 fall 4 cris@cbuuse Properties of the Fourier transform Spatial scaling Linear Amplitude scaling Addition Translation Convolution phase change ℱ {f h} = ℱ {f } ℱ {h} = ℱ {f } ℱ {h} ℱ {f h}
49 Cris Luengo TD396 fall 4 cris@cbuuse Summary of today s lecture Virtually all filtering is a local neighbourhood operation Convolution = linear and shift-invariant filters Many non-linear filters exist also eg mean filter, Gaussian weighted filter kernel can sometimes be decomposed eg median filter, bilateral filter The Fourier transform decomposes a function (image) into trigonometric basis functions (sines & cosines) The Fourier transform is used to analyse frequency components of an image
50 Cris Luengo TD396 fall 4 cris@cbuuse Reading assignment Filtering The Fourier transform Sections 34, 3, 36, 37, 4, 3 Section 4 Exercises: 3, 36, 39, 33, 37, 38, 43
Today s lecture. The Fourier transform. Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm
Today s lecture The Fourier transform What is it? What is it useful for? What are its properties? Sampling, aliasing, interpolation The Fast Fourier Transform (FFT) algorithm Jean Baptiste Joseph Fourier
More informationReading. 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 informationLocal enhancement. Local Enhancement. Local histogram equalized. Histogram equalized. Local Contrast Enhancement. Fig 3.23: Another example
Local enhancement Local Enhancement Median filtering Local Enhancement Sometimes Local Enhancement is Preferred. Malab: BlkProc operation for block processing. Left: original tire image. 0/07/00 Local
More informationLinear Operators and Fourier Transform
Linear Operators and Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013
More informationFiltering and Edge Detection
Filtering and Edge Detection Local Neighborhoods Hard to tell anything from a single pixel Example: you see a reddish pixel. Is this the object s color? Illumination? Noise? The next step in order of complexity
More informationLocal Enhancement. Local enhancement
Local Enhancement Local Enhancement Median filtering (see notes/slides, 3.5.2) HW4 due next Wednesday Required Reading: Sections 3.3, 3.4, 3.5, 3.6, 3.7 Local Enhancement 1 Local enhancement Sometimes
More informationCOMP344 Digital Image Processing Fall 2007 Final Examination
COMP344 Digital Image Processing Fall 2007 Final Examination Time allowed: 2 hours Name Student ID Email Question 1 Question 2 Question 3 Question 4 Question 5 Question 6 Total With model answer HK University
More informationReview Smoothing Spatial Filters Sharpening Spatial Filters. Spatial Filtering. Dr. Praveen Sankaran. Department of ECE NIT Calicut.
Spatial Filtering Dr. Praveen Sankaran Department of ECE NIT Calicut January 7, 203 Outline 2 Linear Nonlinear 3 Spatial Domain Refers to the image plane itself. Direct manipulation of image pixels. Figure:
More informationEmpirical Mean and Variance!
Global Image Properties! Global image properties refer to an image as a whole rather than components. Computation of global image properties is often required for image enhancement, preceding image analysis.!
More informationECE Digital Image Processing and Introduction to Computer Vision. Outline
ECE592-064 Digital mage Processing and ntroduction to Computer Vision Depart. of ECE, NC State University nstructor: Tianfu (Matt) Wu Spring 2017 1. Recap Outline 2. Thinking in the frequency domain Convolution
More informationLecture # 06. Image Processing in Frequency Domain
Digital Image Processing CP-7008 Lecture # 06 Image Processing in Frequency Domain Fall 2011 Outline Fourier Transform Relationship with Image Processing CP-7008: Digital Image Processing Lecture # 6 2
More informationImage Filtering, Edges and Image Representation
Image Filtering, Edges and Image Representation Capturing what s important Req reading: Chapter 7, 9 F&P Adelson, Simoncelli and Freeman (handout online) Opt reading: Horn 7 & 8 FP 8 February 19, 8 A nice
More informationLecture 4 Filtering in the Frequency Domain. Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2016
Lecture 4 Filtering in the Frequency Domain Lin ZHANG, PhD School of Software Engineering Tongji University Spring 2016 Outline Background From Fourier series to Fourier transform Properties of the Fourier
More informationThe Frequency Domain, without tears. Many slides borrowed from Steve Seitz
The Frequency Domain, without tears Many slides borrowed from Steve Seitz Somewhere in Cinque Terre, May 2005 CS194: Image Manipulation & Computational Photography Alexei Efros, UC Berkeley, Fall 2016
More informationLecture 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 informationFiltering in the Frequency Domain
Filtering in the Frequency Domain Outline Fourier Transform Filtering in Fourier Transform Domain 2/20/2014 2 Fourier Series and Fourier Transform: History Jean Baptiste Joseph Fourier, French mathematician
More informationTRACKING 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 informationComputer 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 informationThe Frequency Domain : Computational Photography Alexei Efros, CMU, Fall Many slides borrowed from Steve Seitz
The Frequency Domain 15-463: Computational Photography Alexei Efros, CMU, Fall 2008 Somewhere in Cinque Terre, May 2005 Many slides borrowed from Steve Seitz Salvador Dali Gala Contemplating the Mediterranean
More informationRecap of Monday. Linear filtering. Be aware of details for filter size, extrapolation, cropping
Recap of Monday Linear filtering h[ m, n] k, l f [ k, l] I[ m Not a matrix multiplication Sum over Hadamard product k, n l] Can smooth, sharpen, translate (among many other uses) 1 1 1 1 1 1 1 1 1 Be aware
More informationFILTERING IN THE FREQUENCY DOMAIN
1 FILTERING IN THE FREQUENCY DOMAIN Lecture 4 Spatial Vs Frequency domain 2 Spatial Domain (I) Normal image space Changes in pixel positions correspond to changes in the scene Distances in I correspond
More informationECE Digital Image Processing and Introduction to Computer Vision. Outline
2/9/7 ECE592-064 Digital Image Processing and Introduction to Computer Vision Depart. of ECE, NC State University Instructor: Tianfu (Matt) Wu Spring 207. Recap Outline 2. Sharpening Filtering Illustration
More informationDigital Image Processing COSC 6380/4393
Digital Image Processing COSC 6380/4393 Lecture 11 Oct 3 rd, 2017 Pranav Mantini Slides from Dr. Shishir K Shah, and Frank Liu Review: 2D Discrete Fourier Transform If I is an image of size N then Sin
More informationTheory of signals and images I. Dr. Victor Castaneda
Theory of signals and images I Dr. Victor Castaneda Image as a function Think of an image as a function, f, f: R 2 R I=f(x, y) gives the intensity at position (x, y) The image only is defined over a rectangle,
More informationFiltering, Frequency, and Edges
CS450: Introduction to Computer Vision Filtering, Frequency, and Edges Various slides from previous courses by: D.A. Forsyth (Berkeley / UIUC), I. Kokkinos (Ecole Centrale / UCL). S. Lazebnik (UNC / UIUC),
More informationComputer Vision. Filtering in the Frequency Domain
Computer Vision Filtering in the Frequency Domain Filippo Bergamasco (filippo.bergamasco@unive.it) http://www.dais.unive.it/~bergamasco DAIS, Ca Foscari University of Venice Academic year 2016/2017 Introduction
More informationDiscrete Fourier Transform
Discrete Fourier Transform DD2423 Image Analysis and Computer Vision Mårten Björkman Computational Vision and Active Perception School of Computer Science and Communication November 13, 2013 Mårten Björkman
More informationConvolution and Linear Systems
CS 450: Introduction to Digital Signal and Image Processing Bryan Morse BYU Computer Science Introduction Analyzing Systems Goal: analyze a device that turns one signal into another. Notation: f (t) g(t)
More informationImage preprocessing in spatial domain
Image preprocessing in spatial domain Sharpening, image derivatives, Laplacian, edges Revision: 1.2, dated: May 25, 2007 Tomáš Svoboda Czech Technical University, Faculty of Electrical Engineering Center
More informationTaking derivative by convolution
Taking derivative by convolution Partial derivatives with convolution For 2D function f(x,y), the partial derivative is: For discrete data, we can approximate using finite differences: To implement above
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Filtering in the Frequency Domain http://www.ee.unlv.edu/~b1morris/ecg782/ 2 Outline Background
More informationEdge 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 informationSlow mo guys Saccades. https://youtu.be/fmg9zohesgq?t=4s
Slow mo guys Saccades https://youtu.be/fmg9zohesgq?t=4s Thinking in Frequency Computer Vision James Hays Slides: Hoiem, Efros, and others Recap of Wednesday Linear filtering is dot product at each position
More informationIntroduction 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 informationDigital Image Processing. Image Enhancement: Filtering in the Frequency Domain
Digital Image Processing Image Enhancement: Filtering in the Frequency Domain 2 Contents In this lecture we will look at image enhancement in the frequency domain Jean Baptiste Joseph Fourier The Fourier
More informationLecture 7: Edge Detection
#1 Lecture 7: Edge Detection Saad J Bedros sbedros@umn.edu Review From Last Lecture Definition of an Edge First Order Derivative Approximation as Edge Detector #2 This Lecture Examples of Edge Detection
More informationFourier 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 informationHistogram Processing
Histogram Processing The histogram of a digital image with gray levels in the range [0,L-] is a discrete function h ( r k ) = n k where r k n k = k th gray level = number of pixels in the image having
More informationDISCRETE FOURIER TRANSFORM
DD2423 Image Processing and Computer Vision DISCRETE FOURIER TRANSFORM Mårten Björkman Computer Vision and Active Perception School of Computer Science and Communication November 1, 2012 1 Terminology:
More informationImage Gradients and Gradient Filtering Computer Vision
Image Gradients and Gradient Filtering 16-385 Computer Vision What is an image edge? Recall that an image is a 2D function f(x) edge edge How would you detect an edge? What kinds of filter would you use?
More informationDigital Image Processing COSC 6380/4393
Digital Image Processing COSC 6380/4393 Lecture 13 Oct 2 nd, 2018 Pranav Mantini Slides from Dr. Shishir K Shah, and Frank Liu Review f 0 0 0 1 0 0 0 0 w 1 2 3 2 8 Zero Padding 0 0 0 0 0 0 0 1 0 0 0 0
More informationOutline. 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 informationI Chen Lin, Assistant Professor Dept. of CS, National Chiao Tung University. Computer Vision: 4. Filtering
I Chen Lin, Assistant Professor Dept. of CS, National Chiao Tung University Computer Vision: 4. Filtering Outline Impulse response and convolution. Linear filter and image pyramid. Textbook: David A. Forsyth
More informationCS 4495 Computer Vision. Frequency and Fourier Transforms. Aaron Bobick School of Interactive Computing. Frequency and Fourier Transform
CS 4495 Computer Vision Frequency and Fourier Transforms Aaron Bobick School of Interactive Computing Administrivia Project 1 is (still) on line get started now! Readings for this week: FP Chapter 4 (which
More informationThe Frequency Domain. Many slides borrowed from Steve Seitz
The Frequency Domain Many slides borrowed from Steve Seitz Somewhere in Cinque Terre, May 2005 15-463: Computational Photography Alexei Efros, CMU, Spring 2010 Salvador Dali Gala Contemplating the Mediterranean
More informationIMAGE ENHANCEMENT: FILTERING IN THE FREQUENCY DOMAIN. Francesca Pizzorni Ferrarese
IMAGE ENHANCEMENT: FILTERING IN THE FREQUENCY DOMAIN Francesca Pizzorni Ferrarese Contents In this lecture we will look at image enhancement in the frequency domain Jean Baptiste Joseph Fourier The Fourier
More informationComputer Vision Lecture 3
Computer Vision Lecture 3 Linear Filters 03.11.2015 Bastian Leibe RWTH Aachen http://www.vision.rwth-aachen.de leibe@vision.rwth-aachen.de Demo Haribo Classification Code available on the class website...
More informationWhy does a lower resolution image still make sense to us? What do we lose? Image:
2D FREQUENCY DOMAIN The slides are from several sources through James Hays (Brown); Srinivasa Narasimhan (CMU); Silvio Savarese (U. of Michigan); Bill Freeman and Antonio Torralba (MIT), including their
More information3F1 Random Processes Examples Paper (for all 6 lectures)
3F Random Processes Examples Paper (for all 6 lectures). Three factories make the same electrical component. Factory A supplies half of the total number of components to the central depot, while factories
More informationAnnouncements. Filtering. Image Filtering. Linear Filters. Example: Smoothing by Averaging. Homework 2 is due Apr 26, 11:59 PM Reading:
Announcements Filtering Homework 2 is due Apr 26, :59 PM eading: Chapter 4: Linear Filters CSE 52 Lecture 6 mage Filtering nput Output Filter (From Bill Freeman) Example: Smoothing by Averaging Linear
More informationDiscrete Fourier Transform
Last lecture I introduced the idea that any function defined on x 0,..., N 1 could be written a sum of sines and cosines. There are two different reasons why this is useful. The first is a general one,
More informationGBS765 Electron microscopy
GBS765 Electron microscopy Lecture 1 Waves and Fourier transforms 10/14/14 9:05 AM Some fundamental concepts: Periodicity! If there is some a, for a function f(x), such that f(x) = f(x + na) then function
More informationCITS 4402 Computer Vision
CITS 4402 Computer Vision Prof Ajmal Mian Adj/A/Prof Mehdi Ravanbakhsh, CEO at Mapizy (www.mapizy.com) and InFarm (www.infarm.io) Lecture 04 Greyscale Image Analysis Lecture 03 Summary Images as 2-D signals
More informationColorado School of Mines Image and Multidimensional Signal Processing
Image and Multidimensional Signal Processing Professor William Hoff Department of Electrical Engineering and Computer Science Spatial Filtering Main idea Spatial filtering Define a neighborhood of a pixel
More informationImage Enhancement in the frequency domain. GZ Chapter 4
Image Enhancement in the frequency domain GZ Chapter 4 Contents In this lecture we will look at image enhancement in the frequency domain The Fourier series & the Fourier transform Image Processing in
More informationDigital Image Processing. Filtering in the Frequency Domain
2D Linear Systems 2D Fourier Transform and its Properties The Basics of Filtering in Frequency Domain Image Smoothing Image Sharpening Selective Filtering Implementation Tips 1 General Definition: System
More informationIntroduction to the Discrete Fourier Transform
Introduction to the Discrete ourier Transform Lucas J. van Vliet www.ph.tn.tudelft.nl/~lucas TNW: aculty of Applied Sciences IST: Imaging Science and technology PH: Linear Shift Invariant System A discrete
More informationEdge 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 informationSpatial 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 informationImage Enhancement: Methods. Digital Image Processing. No Explicit definition. Spatial Domain: Frequency Domain:
Image Enhancement: No Explicit definition Methods Spatial Domain: Linear Nonlinear Frequency Domain: Linear Nonlinear 1 Spatial Domain Process,, g x y T f x y 2 For 1 1 neighborhood: Contrast Enhancement/Stretching/Point
More informationEAS 305 Random Processes Viewgraph 1 of 10. Random Processes
EAS 305 Random Processes Viewgraph 1 of 10 Definitions: Random Processes A random process is a family of random variables indexed by a parameter t T, where T is called the index set λ i Experiment outcome
More informationFourier Transform. sin(n# x)), where! = 2" / L and
Fourier Transform Henning Stahlberg Introduction The tools provided by the Fourier transform are helpful for the analysis of 1D signals (time and frequency (or Fourier) domains), as well as 2D/3D signals
More informationG52IVG, School of Computer Science, University of Nottingham
Image Transforms Fourier Transform Basic idea 1 Image Transforms Fourier transform theory Let f(x) be a continuous function of a real variable x. The Fourier transform of f(x) is F ( u) f ( x)exp[ j2πux]
More informationLecture 1 January 5, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture January 5, 26 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long; Edited by E. Candes & E. Bates Outline
More informationCommunication Signals (Haykin Sec. 2.4 and Ziemer Sec Sec. 2.4) KECE321 Communication Systems I
Communication Signals (Haykin Sec..4 and iemer Sec...4-Sec..4) KECE3 Communication Systems I Lecture #3, March, 0 Prof. Young-Chai Ko 년 3 월 일일요일 Review Signal classification Phasor signal and spectra Representation
More informationITK 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 informationSEISMIC WAVE PROPAGATION. Lecture 2: Fourier Analysis
SEISMIC WAVE PROPAGATION Lecture 2: Fourier Analysis Fourier Series & Fourier Transforms Fourier Series Review of trigonometric identities Analysing the square wave Fourier Transform Transforms of some
More informationConvolution Spatial Aliasing Frequency domain filtering fundamentals Applications Image smoothing Image sharpening
Frequency Domain Filtering Correspondence between Spatial and Frequency Filtering Fourier Transform Brief Introduction Sampling Theory 2 D Discrete Fourier Transform Convolution Spatial Aliasing Frequency
More informationDigital 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 informationAnalysis II: Fourier Series
.... Analysis II: Fourier Series Kenichi Maruno Department of Mathematics, The University of Texas - Pan American May 3, 011 K.Maruno (UT-Pan American) Analysis II May 3, 011 1 / 16 Fourier series were
More informationECG782: Multidimensional Digital Signal Processing
Professor Brendan Morris, SEB 3216, brendan.morris@unlv.edu ECG782: Multidimensional Digital Signal Processing Spring 2014 TTh 14:30-15:45 CBC C313 Lecture 05 Image Processing Basics 13/02/04 http://www.ee.unlv.edu/~b1morris/ecg782/
More informationImage Analysis. 3. Fourier Transform
Image Analysis Image Analysis 3. Fourier Transform Lars Schmidt-Thieme Information Systems and Machine Learning Lab (ISMLL) Institute for Business Economics and Information Systems & Institute for Computer
More informationPSET 0 Due Today by 11:59pm Any issues with submissions, post on Piazza. Fall 2018: T-R: Lopata 101. Median Filter / Order Statistics
CSE 559A: Computer Vision OFFICE HOURS This Friday (and this Friday only): Zhihao's Office Hours in Jolley 43 instead of 309. Monday Office Hours: 5:30-6:30pm, Collaboration Space @ Jolley 27. PSET 0 Due
More informationLecture 3: Linear Filters
Lecture 3: Linear Filters Professor Fei Fei Li Stanford Vision Lab 1 What we will learn today? Images as functions Linear systems (filters) Convolution and correlation Discrete Fourier Transform (DFT)
More informationEE 3054: Signals, Systems, and Transforms Summer It is observed of some continuous-time LTI system that the input signal.
EE 34: Signals, Systems, and Transforms Summer 7 Test No notes, closed book. Show your work. Simplify your answers. 3. It is observed of some continuous-time LTI system that the input signal = 3 u(t) produces
More information26. The Fourier Transform in optics
26. The Fourier Transform in optics What is the Fourier Transform? Anharmonic waves The spectrum of a light wave Fourier transform of an exponential The Dirac delta function The Fourier transform of e
More informationImage 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 informationThinking in Frequency
09/05/17 Thinking in Frequency Computational Photography University of Illinois Derek Hoiem Administrative Matlab/linear algebra tutorial tomorrow, planned for 6:30pm Probably 1214 DCL (will send confirmation
More informationLecture 3: Linear Filters
Lecture 3: Linear Filters Professor Fei Fei Li Stanford Vision Lab 1 What we will learn today? Images as functions Linear systems (filters) Convolution and correlation Discrete Fourier Transform (DFT)
More informationf (t) K(t, u) d t. f (t) K 1 (t, u) d u. Integral Transform Inverse Fourier Transform
Integral Transforms Massoud Malek An integral transform maps an equation from its original domain into another domain, where it might be manipulated and solved much more easily than in the original domain.
More informationVisual features: From Fourier to Gabor
Visual features: From Fourier to Gabor Deep Learning Summer School 2015, Montreal Hubel and Wiesel, 1959 from: Natural Image Statistics (Hyvarinen, Hurri, Hoyer; 2009) Alexnet ICA from: Natural Image Statistics
More informationTemplates, Image Pyramids, and Filter Banks
Templates, Image Pyramids, and Filter Banks 09/9/ Computer Vision James Hays, Brown Slides: Hoiem and others Review. Match the spatial domain image to the Fourier magnitude image 2 3 4 5 B A C D E Slide:
More informationLECTURE 12 Sections Introduction to the Fourier series of periodic signals
Signals and Systems I Wednesday, February 11, 29 LECURE 12 Sections 3.1-3.3 Introduction to the Fourier series of periodic signals Chapter 3: Fourier Series of periodic signals 3. Introduction 3.1 Historical
More informationCS711008Z Algorithm Design and Analysis
CS711008Z Algorithm Design and Analysis Lecture 5 FFT and Divide and Conquer Dongbo Bu Institute of Computing Technology Chinese Academy of Sciences, Beijing, China 1 / 56 Outline DFT: evaluate a polynomial
More informationFourier Transform and Frequency Domain
Fourier Transform and Frequency Domain http://www.cs.cmu.edu/~16385/ 16-385 Computer Vision Spring 2018, Lecture 3 (part 2) Overview of today s lecture Some history. Fourier series. Frequency domain. Fourier
More informationFourier Transforms 1D
Fourier Transforms 1D 3D Image Processing Alireza Ghane 1 Overview Recap Intuitions Function representations shift-invariant spaces linear, time-invariant (LTI) systems complex numbers Fourier Transforms
More informationImages have structure at various scales
Images have structure at various scales Frequency Frequency of a signal is how fast it changes Reflects scale of structure A combination of frequencies 0.1 X + 0.3 X + 0.5 X = Fourier transform Can we
More informationLecture 5. The Digital Fourier Transform. (Based, in part, on The Scientist and Engineer's Guide to Digital Signal Processing by Steven Smith)
Lecture 5 The Digital Fourier Transform (Based, in part, on The Scientist and Engineer's Guide to Digital Signal Processing by Steven Smith) 1 -. 8 -. 6 -. 4 -. 2-1 -. 8 -. 6 -. 4 -. 2 -. 2. 4. 6. 8 1
More informationReview of Linear System Theory
Review of Linear System Theory The following is a (very) brief review of linear system theory and Fourier analysis. I work primarily with discrete signals. I assume the reader is familiar with linear algebra
More informationFiltering 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 informationIntroduction to Linear Systems
cfl David J Fleet, 998 Introduction to Linear Systems David Fleet For operator T, input I, and response R = T [I], T satisfies: ffl homogeniety: iff T [ai] = at[i] 8a 2 C ffl additivity: iff T [I + I 2
More informationPhysics 342 Lecture 2. Linear Algebra I. Lecture 2. Physics 342 Quantum Mechanics I
Physics 342 Lecture 2 Linear Algebra I Lecture 2 Physics 342 Quantum Mechanics I Wednesday, January 3th, 28 From separation of variables, we move to linear algebra Roughly speaking, this is the study of
More informationComputational Photography
Computational Photography Si Lu Spring 208 http://web.cecs.pdx.edu/~lusi/cs50/cs50_computati onal_photography.htm 04/0/208 Last Time o Digital Camera History of Camera Controlling Camera o Photography
More informationSignals and Systems Lecture (S2) Orthogonal Functions and Fourier Series March 17, 2008
Signals and Systems Lecture (S) Orthogonal Functions and Fourier Series March 17, 008 Today s Topics 1. Analogy between functions of time and vectors. Fourier series Take Away Periodic complex exponentials
More informationPhysics 342 Lecture 2. Linear Algebra I. Lecture 2. Physics 342 Quantum Mechanics I
Physics 342 Lecture 2 Linear Algebra I Lecture 2 Physics 342 Quantum Mechanics I Wednesday, January 27th, 21 From separation of variables, we move to linear algebra Roughly speaking, this is the study
More informationContinuous-time Fourier Methods
ELEC 321-001 SIGNALS and SYSTEMS Continuous-time Fourier Methods Chapter 6 1 Representing a Signal The convolution method for finding the response of a system to an excitation takes advantage of the linearity
More informationAdditional 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 informationGaussian derivatives
Gaussian derivatives UCU Winter School 2017 James Pritts Czech Tecnical University January 16, 2017 1 Images taken from Noah Snavely s and Robert Collins s course notes Definition An image (grayscale)
More informationFourier Transform and its Applications
Fourier Transform and its Applications Prof. (Dr.) K.R. Chowdhary, Director COE Email: kr.chowdhary@iitj.ac.in webpage: http://www.krchowdhary.com JIET College of Engineering August 18, 2017 kr chowdhary
More informationFundamental Solutions and Green s functions. Simulation Methods in Acoustics
Fundamental Solutions and Green s functions Simulation Methods in Acoustics Definitions Fundamental solution The solution F (x, x 0 ) of the linear PDE L {F (x, x 0 )} = δ(x x 0 ) x R d Is called the fundamental
More information