A sampled image gives us our usual 2D array of pixels
fmn
f
m
n
(Figure 1):
We can filter
fmn
f
m
n
by applying a 2D discrete-space convolution as shown below (where
hmn
h
m
n
is our PSF):
gmn=hmn*fmn=∑k=−∞∞∑l=−∞∞hm−kn−lfkl
g
m
n
h
m
n
f
m
n
k
l
h
m
k
n
l
f
k
l
(1)
We also have discrete-space FTS:
Fuv=∑m=−∞∞∑n=−∞∞fmne−(ium)e−(ivm)
F
u
v
m
n
f
m
n
u
m
v
m
(2)
where
Fuv
F
u
v
is analogous to
DTFT in 1D.
"Convolution in Time" is the same as
"Multiplication in Frequency"
gmn=hmn*fmn
g
m
n
h
m
n
f
m
n
(3)
which, as stated above, is the same as:
Guv=HuvFuv
G
u
v
H
u
v
F
u
v
(4)
To get a better image, we can use the
fftshift command in Matlab to center
the Fourier Transform. The resulting image is shown in
Figure 4: