940 likes | 956 Views
Dive into the world of 2D Fourier transforms, image sampling, and reconstruction. Learn about properties, algorithms, and applications to enhance your computational photography skills.
E N D
Computational Photography:Image Sampling and Reconstruction Jinxiang Chai
Review: 1D Fourier Transform A function f(x) can be represented as a weighted combination of phase-shifted sine waves How to compute F(u)? Inverse Fourier Transform Fourier Transform
Review: Box Function f(x) x |F(u)| u If f(x) is bounded, F(u) is unbounded
Review: Cosine -1 1 If f(x) is even, so is F(u)
Review: Gaussian If f(x) is gaussian, F(u) is also guassian.
Review: Properties Linearity: Time shift: Derivative: Integration: Convolution:
Outline 2D Fourier Transform Nyquist sampling theory Antialiasing Gaussian pyramid
Extension to 2D Fourier Transform: Inverse Fourier transform:
Building Block for 2D Transform Building block: Frequency: Orientation: Oriented wave fields
Building Block for 2D Transform Building block: Frequency: Orientation: Oriented wave fields Higher frequency
Some 2D Transforms From Lehar
Some 2D Transforms Why we have a DC component? From Lehar
Some 2D Transforms Why we have a DC component? From Lehar
Some 2D Transforms Why we have a DC component? From Lehar
Some 2D Transforms Why we have a DC component? From Lehar
Some 2D Transforms Why we have a DC component? - the sum of all pixel values From Lehar
Some 2D Transforms Why we have a DC component? - the sum of all pixel values Oriented stripe in spatial domain = an oriented line in spatial domain From Lehar
2D Fourier Transform Why? - Any relationship between two slopes?
2D Fourier Transform Why? - Any relationship between two slopes? Linearity
2D Fourier Transform Why? Linearity Why is the spectrum bounded?
Online Java Applet Click here.
2D Fourier Transform Pairs Gaussian Gaussian
2D Image Filtering Fourier transform Inverse transform From Lehar
2D Image Filtering Fourier transform Inverse transform Low-pass filter From Lehar
2D Image Filtering Fourier transform Inverse transform Low-pass filter high-pass filter From Lehar
2D Image Filtering Fourier transform Inverse transform Low-pass filter high-pass filter band-pass filter From Lehar
Aliasing Why does this happen?
Aliasing How to reduce it?
f(x) x fs(x) x … … -2T -T 0 T 2T Sampling Analysis Sampling
f(x) x fs(x) x … … -2T -T 0 T 2T Sampling Analysis Sampling Reconstruction
f(x) x fs(x) x … … -2T -T 0 T 2T Sampling Analysis What sampling rate (T) is sufficient to reconstruct the continuous version of the sampled signal? Sampling Reconstruction
Sampling Theory • How many samples are required to represent a given signal without loss of information? • What signals can be reconstructed without loss for a given sampling rate?
fs(x) x … … -2T -T 0 T 2T Sampling Analysis: Spatial Domain f(x) X … … x -2T -T 0 T 2T x ?
fs(x) x … … -2T -T 0 T 2T Sampling Analysis: Spatial Domain f(x) X … … x -2T -T 0 T 2T x ? What happens in Frequency domain?
Review: Dirac Delta and its Transform f(x) x |F(u)| 1 u Fourier transform and inverse Fourier transform are qualitatively the same, so knowing one direction gives you the other
Review: Fourier Transform Properties Linearity: Time shift: Derivative: Integration: Convolution:
Fourier Transform of Dirac Comb T 1/T Moving the spikes closer together in the spatial domain moves them farther apart in the frequency domain!
fs(x) x … … -2T -T 0 T 2T Sampling Analysis: Spatial Domain f(x) X … … x -2T -T 0 T 2T x ? What happens in Frequency domain?
F(u) fmax u -fmax Sampling Analysis: Freq. Domain … -1/T 0 1/T … u
F(u) fmax u -fmax Sampling Analysis: Freq. Domain … -1/T 0 1/T … u How does the convolution result look like?
F(u) fmax u -fmax Sampling Analysis: Freq. Domain … -1/T 0 1/T … u
F(u) fmax u -fmax Sampling Analysis: Freq. Domain … -1/T 0 1/T … u
F(u) fmax u -fmax Sampling Analysis: Freq. Domain … -1/T 0 1/T … u G(0)? G(fmax)? G(u)?
F(u) fmax u -fmax Sampling Analysis: Freq. Domain … -1/T 0 1/T … u G(0) = F(0) G(fmax) = F(fmax) G(u) = F(u)
F(u) fmax u -fmax Sampling Analysis: Freq. Domain How about … -1/T 0 1/T … u Fs(u) fmax u -fmax -1/T 1/T
F(u) fmax u -fmax Sampling Analysis: Freq. Domain How about … -1/T 0 1/T … u Fs(u) fmax u -fmax -1/T 1/T
F(u) fmax u -fmax Sampling Analysis: Freq. Domain … -1/T 0 1/T … u Fs(u) fmax u -fmax -1/T 1/T