Image Analysis. 3. Fourier Transform
|
|
- Harold Wheeler
- 5 years ago
- Views:
Transcription
1 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 Science University of Hildesheim Course on Image Analysis, winter term /66 Image Analysis 1. Fourier Series Representation 2. The Fourier Transform 3. Discrete Signals 4. Discrete Fourier Transform 5. Two-dimensional Fourier Transforms 6. Applications Course on Image Analysis, winter term /66
2 Image Analysis / 1. Fourier Series Representation Periodic Functions f () Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Periodic Functions f () Course on Image Analysis, winter term /66
3 Image Analysis / 1. Fourier Series Representation Periodic Functions A function f : R C is called T -periodic if f( + T ) = f() R Eample: the functions sin and cos are 2π-periodic. Eample: the function sin(ω) is 2π/ω-periodic. Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Periodic Functions / Approimation f () ˆf() = 1 sin() Course on Image Analysis, winter term /66
4 Image Analysis / 1. Fourier Series Representation Periodic Functions / Approimation f () ˆf() = 1 sin() + 2 cos Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Periodic Functions / Approimation f () ˆf() = 1 sin() + 2 cos + 2 sin 2 Course on Image Analysis, winter term /66
5 Image Analysis / 1. Fourier Series Representation Periodic Functions / Approimation f () ˆf() = 1 sin() + 2 cos + 2 sin cos 2 Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Periodic Functions / Approimation f () ˆf() = 1 sin() + 2 cos + 2 sin cos sin 3 Course on Image Analysis, winter term /66
6 Image Analysis / 1. Fourier Series Representation Periodic Functions / Approimation f () ˆf() = 1 sin() + 2 cos + 2 sin cos sin cos 3 Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Periodic Functions / Approimation f () ˆf() = 1 sin()+2 cos +2 sin 2+1 cos 2+3 sin 3+3 cos 3+2 sin 4 Course on Image Analysis, winter term /66
7 Image Analysis / 1. Fourier Series Representation Periodic Functions / Approimation f () ˆf() = 1 sin()+2 cos +2 sin 2+1 cos 2+3 sin 3+3 cos 3+2 sin 4+1 cos 4 Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Fourier Series Representation Theorem (Fourier Series Representation). Any continuous, differentiable and T -periodic function f : R C can be written as f() = a a k cos kω + b k sin kω, ω := 2π T k=1 with coefficients a k, b k C, called a Fourier Series of f. How to compute Fourier coefficients a k, b k for a given function f? Jean Baptiste Joseph Fourier ( ), French Mathematician and Physicist Note. Actually, the class of functions that can be represented as Fourier series is much larger (see, e.g., [Kö90, p. 314]). Course on Image Analysis, winter term /66
8 Image Analysis / 1. Fourier Series Representation Trigonometric Addition Formulas Lemma (Trigonometric Addition Formulas). For all, y R: cos( + y) = cos cos y sin sin y sin( + y) = sin cos y + sin y cos Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation cos a cos b d = sin a sin b d = sin a cos b d = Some Trigonometric Integrals 1 2 ( 1 sin(a + b) + 2 a + b ) sin(a b), if a b a b sin(2a) + 2a, else 4a ( sin(a + b) a + b ) sin(a b), if a b a b sin(2a) 2a, else 4a 1 2 ( cos(a + b) + a + b ) cos(a b), if a b a b sin 2 (a), else 2a Course on Image Analysis, winter term /66
9 Image Analysis / 1. Fourier Series Representation Some Trigonometric Integrals The si formulas easily can be proven by differentiation, e.g., cos a cos b d =? 1 ( ) sin(a + b) sin(a b) + 2 a + b a b for a b. Derivation d d yields: cos a cos b =? 1 ( (a + b) cos(a + b) + 2 a + b = 1 (cos(a + b) + cos(a b)) 2 ) (a b) cos(a b) a b = 1 (cos a cos b sin a sin b + cos a cos( b) sin a sin( b)) 2 = 1 (cos a cos b sin a sin b + cos a cos b + sin a sin b) 2 = cos a cos b Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Trigonometric Orthogonality Relations Let ω R +. The functions {sin kω k N, k > 0} {cos kω k N, k > 0} are pairwise othogonal with respect to f, g := +π/ω π/ω f()g()d i.e., for any two distinct such functions f, g f, g = +π/ω π/ω but f, f = π ω 0 f()g()d = 0 Course on Image Analysis, winter term /66
10 Image Analysis / 1. Fourier Series Representation Trigonometric Orthogonality Relations Proof: let f = cos kω and g = cos lω with k l, then: +π/ω [ ( 1 sin(k + l)ω sin(k l)ω f, g = cos kω cos lω d = + π/ω 2 (k + l)ω (k l)ω = 1 ( ) sin(k + l)π sin(k l)π sin (k + l)π sin (k l)π + 2 (k + l)ω (k l)ω (k + l)ω (k l)ω sin(k + l)π sin(k l)π = + (k + l)ω (k l)ω = 0 but [ sin(2kω) + 2kω f, f = 4kω ] +π/ω π/ω = )] +π/ω π/ω sin(2kπ) + 2kπ sin( 2kπ) + ( 2kπ) 4kω 4kω = π ω Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Fourier Series Representation Theorem (Fourier Series Representation). Any continuous, differentiable and T -periodic function f : R C can be written as f() = a a k cos kω + b k sin kω, k=1 ω := 2π T with coefficients a k, b k C, called a Fourier Series of f. How to compute Fourier coefficients a k, b k for a given function f? a k = 1 +π/ω π/ω b k = 1 π/ω π/ω +π/ω π/ω f() cos kω d f() sin kω d Note. Actually, the class of functions that can be represented as Fourier series is much larger (see, e.g., [Kö90, p. 314]). Course on Image Analysis, winter term /66
11 Image Analysis / 1. Fourier Series Representation Fourier Series Representation Proof. a? k = 1 +π/ω π/ω = 1 π/ω = 1 π/ω π/ω +π/ω π/ω = 1 π π/ω ω a k =a k ( +π/ω π/ω f() cos kω d ( a 0 ) 2 + a l cos lω + b l sin lω cos kω d l=1 a 0 2 cos kω d + l=1 +π/ω + π/ω +π/ω π/ω a l cos lω cos kω d b l sin lω cos kω d ) Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Fourier Representation / Rectangular function Let f be the 2π-periodic rectangular function 1, if ( π, 0) f() = 0, if { π, 0, π} +1, if (0, π) f () The Fourier representation of f is Proof: b k = 2 π a k = 2 π +π π +π π f() = 4 π f() sin k d = 2 +π π 2 sin k d = 4 π f() cos k d = 0 0 k N odd sin k k [ 1 ] π k cos k = { 4 4 (0 1) = πk (1 1) = 0, 4 πk πk, if k odd if k even in general, Fourier representations are infinite as in this eample! Course on Image Analysis, winter term /66
12 Image Analysis / 1. Fourier Series Representation Fourier Representation / Rectangular function f () Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Fourier Representation / Rectangular function / k = 1 f () Course on Image Analysis, winter term /66
13 Image Analysis / 1. Fourier Series Representation Fourier Representation / Rectangular function / k = 5 f () Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Fourier Representation / Rectangular function / k = 21 f () Course on Image Analysis, winter term /66
14 Image Analysis / 1. Fourier Series Representation Eulers Formula cos := n N( 1) n sin := n N( 1) n ep := n N n n! 2n (2n)! 2n+1 (2n + 1)! = 1 2 2! + 4 4! 6 6! + = 3 3! + 5 5! 7 7! + = ! + 3 3! + 4 4! + Lemma (Eulers formula). For C: e i = cos() + i sin(), with i := 1 the imaginary unit Proof: e i = (i) n n! n N = (i) 2n (2n)! + (i) 2n+1 (2n + 1)! n N n N = ( 1) n 2n (2n)! + i ( 1) n 2n+1 (2n + 1)! n N n N = cos() + i sin() Course on Image Analysis, winter term /66 Image Analysis / 1. Fourier Series Representation Comple Fourier Series Representation Theorem (Comple Fourier Series Representation). Any continuous, differentiable and T -periodic function f : R C can be written as f = c k e ikω, ω := 2π T k Z with coefficients c k C, called a Fourier Series of f. The coefficients of the comple Fourier series can be computed via c k = 1 +π/ω f() e ikω d 2π/ω π/ω Course on Image Analysis, winter term /66
15 Image Analysis / 1. Fourier Series Representation Comple Fourier Series Representation Proof. f = k Z c k e ikω with = k Z =c 0 + c k (cos kω + i sin kω) (c k + c k ) cos kω + (c k c k )i sin kω k=1 = a a k cos kω + b k sin kω k=1 and vice versa via: a 0 = 2c 0, a k = c k + c k, b k = i(c k c k ) c 0 = a 0 2, c k = 1 2 (a k ib k ), c k = 1 2 (a k + ib k ) Course on Image Analysis, winter term /66 Image Analysis 1. Fourier Series Representation 2. The Fourier Transform 3. Discrete Signals 4. Discrete Fourier Transform 5. Two-dimensional Fourier Transforms 6. Applications Course on Image Analysis, winter term /66
16 Image Analysis / 2. The Fourier Transform Fourier Transform Let f : R C be a function (that satisfies some regularity conditions). Then F : R C ω 1 2π + f()e iω d eists for each ω and is a continuous function called Fourier Transform of f (aka Fourier spectrum of f). One can show that (if F also satisfies some regularity conditions): f() = 1 2π + F (ω)e iω dω This is called Inverse Fourier Transform. We will write F(f) := F for the Fourier transform of a function f and F 1 (F ) for the inverse Fourier transform of a function F. Course on Image Analysis, winter term /66 Image Analysis / 2. The Fourier Transform Fourier Transforms / Eamples / Gaussian f () F () f() = 1 σ e 2 2σ 2 F () = e σ2 2 2 Course on Image Analysis, winter term /66
17 Image Analysis / 2. The Fourier Transform Fourier Transforms / Eamples / Uniform f () F () f() = δ( [ b, b]) F () = 2b sin(b) 2π Course on Image Analysis, winter term /66 Image Analysis / 2. The Fourier Transform Fourier Transforms / Eamples / Cosine cos () F(cos)() f() = cos(ω) F () = π (δ( ω) + δ( + ω)) 2 Course on Image Analysis, winter term /66
18 Image Analysis / 2. The Fourier Transform Fourier Transforms / Eamples / Sine sin () Im(F(sin)()) f() = sin(ω) F () = i π (δ( ω) δ( + ω)) 2 Course on Image Analysis, winter term /66 Image Analysis / 2. The Fourier Transform Properties function h Fourier transform H property name f F F f( ) inverse af + bg af + bg linearity f g F ()G() convolution 1 f()g() 2π F G multiplication f( a) e ia F () translation e ia f() F ( a) modulation f(/a) a F (a) scaling f F ( ) comple conjugate f() R F ( ) = F () hermitian symmetry Course on Image Analysis, winter term /66
19 Image Analysis 1. Fourier Series Representation 2. The Fourier Transform 3. Discrete Signals 4. Discrete Fourier Transform 5. Two-dimensional Fourier Transforms 6. Applications Course on Image Analysis, winter term /66 Image Analysis / 3. Discrete Signals The symbol Dirac Comb T () := n Z δ( T n) is called Dirac comb (aka impulse train, sampling function, Shah function) with sampling interval T. f() Lemma. The Fourier series of the Dirac comb T is T () = 1 e i2πk T T and its Fourier transform k Z F( T )() = 1 T 2π/T() = 1 T n Z δ( 2π T n) Course on Image Analysis, winter term /66
20 Image Analysis / 3. Discrete Signals Dirac Comb Proof: Obviously T is periodic with period T. Therefore f() = k Z c k e i2πk T with c k = 1 T = 1 T = 1 T +T +T/2 T/2 +T/2 T/2 = 1 T e i2πk T 0 T (y) e i2πk T y dy T (y) e i2πk T y dy δ(y) e i2πk T y dy = 1 T Course on Image Analysis, winter term /66 Image Analysis / 3. Discrete Signals Sampling f() f() Course on Image Analysis, winter term /66
21 Image Analysis / 3. Discrete Signals Sampling f() Sampling a function f at equidistant points T Z can be understood as f sampled () = f() T () Course on Image Analysis, winter term /66 Image Analysis / 3. Discrete Signals Fourier Transform of a Sampled Function Let f : R C be a function and f sampled be a sample of f with sampling period T. Then its Fourier transform F(f sampled ) is 2π/T -periodic and aggregates the Fourier transform F(f) over a 2π/T -periodic grid: F(f sampled )() = n Z F(f)( + n 2π T ) If F(f) vanishes for > π/t, then the Fourier transform F(f) is replicated in a period of the Fourier transform F(f sampled ). Otherwise replicas overlap and the Fourier transform becomes corrupted. This effect is called aliasing. Course on Image Analysis, winter term /66
22 Image Analysis / 3. Discrete Signals Fourier Transform of a Sampled Function The maimal occurring frequency ω ma := ma{ R, F(f)() 0} of the Fourier transform is called its bandwidth. The frequency ω s := 2π T of the sampling function is called its sampling frequency. Then the sampling frequency must be at least twice the bandwith: ω s > 2ω ma Course on Image Analysis, winter term /66 Image Analysis / 3. Discrete Signals Fourier Transform of a Sampled Function F(f) ω ma ω F(f_sampled1) ω ma ω ωs1 F(f_sampled2) ω ma ωs2 ω (cf. [BB08, p. 330]) Course on Image Analysis, winter term /66
23 Image Analysis / 3. Discrete Signals Fourier Transform of a Sampled Function Proof. F(f sampled )() =F(f T )() =F(f) F( T )() =F(f) 1 T 2π/T() =F(f) 1 δ( + n 2π T T ) n Z = 1 F(f) δ( + n 2π T T ) n Z = 1 F(f)( + n 2π T T ) n Z If < π/t and F(f)(y) vanishes for y > π/t, then F(f sampled )() = 1 T F(f)() Course on Image Analysis, winter term /66 Image Analysis / 3. Discrete Signals The Fourier Transform of a Discrete Function Let f be a discrete function: Then its Fourier transform is: F(f)(ω) = 1 2π f() = n Z y n δ( n), y n R = 1 2π = 1 2π n Z n Z f()e iω d y n δ( n) e iω d = 1 y n e iωn 2π n Z = 1 f(n) e iωn 2π n Z y n e iω δ( n)d i.e., a periodic function or equivalently: a function defined on an interval. Course on Image Analysis, winter term /66
24 Image Analysis 1. Fourier Series Representation 2. The Fourier Transform 3. Discrete Signals 4. Discrete Fourier Transform 5. Two-dimensional Fourier Transforms 6. Applications Course on Image Analysis, winter term /66 Image Analysis / 4. Discrete Fourier Transform Definition Let f : {0, 1,..., N 1} C be a finite discrete function, then F : {0, 1,..., N 1} C ω 1 N 1 N =0 = 1 N 1 N f()(cos(2π ω N ) i sin(2πω N )) =0 i2π ω f()e N is called discrete Fourier transform of f, denoted DFT(f). Then f() = 1 N 1 N ω=0 = 1 N 1 N ω=0 F (ω)(cos(2π ω N ) + i sin(2πω N )) i2π ω F (ω)e N This is called inverse discrete Fourier transform of f, denoted DFT 1 (F ).. Course on Image Analysis, winter term /66
25 Image Analysis / 4. Discrete Fourier Transform Eample DFT(f)(0) = =0 f() ω F (ω) i i i DFT i i i i i i i i i i i i DFT i i i i i f()(cos(2π 0 10 ) i sin(2π0 10 )) = 1 10 (cf. [BB08, p. 333]) 9 =0 f() = = Course on Image Analysis, winter term /66 Image Analysis / 4. Discrete Fourier Transform Eample DFT(f)(1) = 1 10 = =0 9 =0 f() ω F (ω) i i i DFT i i i i i i i i i i i i DFT i i i i i (cf. [BB08, p. 333]) f()(cos(2π 1 10 ) i sin(2π1 10 )) f() cos(π 5 ) i =0 f() sin(π ) = i 5 Course on Image Analysis, winter term /66
26 Image Analysis / 4. Discrete Fourier Transform Discrete Fourier Transform / Algorithm (naive) For C, denote R() its real part and I() its imaginary part, i.e., = R() + ii(), R(), I() R (for R one often also uses Re, for I also Im). To compute the discrete Fourier transform, one computes dft(f)(ω) for ω = 0,..., N 1 via: DFT(f)(ω) = 1 N 1 N =0 f()(cos(2π ω N ) i sin(2πω N )) = 1 N 1 (R(f()) + ii(f()))(cos(2π ω N N ) i sin(2πω N )) =0 = 1 N 1 N =0 N i N R(f()) cos(2π ω ) + I(f()) sin(2πω N N ) =0 R(f()) sin(2π ω ) + I(f()) cos(2πω N N ) Course on Image Analysis, winter term /66 Image Analysis / 4. Discrete Fourier Transform Discrete Fourier Transform / Algorithm (naive) 1 dft-naive(sequence f = (f() 0, f() 1 ) =0,...,N 1 ) : 2 N := length(f) 3 F := (F () 0, F () 1 ) =0,...,N 1 = (0, 0) =0,...,N 1 4 for ω := 0,..., N 1 do 5 c := (c 0, c 1 ) := (0, 0) 6 for := 0,..., N 1 do 7 c 0 := c 0 + f() 0 cos(2πω/n) + f() 1 sin(2πω/n) 8 c 1 := c 1 f() 0 sin(2πω/n) + f() 1 cos(2πω/n) 9 od 10 c 0 := c 0 / N 11 c 1 := c 1 / N 12 F (ω) := c 13 od 14 return F Course on Image Analysis, winter term /66
27 Image Analysis / 4. Discrete Fourier Transform Discrete Fourier Transform / Algorithm (naive) When computing the values of the discrete Fourier transform for different arguments ω, the cosine and sine functions repeatedly are callled with the same arguments: dft(f)(1): cos(2π 0/10) cos(2π 1/10) cos(2π 2/10) cos(2π 3/10) cos(2π 4/10) cos(2π 5/10) cos(2π 6/10) cos(2π 7/10) cos(2π 8/10) cos(2π 9/10) dft(f)(2): cos(2π 0/10) cos(2π 2/10) cos(2π 4/10) cos(2π 6/10) cos(2π 8/10) cos(2π 10/10) = cos(2π 0/10) cos(2π 12/10) = cos(2π 2/10) cos(2π 14/10) = cos(2π 4/10) cos(2π 16/10) = cos(2π 6/10) cos(2π 18/10) = cos(2π 8/10) Caching the epensive sine and cosine computations accelerates the algorithm! Course on Image Analysis, winter term /66 Image Analysis / 4. Discrete Fourier Transform Discrete Fourier Transform / Algorithm (naive, cached) 1 dft-naive-cached(sequence f = (f() 0, f() 1 ) =0,...,N 1 ) : 2 N := length(f) 3 for ω := 0,..., N 1 do 4 C(ω) := cos(2πω/n) 5 S(ω) := sin(2πω/n) 6 od 7 F := (F () 0, F () 1 ) =0,...,N 1 = (0, 0) =0,...,N 1 8 for ω := 0,..., N 1 do 9 c := (c 0, c 1 ) := (0, 0) 10 for := 0,..., N 1 do 11 c 0 := c 0 + f() 0 C(ω mod N) + f() 1 S(ω mod N) 12 c 1 := c 1 f() 0 S(ω mod N) + f() 1 C(ω mod N) 13 od 14 c 0 := c 0 / N 15 c 1 := c 1 / N 16 F (ω) := c 17 od 18 return F Course on Image Analysis, winter term /66
28 Image Analysis / 4. Discrete Fourier Transform Fast Fourier Transform (Gauss ca. 1805; Cooley/Tukey 1965) The naive algorithm for DFT still has compleity O(N 2 ). The Fast Fourier Transform algorithm is based on a decomposition of the DFT for sequences f of even length N: DFT(f)(ω) =DFT(f even )(ω mod N/2) + e i2πω/n DFT(f odd )(ω mod N/2) for ω = 0,..., N 1 and where f even () :=f(2), = 0,..., N/2 f odd () :=f(2 + 1) Course on Image Analysis, winter term /66 Image Analysis / 4. Discrete Fourier Transform Fast Fourier Transform / Proof Proof. N 1 i2π ω DFT(f)(ω) = f()e N = = = =0 N 1 =0 even N/2 1 =0 N/2 1 =0 N 1 i2π ω f()e N + i2π ω2 f(2)e N f even i2π ω ()e =0 odd N/2 1 + =0 i2π ω f()e N ω(2+1) i2π f(2 + 1)e N N/2 1 N/2 + e i2πω/n =0 =DFT(f even )(ω) + e i2πω/n DFT(f odd )(ω) f odd i2π ω ()e N/2 or more eaclty, as DFT(f even )(ω) is only defined for ω < N/2: =DFT(f even )(ω mod N/2) + e i2πω/n DFT(f odd )(ω mod N/2) Course on Image Analysis, winter term /66
29 Image Analysis / 4. Discrete Fourier Transform Fast Fourier Transform / Real version DFT(f)(ω) =DFT(f even )(ω mod N/2) + e i2πω/n DFT(f odd )(ω mod N/2) =DFT(f even )(ω mod N/2) and thus + (cos 2πω/N i sin 2πω/N)DFT(f odd )(ω mod N/2) R(DFT(f)(ω)) =R(DFT(f even )(ω mod N/2)) + cos 2πω/N R(DFT(f odd )(ω mod N/2)) + sin 2πω/N I(DFT(f odd )(ω mod N/2)) I(DFT(f)(ω)) =I(DFT(f even )(ω mod N/2)) + cos 2πω/N I(DFT(f odd )(ω mod N/2)) sin 2πω/N R(DFT(f odd )(ω mod N/2)) Course on Image Analysis, winter term /66 Image Analysis / 4. Discrete Fourier Transform Fast Fourier Transform / Algorithm 1 fft(sequence f = (f() 0, f() 1 ) =0,...,N 1 ) : 2 N := length(f) 3 if N is even 4 F := (F () 0, F () 1 ) =0,...,N 1 = (0, 0) =0,...,N 1 5 A := fft((f()) =0,2,4,...,N 2 ) 6 B := fft((f()) =1,3,5,,...,N 1 ) 7 for ω := 0,..., N 1 do 8 a := A(ω mod N/2) 9 b := B(ω mod N/2) 10 F (ω) 0 := a 0 + cos 2πω/N b 0 + sin 2πω/N b 1 11 F (ω) 1 := a 1 + cos 2πω/N b 1 sin 2πω/N b 0 12 od 13 return 14 else 15 F := dft-naive-cached(f) 16 fi 17 return F Course on Image Analysis, winter term /66
30 Image Analysis / 4. Discrete Fourier Transform Fast Fourier Transform / Outlook The computation of cos 2πω/N and sin 2πω/N also can be done recursively using the addition formulas. In this way, FFT best is applied to sequences of length 2 n (called radi-2 case). The FFT decomposition works with any factorization N = N 1 N 2 in a similar way, and thus also for sequences of length other than 2 n. FFT has compleity O(N log N) (if N is a power of 2). An early eperiment from 1969 reports a runtime of 13 1/2 hours for computing the DFT of a sequence of length 2048 by the naive method and 2.4 seconds using FFT. In practice, FFT is implemented in a linearized version avoiding eplicit recursions (see [CLRS03, p. 839]). Course on Image Analysis, winter term /66 Image Analysis / 4. Discrete Fourier Transform Types of Fourier Transforms type of function f sup f sup F(f) type of Fourier decomp. general (integrable) function R R (general) Fourier decomposition periodic function, function on interval interval I R Z Fourier series general (integrable) discrete function (sum of Diracs) Z interval I R discrete-time Fourier transform periodic discrete function, finite discrete function (finite sum of Diracs) finite I Z I discrete Fourier transform Course on Image Analysis, winter term /66
31 Image Analysis 1. Fourier Series Representation 2. The Fourier Transform 3. Discrete Signals 4. Discrete Fourier Transform 5. Two-dimensional Fourier Transforms 6. Applications Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms General Fourier Transform in 2D For two-dimensional functions f : R R C Fourier Transforms, Fourier Series and Discrete Fourier Transforms can be defined analogously. Let f : R R C be a function (that satisfies some regularity conditions). Then F : R R C (ω 1, ω 2 ) f(, y) e iω1 e iω2y dy d 2π eists for each (ω 1, ω 2 ) and is a continuous function called Fourier Transform of f (aka Fourier spectrum of f). One can show that (if F also satisfies some regularity conditions): f(, y) = F (ω 1, ω 2 )e iω1 e iω2y dω 1 dω 2 2π This is called Inverse Fourier Transform. We will write F(f) := F for the Fourier transform of a function f and F 1 (F ) for the inverse Fourier transform of a function F. Course on Image Analysis, winter term /66
32 Image Analysis / 5. Two-dimensional Fourier Transforms Bases in 2D Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms Fourier Series in 2D If f : R R C is (T1, T2)-periodic, i.e., f (, y) = f ( + T1, y + T2), y R then f can already be reconstructed from Z-many Fourier coefficients: 1 XX i 2π n i 2π my f (, y) = F (n, m) e T1 e T2 2π n Z m Z with 1 F (n, m) := 2π Z +T1 /2 Z +T2 /2 T1 /2 T2 /2 f (, y) e inω1 e imω2y dy d Course on Image Analysis, winter term /66
33 Image Analysis / 5. Two-dimensional Fourier Transforms If f is discrete, i.e., f(, y) = n Z Discrete Fourier Transform in 2D y n,m δ( T 1 n, y T 2 m), m Z then its Fourier transform is periodic: with f(, y) = 1 2π + + F (ω 1, ω 2 ) := 1 2π n Z m Z y n,m R F (ω 1, ω 2 ) e iω 1 e iω 2y dω 1 dω 2 f(n, m) e iω1n e iω 2m Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms Discrete Fourier Transform in 2D And finally, if f is discrete and finite, i.e., f(, y) = N 1 n=0 M 1 m=0 y n,m δ( n, y m), y n,m R then its Fourier transform is periodic and made from finitely many components: with f(, y) = F (ω 1, ω 2 ) := N 1 1 NM n=0 N 1 1 NM n=0 M 1 m=0 M 1 m=0 F (ω 1, ω 2 ) e iω 1 e iω 2y f(n, m) e iω 1n e iω 2m Course on Image Analysis, winter term /66
34 Image Analysis / 5. Two-dimensional Fourier Transforms Discrete Fourier Transform in 2D / Algorithm (naive, cached) 1 dft-2d-naive-cached(array f = (f() 0, f() 1 ) =0,...,N 1,y=0,...,M 1 ) : 2 for ω := 0,..., N 1 do 3 C 1 (ω) := cos(2πω/n) 4 S 1 (ω) := sin(2πω/n) 5 od 6 for ω := 0,..., M 1 do 7 C 2 (ω) := cos(2πω/m) 8 S 2 (ω) := sin(2πω/m) 9 od 10 F := (F () 0, F () 1 ) =0,...,N 1,y=0,...,M 1 = (0, 0) =0,...,N 1,y=0,...,M 1 11 for ω 1 := 0,..., N 1 do 12 for ω 2 := 0,..., M 1 do 13 c := (c 0, c 1 ) := (0, 0) 14 for := 0,..., N 1 do 15 for y := 0,..., M 1 do 16 C := C 1 (ω 1 mod N) C 2 (ω 2 y mod M) S 1 (ω 1 mod N) S 2 (ω 2 y mod M) 17 S := S 1 (ω 1 mod N) C 2 (ω 2 y mod M) + C 1 (ω 1 mod N) S 2 (ω 2 y mod M) 18 c 0 := c 0 + f() 0 C f() 1 S 19 c 1 := c 1 f() 0 S + f() 1 C 20 od 21 od 22 c 0 := c 0 / NM 23 c 1 := c 1 / NM 24 F (ω) := c 25 od 26 od 27 return F Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms Discrete Fourier Transform in 2D / Separability F (ω 1, ω 2 ) := 1 N 1 NM n=0 n=0 M 1 m=0 f(n, m) e iω 1n e iω 2m ( ) = 1 N 1 M 1 1 f(n, m) e iω 2m N M m=0 e iω 1n = 1 N 1 DFT((f(n, m)) m=1,...,m )(ω 2 ) e iω 1n N n=0 Course on Image Analysis, winter term /66
35 Image Analysis / 5. Two-dimensional Fourier Transforms Discrete Fourier Transform in 2D / FFT 1 fft-2d(array f = (f() 0, f() 1 ) =0,...,N 1,y=0,...,M 1 ) : 2 G := (G() 0, G() 1 ) =0,...,N 1,y=0,...,M 1 = (0, 0) =0,...,N 1,y=0,...,M 1 3 for ω 1 := 0,..., N 1 do 4 G(ω 1,.) := fft(f(ω 1, y) y=0,...,m 1 ) 5 od 6 F := (F () 0, F () 1 ) =0,...,N 1,y=0,...,M 1 = (0, 0) =0,...,N 1,y=0,...,M 1 7 for ω 2 := 0,..., M 1 do 8 F (., ω 2 ) := fft(g(, ω 2 ) =0,...,N 1 ) 9 od 10 return F Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms Power Spectrum The fourier spectrum of is a discrete N M gray-scale image f, i.e., with one channel, a discrete N M image F with comple intensity values, i.e., two channels. For visualization one usually shows the power spectrum defined as: F power (f)() := (RF(f)()) 2 + (IF(f)()) 2 The power spectrum measures the absolute value of the comple amplitude. The complementary information θ called phase is not shown. Course on Image Analysis, winter term /66
36 Image Analysis / 5. Two-dimensional Fourier Transforms What does Power Spectrum mean? Comple Coordinates Cartesian coordinates: Polar coordinates: Im r θ Re = (R, I) R R R =r cos θ I =r sin θ = (r, θ) R + 0 [0, 2π) r = (R) 2 + (I) 2 arctan(r/i), if > 0, y > 0 θ = arctan(r/i) + π, if y < 0 arctan(r/i) + 2π, if > 0, y < 0 = R + ii = r cos θ + ir sin θ = re iθ Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66
37 Image Analysis / 5. Two-dimensional Fourier Transforms Displaying Fourier Power Spectra For displaying spectra, some further conventions are used: As the scale of many power spectra is dominated by a few large values, one usually plots log F power (f)() or (F power (f)()) 1 2 instead of the raw power spectrum values F power (f)(). Usually, the centered spectrum is shown, i.e., the intensities for and { N/2, N/2 + 1,..., 1, 0, 1,..., N/2 1, N/2} y { M/2, M/2 + 1,..., 1, 0, 1,..., M/2 1, M/2} instead of the intensities for 0, 1,..., N 1 and 0, 1,..., M 1. Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms Centered Spectrum A B C D Course on Image Analysis, winter term /66
38 Image Analysis / 5. Two-dimensional Fourier Transforms Centered Spectrum A B C A D A C B C B D A D B C D Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms Centered Spectrum original spectrum: A C centered spectrum: B D C D B A Course on Image Analysis, winter term /66
39 Image Analysis / 5. Two-dimensional Fourier Transforms Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms Symmetry of Fouier Power Spectra for Real Images Usually images are real, i.e., f() R (not C). For real functions, we know that F(f)( ) = F (f)() is hermitian (with := R i I). As C has the same radius as : r() = (R) 2 + (I) 2 =? r( ) = r(r ii) = (R) 2 + ( I) 2 for real functions F(f) and F (f) have the same radius and thus F power (f)( ) = F power (f)() i.e., the power spectrum is symmetric around the origin. Course on Image Analysis, winter term /66
40 Image Analysis / 5. Two-dimensional Fourier Transforms Symmetry of Fouier Power Spectra for Real Images image f: power spectrum F power (f): Spectra of real images are symmetric around the origin (red circle). So storing just half of the power spectrum is sufficient (e.g., above green line any line through the origin will do). Course on Image Analysis, winter term /66 Image Analysis / 5. Two-dimensional Fourier Transforms Another Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66
41 Image Analysis 1. Fourier Series Representation 2. The Fourier Transform 3. Discrete Signals 4. Discrete Fourier Transform 5. Two-dimensional Fourier Transforms 6. Applications Course on Image Analysis, winter term /66 Image Analysis / 6. Applications Low Pass Filters High frequencies are responsible for sharp edges, low frequencies for constant and slowly changing areas. Low pass filters retain only low frequencies ω ω ma, i.e., filter out high frequencies. and thus smooth / blur an image. Course on Image Analysis, winter term /66
42 Image Analysis / 6. Applications Low Pass Filters / Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66 Image Analysis / 6. Applications Low Pass Filters / Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66
43 Image Analysis / 6. Applications Low Pass Filters / Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66 Image Analysis / 6. Applications High Pass Filters High frequencies are responsible for sharp edges, low frequencies for constant and slowly changing areas. High pass filters retain only high frequencies ω ω min, i.e., filter out low frequencies. and thus sharpen an image and detect edges. Course on Image Analysis, winter term /66
44 Image Analysis / 6. Applications High Pass Filters / Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66 Image Analysis / 6. Applications High Pass Filters / Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66
45 Image Analysis / 6. Applications High Pass Filters / Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66 Image Analysis / 6. Applications Band Pass Filters High frequencies are responsible for sharp edges, low frequencies for constant and slowly changing areas. Band pass filters retain only frequencies ω [ω min, ω ma ] in a given interval (the frequency band), i.e., filter out low and high frequencies. and thus detect edges. Course on Image Analysis, winter term /66
46 Image Analysis / 6. Applications Band Pass Filters / Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66 Image Analysis / 6. Applications Band Pass Filters / Eample image f: power spectrum F power (f): Course on Image Analysis, winter term /66
47 Image Analysis / 6. Applications Reducing Periodic Noise image f: power spectrum F power (f): Course on Image Analysis, winter term /66 Image Analysis / 6. Applications Reducing Periodic Noise image f: power spectrum F power (f): Course on Image Analysis, winter term /66
48 Image Analysis / 6. Applications Reducing Periodic Noise image f: power spectrum F power (f): Periodic noise can be reduced by filtering out the frequencies belonging to the periodic noise pattern. This also can be understood as a simple method for inpainting. Course on Image Analysis, winter term /66 Image Analysis / 6. Applications Reducing Salt and Pepper Noise image f: power spectrum F power (f): Course on Image Analysis, winter term /66
49 Image Analysis / 6. Applications Reducing Salt and Pepper Noise image f: power spectrum F power (f): Reducing non-periodic noise patterns via frequency filters is difficult. Course on Image Analysis, winter term /66 Image Analysis / 6. Applications Deconvolution via Fourier Transform Assume, an image f has been corrupted by a convolution with a kernel k (e.g., blurred): g = k f If the kernel k is known, one can undo the convolution using the Fourier transform: F(g) = F(k f) = F(k) F(f) F(f) = F(g) F(k) ( ) F(g) f = F 1 F(f) = F 1 F(k) Course on Image Analysis, winter term /66
50 Image Analysis / 6. Applications Schedule Schedule until Christmas: net Tue., 9.12., no lecture. net Wed., 10.12, lecture. Tue., , no lecture. Wed., 17.12, lecture. Course on Image Analysis, winter term /66 Image Analysis / 6. Applications References [BB08] Wilhelm Burger and Mark J. Burge. Digital Image Processing, An Algorithmic Introduction Using Java. Springer, [CLRS03] Thomas H. Corman, Charles E. Leiserson, Ronald R. Rivest, and Clifford Stein. An Introduction to Algorithms. The MIT Press, [Kö90] Konrad Königsberger. Analysis 1. Springer, Course on Image Analysis, winter term /66
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 informationFourier Series Example
Fourier Series Example Let us compute the Fourier series for the function on the interval [ π,π]. f(x) = x f is an odd function, so the a n are zero, and thus the Fourier series will be of the form f(x)
More informationFourier Transform 2D
Image Processing - Lesson 8 Fourier Transform 2D Discrete Fourier Transform - 2D Continues Fourier Transform - 2D Fourier Properties Convolution Theorem Eamples = + + + The 2D Discrete Fourier Transform
More informationOutline. Math Partial Differential Equations. Fourier Transforms for PDEs. Joseph M. Mahaffy,
Outline Math 53 - Partial Differential Equations s for PDEs Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu Department of Mathematics and Statistics Dynamical Systems Group Computational Sciences Research Center
More informationThe Fourier Transform
The Fourier Transorm Fourier Series Fourier Transorm The Basic Theorems and Applications Sampling Bracewell, R. The Fourier Transorm and Its Applications, 3rd ed. New York: McGraw-Hill, 2. Eric W. Weisstein.
More informationToday 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 information5. THE CLASSES OF FOURIER TRANSFORMS
5. THE CLASSES OF FOURIER TRANSFORMS There are four classes of Fourier transform, which are represented in the following table. So far, we have concentrated on the discrete Fourier transform. Table 1.
More informationECE 425. Image Science and Engineering
ECE 425 Image Science and Engineering Spring Semester 2000 Course Notes Robert A. Schowengerdt schowengerdt@ece.arizona.edu (520) 62-2706 (voice), (520) 62-8076 (fa) ECE402 DEFINITIONS 2 Image science
More informationFOURIER TRANSFORM METHODS David Sandwell, January, 2013
1 FOURIER TRANSFORM METHODS David Sandwell, January, 2013 1. Fourier Transforms Fourier analysis is a fundamental tool used in all areas of science and engineering. The fast fourier transform (FFT) algorithm
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 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 information3. 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 informationReview 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 informationContents. Signals as functions (1D, 2D)
Fourier Transform The idea A signal can be interpreted as en electromagnetic wave. This consists of lights of different color, or frequency, that can be split apart usign an optic prism. Each component
More informationECE 301 Fall 2011 Division 1 Homework 10 Solutions. { 1, for 0.5 t 0.5 x(t) = 0, for 0.5 < t 1
ECE 3 Fall Division Homework Solutions Problem. Reconstruction of a continuous-time signal from its samples. Let x be a periodic continuous-time signal with period, such that {, for.5 t.5 x(t) =, for.5
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 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 informationLecture 34. Fourier Transforms
Lecture 34 Fourier Transforms In this section, we introduce the Fourier transform, a method of analyzing the frequency content of functions that are no longer τ-periodic, but which are defined over the
More informationContents. Signals as functions (1D, 2D)
Fourier Transform The idea A signal can be interpreted as en electromagnetic wave. This consists of lights of different color, or frequency, that can be split apart usign an optic prism. Each component
More information1-D MATH REVIEW CONTINUOUS 1-D FUNCTIONS. Kronecker delta function and its relatives. x 0 = 0
-D MATH REVIEW CONTINUOUS -D FUNCTIONS Kronecker delta function and its relatives delta function δ ( 0 ) 0 = 0 NOTE: The delta function s amplitude is infinite and its area is. The amplitude is shown as
More informationLecture 7 January 26, 2016
MATH 262/CME 372: Applied Fourier Analysis and Winter 26 Elements of Modern Signal Processing Lecture 7 January 26, 26 Prof Emmanuel Candes Scribe: Carlos A Sing-Long, Edited by E Bates Outline Agenda:
More informationGabor filters. Konstantinos G. Derpanis York University. Version 1.3. April 23, 2007
Gabor filters Konstantinos G. Derpanis York University kosta@cs.yorku.ca Version.3 April 3, 7 In this note the Gabor filter is reviewed. The Gabor filter was originally introduced by Dennis Gabor (Gabor,
More informationFigure 3.1 Effect on frequency spectrum of increasing period T 0. Consider the amplitude spectrum of a periodic waveform as shown in Figure 3.2.
3. Fourier ransorm From Fourier Series to Fourier ransorm [, 2] In communication systems, we oten deal with non-periodic signals. An extension o the time-requency relationship to a non-periodic signal
More informationFourier Series and Fourier Transforms
Fourier Series and Fourier Transforms EECS2 (6.082), MIT Fall 2006 Lectures 2 and 3 Fourier Series From your differential equations course, 18.03, you know Fourier s expression representing a T -periodic
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 informationDigital Image Processing
Digital Image Processing Part 3: Fourier Transform and Filtering in the Frequency Domain AASS Learning Systems Lab, Dep. Teknik Room T109 (Fr, 11-1 o'clock) achim.lilienthal@oru.se Course Book Chapter
More informationStatistical Geometry Processing Winter Semester 2011/2012
Statistical Geometry Processing Winter Semester 2011/2012 Linear Algebra, Function Spaces & Inverse Problems Vector and Function Spaces 3 Vectors vectors are arrows in space classically: 2 or 3 dim. Euclidian
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 informationa k cos kω 0 t + b k sin kω 0 t (1) k=1
MOAC worksheet Fourier series, Fourier transform, & Sampling Working through the following exercises you will glean a quick overview/review of a few essential ideas that you will need in the moac course.
More informationTHE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3
THE FOURIER TRANSFORM (Fourier series for a function whose period is very, very long) Reading: Main 11.3 Any periodic function f(t) can be written as a Fourier Series a 0 2 + a n cos( nωt) + b n sin n
More informationChapter 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 informationContents. Signals as functions (1D, 2D)
Fourier Transform The idea A signal can be interpreted as en electromagnetic wave. This consists of lights of different color, or frequency, that can be split apart usign an optic prism. Each component
More informationPS403 - Digital Signal processing
PS403 - Digital Signal processing III. DSP - Digital Fourier Series and Transforms Key Text: Digital Signal Processing with Computer Applications (2 nd Ed.) Paul A Lynn and Wolfgang Fuerst, (Publisher:
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 informationFundamentals of the Discrete Fourier Transform
Seminar presentation at the Politecnico di Milano, Como, November 12, 2012 Fundamentals of the Discrete Fourier Transform Michael G. Sideris sideris@ucalgary.ca Department of Geomatics Engineering University
More informationPoles and Zeros in z-plane
M58 Mixed Signal Processors page of 6 Poles and Zeros in z-plane z-plane Response of discrete-time system (i.e. digital filter at a particular frequency ω is determined by the distance between its poles
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 informationSignal Processing COS 323
Signal Processing COS 323 Digital Signals D: functions of space or time e.g., sound 2D: often functions of 2 spatial dimensions e.g. images 3D: functions of 3 spatial dimensions CAT, MRI scans or 2 space,
More informationEE16B - Spring 17 - Lecture 11B Notes 1
EE6B - Spring 7 - Lecture B Notes Murat Arcak 6 April 207 Licensed under a Creative Commons Attribution-NonCommercial-ShareAlike 4.0 International License. Interpolation with Basis Functions Recall that
More informationMARN 5898 Fourier Analysis.
MAR 5898 Fourier Analysis. Dmitriy Leykekhman Spring 2010 Goals Fourier Series. Discrete Fourier Transforms. D. Leykekhman - MAR 5898 Parameter estimation in marine sciences Linear Least Squares 1 Complex
More informationSound & Vibration Magazine March, Fundamentals of the Discrete Fourier Transform
Fundamentals of the Discrete Fourier Transform Mark H. Richardson Hewlett Packard Corporation Santa Clara, California The Fourier transform is a mathematical procedure that was discovered by a French mathematician
More informationFourier transforms. R. C. Daileda. Partial Differential Equations April 17, Trinity University
The Fourier Transform R. C. Trinity University Partial Differential Equations April 17, 214 The Fourier series representation For periodic functions Recall: If f is a 2p-periodic (piecewise smooth) function,
More informationQueens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 2018 Instructor: Dr. Sateesh Mane.
Queens College, CUNY, Department of Computer Science Numerical Methods CSCI 361 / 761 Spring 018 Instructor: Dr. Sateesh Mane c Sateesh R. Mane 018 30 Lecture 30 April 7, 018 Fourier Series In this lecture,
More informationFourier Series. Fourier Transform
Math Methods I Lia Vas Fourier Series. Fourier ransform Fourier Series. Recall that a function differentiable any number of times at x = a can be represented as a power series n= a n (x a) n where the
More information1 From Fourier Series to Fourier transform
Differential Equations 2 Fall 206 The Fourier Transform From Fourier Series to Fourier transform Recall: If fx is defined and piecewise smooth in the interval [, ] then we can write where fx = a n = b
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 informationTopic 3: Fourier Series (FS)
ELEC264: Signals And Systems Topic 3: Fourier Series (FS) o o o o Introduction to frequency analysis of signals CT FS Fourier series of CT periodic signals Signal Symmetry and CT Fourier Series Properties
More informationDiscrete Fourier Transform
Discrete Fourier Transform Virtually all practical signals have finite length (e.g., sensor data, audio records, digital images, stock values, etc). Rather than considering such signals to be zero-padded
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 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 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
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 information2 Fourier Transforms and Sampling
2 Fourier ransforms and Sampling 2.1 he Fourier ransform he Fourier ransform is an integral operator that transforms a continuous function into a continuous function H(ω) =F t ω [h(t)] := h(t)e iωt dt
More informationDiscrete Time Fourier Transform (DTFT)
Discrete Time Fourier Transform (DTFT) 1 Discrete Time Fourier Transform (DTFT) The DTFT is the Fourier transform of choice for analyzing infinite-length signals and systems Useful for conceptual, pencil-and-paper
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 informationTopic 4 Notes Jeremy Orloff
Topic 4 Notes Jeremy Orloff 4 Complex numbers and exponentials 4.1 Goals 1. Do arithmetic with complex numbers.. Define and compute: magnitude, argument and complex conjugate of a complex number. 3. Euler
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 informationFrom Fourier to Wavelets in 60 Slides
From Fourier to Wavelets in 60 Slides Bernhard G. Bodmann Math Department, UH September 20, 2008 B. G. Bodmann (UH Math) From Fourier to Wavelets in 60 Slides September 20, 2008 1 / 62 Outline 1 From Fourier
More informationDIFFERENTIATING UNDER THE INTEGRAL SIGN
DIFFEENTIATING UNDE THE INTEGAL SIGN KEITH CONAD I had learned to do integrals by various methods shown in a book that my high school physics teacher Mr. Bader had given me. [It] showed how to differentiate
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 informationLecture 6 January 21, 2016
MATH 6/CME 37: Applied Fourier Analysis and Winter 06 Elements of Modern Signal Processing Lecture 6 January, 06 Prof. Emmanuel Candes Scribe: Carlos A. Sing-Long, Edited by E. Bates Outline Agenda: Fourier
More informationFourier Transforms For additional information, see the classic book The Fourier Transform and its Applications by Ronald N. Bracewell (which is on the shelves of most radio astronomers) and the Wikipedia
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 informationUnstable Oscillations!
Unstable Oscillations X( t ) = [ A 0 + A( t ) ] sin( ω t + Φ 0 + Φ( t ) ) Amplitude modulation: A( t ) Phase modulation: Φ( t ) S(ω) S(ω) Special case: C(ω) Unstable oscillation has a broader periodogram
More informationFourier transform. Stefano Ferrari. Università degli Studi di Milano Methods for Image Processing. academic year
Fourier transform Stefano Ferrari Università degli Studi di Milano stefano.ferrari@unimi.it Methods for Image Processing academic year 27 28 Function transforms Sometimes, operating on a class of functions
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 informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 12 Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign Copyright c 2002. Reproduction permitted for noncommercial,
More informationAlgorithms of Scientific Computing
Algorithms of Scientific Computing Discrete Fourier Transform (DFT) Michael Bader Technical University of Munich Summer 2018 Fast Fourier Transform Outline Discrete Fourier transform Fast Fourier transform
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 informationFourier Series. (Com S 477/577 Notes) Yan-Bin Jia. Nov 29, 2016
Fourier Series (Com S 477/577 otes) Yan-Bin Jia ov 9, 016 1 Introduction Many functions in nature are periodic, that is, f(x+τ) = f(x), for some fixed τ, which is called the period of f. Though function
More informationWave Phenomena Physics 15c
Wave Phenomena Physics 5c Lecture Fourier Analysis (H&L Sections 3. 4) (Georgi Chapter ) What We Did Last ime Studied reflection of mechanical waves Similar to reflection of electromagnetic waves Mechanical
More informationComputational Methods CMSC/AMSC/MAPL 460
Computational Methods CMSC/AMSC/MAPL 460 Fourier transform Balaji Vasan Srinivasan Dept of Computer Science Several slides from Prof Healy s course at UMD Last time: Fourier analysis F(t) = A 0 /2 + A
More informationRepresenting a Signal
The Fourier Series Representing a Signal The convolution method for finding the response of a system to an excitation takes advantage of the linearity and timeinvariance of the system and represents the
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 informationENSC327 Communications Systems 2: Fourier Representations. School of Engineering Science Simon Fraser University
ENSC37 Communications Systems : Fourier Representations School o Engineering Science Simon Fraser University Outline Chap..5: Signal Classiications Fourier Transorm Dirac Delta Function Unit Impulse Fourier
More informationThe O () notation. Definition: Let f(n), g(n) be functions of the natural (or real)
The O () notation When analyzing the runtime of an algorithm, we want to consider the time required for large n. We also want to ignore constant factors (which often stem from tricks and do not indicate
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 informationFourier Transform and Frequency Domain
Fourier Transform and Frequency Domain http://graphics.cs.cmu.edu/courses/15-463 15-463, 15-663, 15-862 Computational Photography Fall 2017, Lecture 6 Course announcements Last call for responses to Doodle
More information4.5 The Open and Inverse Mapping Theorem
4.5 The Open and Inverse Mapping Theorem Theorem 4.5.1. [Open Mapping Theorem] Let f be analytic in an open set U such that f is not constant in any open disc. Then f(u) = {f() : U} is open. Lemma 4.5.1.
More informationHarmonic Analysis Homework 5
Harmonic Analysis Homework 5 Bruno Poggi Department of Mathematics, University of Minnesota November 4, 6 Notation Throughout, B, r is the ball of radius r with center in the understood metric space usually
More informationSignal and systems. Linear Systems. Luigi Palopoli. Signal and systems p. 1/5
Signal and systems p. 1/5 Signal and systems Linear Systems Luigi Palopoli palopoli@dit.unitn.it Wrap-Up Signal and systems p. 2/5 Signal and systems p. 3/5 Fourier Series We have see that is a signal
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 information6 From Analog to Discrete Signals
6 From Analog to Discrete Signals You don t have to be a mathematician to have a feel for numbers. - John Forbes ash (1928 - As you may recall, a goal of this course has been to introduce some of the tools
More informationLAB 6: FIR Filter Design Summer 2011
University of Illinois at Urbana-Champaign Department of Electrical and Computer Engineering ECE 311: Digital Signal Processing Lab Chandra Radhakrishnan Peter Kairouz LAB 6: FIR Filter Design Summer 011
More informationLAST Name FOURIER FIRST Name Jean Baptiste Joseph Lab Time 365/24/7
EECS 20N: Structure and Interpretation of Signals and Systems Final Exam Sol. Department of Electrical Engineering and Computer Sciences 13 December 2005 UNIVERSITY OF CALIFORNIA BERKELEY LAST Name FOURIER
More informationExpectation Maximization Deconvolution Algorithm
Epectation Maimization Deconvolution Algorithm Miaomiao ZHANG March 30, 2011 Abstract In this paper, we use a general mathematical and eperimental methodology to analyze image deconvolution. The main procedure
More informationVII. Wavelet Bases. One can construct wavelets ψ such that the dilated and translated family. 2 j
VII Wavelet Bases One can construct wavelets ψ such that the dilated and translated family { ψ j,n (t) = ( t j )} ψ n j j (j,n) Z is an orthonormal basis of L (R). Behind this simple statement lie very
More informationChapter 4 The Fourier Series and Fourier Transform
Chapter 4 The Fourier Series and Fourier Transform Representation of Signals in Terms of Frequency Components Consider the CT signal defined by N xt () = Acos( ω t+ θ ), t k = 1 k k k The frequencies `present
More informationE2.5 Signals & Linear Systems. Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & 2)
E.5 Signals & Linear Systems Tutorial Sheet 1 Introduction to Signals & Systems (Lectures 1 & ) 1. Sketch each of the following continuous-time signals, specify if the signal is periodic/non-periodic,
More informationMEDE2500 Tutorial Nov-7
(updated 2016-Nov-4,7:40pm) MEDE2500 (2016-2017) Tutorial 3 MEDE2500 Tutorial 3 2016-Nov-7 Content 1. The Dirac Delta Function, singularity functions, even and odd functions 2. The sampling process and
More informationMathematics for Chemists 2 Lecture 14: Fourier analysis. Fourier series, Fourier transform, DFT/FFT
Mathematics for Chemists 2 Lecture 14: Fourier analysis Fourier series, Fourier transform, DFT/FFT Johannes Kepler University Summer semester 2012 Lecturer: David Sevilla Fourier analysis 1/25 Remembering
More informationCMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation
CMPT 318: Lecture 5 Complex Exponentials, Spectrum Representation Tamara Smyth, tamaras@cs.sfu.ca School of Computing Science, Simon Fraser University January 23, 2006 1 Exponentials The exponential is
More informationTutorial 9 The Discrete Fourier Transform (DFT) SIPC , Spring 2017 Technion, CS Department
Tutorial 9 The Discrete Fourier Transform (DFT) SIPC 236327, Spring 2017 Technion, CS Department The DFT Matrix The DFT matrix of size M M is defined as DFT = 1 M W 0 0 W 0 W 0 W where W = e i2π M i =
More informationsech(ax) = 1 ch(x) = 2 e x + e x (10) cos(a) cos(b) = 2 sin( 1 2 (a + b)) sin( 1 2
Math 60 43 6 Fourier transform Here are some formulae that you may want to use: F(f)(ω) def = f(x)e iωx dx, F (f)(x) = f(ω)e iωx dω, (3) F(f g) = F(f)F(g), (4) F(e α x ) = α π ω + α, F( α x + α )(ω) =
More informationReview of Discrete-Time System
Review of Discrete-Time System Electrical & Computer Engineering University of Maryland, College Park Acknowledgment: ENEE630 slides were based on class notes developed by Profs. K.J. Ray Liu and Min Wu.
More informationSolution Set #1
05-738-0083 Solution Set # i. Evaluate δ [SINC []] = δ sin[π] π Recall relation for Dirac delta function wit functional argument (Fourier Metods, δ [g []] = δ [ 0] (6.38) dg =0 Te SINC function as zeros
More informationInformation and Communications Security: Encryption and Information Hiding
Short Course on Information and Communications Security: Encryption and Information Hiding Tuesday, 10 March Friday, 13 March, 2015 Lecture 5: Signal Analysis Contents The complex exponential The complex
More informationWave Phenomena Physics 15c. Lecture 10 Fourier Transform
Wave Phenomena Physics 15c Lecture 10 Fourier ransform What We Did Last ime Reflection of mechanical waves Similar to reflection of electromagnetic waves Mechanical impedance is defined by For transverse/longitudinal
More informationLecture 6: Discrete Fourier Transform
Lecture 6: Discrete Fourier Transform In the previous lecture we introduced the discrete Fourier transform as given either by summations or as a matrix vector product The discrete Fourier transform of
More informationLOPE3202: Communication Systems 10/18/2017 2
By Lecturer Ahmed Wael Academic Year 2017-2018 LOPE3202: Communication Systems 10/18/2017 We need tools to build any communication system. Mathematics is our premium tool to do work with signals and systems.
More information