4.3 平滑的频率域滤波器
- 边缘和其他尖锐变化(如噪声)在图像的灰度级中主要处于傅立叶变换的高频部分。因此,平滑(模糊)可以通过衰减指定图像傅立叶变换中高频成分的范围来实现。
G ( u , v ) = H ( u , v ) F ( u , v ) G(u,v) = H(u,v)F(u,v) G(u,v)=H(u,v)F(u,v)
目标是选择一个滤波器变换函数 H ( u , v ) H(u,v) H(u,v)以通过衰减 F ( u , v ) F(u,v) F(u,v)的高频成分产生 G ( u , v ) G(u,v) G(u,v)。
三种滤波器:
- 理想低通滤波器
- 巴特沃思低通滤波器
- 高斯低通滤波器
理想低通滤波器
- 最简单的低通滤波器是“截断”傅立叶变换中所有高频成分,这些成分处在距变换原点的距离比指定距离 D 0 D_0 D0要远得多的位置。这种滤波器称为二维理想低通滤波器,其变换函数为:
其中,
D
0
D_0
D0(截止频率)是指定的非负数值,
D
(
u
,
v
)
D(u,v)
D(u,v)是
(
u
,
v
)
(u,v)
(u,v)点距离频率矩形原点的距离。由于变换被中心化了,如果要研究的图像尺寸为
M
×
N
M \times N
M×N,从点
(
u
,
v
)
(u,v)
(u,v)到傅立叶变换中心(原点)的距离如下所示:
D ( u , v ) = [ ( u − M / 2 ) 2 + ( v − N / 2 ) 2 ] 1 / 2 D(u,v) = [(u - M/2)^2 + (v - N/2)^2]^{1/2} D(u,v)=[(u−M/2)2+(v−N/2)2]1/2
- 为了研究理想低通滤波器的特性,建立一组具有不同截止频率的理想滤波器。建立方法是衡量截止频率以内的功率与整个图像的功率的比值。
总功率是将每个点
(
u
,
v
)
(u,v)
(u,v)的功率谱成分相加得到的,其中
u
=
0
,
1
,
2
,
.
.
.
,
M
−
1
,
v
=
0
,
1
,
1
,
.
.
.
,
N
−
1
u=0,1,2,...,M-1,v=0,1,1,...,N-1
u=0,1,2,...,M−1,v=0,1,1,...,N−1,即:
P
(
r
)
=
∑
u
=
0
M
−
1
∑
v
=
0
N
−
1
P
(
u
,
v
)
P(r) = \sum_{u=0}^{M-1}\sum_{v=0}^{N-1}P(u,v)
P(r)=∑u=0M−1∑v=0N−1P(u,v)
如果变换被中心化,原点在频率矩形的中心,半径为
r
r
r的圆包含
a
%
a\%
a%的功率,其中:
a
=
100
[
∑
u
∑
v
P
(
u
,
v
)
/
P
r
]
a = 100[\sum_{u}\sum_{v}P(u,v)/P_r]
a=100[∑u∑vP(u,v)/Pr]
理想滤波器的截止频率越低,那么中心成分的半径就越大,相当于在做空间域的线性平滑滤波时使用的模板越大,就意味着图像更加模糊。
巴特沃思低通滤波器
- 由第三个图可知, n n n越大,图像越陡峭。由此猜测,当 n n n比较大的时候,巴特沃思滤波器类似于理想低通滤波器的特性。
高斯低通滤波器
完全没有振铃
二维高斯低通滤波器(GLPF)的形式由下式给出:
H
(
u
,
v
)
=
e
−
D
2
(
u
,
v
)
/
2
D
0
2
H(u,v) = e^{-D^2(u,v)/2D_0^2}
H(u,v)=e−D2(u,v)/2D02
其中,
D
0
D_0
D0是截止频率。当
D
(
u
,
v
)
=
D
0
D(u,v) = D_0
D(u,v)=D0时,滤波器下降到它最大值的0.607处。
4.4 频率域锐化滤波器
- 灰度级的边缘和其他地方的急剧变与高频成分有关,图像的锐化能够在频率域用高通滤波器处理实现,而衰减低频成分并不会扰乱傅立叶变换的高频信息。