1 / 61

ENG4BF3 Medical Image Processing

ENG4BF3 Medical Image Processing. Image Enhancement in Frequency Domain. Image Enhancement. Original image. Enhanced image. Enhancement: to process an image for more suitable output for a specific application. Image Enhancement. Image enhancement techniques: Spatial domain methods

ebliss
Download Presentation

ENG4BF3 Medical Image Processing

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. ENG4BF3Medical Image Processing Image Enhancement in Frequency Domain

  2. Image Enhancement Original image Enhanced image Enhancement: to process an image for more suitable output for a specific application.

  3. Image Enhancement • Image enhancement techniques: • Spatial domain methods • Frequency domain methods • Spatial (time) domain techniques are techniques that operate directly on pixels. • Frequency domain techniques are based on modifying the Fourier transform of an image.

  4. Fourier Transform: a review • Basic ideas: • A periodic function can be represented by the sum of sines/cosines functions of different frequencies, multiplied by a different coefficient. • Non-periodic functions can also be represented as the integral of sines/cosines multiplied by weighing function.

  5. Joseph Fourier(1768-1830) Fourier was obsessed with the physics of heat and developed the Fourier transform theory to model heat-flow problems.

  6. Fourier transformbasis functions Approximating a square wave as the sum of sine waves.

  7. E(-x) = E(x) O(-x) = -O(x) Any function can be written as thesum of an even and an odd function

  8. Fourier Cosine Series Because cos(mt) is an even function, we can write an even function, f(t), as: where series Fm is computed as Here we suppose f(t) is over the interval (–π,π).

  9. Fourier Sine Series Because sin(mt) is an odd function, we can write any odd function, f(t), as: where the series F’m is computed as

  10. Fourier Series So if f(t) is a general function, neither even nor odd, it can be written: Even component Odd component where the Fourier series is

  11. The Fourier Transform Let F(m) incorporates both cosine and sine series coefficients, with the sine series distinguished by making it the imaginary component: Let’s now allow f(t) range from – to , we rewrite: F(u) is called the Fourier Transform of f(t). We say that f(t) lives in the “time domain,” and F(u) lives in the “frequency domain.” u is called the frequency variable.

  12. The Inverse Fourier Transform Given F(u), f(t) can be obtained by the inverse Fourier transform We go from f(t) to F(u) by Fourier Transform Inverse Fourier Transform

  13. 2-D Fourier Transform and the inverse Fourier transform Fourier transform for f(x,y) with two variables

  14. Discrete Fourier Transform (DFT) • A continuous function f(x) is discretized as:

  15. Discrete Fourier Transform (DFT) Let x denote the discrete values (x=0,1,2,…,M-1), i.e. then

  16. Discrete Fourier Transform (DFT) • The discrete Fourier transform pair that applies to sampled functions is given by: u=0,1,2,…,M-1 and x=0,1,2,…,M-1

  17. 2-D Discrete Fourier Transform • In 2-D case, the DFT pair is: u=0,1,2,…,M-1 and v=0,1,2,…,N-1 and: x=0,1,2,…,M-1 and y=0,1,2,…,N-1

  18. Polar Coordinate Representation of FT • The Fourier transform of a real function is generally complex and we use polar coordinates: Polar coordinate Magnitude: Phase:

  19. Fourier Transform: shift • It is common to multiply input image by (-1)x+y prior to computing the FT. This shift the center of the FT to (M/2,N/2). Shift

  20. Symmetry of FT • For real image f(x,y), FT is conjugate symmetric: • The magnitude of FT is symmetric:

  21. FT IFT

  22. IFT IFT

  23. The central part of FT, i.e. the low frequency components are responsible for the general gray-level appearance of an image. The high frequency components of FT are responsible for the detail information of an image.

  24. v u Image Frequency Domain (log magnitude) Detail General appearance

  25. 5 % 10 % 20 % 50 %

  26. Frequency Domain Filtering

  27. Frequency Domain Filtering • Edges and sharp transitions (e.g., noise) in an image contribute significantly to high-frequency content of FT. • Low frequency contents in the FT are responsible to the general appearance of the image over smooth areas. • Blurring (smoothing) is achieved by attenuating range of high frequency components of FT.

  28. Convolution Theorem • f(x,y) is the input image • g(x,y) is the filtered • h(x,y): impulse response Multiplication in Frequency Domain G(u,v)=F(u,v)●H(u,v) Convolution in Time Domain g(x,y)=h(x,y)*f(x,y) • Filtering in Frequency Domain with H(u,v) is equivalent to filtering in Spatial Domain with f(x,y).

  29. Gaussian lowpass filter Gaussian highpass filter Examples of Filters Frequency domain Spatial domain

  30. Ideal low-pass filter (ILPF) (M/2,N/2): center in frequency domain D0 is called the cutoff frequency.

  31. Shape of ILPF Frequency domain Spatial domain

  32. FT ringing and blurring Ideal in frequency domain means non-ideal in spatial domain, vice versa.

  33. Butterworth Lowpass Filters (BLPF) • Smooth transfer function, no sharp discontinuity, no clear cutoff frequency.

  34. Butterworth Lowpass Filters (BLPF) n=1 n=2 n=5 n=20

  35. No serious ringing artifacts

  36. Gaussian Lowpass Filters (GLPF) • Smooth transfer function, smooth impulse response, no ringing

  37. Gaussian lowpass filter GLPF Frequency domain Spatial domain

  38. No ringing artifacts

  39. Examples of Lowpass Filtering

  40. Examples of Lowpass Filtering Low-pass filter H(u,v) Original image and its FT Filtered image and its FT

  41. Sharpening High-pass Filters • Hhp(u,v)=1-Hlp(u,v) • Ideal: • Butterworth: • Gaussian:

  42. High-pass Filters

  43. Ideal High-pass Filtering ringing artifacts

  44. Butterworth High-pass Filtering

  45. Gaussian High-pass Filtering

  46. Gaussian High-pass Filtering Gaussian filter H(u,v) Original image Filtered image and its FT

  47. ¶ 2 2 f ( x , y ) f ( x , y ) Á + = - + 2 2 [ ] ( u v ) F ( u , v ) = - + 2 2 H ( u , v ) ( u v ) ¶ ¶ 2 2 x y 1 Laplacian in Frequency Domain Frequency domain Spatial domain Laplacian operator

  48. = - Ñ 2 g ( x , y ) f ( x , y ) f ( x , y ) = + + 2 2 G ( u , v ) F ( u , v ) ( u v ) F ( u , v ) = + + = - 2 2 H ( u , v ) 1 ( u v ) 1 H ( u , v ) 2 1 Subtract Laplacian from the Original Image to Enhance It enhanced image Original image Laplacian output Spatial domain Frequency domain new operator Laplacian

More Related