1 / 57

Fourier Transform (Chapter 4)

Fourier Transform (Chapter 4). CS474/674 – Prof. Bebis. Mathematical Background: Complex Numbers. A complex number x is of the form: α : real part , b: imaginary part Addition: Multiplication:. Mathematical Background: Complex Numbers (cont’d).

jana
Download Presentation

Fourier Transform (Chapter 4)

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. Fourier Transform (Chapter 4) CS474/674 – Prof. Bebis

  2. Mathematical Background:Complex Numbers • A complex number x is of the form: α: real part, b: imaginary part • Addition: • Multiplication:

  3. Mathematical Background:Complex Numbers (cont’d) • Magnitude-Phase (i.e.,vector) representation Magnitude: Phase: φ Magnitude-Phase notation:

  4. Mathematical Background:Complex Numbers (cont’d) • Multiplication using magnitude-phase representation • Complex conjugate • Properties

  5. Mathematical Background:Complex Numbers (cont’d) • Euler’s formula • Properties j

  6. Mathematical Background:Sine and Cosine Functions • Periodic functions • General form of sine and cosine functions:

  7. Mathematical Background:Sine and Cosine Functions Special case: A=1, b=0, α=1 3π/2 π π/2 π 3π/2 π/2

  8. Mathematical Background:Sine and Cosine Functions (cont’d) • Shifting or translating the sine function by a const b Note: cosine is a shifted sine function:

  9. Mathematical Background:Sine and Cosine Functions (cont’d) • Changing the amplitude A

  10. Mathematical Background:Sine and Cosine Functions (cont’d) • Changing the period T=2π/|α| consider A=1, b=0: y=cos(αt) α =4 period 2π/4=π/2 shorter period higher frequency (i.e., oscillates faster) Frequency is defined as f=1/T Alternative notation: sin(αt)=sin(2πt/T)=sin(2πft)

  11. Basis Functions • Given a vector space of functions, S, then if any f(t) ϵ S can be expressed as the set of functions φk(t) are called the expansion set of S. • If the expansion is unique, the set φk(t) is a basis.

  12. Image Transforms • Many times, image processing tasks are best performed in a domain other than the spatial domain. • Key steps: (1) Transform the image (2) Carry the task(s) in the transformed domain. (3) Apply inverse transform to return to the spatial domain.

  13. Transformation Kernels forward transformation kernel • Forward Transformation • Inverse Transformation inverse transformation kernel

  14. Kernel Properties • A kernel is said to be separable if: • A kernel is said to be symmetric if:

  15. Notation • Continuous Fourier Transform (FT) • Discrete Fourier Transform (DFT) • Fast Fourier Transform (FFT)

  16. Fourier Series Theorem • Any periodic function f(t) can be expressed as a weighted sum (infinite) of sine and cosine functions of varying frequency: is called the “fundamental frequency”

  17. Fourier Series (cont’d) α1 α2 α3

  18. Continuous Fourier Transform (FT) • Transforms a signal (i.e., function) from the spatial(x) domain to the frequency (u) domain. (IFT) where

  19. Why is FT Useful? • Easier to remove undesirable frequencies. • Faster perform certain operations in the frequency domain than in the spatial domain.

  20. Example: Removing undesirable frequencies frequencies noisy signal remove high frequencies reconstructed signal To remove certain frequencies, set their corresponding F(u) coefficients to zero!

  21. How do frequencies show up in an image? • Low frequencies correspond to slowly varying information (e.g., continuous surface). • High frequencies correspond to quickly varying information (e.g., edges) Original Image Low-passed

  22. Example of noise reduction using FT

  23. Frequency Filtering Steps 1. Take the FT of f(x): 2. Remove undesired frequencies: 3. Convert back to a signal: We’ll talk more about this later .....

  24. Definitions • F(u) is a complex function: • Magnitude of FT (spectrum): • Phase of FT: • Magnitude-Phase representation: • Power of f(x): P(u)=|F(u)|2=

  25. Example: rectangular pulse magnitude rect(x) function sinc(x)=sin(x)/x

  26. Example: impulse or “delta” function • Definition of delta function: • Properties:

  27. 1 x u Example: impulse or “delta” function (cont’d) • FT of delta function:

  28. Example: spatial/frequency shifts Special Cases:

  29. Example: sine and cosine functions • FT of the cosine function cos(2πu0x) F(u) 1/2

  30. Example: sine and cosine functions (cont’d) • FT of the sine function -jF(u) sin(2πu0x)

  31. Extending FT in 2D • Forward FT • Inverse FT

  32. Example: 2D rectangle function • FT of 2D rectangle function 2D sinc()

  33. Discrete Fourier Transform (DFT)

  34. Discrete Fourier Transform (DFT) (cont’d) • Forward DFT • Inverse DFT 1/NΔx

  35. Example

  36. Extending DFT to 2D • Assume that f(x,y) is M x N. • Forward DFT • Inverse DFT:

  37. Extending DFT to 2D (cont’d) • Special case: f(x,y) is N x N. • Forward DFT • Inverse DFT u,v = 0,1,2, …, N-1 x,y = 0,1,2, …, N-1

  38. Extending DFT to 2D (cont’d) 2D cos/sin functions

  39. Visualizing DFT • Typically, we visualize |F(u,v)| • The dynamic range of |F(u,v)| is typically very large • Apply streching: (c is const) |D(u,v)| |F(u,v)| original image before stretching after stretching

  40. DFT Properties: (1) Separability • The 2D DFT can be computed using 1D transforms only: Forward DFT: kernel is separable:

  41. DFT Properties: (1) Separability (cont’d) • Rewrite F(u,v) as follows: • Let’s set: • Then:

  42. ) DFT Properties: (1) Separability (cont’d) • How can we compute F(x,v)? • How can we compute F(u,v)? N x DFT of rows of f(x,y) DFT of cols of F(x,v)

  43. DFT Properties: (1) Separability (cont’d)

  44. DFT Properties: (2) Periodicity • The DFT and its inverse are periodic with period N

  45. DFT Properties: (3) Symmetry • If f(x,y) is real, then (see Table 4.1 for more properties)

  46. More symmetry properties

  47. f(x,y) F(u,v) ) N DFT Properties: (4) Translation • Translation in spatial domain: • Translation in frequency domain:

  48. DFT Properties: (4) Translation (cont’d) • Warning: to show a full period, we need to translate the origin of the transform at u=N/2 (or at (N/2,N/2) in 2D) |F(u)| |F(u-N/2)|

  49. ) ) N N DFT Properties: (4) Translation (cont’d) • To move F(u,v) at (N/2, N/2), take

  50. DFT Properties: (4) Translation (cont’d) no translation after translation

More Related