VTUlive.

com 1

Image Processing 10EC763

Unit-3
UNITARY TRANSFORMS
One dimensional signals

For a one dimensional sequence { f ( x), 0 x N 1} represented as a vector

T
f f (0) f (1)  f ( N 1) of size N , a transformation may be written as
N 1
g T f g (u) T (u, x) f ( x), 0 u N 1
x 0

where g (u) is the transform (or transformation) of f (x) , and T (u, x) is the so
called

Forward transformation kernel. Similarly, the inverse transform is the

relation
N 1
f ( x) I ( x, u) g (u), 0 x N 1
u 0
or written in a matrix form
1
f I g T g
where I ( x, u) is the so called inverse transformation kernel.

If
1 T
I T T

the matrix T is called unitary, and the transformation is called unitary as

well. It can be proven (how?) that the columns (or rows) of an N N unitary
matrix are orthonormal and therefore, form a complete set of basis vectors in
the N dimensional vector space.
In that case
T N 1
f T g f ( x) T (u, x) g (u )
u 0
T T
The columns of T , that is, the vectors T u T (u,0) T (u,1) T (u, N 1)
are called the basis vectors of T .

Two dimensional signals (images)

As a one dimensional signal can be represented by an orthonormal set of

basis vectors, an image can also be expanded in terms of a discrete set of
basis arrays called basis images through a two dimensional (image)
transform.

Dept of ECE/SJBIT Page 19

VTUlive.com 2

Image Processing 10EC763

For an N N image f ( x, y ) the forward and inverse transforms are given

below
N 1N 1
g (u, v) T (u, v, x, y) f ( x, y)
x 0 y 0
N 1N 1
f ( x, y) I ( x, y, u, v) g (u, v)
u 0v 0
where, again, T (u, v, x, y) and I ( x, y, u, v) are called the forward and inverse
transformation kernels, respectively.

The forward kernel is said to be separable if

T (u, v, x, y) T1 (u, x)T2 (v, y)

It is said to be symmetric if T1 is functionally equal to T2 such that

T (u, v, x, y) T1 (u, x)T1 (v, y)

The same comments are valid for the inverse kernel.

If the kernel T (u, v, x, y) of an image transform is separable and symmetric,
N 1N 1 N 1N 1
then the transform g (u, v) T (u, v, x, y) f ( x, y) T1 (u, x)T1 (v, y) f ( x, y)
x 0 y 0 x 0 y 0

can be written in matrix form as follows

T
g T1 f T1
where f is the original image of size N N , and T 1 is an N N
transformation matrix with elements tij T1 (i, j ) . If, in addition, T 1 is a
unitary matrix then the transform is called separable unitary and the original
image is recovered through the relationship

T
f T1 g T1

Fundamental properties of unitary transforms

The property of energy preservation

In the unitary transformation

g T f
1 T
it is easily proven (try the proof by using the relation T T ) that
2 2
g f

Thus, a unitary transformation preserves the signal energy. This property is

called energy preservation property.

Dept of ECE/SJBIT Page 20

VTUlive.com 3

Image Processing 10EC763

This means that every unitary transformation is simply a rotation of the

vector f in the N - dimensional vector space.

For the 2-D case the energy preservation property is written as

N 1N 1 2 N 1N 1 2
f ( x, y ) g( u, v )
x 0y 0 u 0v 0

The property of energy compaction

Most unitary transforms pack a large fraction of the energy of the image into
relatively few of the transform coefficients. This means that relatively few of
the transform coefficients have significant values and these are the coefficients
that are close to the origin (small index coefficients).
This property is very useful for compression purposes. (Why?)

THE TWO DIMENSIONAL FOURIER TRANSFORM

Continuous space and continuous frequency

The Fourier transform is extended to a function f ( x, y ) of two variables. If

f ( x, y ) is continuous and integrable and F(u, v ) is integrable, the following
Fourier transform pair exists:
j 2 ( ux vy)
F (u , v) f ( x, y )e dxdy

1
f ( x, y ) F (u , v)e j 2 ( ux vy)
dudv
(2 ) 2
In general F(u, v ) is a complex-valued function of two real frequency
variables u, v and hence, it can be written as:
F (u, v) R(u, v) jI (u, v)
The amplitude spectrum, phase spectrum and power spectrum, respectively,
are defined as follows.
F (u, v ) R2 (u, v ) I 2 (u, v )
1 I (u , v )
(u , v ) tan
R (u , v )
2
P(u, v ) F(u, v ) R2 (u, v ) I 2 (u, v )

2.2 Discrete space and continuous frequency

Dept of ECE/SJBIT Page 21

VTUlive.com 4

Image Processing 10EC763

For the case of a discrete sequence f ( x, y ) of infinite duration we can define

the 2-D discrete space Fourier transform pair as follows

j ( xu vy)
F ( u, v ) f ( x, y ) e
x y

1
f ( x, y ) F (u, v )e j ( xu vy )
dudv
(2 ) 2 u v

F(u, v ) is again a complex-valued function of two real frequency variables

u, v and it is

periodic with a period 2 2 , that is to say

F(u, v )
F(u 2 , v ) F(u, v 2 )
The Fourier transform of f ( x, y ) is said to converge uniformly when F(u, v )
is finite and
N1 N2
j ( xu vy )
lim lim f ( x, y )e F (u, v ) for all u, v .
N1 N2 x N1 y N2

When the Fourier transform of f ( x, y ) converges uniformly, F(u, v ) is an

analytic function and is infinitely differentiable with respect to u and v .

Discrete space and discrete frequency: The two dimensional Discrete Fourier
Transform (2-D DFT)

If f ( x, y ) is an M N array, such as that obtained by sampling a continuous

function of two dimensions at dimensions M and N on a rectangular grid,
then its two dimensional Discrete Fourier transform (DFT) is the array given
by

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

u 0,, M 1 , v 0,, N 1

and the inverse DFT (IDFT) is

M 1N 1
f ( x, y ) F ( u, v )e j 2 ( ux / M vy / N )

u 0 v 0

When images are sampled in a square array, M N and

1 N 1N 1
j 2 ( ux vy )/ N
F ( u, v ) f ( x, y )e
N x 0y 0

Dept of ECE/SJBIT Page 22

VTUlive.com 5

Image Processing 10EC763

1 N 1N 1
f ( x, y ) F ( u, v )e j 2 ( ux vy )/ N
N u 0v 0

It is straightforward to prove that the two dimensional Discrete Fourier

Transform is separable, symmetric and unitary.

Properties of the 2-D DFT

Most of them are straightforward extensions of the properties of the 1-D

Fourier Transform. Advise any introductory book on Image Processing.

The importance of the phase in 2-D DFT. Image reconstruction from amplitude
or phase only.

The Fourier transform of a sequence is, in general, complex-valued, and the

unique representation of a sequence in the Fourier transform domain requires
both the phase and the magnitude of the Fourier transform. In various
contexts it is often desirable to reconstruct a signal from only partial domain
information. Consider a 2-D sequence f ( x, y ) with Fourier transform
F ( u, v ) f ( x, y ) so that

j f (u , v )
F (u, v) { f ( x, y} F (u, v) e

It has been observed that a straightforward signal synthesis from the Fourier
transform phase f (u , v ) alone, often captures most of the intelligibility of the

original image f ( x, y ) (why?). A straightforward synthesis from the Fourier

transform magnitude F (u , v) alone, however, does not generally capture the

original signal’s intelligibility. The above observation is valid for a large
number of signals (or images). To illustrate this, we can synthesise the phase-
only signal f p ( x, y ) and the magnitude-only signal f m ( x, y) by

1 j f (u , v )
f p ( x, y) 1e

1
f m ( x, y ) F (u , v) e j 0

and observe the two results (Try this exercise in MATLAB).

Dept of ECE/SJBIT Page 23

VTUlive.com 6

Image Processing 10EC763

An experiment which more dramatically illustrates the observation that

phase-only signal synthesis captures more of the signal intelligibility than
magnitude-only synthesis, can be performed as follows.
Consider two images f ( x, y ) and g ( x, y) . From these two images, we
synthesise two other images f1 ( x, y ) and g1 ( x, y) by mixing the amplitudes
and phases of the original images as follows:
1 j f (u , v )
f1 ( x, y) G(u, v) e
1 j (u , v )
g1 ( x, y) F (u, v) e g
In this experiment f1 ( x, y ) captures the intelligibility of f ( x, y ) , while
g1 ( x, y) captures the intelligibility of g ( x, y) (Try this exercise in
MATLAB).

THE DISCRETE COSINE TRANSFORM (DCT)

One dimensional signals

This is a transform that is similar to the Fourier transform in the sense that the
new independent variable represents again frequency. The DCT is defined
below.
N 1 (2 x 1)u
C (u ) a(u ) f ( x) cos , u 0,1,, N 1
x 0 2N
with a(u ) a parameter that is defined below.
1/ N u 0
a (u )
2/ N u 1,, N 1
The inverse DCT (IDCT) is defined below.
N 1 (2 x 1)u
f ( x) a(u )C (u ) cos
u 0 2N

Two dimensional signals (images)

For 2-D signals it is defined as

N 1N 1 (2 x 1)u (2 y 1)v
C (u, v) a(u )a(v) f ( x, y ) cos cos
x 0 y 0 2N 2N
N 1N 1 (2 x 1)u (2 y 1)v
f ( x, y ) a(u )a(v)C (u, v) cos cos
u 0v 0 2N 2N
a(u ) is defined as above and u, v 0,1, , N 1

Properties of the DCT transform

Dept of ECE/SJBIT Page 24
VTUlive.com 7

Image Processing 10EC763

 The DCT is a real transform. This property makes it attractive in comparison

to the Fourier transform.
 The DCT has excellent energy compaction properties. For that reason it is
widely used in image compression standards (as for example JPEG
standards).
 There are fast algorithms to compute the DCT, similar to the FFT for
computing the DFT.

Recommended Questions

1. Define two-dimensional DFT. Explain the following properties of 2-DFT.

i) Translation ii) rotation iii) distributivity and scaling iv) separability
2. What are basis vectors?
3. Derive the expression for 2D circular convolution theorem.
4. Define two – dimensional unitary transform. Check whether the unitary DFT
matrix is unitary or not for N = 4.
5. For the 2 X 2 transform A and the image U
1 1 1 2
A = 1/1 1 -1 and U = 8 4

Calculate the transformed image V and the basis images.

6. Consider the image segment shown in fig
i) Let V = {0, 1}. Compute the lengths of shortest 4 - , 8 – and m – paths between p
and q.
ii) Repeat for V = {1, 2}.
3 1 2 1 (q)
3 2 0 2
( p) 1 2 1 1

Dept of ECE/SJBIT Page 25