282 56
538 57
1088 58
2094 59
4052. 4052 60
8056. 61
15366 62
28743 63
49190. 64
65536. 65
Fourier transform magnitude 66
Masking out the fundamental and harmonics from periodic pillars 67
Name as many functions as you can that retain that same functional form in the transform domain 68
69 Forsyth&Ponce
Discrete-time, continuous frequency Fourier transform Oppenheim, Schafer and Buck, Discrete-time signal processing, Prentice Hall, 1999 70
Discrete-time, continuous frequency Fourier transform pairs Oppenheim, Schafer and Buck, Discrete-time signal processing, Prentice Hall, 1999 71
Bracewell’s pictorial dictionary of Fourier transform pairs 73 Bracewell, The Fourier Transform and its Applications, McGraw Hill 1978
Bracewell’s pictorial dictionary of Fourier transform pairs 74 Bracewell, The Fourier Transform and its Applications, McGraw Hill 1978
Bracewell’s pictorial dictionary of Fourier transform pairs 75 Bracewell, The Fourier Transform and its Applications, McGraw Hill 1978
Bracewell’s pictorial dictionary of Fourier transform pairs 76 Bracewell, The Fourier Transform and its Applications, McGraw Hill 1978
Why is the Fourier domain particularly useful? • It tells us the effect of linear convolutions. 77
Fourier transform of convolution Consider a (circular) convolution of g and h f � g � h 78
Fourier transform of convolution f � g � h Take DFT of both sides � � F m n DFT g h [ , ] � � 79
Fourier transform of convolution f � g � h � � F [ m , n ] � DFT g � h Write the DFT and convolution explicitly � � um vn � � M 1 N 1 � à � � i � � M N � � ��� � F m , n g [ u k , v l ] h [ k , l ] e � � � � � u 0 v 0 k , l 80
Fourier transform of convolution f � g � h � � F [ m , n ] � DFT g � h � � um vn � � M 1 N 1 � à i � � � � M N � � ��� � F m , n g [ u k , v l ] h [ k , l ] e � � � � � u 0 v 0 k , l Move the exponent in � � um vn � � M 1 N 1 � à � � i � � M N � g [ u k , v l ] e h [ k , l ] ��� � � � � � u 0 v 0 k , l 81
Fourier transform of convolution f � g � h � � F [ m , n ] � DFT g � h � � um vn � � M 1 N 1 � à i � � � � M N � � ��� � F m , n g [ u k , v l ] h [ k , l ] e � � � � � u 0 v 0 k , l � � um vn � � M 1 N 1 � à i � � � M N � � g [ u k , v l ] e h [ k , l ] ��� � � � � � u 0 v 0 k , l Change variables in the sum � à � � à � � � k m l n � � � � � � M k 1 N l 1 � à i � � � � M N � g [ , ] e h [ k , l ] � � � � � � � � � � k l k , l � � 82
Fourier transform of convolution f � g � h � � F [ m , n ] � DFT g � h � � um vn � � M 1 N 1 � à i � � � � M N � � ��� � F m , n g [ u k , v l ] h [ k , l ] e � � � � � u 0 v 0 k , l � � um vn � � M 1 N 1 � à i � � � M N � � g [ u k , v l ] e h [ k , l ] ��� � � � � � u 0 v 0 k , l � à � � à � � � k m l n � � � � � � M k 1 N l 1 � à i � � � � M N � g e h k l � � � [ , ] [ , ] � � � � � � � k l k , l � � Perform the DFT (circular boundary conditions) � � km ln � à � � i � � M N � � � G m , n e h [ k , l ] � � k , l 83
Fourier transform of convolution f � g � h � � F [ m , n ] � DFT g � h � � um vn � � M 1 N 1 � à i � � � � M N � � ��� � F m , n g [ u k , v l ] h [ k , l ] e � � � � � u 0 v 0 k , l � � um vn � � M 1 N 1 � à i � � � M N � � g [ u k , v l ] e h [ k , l ] ��� � � � � � u 0 v 0 k , l � à � � à � � � k m l n � � � � � � M k 1 N l 1 � à i � � � � M N � g e h k l � � � [ , ] [ , ] � � � � � � � k l k , l � � � km ln � � à i � � � � M N � � � � G m , n e h [ k , l ] � k , l Perform the other DFT (circular boundary conditions) � � � � G m n H m n � , , 84
Analysis of our simple filters 85
Analysis of our simple filters coefficient 1.0 0 Pixel offset original Filtered (no change) � � km ln � � M 1 N 1 � à i � � � � M N � � �� F [ m , n ] f [ k , l ] e � � k 0 l 0 1.0 constant 1 � 0 86
Analysis of our simple filters coefficient 1.0 0 Pixel offset original shifted � � km ln � � M 1 N 1 � à � i � � � M N � � �� F [ m , n ] f [ k � , l ] e � Constant � � k 0 l 0 magnitude, � m 1.0 � � i linearly shifted M e � phase 0 87
Analysis of our simple filters coefficient 0.3 0 Pixel offset original blurred � km ln � � � M 1 N 1 � à i � � � � M N � � �� F m n f k l e [ , ] [ , ] � � k 0 l 0 Low-pass filter � � � � 1 m � 1.0 � � à � � � 1 2 cos � � M 3 � � � � 0 88
Analysis of our simple filters 2.0 0.33 0 0 sharpened original � km ln � � � M 1 N 1 � à i � � � � M N � � �� F m n f k l e [ , ] [ , ] high-pass filter � � k 0 l 0 2.3 � � � � 1 m � � � � à � � � 2 1 2 cos � � 1.0 M 3 � � � � 0 89
Sampling and aliasing 90
Sampling in 1D takes a continuous function and replaces it with a vector of values, consisting of the function’s values at a set of sample points. We’ll assume that these sample points are on a regular grid, and can place one at each integer for convenience. 91
Sampling in 2D does the same thing, only in 2D. We’ll assume that these sample points are on a regular grid, and can place one at each integer point for convenience. 92
A continuous model for a sampled function • We want to be able to approximate integrals sensibly • Leads to – the delta function – model on right � � � � � � � Sample 2D f ( x , y ) f ( x , y ) � ( x � i , y � j ) i ��� i ��� � � � f ( x , y ) � � � ( x � i , y � j ) i ��� i ��� 93
The Fourier transform of a sampled signal � � �� �� � � F f ( x , y ) � ( x � i , y � j ) � F Sample 2D f ( x , y ) � � � � �� �� �� �� i ��� i ��� � � �� �� � F f ( x , y ) � ( x � i , y � j ) � � ** F �� � � �� �� �� i ��� i ��� � � � F u � i , v � j � � � � i ��� j ��� 94
95
96
Aliasing • Can’t shrink an image by taking every second pixel • If we do, characteristic errors appear – In the next few slides – Typically, small phenomena look bigger; fast phenomena can look slower – Common phenomenon • Wagon wheels rolling the wrong way in movies • Checkerboards misrepresented in ray tracing • Striped shirts look funny on colour television 97
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. 98
Constructing a pyramid by taking every second pixel leads to layers that badly misrepresent the top layer 99
Smoothing as low-pass filtering • The message of the FT is • A filter whose FT is a that high frequencies lead box is bad, because the to trouble with sampling. filter kernel has • Solution: suppress high infinite support frequencies before • Common solution: use sampling a Gaussian – multiply the FT of the signal with something – multiplying FT by that suppresses high Gaussian is equivalent frequencies to convolving image with Gaussian. – or convolve with a low-pass filter 100
Sampling without smoothing. Top row shows the images, sampled at every second pixel to get the next; bottom row shows the magnitude spectrum of these images. 101
Recommend
More recommend