0 前言
本文记录了笔者数字图像处理作业的一道二维离散傅里叶变换的平移性质的证明题,希望能够帮助到读者。
1 证明
f
(
x
,
y
)
e
j
2
π
(
u
0
x
M
+
v
0
y
N
)
<
=
>
F
(
u
−
u
0
,
v
−
v
0
)
f
(
x
−
x
0
,
y
−
y
0
)
<
=
>
F
(
u
,
v
)
e
−
j
2
π
(
u
0
x
M
+
u
0
y
N
)
\begin{aligned}f(x,y)e^{j2\pi(\frac{u_0x}M+\frac{v_0y}N)}&<=>F(u-u_0,v-v_0)\\f(x-x_0,y-y_0)&<=>F(u,v)e^{-j2\pi(\frac{u_0x}M+\frac{u_0y}N)}\end{aligned}
f(x,y)ej2π(Mu0x+Nv0y)f(x−x0,y−y0)<=>F(u−u0,v−v0)<=>F(u,v)e−j2π(Mu0x+Nu0y)
也就是说,用指数项乘以
f
(
x
,
y
)
f(x,y)
f(x,y) 将使DFT的原点移到
(
u
0
,
v
0
)
(u_0,v_0)
(u0,v0);反之,用负指数乘以
F
(
u
,
v
)
F(u,v)
F(u,v) 将使
f
(
x
,
y
)
f(x,y)
f(x,y) 的原点平移到点
(
x
0
,
y
0
)
(x_0,y_0)
(x0,y0)。即图像平移不影响幅度谱,但影响相位谱。
1.1 空间位移
f ( x − x 0 , y − y 0 ) < = > F ( u ′ , v ′ ) f(x-x_0,y-y_0)<=>F(u^{\prime},v^{\prime}) f(x−x0,y−y0)<=>F(u′,v′)
F ( u ′ , v ′ ) = 1 M N ∑ x = 0 M − 1 ∑ y = 0 N − 1 f ( x − x 0 , y − y 0 ) e − j 2 π ( u x M + v y N ) F(u^{\prime},v^{\prime})=\frac1{\sqrt{MN}}\sum_{x=0}^{M-1}\sum_{y=0}^{N-1}f(x-x_0,y-y_0)e^{-j2\pi(\frac{\mathrm{u}x}M+\frac{\mathrm{v}y}N)} F(u′,v′)=MN1x=0∑M−1y=0∑N−1f(x−x0,y−y0)e−j2π(Mux+Nvy)
令
x
−
x
0
=
a
,
y
−
y
0
=
b
x-x_0 = a,y-y_0=b
x−x0=a,y−y0=b
则
F
(
u
′
,
v
′
)
=
1
M
N
∑
x
=
0
M
−
1
∑
y
=
0
N
−
1
f
(
a
,
b
)
e
−
j
2
π
(
u
(
a
+
x
0
)
M
+
ν
(
b
+
y
0
)
N
)
=
e
−
j
2
π
(
u
x
0
M
+
v
y
0
N
)
⋅
1
M
N
∑
x
=
0
M
−
1
∑
y
=
0
N
−
1
f
(
a
,
b
)
e
−
j
2
π
(
u
a
M
+
v
b
N
)
=
e
−
j
2
π
(
u
x
0
M
+
v
y
0
N
)
F
(
u
,
v
)
\begin{aligned} F(u',v')& \begin{aligned}=\frac{1}{\sqrt{MN}}\sum_{x=0}^{M-1}\sum_{y=0}^{N-1}f(a,b)e^{-j2\pi(\frac{u(a+x_0)}M+\frac{\nu(b+y_0)}N)}\end{aligned} \\ &=e^{-j2\pi(\frac{ux_0}M+\frac{vy_0}N)}\cdot\frac1{\sqrt{MN}}\sum_{x=0}^{M-1}\sum_{y=0}^{N-1}f(a,b)e^{-j2\pi(\frac{ua}M+\frac{vb}N)} \\ &=e^{-j2\pi(\frac{ux_{0}}{M}+\frac{vy_{0}}{N})}F(u,v) \end{aligned}
F(u′,v′)=MN1x=0∑M−1y=0∑N−1f(a,b)e−j2π(Mu(a+x0)+Nν(b+y0))=e−j2π(Mux0+Nvy0)⋅MN1x=0∑M−1y=0∑N−1f(a,b)e−j2π(Mua+Nvb)=e−j2π(Mux0+Nvy0)F(u,v)
即
f
(
x
−
x
0
,
y
−
y
0
)
<
=
>
F
(
u
,
v
)
⋅
e
−
j
2
π
(
u
0
x
M
+
u
0
y
N
)
f(x-x_0,y-y_0)<=>F(u,v)\cdot e^{-j2\pi(\frac{u_0x}M+\frac{u_0y}N)}
f(x−x0,y−y0)<=>F(u,v)⋅e−j2π(Mu0x+Nu0y)
1.2 频率位移
F
(
u
−
u
0
,
v
−
v
0
)
<
=
>
f
(
x
′
,
y
′
)
F(u-u_0,v-v_0)<=>f(x^{\prime},y^{\prime})
F(u−u0,v−v0)<=>f(x′,y′)
f
(
x
′
,
y
′
)
=
1
M
N
∑
x
=
0
M
−
1
∑
y
=
0
N
−
1
F
(
u
−
u
0
,
v
−
v
0
)
e
j
2
π
(
u
x
M
+
v
y
N
)
f(x^{\prime},y^{\prime})=\frac1{\sqrt{MN}}\sum_{x=0}^{M-1}\sum_{y=0}^{N-1}F(u-u_0,v-v_0)e^{j2\pi(\frac{\mathrm{u}x}M+\frac{\mathrm{v}y}N)}
f(x′,y′)=MN1x=0∑M−1y=0∑N−1F(u−u0,v−v0)ej2π(Mux+Nvy)
令
u
−
u
0
=
a
,
v
−
v
0
=
b
u-u_0=a,v-v_0=b
u−u0=a,v−v0=b
则
f
(
x
′
,
y
′
)
=
1
M
N
∑
x
=
0
M
−
1
∑
y
=
0
N
−
1
f
(
a
,
b
)
e
j
2
π
(
(
a
+
u
0
)
x
M
+
(
b
+
v
0
)
y
N
)
=
e
j
2
π
(
u
x
0
M
+
ν
y
0
N
)
⋅
1
M
N
∑
x
=
0
M
−
1
∑
y
=
0
N
−
1
f
(
a
,
b
)
e
j
2
π
(
a
x
M
+
b
y
N
)
=
e
j
2
π
(
u
x
0
M
+
v
y
0
N
)
f
(
x
,
y
)
\begin{aligned} f(x^{\prime},y^{\prime})& =\frac1{\sqrt{MN}}\sum_{x=0}^{M-1}\sum_{y=0}^{N-1}f(a,b)e^{j2\pi(\frac{(a+u_0)x}M+\frac{(b+v_0)y}N)} \\ &=e^{j2\pi(\frac{ux_0}M+\frac{\nu y_0}N)}\cdot\frac1{\sqrt{MN}}\sum_{x=0}^{M-1}\sum_{y=0}^{N-1}f(a,b)e^{j2\pi(\frac{ax}M+\frac{by}N)} \\ &=e^{j2\pi(\frac{ux_{0}}{M}+\frac{vy_{0}}{N})}f(x,y) \end{aligned}
f(x′,y′)=MN1x=0∑M−1y=0∑N−1f(a,b)ej2π(M(a+u0)x+N(b+v0)y)=ej2π(Mux0+Nνy0)⋅MN1x=0∑M−1y=0∑N−1f(a,b)ej2π(Max+Nby)=ej2π(Mux0+Nvy0)f(x,y)
即
F
(
u
−
u
0
,
v
−
v
0
)
<
=
>
f
(
x
,
y
)
⋅
e
j
2
π
(
u
0
x
M
+
u
0
y
N
)
F(u-u_0,v-v_0)<=>f(x,y)\cdot e^{j2\pi(\frac{u_0x}M+\frac{u_0y}N)}
F(u−u0,v−v0)<=>f(x,y)⋅ej2π(Mu0x+Nu0y)
综上,得证!
2 总结
以上就是笔者要分享的内容了,希望可以帮助到您,倘若大佬发现问题,直接评论区斧正即可,我们下期再见!🌞🌞🌞