CITS 4402 Computer Vision

Size: px
Start display at page:

Download "CITS 4402 Computer Vision"

Transcription

1 CITS 4402 Computer Vision Prof Ajmal Mian Adj/A/Prof Mehdi Ravanbakhsh, CEO at Mapizy ( and InFarm ( Lecture 04 Greyscale Image Analysis

2 Lecture 03 Summary Images as 2-D signals Linear and non-linear filters Fourier series Fourier Transform Discrete Fourier Transform 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 2

3 Image Enhancement Aim: Make images easier to interpret for the human eye Generate better input for other image processing techniques There are two main categories of techniques Spatial domain methods which operate directly on pixels Frequency domain methods which operate on the Fourier transform of the image 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 3

4 Single Pixel Manipulation The value of g(x, y) depends directly on the value of f x, y This is a greyscale transformation or simply mapping Common mapping functions: Identity transformation Thresholding 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 4

5 Common Mapping Functions lightening increase contrast 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 5

6 Common Mapping Functions enhance contrast in dark regions enhance contrast in light regions 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 6

7 (Hawkes Bay, NZ) 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 7

8 Sample MATLAB Code The Matlab code for enhancing the Hawkes Bay image on the previous slide is extremely simple! im = imread('unequalized_hawkes_bay_nz.jpg'); Imshow(im); % after obesrving the intensity histogram of the input image... f = 0:255; % input intensity g(1:100) = linspace(0,20,100); g(101:200) = linspace(21,240,100); g(201:256) = linspace(241,255,56); figure, plot(f,g,'b-'); % generate output image newim = g(im+1); newim = uint8(newim); % convert to uint8 figure, imshow(newim); 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 8

9 Histogram Equalization stretch out the histogram to produce a more uniform distribution Image data squashed into a small range of grey values Ideally we want the image data to spread uniformly over all grey values 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 9

10 Histogram Equalization For digital images, we have a discrete formulation. Let n k number of pixels with grey level k N total number of pixels Then the probability of obtaining grey level k in input image f is: P f f k = n k N The transformation is: g k = T f k = Note: the values of g k will have to be scaled up by 255 and rounded to the nearest integer. This discretization of g k means that the transformed image will not have a perfectly uniform distribution Indeed, the total number of distinct grey levels is reduced! k i=0 ni N 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 10

11 Histogram Equalization (cont.) An example: suppose that we have a image with 8 grey levels: f i n i n i /N g i (ideally) /7=0 1/ /7=0.14 2/ /7=0.28 3/ /7=0.42 4/ /7=0.57 5/ /7=0.71 6/ /7= /7=1 We want the cumulative distribution N=4096 of the output image to look as much like a 45 line as possible. 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 11

12 Histogram Equalization (cont.) f i n i n i /N g i (ideally) /7=0 1/ /7=0.14 2/ /7=0.28 3/ /7=0.42 4/ /7=0.57 5/ /7=0.71 6/ /7= /7=1 Thus intensity mapping is: InputIntensity OutputIntensity 0 1/7 1/7 3/7 2/7 5/7 3/7 6/7 4/7 6/7 5/7 1 6/ Input image has 8 grey levels Output image has only 5 grey levels 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 12

13 Histogram Equalization (cont.) f i n i n i /N g i n i of g i / / / / / / Thus intensity mapping is: InputIntensity OutputIntensity 0 1/7 1/7 3/7 2/7 5/7 3/7 6/7 4/7 6/7 5/7 1 6/ Input image has 8 grey levels Output image has only 5 grey levels 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 13

14 Histogram Equalization (cont.) f i n i n i /N g i n i of g i / / / / / / /19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 14

15 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 15

16 A very noisy image 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 16

17 Using Neighborhood Pixels So far we modified values of single pixels What if we want to take neighborhood into account? A common application of linear filtering is image smoothing using an averaging filter or averaging mask or averaging kernel Each point in the smoothed image g(x, y) is obtained from the average pixel value in a neighbourhood of (x, y) in the input image The averaging filter is also known as the box filter. Each pixel under the mask is multiplied by 1/9, summed, and the result is placed in the output image The mask is successively moved across the image. That is, we convolve the image with the mask /19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 17

18 Example of 2-D Convolution Convolution with box filters of size 1,5,15,3,9,35 (reading order) 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 18

19 Image filtering g.,. g.,. = 1 9 f [.,.] h[.,.] h[ m, n] g[ k, l] f [ m k, n l] k, l 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 19 Credit: S. Seitz

20 Image filtering g.,. g.,. = 1 9 f [.,.] h[.,.] h[ m, n] g[ k, l] f [ m k, n l] k, l 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 20 Credit: S. Seitz

21 Image filtering g.,. g.,. = 1 9 f [.,.] h[.,.] h[ m, n] g[ k, l] f [ m k, n l] k, l 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 21 Credit: S. Seitz

22 Image filtering g.,. g.,. = 1 9 f [.,.] h[.,.] h[ m, n] g[ k, l] f [ m k, n l] k, l 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 22 Credit: S. Seitz

23 Image filtering g.,. g.,. = 1 9 f [.,.] h[.,.] h[ m, n] g[ k, l] f [ m k, n l] k, l 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 23 Credit: S. Seitz

24 Image filtering g.,. g.,. = 1 9 f [.,.] h[.,.] ? h[ m, n] g[ k, l] f [ m k, n l] k, l 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 24 Credit: S. Seitz

25 Image filtering g.,. g.,. = 1 9 f [.,.] h[.,.] ? h[ m, n] g[ k, l] f [ m k, n l] k, l 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 25 Credit: S. Seitz

26 Image filtering g.,. g.,. = 1 9 f [.,.] h[.,.] h[ m, n] g[ k, l] f [ m k, n l] k, l 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 26 Credit: S. Seitz

27 Practice with Linear Filters ? Original 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 27

28 Practice with Linear Filters Original Filtered (no change) 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 28

29 Practice with Linear Filters ? Original 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 29

30 Practice with Linear Filters Original Shifted left By 1 pixel 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 30

31 Practice with Linear Filters ? Original (Note that filter sums to 1) 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 31

32 Practice with Linear Filters Original Sharpening filter - Accentuates differences with local average 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 32

33 Sharpening 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 33 Source: D. Lowe

34 Other Averaging Filters One expects the value of a pixel to be more closely related to the values of pixels close to it than to those further away Accordingly it is usual to weight the pixels near the centre of the mask more strongly than those at the edge 1? /19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 34

35 Other Averaging Filters One expects the value of a pixel to be more closely related to the values of pixels close to it than to those further away Accordingly it is usual to weight the pixels near the centre of the mask more strongly than those at the edge /19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 35

36 Important filter: Gaussian Gaussian Filter x 5, = 1 σ 2 is also known as the width of the kernel Use Matlab s fspecial function to create a Gaussian filter. 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis Slide credit: The University Christopher of Western Australia Rasmussen 36

37 Smoothing with box filter Example Box Filter Smoothing 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 37

38 Smoothing with Gaussian filter Example Gaussian Smoothing 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 38

39 Key Properties of Linear Filters Linearity: filter(f 1 + f 2 ) = filter(f 1 ) + filter(f 2 ) Shift invariance: same behavior regardless of pixel location filter(shift(f)) = shift(filter(f)) Any linear, shift-invariant operator can be represented as a convolution 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 39 Source: S. Lazebnik

40 Key Properties of Linear Filters Commutative: a * b = b * a Conceptually no difference between filter and signal Associative: a * (b * c) = (a * b) * c Often apply several filters one after another: (((a * b 1 ) * b 2 ) * b 3 ) This is equivalent to applying one filter: a * (b 1 * b 2 * b 3 ) Distributes over addition: a * (b + c) = (a * b) + (a * c) Scalars factor out: ka * b = a * kb = k (a * b) Identity: unit impulse e = [0, 0, 1, 0, 0], a * e = a 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 40 Source: S. Lazebnik

41 Gaussian Filters Linear filters Remove high-frequency components from the image (low-pass filter) Images become more smooth Convolution of a Gaussian with a Gaussian is another Gaussian So can smooth with small-width kernel, repeat, and get same result as larger-width kernel would have Convolving two times with Gaussian kernel of width σ is same as convolving once with kernel of width σ 2 Separable kernels Factors into product of two 1D Gaussians 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 41

42 Separability of the Gaussian filter Gaussian Filters 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 42

43 Practical Matters What about near the edge? The filter window falls off the edge of the image Need to extrapolate Methods: clip filter (black) wrap around copy edge reflect across edge 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 43

44 Q? Practical Matters Methods (MATLAB): clip filter (black): imfilter(f, g, 0) wrap around: imfilter(f, g, circular ) copy edge: imfilter(f, g, replicate ) reflect across edge: imfilter(f, g, symmetric ) 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 44

45 Practical Matters MATLAB: filter2(g, f, shape) shape = full : output size is sum of sizes of f and g shape = same : output size is same as f shape = valid : output size is difference of sizes of f and g g g g g g g f f f g g g g g full same valid g 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 45

46 Low-Pass Filtering Removing all high spatial frequencies from a signal to retain only low spatial frequencies is called low-pass filtering.? Old Spectrum New Spectrum Low-Pass Filtered Image 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 46

47 High-Pass Filtering Removing all low spatial frequencies from a signal to retain only high spatial frequencies is called high-pass filtering.? Old Spectrum New Spectrum High-Pass Filtered Image 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 47

48 Low-Pass Filtering An Example 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 48

49 Low-Pass Filtering An Example 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 49

50 High-Pass Filtering An Example 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 50

51 High-Pass Filtering An Example 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 51

52 Filtering in the Spatial Domain Low-pass filtering -> convolve the image with a box / Gaussian filter High-pass filtering ->? 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 52

53 Filtering in the Spatial Domain Low-pass filtering -> convolve the image with a box / Gaussian filter High-pass filtering -> -1/9-1/9-1/9-1/9 8/9-1/9-1/9-1/9-1/9 Since the sum of the weights is 0, the resulting signal will have a 0 DC value (i.e., the average value or the coefficient of the zero frequency component). To display the image you will need to take the absolute values. 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 53

54 Non-Linear Filtering Neighbourhood averaging or Gaussian smoothing will tend to blur edges because the high frequency components are attenuated An alternative approach is to use median filtering. Here we set the grey level to be the median of the pixel values in the neighbourhood Example: pixel values in 3 3 neighbourhood median value Sort the values Pixels with outlying values are forced to become more like their neighbours 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 54

55 Median Filtering An Example Median filter removes outliers Median filter smooths the image without blurring the edges 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 55

56 High-Boost Filtering Here we take the original image and boost the high frequency components. Can think of HighPass = Original LowPass. Thus HighBoost = b*original LowPass = (b-1)*original + Original LowPass = (b-1)*original + HighPass b is the boosting factor. When b=1, HighBoost = HighPass High-boost filtering useful for emphasizing high frequencies while retaining some low frequency components of the original image to aid in the interpretation of the image 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 56

57 High Boost Filtering (cont.) How can we perform high-boost filtering in the spatial domain? -1/9-1/9-1/9-1/9 w/9-1/9-1/9-1/9-1/9 where w = 9*b - 1 Original High Boosted (intermediate result) High Boosted (contrast stretching) 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 57

58 Homomorphic Filtering The brightness of an image point f(x, y) is a function of the illumination at that point and the reflectance of the object at that point, i.e., f x, y = i x, y r(x, y) illumination 0 i x, y < reflectance 0 r x, y 1 It is the reflectance that tells us information about the scene. We want to reduce the influence of illumination. Assumptions: 1. Illumination variations vary with low spatial frequency 2. Features of interest are the result of different reflectance properties of objects across the image and these vary with high(er) spatial frequency 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 58

59 Homomorphic Filtering (cont.) Let z x, y = log f x, y = log i x, y + log r x, y In the frequency domain, we have Z ω, ν = I ω, ν + R ω, ν Fourier transform of log(i x, y ) Fourier transform of log(r x, y ) If we apply a high boost filter, then we can suppress the low frequency illumination components and enhance the reflectance components Z ω, ν represents the Fourier Transform of the sum of two images: a low frequency illumination image a high frequency reflectance image 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 59

60 Homomorphic Filtering (cont.) Apply the high-boost filter H ω, ν : S ω, ν = Z ω, ν. H ω, ν Take the inverse FFT: s x, y = F 1 (S ω, ν ) Finally, exponentiate the result to account for taking the log of the original image: g x, y = exp(s x, y ) High-boost f(x, y) log FFT FFT filtering -1 exp g(x, y) 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 60

61 Homomorphic Filtering An Example Original image Filtered image 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 61

62 Smoothing and Sub-sampling In many Computer Vision applications, sub-sampling is often needed, e.g., to build an image pyramid, or simply to reduce the resolution for efficient storage, transmission, processing, 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 62

63 Throw away every other row and column to create a 1/2 size image 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 63

64 Sub-sampling Issues Sub-sampling may be dangerous. Why? Resample the checkerboard by taking one sample at each circle. In the case of the top left board, new representation is reasonable. Top right also yields a reasonable representation. Bottom left is all black (dubious) and bottom right has checks that are too big. 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 64

65 Aliasing in Videos 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 65

66 Aliasing in Graphics 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 66

67 Sampling and Aliasing Top row shows the images, sampled at every second pixel to get the next; bottom row shows the magnitude of frequency spectrum of these images. 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 67

68 Anti-aliasing 1. Sample more often 2. Get rid of all frequencies that are greater than half the new sampling frequency (Nyquist frequency) Will lose information But it s better than aliasing Apply a smoothing filter 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 68

69 Sampling and Aliasing Sampling with smoothing. Top row shows the images. We get the next image by smoothing the image with a Gaussian with sigma 1 pixel, then sampling at every second pixel to get the next; bottom row shows the magnitude spectrum of these images. 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 69

70 Subsampling without Pre-filtering 1/2 1/4 (2x zoom) 1/8 (4x zoom) 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 70

71 Subsampling with Gaussian Pre-filtering Gaussian 1/2 G 1/4 G 1/8 Slide by Steve Seitz 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 71

72 Image Pyramids Key component of many high level computer vision tasks How to create an image pyramid? Represent each image as a layer 1. Convolve layer i with a Gaussian filter 2. Subsample layer i by a factor of two (remove all even-numbered rows and columns) to get layer i + 1 Repeat steps 1 and 2 until a stopping criteria is satisfied 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 72

73 Image Pyramids 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 73

74 Template Matching Finding object of known shape and appearance in an image To identify the object, we have to compare the template image against the source image by sliding it. At each location, we find the matching score between the image and the template 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 74

75 Template Matching Scores Sum of Squared Difference (SSD) in pixel values (Normalized) Correlation coefficient (Normalized) Cross-correlation Template Source Image 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 75

76 Template Matching Scores SSD Normalized Correlation Coefficient Normalized Cross-Correlation 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 76

77 Summary Single Pixel Operations Histogram Equalization Filtering in the Spatial Domain Subsampling and Anti-aliasing Template Matching Acknowledgements: The slides are based on previous lectures by A/Prof Du Huynh and Prof Peter Koveski. Other material has been taken from Wikipedia, computer vision textbook by Forsyth & Ponce, and OpenCV documentation. 3/19/2018 Computer Vision - Lecture 04 - Greyscale Image Analysis The University of Western Australia 77

Computer Vision Lecture 3

Computer 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 information

Taking derivative by convolution

Taking 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 information

Lecture 3: Linear Filters

Lecture 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 information

Lecture 3: Linear Filters

Lecture 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 information

Templates, Image Pyramids, and Filter Banks

Templates, 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 information

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

Computer Vision Lecture 3

Computer Vision Lecture 3 Demo Haribo Classification Computer Vision Lecture 3 Linear Filters 23..24 Bastian Leibe RWTH Aachen http://www.vision.rwth-aachen.de leibe@vision.rwth-aachen.de Code available on the class website...

More information

CITS 4402 Computer Vision

CITS 4402 Computer Vision CITS 4402 Computer Vision A/Prof Ajmal Mian Adj/A/Prof Mehdi Ravanbakhsh Lecture 06 Object Recognition Objectives To understand the concept of image based object recognition To learn how to match images

More information

ECE Digital Image Processing and Introduction to Computer Vision. Outline

ECE 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 information

CEE598 - Visual Sensing for Civil Infrastructure Eng. & Mgmt.

CEE598 - Visual Sensing for Civil Infrastructure Eng. & Mgmt. CEE598 - Visual Sensing for Civil Infrastructure Eng. & Mgmt. Session 8- Linear Filters From Spatial Domain to Frequency Domain Mani Golparvar-Fard Department of Civil and Environmental Engineering 329D,

More information

Digital Image Processing. Filtering in the Frequency Domain

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

Thinking in Frequency

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

CSE 473/573 Computer Vision and Image Processing (CVIP)

CSE 473/573 Computer Vision and Image Processing (CVIP) CSE 473/573 Computer Vision and Image Processing (CVIP) Ifeoma Nwogu inwogu@buffalo.edu Lecture 11 Local Features 1 Schedule Last class We started local features Today More on local features Readings for

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

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

ECG782: Multidimensional Digital Signal Processing

ECG782: 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 information

ECE Digital Image Processing and Introduction to Computer Vision

ECE Digital Image Processing and Introduction to Computer Vision ECE592-064 Digital Image Processing and Introduction to Computer Vision Depart. of ECE, NC State University Instructor: Tianfu (Matt) Wu Spring 2017 Outline Recap, image degradation / restoration Template

More information

Frequency Filtering CSC 767

Frequency Filtering CSC 767 Frequency Filtering CSC 767 Outline Fourier transform and frequency domain Frequency view of filtering Hybrid images Sampling Slide: Hoiem Why does the Gaussian give a nice smooth image, but the square

More information

Subsampling and image pyramids

Subsampling and image pyramids Subsampling and image pyramids http://www.cs.cmu.edu/~16385/ 16-385 Computer Vision Spring 2018, Lecture 3 Course announcements Homework 0 and homework 1 will be posted tonight. - Homework 0 is not required

More information

Multiscale Image Transforms

Multiscale Image Transforms Multiscale Image Transforms Goal: Develop filter-based representations to decompose images into component parts, to extract features/structures of interest, and to attenuate noise. Motivation: extract

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

Announcements. Filtering. Image Filtering. Linear Filters. Example: Smoothing by Averaging. Homework 2 is due Apr 26, 11:59 PM Reading:

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

Recap of Monday. Linear filtering. Be aware of details for filter size, extrapolation, cropping

Recap 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 information

COMP344 Digital Image Processing Fall 2007 Final Examination

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

Advanced Features. Advanced Features: Topics. Jana Kosecka. Slides from: S. Thurn, D. Lowe, Forsyth and Ponce. Advanced features and feature matching

Advanced Features. Advanced Features: Topics. Jana Kosecka. Slides from: S. Thurn, D. Lowe, Forsyth and Ponce. Advanced features and feature matching Advanced Features Jana Kosecka Slides from: S. Thurn, D. Lowe, Forsyth and Ponce Advanced Features: Topics Advanced features and feature matching Template matching SIFT features Haar features 2 1 Features

More information

Linear Operators and Fourier Transform

Linear 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 information

Filtering and Edge Detection

Filtering 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 information

Lecture 7: Edge Detection

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

Machine vision. Summary # 4. The mask for Laplacian is given

Machine vision. Summary # 4. The mask for Laplacian is given 1 Machine vision Summary # 4 The mask for Laplacian is given L = 0 1 0 1 4 1 (6) 0 1 0 Another Laplacian mask that gives more importance to the center element is L = 1 1 1 1 8 1 (7) 1 1 1 Note that the

More information

Histogram Processing

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

Computer Vision & Digital Image Processing. Periodicity of the Fourier transform

Computer Vision & Digital Image Processing. Periodicity of the Fourier transform Computer Vision & Digital Image Processing Fourier Transform Properties, the Laplacian, Convolution and Correlation Dr. D. J. Jackson Lecture 9- Periodicity of the Fourier transform The discrete Fourier

More information

Image Enhancement in the frequency domain. GZ Chapter 4

Image 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 information

Machine vision, spring 2018 Summary 4

Machine vision, spring 2018 Summary 4 Machine vision Summary # 4 The mask for Laplacian is given L = 4 (6) Another Laplacian mask that gives more importance to the center element is given by L = 8 (7) Note that the sum of the elements in the

More information

Today s lecture. Local neighbourhood processing. The convolution. Removing uncorrelated noise from an image The Fourier transform

Today s lecture. Local neighbourhood processing. The convolution. Removing uncorrelated noise from an image The Fourier transform 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

More information

Computational Photography

Computational 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 information

Colorado School of Mines Image and Multidimensional Signal Processing

Colorado 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 information

Image Enhancement: Methods. Digital Image Processing. No Explicit definition. Spatial Domain: Frequency Domain:

Image 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 information

Vlad Estivill-Castro (2016) Robots for People --- A project for intelligent integrated systems

Vlad Estivill-Castro (2016) Robots for People --- A project for intelligent integrated systems 1 Vlad Estivill-Castro (2016) Robots for People --- A project for intelligent integrated systems V. Estivill-Castro 2 Perception Concepts Vision Chapter 4 (textbook) Sections 4.3 to 4.5 What is the course

More information

Image Processing /6.865 Frédo Durand A bunch of slides by Bill Freeman (MIT) & Alyosha Efros (CMU)

Image Processing /6.865 Frédo Durand A bunch of slides by Bill Freeman (MIT) & Alyosha Efros (CMU) Image Processing 6.815/6.865 Frédo Durand A bunch of slides by Bill Freeman (MIT) & Alyosha Efros (CMU) define cumulative histogram work on hist eq proof rearrange Fourier order discuss complex exponentials

More information

Digital Image Processing. Chapter 4: Image Enhancement in the Frequency Domain

Digital Image Processing. Chapter 4: Image Enhancement in the Frequency Domain Digital Image Processing Chapter 4: Image Enhancement in the Frequency Domain Image Enhancement in Frequency Domain Objective: To understand the Fourier Transform and frequency domain and how to apply

More information

Review Smoothing Spatial Filters Sharpening Spatial Filters. Spatial Filtering. Dr. Praveen Sankaran. Department of ECE NIT Calicut.

Review 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 information

G52IVG, School of Computer Science, University of Nottingham

G52IVG, 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 information

Filtering, Frequency, and Edges

Filtering, 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 information

Empirical Mean and Variance!

Empirical 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 information

Edges and Scale. Image Features. Detecting edges. Origin of Edges. Solution: smooth first. Effects of noise

Edges and Scale. Image Features. Detecting edges. Origin of Edges. Solution: smooth first. Effects of noise Edges and Scale Image Features From Sandlot Science Slides revised from S. Seitz, R. Szeliski, S. Lazebnik, etc. Origin of Edges surface normal discontinuity depth discontinuity surface color discontinuity

More information

The Frequency Domain, without tears. Many slides borrowed from Steve Seitz

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

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

EXAMINATION QUESTION PAPER

EXAMINATION QUESTION PAPER Faculty of Science and Technology EXAMINATION QUESTION PAPER Exam in: FYS-2010 Digital Image Processing Date: Monday 26 September 2016 Time: 09.00 13.00 Place: Approved aids: Administrasjonsbygget, Aud.Max.

More information

ECG782: Multidimensional Digital Signal Processing

ECG782: 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 information

Chapter 4: Filtering in the Frequency Domain. Fourier Analysis R. C. Gonzalez & R. E. Woods

Chapter 4: Filtering in the Frequency Domain. Fourier Analysis R. C. Gonzalez & R. E. Woods Fourier Analysis 1992 2008 R. C. Gonzalez & R. E. Woods Properties of δ (t) and (x) δ : f t) δ ( t t ) dt = f ( ) f x) δ ( x x ) = f ( ) ( 0 t0 x= ( 0 x0 1992 2008 R. C. Gonzalez & R. E. Woods Sampling

More information

Fourier Transforms 1D

Fourier 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 information

Chapter 4 Image Enhancement in the Frequency Domain

Chapter 4 Image Enhancement in the Frequency Domain Chapter 4 Image Enhancement in the Frequency Domain Yinghua He School of Computer Science and Technology Tianjin University Background Introduction to the Fourier Transform and the Frequency Domain Smoothing

More information

Edge Detection PSY 5018H: Math Models Hum Behavior, Prof. Paul Schrater, Spring 2005

Edge Detection PSY 5018H: Math Models Hum Behavior, Prof. Paul Schrater, Spring 2005 Edge Detection PSY 5018H: Math Models Hum Behavior, Prof. Paul Schrater, Spring 2005 Gradients and edges Points of sharp change in an image are interesting: change in reflectance change in object change

More information

Harris Corner Detector

Harris Corner Detector Multimedia Computing: Algorithms, Systems, and Applications: Feature Extraction By Dr. Yu Cao Department of Computer Science The University of Massachusetts Lowell Lowell, MA 01854, USA Part of the slides

More information

Digital Image Processing. Image Enhancement: Filtering in the Frequency Domain

Digital 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 information

Fourier series: Any periodic signals can be viewed as weighted sum. different frequencies. view frequency as an

Fourier series: Any periodic signals can be viewed as weighted sum. different frequencies. view frequency as an Image Enhancement in the Frequency Domain Fourier series: Any periodic signals can be viewed as weighted sum of sinusoidal signals with different frequencies Frequency Domain: view frequency as an independent

More information

Introduction to Computer Vision

Introduction to Computer Vision Introduction to Computer Vision Michael J. Black Sept 2009 Lecture 8: Pyramids and image derivatives Goals Images as functions Derivatives of images Edges and gradients Laplacian pyramids Code for lecture

More information

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

Lecture 6: Edge Detection. CAP 5415: Computer Vision Fall 2008

Lecture 6: Edge Detection. CAP 5415: Computer Vision Fall 2008 Lecture 6: Edge Detection CAP 5415: Computer Vision Fall 2008 Announcements PS 2 is available Please read it by Thursday During Thursday lecture, I will be going over it in some detail Monday - Computer

More information

IMAGE ENHANCEMENT II (CONVOLUTION)

IMAGE ENHANCEMENT II (CONVOLUTION) MOTIVATION Recorded images often exhibit problems such as: blurry noisy Image enhancement aims to improve visual quality Cosmetic processing Usually empirical techniques, with ad hoc parameters ( whatever

More information

Introduction to Computer Vision

Introduction to Computer Vision Introduction to Computer Vision Michael J. Black Oct. 2009 Lecture 10: Images as vectors. Appearance-based models. News Assignment 1 parts 3&4 extension. Due tomorrow, Tuesday, 10/6 at 11am. Goals Images

More information

Sampling in 1D ( ) Continuous time signal f(t) Discrete time signal. f(t) comb

Sampling in 1D ( ) Continuous time signal f(t) Discrete time signal. f(t) comb Sampling in 2D 1 Sampling in 1D Continuous time signal f(t) Discrete time signal t ( ) f [ k] = f( kt ) = f( t) δ t kt s k s f(t) comb k 2 Nyquist theorem (1D) At least 2 sample/period are needed to represent

More information

Discrete Fourier Transform

Discrete 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 information

Frequency2: Sampling and Aliasing

Frequency2: Sampling and Aliasing CS 4495 Computer Vision Frequency2: Sampling and Aliasing Aaron Bobick School of Interactive Computing Administrivia Project 1 is due tonight. Submit what you have at the deadline. Next problem set stereo

More information

Images have structure at various scales

Images 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 information

Digital Image Processing COSC 6380/4393

Digital 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 information

Blob Detection CSC 767

Blob Detection CSC 767 Blob Detection CSC 767 Blob detection Slides: S. Lazebnik Feature detection with scale selection We want to extract features with characteristic scale that is covariant with the image transformation Blob

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

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

Intensity Transformations and Spatial Filtering: WHICH ONE LOOKS BETTER? Intensity Transformations and Spatial Filtering: WHICH ONE LOOKS BETTER?

Intensity Transformations and Spatial Filtering: WHICH ONE LOOKS BETTER? Intensity Transformations and Spatial Filtering: WHICH ONE LOOKS BETTER? : WHICH ONE LOOKS BETTER? 3.1 : WHICH ONE LOOKS BETTER? 3.2 1 Goal: Image enhancement seeks to improve the visual appearance of an image, or convert it to a form suited for analysis by a human or a machine.

More information

Chapter 16. Local Operations

Chapter 16. Local Operations Chapter 16 Local Operations g[x, y] =O{f[x ± x, y ± y]} In many common image processing operations, the output pixel is a weighted combination of the gray values of pixels in the neighborhood of the input

More information

Corners, Blobs & Descriptors. With slides from S. Lazebnik & S. Seitz, D. Lowe, A. Efros

Corners, Blobs & Descriptors. With slides from S. Lazebnik & S. Seitz, D. Lowe, A. Efros Corners, Blobs & Descriptors With slides from S. Lazebnik & S. Seitz, D. Lowe, A. Efros Motivation: Build a Panorama M. Brown and D. G. Lowe. Recognising Panoramas. ICCV 2003 How do we build panorama?

More information

Why does a lower resolution image still make sense to us? What do we lose? Image:

Why 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 information

Digital Image Processing COSC 6380/4393

Digital 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 information

Generalized Laplacian as Focus Measure

Generalized Laplacian as Focus Measure Generalized Laplacian as Focus Measure Muhammad Riaz 1, Seungjin Park, Muhammad Bilal Ahmad 1, Waqas Rasheed 1, and Jongan Park 1 1 School of Information & Communications Engineering, Chosun University,

More information

Slow mo guys Saccades. https://youtu.be/fmg9zohesgq?t=4s

Slow 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 information

ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin

ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM. Dr. Lim Chee Chin ENT 315 Medical Signal Processing CHAPTER 2 DISCRETE FOURIER TRANSFORM Dr. Lim Chee Chin Outline Introduction Discrete Fourier Series Properties of Discrete Fourier Series Time domain aliasing due to frequency

More information

Multiresolution schemes

Multiresolution schemes Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Multi-dimensional signal processing Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione

More information

DISCRETE FOURIER TRANSFORM

DISCRETE 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 information

IMAGE ENHANCEMENT: FILTERING IN THE FREQUENCY DOMAIN. Francesca Pizzorni Ferrarese

IMAGE 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 information

Multiresolution schemes

Multiresolution schemes Multiresolution schemes Fondamenti di elaborazione del segnale multi-dimensionale Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Elaborazione dei Segnali Multi-dimensionali e

More information

Screen-space processing Further Graphics

Screen-space processing Further Graphics Screen-space processing Rafał Mantiuk Computer Laboratory, University of Cambridge Cornell Box and tone-mapping Rendering Photograph 2 Real-world scenes are more challenging } The match could not be achieved

More information

Convolution Spatial Aliasing Frequency domain filtering fundamentals Applications Image smoothing Image sharpening

Convolution 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 information

Review of Linear Systems Theory

Review of Linear Systems Theory Review of Linear Systems Theory The following is a (very) brief review of linear systems theory, convolution, and Fourier analysis. I work primarily with discrete signals, but each result developed in

More information

Review of Linear System Theory

Review 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 information

FILTERING IN THE FREQUENCY DOMAIN

FILTERING 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 information

Computer Vision. Filtering in the Frequency Domain

Computer 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 information

Spatially adaptive alpha-rooting in BM3D sharpening

Spatially adaptive alpha-rooting in BM3D sharpening Spatially adaptive alpha-rooting in BM3D sharpening Markku Mäkitalo and Alessandro Foi Department of Signal Processing, Tampere University of Technology, P.O. Box FIN-553, 33101, Tampere, Finland e-mail:

More information

Lecture 8: Interest Point Detection. Saad J Bedros

Lecture 8: Interest Point Detection. Saad J Bedros #1 Lecture 8: Interest Point Detection Saad J Bedros sbedros@umn.edu Review of Edge Detectors #2 Today s Lecture Interest Points Detection What do we mean with Interest Point Detection in an Image Goal:

More information

Edge Detection. Computer Vision P. Schrater Spring 2003

Edge Detection. Computer Vision P. Schrater Spring 2003 Edge Detection Computer Vision P. Schrater Spring 2003 Simplest Model: (Canny) Edge(x) = a U(x) + n(x) U(x)? x=0 Convolve image with U and find points with high magnitude. Choose value by comparing with

More information

9. Image filtering in the spatial and frequency domains

9. Image filtering in the spatial and frequency domains Image Processing - Laboratory 9: Image filtering in the spatial and frequency domains 9. Image filtering in the spatial and frequency domains 9.. Introduction In this laboratory the convolution operator

More information

LoG Blob Finding and Scale. Scale Selection. Blobs (and scale selection) Achieving scale covariance. Blob detection in 2D. Blob detection in 2D

LoG Blob Finding and Scale. Scale Selection. Blobs (and scale selection) Achieving scale covariance. Blob detection in 2D. Blob detection in 2D Achieving scale covariance Blobs (and scale selection) Goal: independently detect corresponding regions in scaled versions of the same image Need scale selection mechanism for finding characteristic region

More information

Quality Improves with More Rays

Quality Improves with More Rays Recap Quality Improves with More Rays Area Area 1 shadow ray 16 shadow rays CS348b Lecture 8 Pat Hanrahan / Matt Pharr, Spring 2018 pixelsamples = 1 jaggies pixelsamples = 16 anti-aliased Sampling and

More information