0% found this document useful (0 votes)
38 views

Chapter 5

The document discusses various image operations that can be performed in different domains: 1. Gray value (histogram) domain - Operations like histogram equalization and specification. 2. Spatial (image) domain - Filters like average, median, and gradient filters that operate directly on pixel values. 3. Frequency (Fourier) domain - Operations performed on the Fourier transform of an image, such as emphasis or de-emphasis of certain frequency components.

Uploaded by

Sunny
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
38 views

Chapter 5

The document discusses various image operations that can be performed in different domains: 1. Gray value (histogram) domain - Operations like histogram equalization and specification. 2. Spatial (image) domain - Filters like average, median, and gradient filters that operate directly on pixel values. 3. Frequency (Fourier) domain - Operations performed on the Fourier transform of an image, such as emphasis or de-emphasis of certain frequency components.

Uploaded by

Sunny
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 49

Image Operations in Different Domains

1) Gray value (histogram) domain


- Histogram equalization, specification, etc...
2) Spatial (image) domain
- Average filter, median filter, gradient, laplacian, etc…
3) Frequency (Fourier) domain

3 X 3 Average 5 X 5 Average
3500

3000
+
Original histogram
3500

3000
=
Equalized histogram
2500 2500

Noisy image
2000 2000

(Salt & Pepper noise) Original image f


1500 1500 Gradient magnitude
f Sharpened
f x 2  f yImage
2
1000 1000

500
Blurry Image Laplacian 500

0 0
0 50 100 150 200 250 0 50 100 150 200 250

7 X 7 Average Median
Image Enhancement
 Purposes: To make an image better appealing and
easier to deal with than the original image
 Three categories:
1. Spatial domain methods: operate on the images
itself,
 Point processing, e.g., image averaging; logic
operation; contrast stretching ...
 Mask processing, e.g., filtering or mask
operation, (blurring, median
Image Enhancement …
2. Frequency domain methods:
work on the Fourier transformed output of the image,
examples: from the convolution theory
 g(x,y) = f(x,y) h(x,y)
 => G(u,v) = F(u,v) • H(u,v)
 => certain properties of F(u,v) can be
 emphasized into G(u,v)

 => spatial domain g(x,y) = F-1{G(u,v)}

3. Combination of the above two categories


Why Do Transforms?
• Fast computation
– E.g., convolution vs. multiplication for filter with wide
support
• Conceptual insights for various image processing
– E.g., spatial frequency info. (smooth, moderate change,
fast change, etc.)
• Obtain transformed data as measurement
– E.g., blurred images, radiology images (medical and
astrophysics)
– Often need inverse transform
– May need to get assistance from other transforms
• For efficient storage and transmission
– Pick a few “representatives” (basis)
– Just store/send the “contribution” from each basis
Listening to an orchestra, you can distinguish
between different instruments, although the sound is
a SINGLE FUNCTION !

Flute

Tabla

Sitar
Fourier Series
 Any function that
periodically repeats
itself can be expressed
as the sum of sines
and/or cosines of
different frequencies,
each multiplied by a
different coefficients. This
sum is called a Fourier
series.
Joseph Fourier (1768-1830)
 Had crazy idea (1807):
 Any periodic function can be
rewritten as a weighted sum of
Sines or Cosines of different
frequencies.
 Don’t believe it?
 Neither did Lagrange,
Laplace, Poisson and other
big wigs
 Not translated into English
until 1878!
 But it’s true!
 called Fourier Series
 Possibly the greatest tool
used in Engineering
Periodic Signals
Fourier Series

y(x) = A sin(fx + p)
Any mathematical function that periodically repeats itself can be
expressed as a sum of sines &/or cosines with different amplitudes A,
frequencies f, and phases p.

It does not matter how complicated the function is – as long as it is


periodic (& meets certain mathematical conditions), it can be
represented by such a sum.

Even functions that are not periodic (but whose area under
the curve is finite) can be expressed as the integral of
sines &/or cosines multiplied by a weighting function.
A sum of sines and cosines
=

3 sin(x) A

+ 1 sin(3x) B A+B

+ 0.8 sin(5x) C
A+B+C

+ 0.4 sin(7x) D
A+B+C+D
Amplitude and Phase
 The spectrum is the set of waves representing a
signal as frequency components. It specifies for
each frequency:
 The amplitude (related to the energy)
 The phase (its „position‟ relative to other frequencies)

Amplitude Phase
Fourier Transform
 We want to understand the frequency w of our signal. So, let‟s
reparametrize the signal by w instead of x:

Fourier
f(x) F()
Transform

• For every  from 0 to inf, F() holds the amplitude A and phase of
the corresponding sine Asin(x   
– How can F hold both? Complex number trick!
F ( )  R( )  iI ( )
I ( )
A   R( )  I ( )
2 2
  tan 1

R( )
Inverse Fourier
F() f(x)
Transform
Time and Frequency

 example : g(t) = sin(2pf t) + (1/3)sin(2p(3f) t)


Time and Frequency
 example : g(t) = sin(2pf t) + (1/3)sin(2p(3f) t)

= +
Frequency Spectra
 example : g(t) = sin(2pf t) + (1/3)sin(2p(3f) t)

= +
The Continuous Fourier Transform
 Represent the signal as an infinite weighted
sum of an infinite number of sinusoids

 f (t )  F (u )   f (t )e(  j 2ut ) dt


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.
The Inverse Fourier Transform
 Given F(u), f(t) can be obtained by the inverse
Fourier transform (IFT)


1F (u )  f (t )   F (u ) e ( j 2ut )
du

AW Fourier
Transform

-1/W 1/W
The Fourier spectrum tells us:
Fourier
 The pulse is made up of infinite Spectrum
frequencies.
 The frequency component at 0 has the
maximum amplitude.
 At higher frequencies, the amplitude goes
on decreasing.
 If we add up all these frequencies (u), we
get back the original signal.
Discrete Fourier Transform (DFT)
 A continuous function f(x) is discretized as:
{ f ( x0 ), f ( x0  x), f ( x0  2x),..., f ( x0  ( M  1)x)}
Discrete Fourier Transform (DFT)
Let x denote the discrete values (x=0,1,2,…,M-1),
i.e.
f ( x)  f ( x0  xx)
then
{ f ( x0 ), f ( x0  x), f ( x0  2x),..., f ( x0  ( M  1)x)}

{ f (0), f (1), f (2),..., f ( M  1)}


Discrete Fourier Transform (DFT)
 The discrete Fourier transform pair that applies to
sampled functions is given by:

M 1
1
F (u ) 
M
 f ( x )e
x 0
(  j 2ux / M )
u=0,1,2,…,M-1

and
M 1
f ( x)   F (u )e ( j 2ux / M )
x=0,1,2,…,M-1
u 0
Example: A simple one-dimensional DFT
Polar Coordinate Representation of FT

 The Fourier transform of a real function is generally


complex and we use polar coordinates:
F (u, v)  R(u, v)  j  I (u, v)
Polar coordinate
F (u, v)  F (u, v) exp( j (u , v))
Magnitude: F (u, v)  [ R 2 (u, v)  I 2 (u, v)]1/2
 I (u, v) 
Phase:  (u, v)  tan 
1

 R (u , v ) 
P(u, v)  F (u, v)  R 2 (u, v)  I 2 (u, v)
2
Power spectrum:
How to compute F(u)?
M 1
1
F (u ) 
M

x 0
f ( x )e (  j 2ux / M ) u =[0,1,2, …, M-1]

1. Substitute u = 0 in the exponential term & then


sum for all values of x.
2. Next substitute u = 1 in the exponential & repeat
the summation over all values of x.
3. Repeat this process for all M values of u in order
to obtain the complete Fourier transform.
DFT: Example
 Find the DFT of f(x) = {0, 1, 2, 1}
 DFT Matrix Method

Twiddle Factor
A new factor defined as
DFT: Example
 The DFT now reduces to

To solve this equation, form a square matrix WM of size


MxM

The equation then reduces to

F=W.f
Here WM will be a 4x4 matrix
2-D Discrete Fourier Transform …

What do we mean by frequency in an


image?

Low frequency region

High frequency region

Low frequency region


Why is FT Useful?
• Remove undesirable frequencies from a signal.
• Easier and faster to perform certain operations in
the frequency domain than in the spatial domain.
 g(x,y)= f(x,y) h(x,y)
 => G(u,v) = F(u,v) • H(u,v)
 => certain properties of F(u,v) can be
 emphasized into G(u,v)

 => spatial domain g(x,y) = F-1{G(u,v)}


Properties of DFT

 The Separability Property


 A 2-D DFT can be separated into
two 1-D DFT
Properties of DFT …

We can rearrange the terms as follows


(0,0) y (0,0) v
Row
Transformation
f(x,y) F(x,v)

x x

(0,0) v

F(u,v) Column
Transformation
u
Properties of DFT …

 The Translation Property


If f(x,y) is multiplied by an
exponential, the original Fourier
Transform, F(u,v), gets shifted in
frequency by f(u-u0, v-v0)
= F(u-u0, v-v0)

Thus, if

Then,
What happens when
u0 = M/2 & v0 = N/2

Why?
For an M x N image, what are the
coordinates of the centre pixel?

Since (x + y) is an integer
So sin(x+y) = 0 & cos(x+y) is
1 if (x+y) is even & -1 if it is odd
F(0, 0) is at the
origin

F(M/2, N/2) is
at the origin
A common
practice …

Multiply the i/p image function by (-1)x + y prior to


computing the Fourier transform

This shifts the origin of F(u, v) to frequency coordinates


(M/2, N/2),which is the centre of the M x N area occupied
by the 2-DFT
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).
 f ( x, y )  F (u, v)
 f ( x, y )(1) x  y   F (u  M / 2, v  N / 2)

Shift
The Average Value Property

What is the average of all the intensity values of


the image?
The Average Value Property …

The value of transform at (u, v) = (0, 0)

Which is the average of f(x, y)


f(x, y) is an image
The value of the Fourier transform at origin is equal to
the average gray-level of the image
2-D Discrete Fourier Transform
 In 2-D case, the DFT pair is:

1 M 1 N 1
F (u, v)  
MN x 0 y 0
f ( x, y ) e (  j 2 ( ux / M  vy / N ))

u=0,1,2,…,M-1 and v=0,1,2,…,N-1


and:
M 1 N 1
f ( x, y )   F (u, v)e ( j 2 (ux / M  vy / N ))
u 0 v 0

x=0,1,2,…,M-1 and y=0,1,2,…,N-1


Problems
Q) Find the DFT of the following image.

0 1 2 1
1 2 3 2
2 3 4 3
1 2 3 2

Hint:
First do row-wise transformation & then column-wise.
Problems …
 Find the DFT of f(x)= {0,1,2,1}
Assignment …
Any Questions ?

You might also like