Circular convolution for 1D signals (Figure 1).
yn=∑k=0N-1hn-kmodNxk
y
n
k
0
N
1
h
n
k
N
x
k
(1)
for
0≤n≤N-1
0
n
N
1
.
- Flip hh
around
- Shift this function
- Multiply by xx
- Add up to get
yn
y
n
- Repeat for each nn
The domain has two dimensions.
Images
(
N×N
N
N
box of numbers), Figure 2.
Figure 3.
When
ℋℋ is
LSI, it is completely determined by its
impulse
response:
h∈ℝN2
h
N
2
hmn
h
m
n
Compute output yy via 2D circular convolution (Figure 4).
2D circular convolution of
x
⊛
N
h
x
⊛
N
h
(each
N×N
N
N
image), Figure 5.
ymn=∑k=0N-1∑l=0N-1hm-lmodNn-kmodNxlk
y
m
n
k
0
N
1
l
0
N
1
h
m
l
N
n
k
N
x
l
k
(2)
Same procedure as 1D:
flip; shift; multiply; add up; repeat.
All vector space theory goes
through to 2D images and general
dd-dimensional functions.
∥x∥pp=∑m=0N-1∑n=0N-1|xmn|p
p
x
p
m
0
N
1
n
0
N
1
x
m
n
p
(3)
<x,y>=∑m=0N-1∑n=0N-1xmnymn¯
x
y
m
0
N
1
n
0
N
1
x
m
n
y
m
n
(4)
<x,y>≤∥x∥2∥y∥2
x
y
2
x
2
y
(5)
(where
Equation 5 is CSI). There exist 2D ONB's
etc...
You will learn more on this in Elec 439.
Much more prevalant
use of nonlinear filters on images.