2D Discrete Fourier Transform (DFT)

Outline
• Circular and linear convolutions
• 2D DFT
• 2D DCT
• Properties
• Other formulations
• Examples

2
 Circular convolution
• Finite length signals (N0 samples) → circular or periodic convolution
N 0 −1
– the summation is over 1 period
– the result is a N0 period sequence
c[k ] = f [k ] ⊗ g[k ] = ∑
n=0
f [ n] g[ k − n]

• The circular convolution is equivalent to the linear convolution of the

zero-padded equal length sequences

f [m]* g[m] ⇔ F [k ]G[k ]

f [ m] g [ m] f [m]* g[m]

* =
m m m
Length=P Length=Q Length=P+Q-1

For the convolution property to hold, M must be greater than or equal to P+Q-1.
3
 Convolution
f [m]* g[m] ⇔ F [k ]G[k ]
• Zero padding

f [ m] g [ m] f [m]* g[m]

* =
m m m

4-point DFT
(M=4)

F [k ] G[k ] F [k ]G[k ]

4
 In words
• Given 2 sequences of length N and M, let y[k] be their linear convolution
+∞
y[k ] = f [k ] ∗ h[k ] = ∑
n =−∞
f [n]h[ k − n]

• y[k] is also equal to the circular convolution of the two suitably zero padded
sequences making them consist of the same number of samples
N 0 −1
c[k ] = f [k ] ⊗ h[k ] = ∑
n=0
f [ n]h[ k − n]
N 0 = N f +N h − 1: length of the zero-padded seq

• In this way, the linear convolution between two sequences having a different length
(filtering) can be computed by the DFT (which rests on the circular convolution)
– The procedure is the following
• Pad f[n] with Nh-1 zeros and h[n] with Nf-1 zeros
• Find Y[r] as the product of F[r] and H[r] (which are the DFTs of the corresponding zero-padded
signals)
• Find the inverse DFT of Y[r]

• Allows to perform linear filtering using DFT

5
 2D Discrete Fourier Transform
• Fourier transform of a 2D signal defined over a discrete finite 2D grid
of size MxN
or equivalently
• Fourier transform of a 2D set of samples forming a bidimensional
sequence
• As in the 1D case, 2D-DFT, though a self-consistent transform, can
be considered as a mean of calculating the transform of a 2D
sampled signal defined over a discrete grid.
• The signal is periodized along both dimensions and the 2D-DFT can
be regarded as a sampled version of the 2D DTFT

6
 2D Discrete Fourier Transform (2D DFT)

• 2D Fourier (discrete time) Transform (DTFT) [Gonzalez]

∞ ∞
F (u , v) = ∑∑
m =−∞ n =−∞
f [m, n]e − j 2π (um + vn ) a-periodic signal
periodic transform

• 2D Discrete Fourier Transform (DFT)

M −1 N −1 ⎛ k l ⎞
1 − j 2π ⎜ m + n ⎟
∑∑
periodized signal
F [k , l ] = f [m, n]e ⎝M N ⎠
periodic and sampled
MN m=0 n =0
transform

2D DFT can be regarded as a sampled version of 2D DTFT.

7
 2D DFT: Periodicity

• A [M,N] point DFT is periodic with period [M,N]

– Proof

M −1 N −1 ⎛ k l ⎞
1 − j 2π ⎜ m + n ⎟
F [k , l ] =
MN
∑ ∑ f [m, n]e
m=0 n =0
⎝M N ⎠

M −1 N −1 ⎛ k +M l+N ⎞
1 − j 2π ⎜ m+ n⎟
F [k + M , l + N ] =
MN
∑ ∑ f [m, n]e
m=0 n =0
⎝ M N ⎠
1
M −1 N −1 ⎛ k l ⎞ ⎛M N ⎞
1 − j 2π ⎜ m + n ⎟ − j 2π ⎜ m + n ⎟
=
MN
∑ ∑ f [m, n]e
m=0 n=0
⎝M N ⎠
e ⎝M N ⎠

= F [k , l ]

(In what follows: spatial coordinates=k,l, frequency coordinates: u,v)

8
 2D DFT: Periodicity
• Periodicity
F [u , v] = F [u + mM , v] = F [u , v + nN ] = F [u + mM , v + nN ]

f [ k , l ] = f [ k + mM , l ] = f [k , l + nN ] = f [k + mM , l + nN ]

• This has important consequences on the implementation and energy

compaction property
– 1D The two inverted periods meet here
∗ f[u]
F [ N − u ] = F [u ]

f[k] real→F[u] is symmetric

M/2 samples are enough

M/2 M u
0

9
 Periodicity: 1D
f [ k ] ↔ F [u ]
u0 k
j 2π changing the sign of every other
f [ k ]e M
↔ F [u − u0 ] sample puts F[0] at the center of the
u k Mk interval [0,M]
M j 2π 0 j 2π
u0 = →e M
= e 2 M = e jπ k = (−1) k
2
M
(−1) k f [k ] ↔ F [u − ]
2
f[u] The two inverted periods meet here

M/2 M u
0
It is more practical to have one complete period positioned in [0, M-1]

10
 Periodicity: 2D

DFT periods

N/2 MxN values

-M/2

(0,0) M/2
4 inverted
periods meet F[u,v]
here
-N/2

11
 Periodicity: 2D data contain one centered
F[u,v] complete period
N-1
j 2π (
u0 k v0 l
+ )
DFT periods
f [ k , l ]e M N
↔ F [u − u0 , v − v0 ]
M N
u0 = , v0 = → MxN values
2 2 N/2
⎡ M N⎤
(−1) k + l f [ k ] ↔ F ⎢u − , v − ⎥
⎣ 2 2⎦
M/2 M-1
(0,0)

4 inverted
periods meet
here

12
 Periodicity: 2D
F[u,v]
N-1

M/2

(0,0) M-1

N/2
4 inverted
periods meet
here

13
 Periodicity in spatial domain

• [M,N] point inverse DFT is periodic with period [M,N]

M −1 N −1 ⎛ k l ⎞
j 2π ⎜ m + n ⎟
f [m, n] = ∑ ∑ F [k , l ]e ⎝M N ⎠

k =0 l =0

M −1 N −1 ⎛ k l ⎞
j 2π ⎜ ( m + M ) + ( n + N ) ⎟
f [m + M , n + N ] = ∑ ∑ F [k , l ]e ⎝M N ⎠

k =0 l =0 1
M −1 N −1 ⎛ k l ⎞ ⎛ k l ⎞
j 2π ⎜ m + n ⎟ j 2π ⎜ M + N ⎟
= ∑ ∑ F [k , l ]e ⎝M N ⎠
e ⎝M N ⎠

k =0 l =0

= f [m, n]

14
 Angle and phase spectra
F [ u , v ] = F [ u , v ] e j Φ [u , v ]

F [u , v ] = ⎡ Re { F [u , v ]} + Im { F [u , v ]} ⎤
2 2 1/ 2
modulus (amplitude spectrum)
⎣ ⎦
⎡ Im { F [u , v ]} ⎤
Φ [u , v ] = arctan ⎢ ⎥ phase
⎣ { [ ]} ⎦
Re F u , v
P[u , v] = F [u , v ]
2 power spectrum

For a real function

F [−u , −v] = F ∗ [u , v] conjugate symmetric with respect to the origin

F [−u , −v] = F [u , v]
Φ[−u , −v] = −Φ[u , v]

15
 Translation and rotation

⎛m n ⎞
j 2π ⎜ k + l ⎟
f [ k , l ]e ⎝M N ⎠
↔ F [u − m, v − l ]
⎛m n ⎞
f [ k − m, l − n ] ↔ F [u , v ]
− j 2π ⎜ k + l ⎟
⎝M N ⎠

⎧k = r cos ϑ ⎧u = ω cos ϕ
⎨ ⎨
⎩l = r sin ϑ ⎩l = ω sin ϕ
f [ r , ϑ + ϑ0 ] ↔ F [ω , ϕ + ϑ0 ]

Rotations in spatial domain correspond equal rotations in Fourier domain

16
 mean value

N −1 M −1
1
F [ 0, 0] =
NM
∑ ∑ f [ n, m ]
n=0 m=0
DC coefficient

17
 Separability
• The discrete two-dimensional Fourier transform of an image array is
defined in series form as
M −1 N −1 ⎛ k l ⎞
1 − j 2π ⎜ m + n ⎟
F [k , l ] =
MN
∑ ∑ f [m, n]e
m=0 n =0
⎝M N ⎠

• inverse transform
M −1 N −1 ⎛ k l ⎞
j 2π ⎜ m + n ⎟
f [m, n] = ∑ ∑ F [k , l ]e ⎝M N ⎠

k =0 l =0
• Because the transform kernels are separable and symmetric, the two
dimensional transforms can be computed as sequential row and column
one-dimensional transforms.
• The basis functions of the transform are complex exponentials that may be
decomposed into sine and cosine components.

18
2D DFT: summary

19
2D DFT: summary

20
2D DFT: summary

21
2D DFT: summary

22
other formulations
 2D Discrete Fourier Transform

• 2D Discrete Fourier Transform (DFT)

M −1 N −1 ⎛ k l ⎞
1 − j 2π ⎜ m + n ⎟
F [k , l ] =
MN
∑ ∑ f [m, n]e
m=0 n =0
⎝M N ⎠

where l = 0,1,..., N − 1
k = 0,1,..., M − 1

• Inverse DFT
M −1 N −1 ⎛ k l ⎞
j 2π ⎜ m + n ⎟
f [m, n] = ∑ ∑ F [k , l ]e ⎝M N ⎠

k =0 l =0

24
 2D Discrete Fourier Transform

• It is also possible to define DFT as follows

M −1 N −1 ⎛ k l ⎞
1 − j 2π ⎜ m + n ⎟
F [k , l ] =
MN
∑ ∑ f [m, n]e
m=0 n=0
⎝M N ⎠

where k = 0,1,..., M − 1
l = 0,1,..., N − 1

• Inverse DFT
M −1 N −1 ⎛ k l ⎞
1 j 2π ⎜ m + n ⎟
f [m, n] =
MN
∑ ∑ F [k , l ]e
k =0 l =0
⎝M N ⎠

25
 2D Discrete Fourier Transform

• Or, as follows

M −1 N −1 ⎛ k l ⎞
− j 2π ⎜ m + n ⎟
F [k , l ] = ∑ ∑ f [m, n]e ⎝M N ⎠

m=0 n =0

where k = 0,1,..., M − 1 and l = 0,1,..., N − 1

• Inverse DFT
M −1 N −1 ⎛ k l ⎞
1 j 2π ⎜ m + n ⎟
f [m, n] =
MN
∑ ∑ F [k , l ]e
k =0 l =0
⎝M N ⎠

26
 2D DFT
• The discrete two-dimensional Fourier transform of an image array is
defined in series form as

• inverse transform

27
 2D DCT

Discrete Cosine Transform

 2D DCT
• based on most common form for 1D DCT

u,x=0,1,…, N-1

“mean” value

29
 1D basis functions
Figure 1

Cosine basis functions are orthogonal

30
 2D DCT
• Corresponding 2D formulation

direct

u,v=0,1,…., N-1

inverse

31
 2D basis functions
• The 2-D basis functions can be generated by multiplying the
horizontally oriented 1-D basis functions (shown in Figure 1) with
vertically oriented set of the same functions.
• The basis functions for N = 8 are shown in Figure 2.
– The basis functions exhibit a progressive increase in frequency both in
the vertical and horizontal direction.
– The top left basis function assumes a constant value and is referred to
as the DC coefficient.

32
 2D DCT basis functions
Figure 2

33
 Separability

The inverse of a multi-dimensional DCT is just a separable product of the inverse(s) of the
corresponding one-dimensional DCT , e.g. the one-dimensional inverses applied along one
dimension at a time

34
 Separability
• Symmetry
– Another look at the row and column operations reveals that these
operations are functionally identical. Such a transformation is called a
symmetric transformation.
– A separable and symmetric transform can be expressed in the form

T = AfA
– where A is a NxN symmetric transformation matrix which entries a(i,j)
are given by

• This is an extremely useful property since it implies that the transformation

matrix can be pre computed offline and then applied to the image thereby
providing orders of magnitude improvement in computation efficiency.

35
 Computational efficiency
• Computational efficiency
– Inverse transform
– DCT basis functions are orthogonal. Thus, the inverse transformation
matrix of A is equal to its transpose i.e. A-1= AT. This property renders
some reduction in the pre-computation complexity.

36
 Block-based implementation
Basis function
Block-based transform

Block size
N=M=8

The source data (8x8) is transformed to a

linear combination of these 64 frequency
squares.

37
Energy compaction

38
Energy compaction

39
 Appendix
• Eulero’s formula

40