首页 理论教育频域增强与同态滤波的视觉测量技术

频域增强与同态滤波的视觉测量技术

【摘要】:巴特沃思高通滤波具有截止频率D0的n阶巴特沃思高通滤波器的传递函数为式中,D(u,v)仍由式或式定义。图3.24 各种高通滤波器的特性曲线2.同态滤波同态滤波一种增强图像对比度的频域处理方法,尤其适合于处理由于光照不均导致的图像不清。同态滤波以式所给的图像模型为基础,在频域中同时进行图像对比度增强和压缩图像亮度范围。图3.28 同态滤波增强

频域增强通过对图像做傅里叶变换,将图像变换到频域,然后有选择地增强某些频段,再进行傅里叶反变换得到增强图像。小波域增强通过对图像做小波变换,然后增强小波分解的高频分量,再进行小波重构,即得到增强图像。流程如下:

原始图像→变换域→增强→反变换→增强图像

1.频域高通滤波锐化

频域锐化主要是通过频域高通滤波进行的。从频谱的角度看,图像模糊的实质是其高频部分被衰减,因此可以用高频滤波来使图像更加清晰。高频滤波器衰减傅里叶变换中的低频部分,而高频部分则无损地通过。与低通滤波器相对应,在高通滤波器中常用的传递函数Huv)有以下几种。

978-7-111-34687-6-Chapter03-102.jpg

图3.23 直方图规定化计算实例

(1)理想高通滤波器(Ideal High Pass Filter,IHPF)

一个二维理想高通滤波器的传递函数Huv)的表达式为

978-7-111-34687-6-Chapter03-103.jpg

式中,D0截止频率Duv)为点(uv)到频率平面原点的距离。理想高通滤波器传递函数的特性曲线如图3.24a所示。可见,这个滤波器与理想低通滤波器的特性正好完全相反,它把半径为D0的圆内所有频率完全衰减,对圆外的所有频率则无损通过。同理,理想高通滤波器在物理上也是不能实现的,只能用计算机进行模拟。

(2)巴特沃思高通滤波

具有截止频率D0n阶巴特沃思高通滤波器的传递函数为

978-7-111-34687-6-Chapter03-104.jpg

式中,Duv)仍由式(3.44)或式(3.45)定义。巴特沃思高通滤波器传递函数的特性曲线如图3.24b所示。

Duv=D0时,Huv)下降到最大值的1/2。如同巴特沃思低通滤波器的情形一样,如果截止频率的Huv)为最大值的978-7-111-34687-6-Chapter03-105.jpg,则式(3.73)可以化简为

978-7-111-34687-6-Chapter03-106.jpg

(3)指数高通滤波器(Exponent High Pass Filter,EHPF)具有截止频率D0的指数高通滤波器其传递函数表达式为

978-7-111-34687-6-Chapter03-107.jpg

式中,Duv)可由式(3.44)或式(3.45)表示;n控制着指数函数的增长率。指数高通滤波器传递函数的特性曲线如图3.24c所示。

Duv=D0时,Huv)下降到最大值的1/2。如同指数低通滤波器的情形一样,如果截止频率的Huv)为最大值的978-7-111-34687-6-Chapter03-108.jpg,则式(3.75)可以化简为

978-7-111-34687-6-Chapter03-109.jpg

(4)梯形滤波器(Trapezia High Pass Filter,THPF)梯形高通滤波器的传递函数表达式为

978-7-111-34687-6-Chapter03-110.jpg

式中,Duv)仍由式(3.44)或式(3.45)定义;D0为截止频率。图3.24d所示为典型的梯形高通滤波器传递函数的特性曲线。

978-7-111-34687-6-Chapter03-111.jpg

图3.24 各种高通滤波器的特性曲线

2.同态滤波

同态滤波一种增强图像对比度的频域处理方法,尤其适合于处理由于光照不均导致的图像不清。它与图像成像模型有密切关系。因此在介绍同态滤波之前,先简单介绍一下图像的成像原理。

人之所以能够看到某个物体,是因为有光照射在物体上,物体将入射光沿着一定方向进行反射或透射。反射/透射光进入人眼后在视网膜上产生感知信号。该感知信号沿着视神经传送到大脑后形成物体的像。按照图像的成像原理,可以将图像函数fxy)分为两部分:入射到可见场景上的光ixy)和场景中物体对入射光的反射/透视。前者也称为照度分量、照射分量,后者用反射率或透射率rxy)表示,也称为反射分量。因此可以将一幅图像进行如下简单建模:(www.chuimin.cn)

978-7-111-34687-6-Chapter03-112.jpg

式中,ixy)的值由光源决定;反射率rxy)由场景中的物体的属性决定,取值范围在全吸收(0)和全反射(1)之间。

在一幅图像中,反射分量决定图像的内容,随图像细节的不同在空间上作快速变化。而照度分量在空间则通常具有缓慢变化的性质。所以照度分量的频谱分布在低频区域,而反射分量的频谱则大部分分布在空间的高频区域。为了增强图像的细节表现,必须改善照度分量和反射分量的分布。

同态滤波以式(3.78)所给的图像模型为基础,在频域中同时进行图像对比度增强和压缩图像亮度范围。其基本思想是,将图像模型的两个分量分开,在频域中选择适当的滤波器参数Huv),在抑制图像低频信号ixy)的同时加强其高频信息rxy)。最后将频域信号反变换到空间域即可得到预期的效果。

为了把照度分量和反射分量分开处理,一种方法就是对式(3.78)进行对数变换,这样就可以把两者混合在一起的乘法运算转变成两者分离的加法运算。然后在频域中抑制低频照度分量并加强反射分量,最后再将处理结果反变换到空间域中,从而得到预期的效果。同态滤波的基本流程如图3.26所示。

978-7-111-34687-6-Chapter03-113.jpg

图3.25 图像模型示意图

978-7-111-34687-6-Chapter03-114.jpg

图3.26 同态滤波的流程

1)首先对式(3.78)的成像模型两边取对数,即

lnfxy=lnixy+lnrxy) (3.79)

2)对式(3.79)进行正向傅里叶变换,即

Fuv=Iuv+Ruv) (3.80)

3)选择一个适当的频域函数Huv)进行滤波,即

HuvFuv=HuvIuv+HuvRuv) (3.81)

4)根据需要对滤波后的数据进行处理。

5)将处理后的频域结果经过傅里叶反变换转换到对数空间域里,则有

hfxy=hixy+hrxy) (3.82)

978-7-111-34687-6-Chapter03-115.jpg

图3.27 同态滤波器圆对称滤波函数的剖面图,Duv)表示点(uv)离原点的距离

6)将式(3.82)进行指数变换,得到数字图像所在的空间,即

978-7-111-34687-6-Chapter03-116.jpg

可见,Huv)的选择决定了同态滤波的效果。为了能够在抑制低频的同时加强高频信息,Huv)通常在低频部分小于1,而在高频部分大于1。图3.27给出了一个满足上述要求的滤波函数剖面图,将它绕纵轴旋转360°即得到完整的二维Huv)。

【例3.10】 用同态滤波增强

在图3.28a所示的原始图像中,掩体内部的细节被遮蔽。图3.28b所示为用同态滤波器处理的结果,其中滤波函数γL=0.5,γH=2.0,减小了亮度的动态范围,增强了对比度可以使掩体内的细节显示出来,并且平衡墙外的灰度。增强后的图像更加锐化。

978-7-111-34687-6-Chapter03-117.jpg

图3.28 同态滤波增强