频域图像增强课件

第四章：频域图像增强频域滤波基础频域平滑滤波频域锐化滤波14.1频域滤波基础任何周期函数可以表示为不同频率的正弦与余弦的加权和形式。

复杂函数可以用由简单的正弦和余弦函数表示。

下面的函数是上面四个函数的和。24.1频域滤波基础傅立叶级数：3单变量连续函数f(x)的傅立叶变换F(u)定义为：

二维连续函数f(x，y)的傅立叶变换F(u，v)定义为：4.1频域滤波基础5一幅图像，在原点的傅立叶变换即等于图像的平均灰度级。二维离散傅立叶变换与反变换4.1频域滤波基础6用极坐标表示F(u)比较方便：R(u)和I(u)分别为F(u)的实部和虚部4.1频域滤波基础74.2傅立叶变换性质

二维傅立叶变换的基本性质：平移可以用于中心化变换,u和v的范围分别为[0,M-1]和[0,N-1],变换后的中心变为u=(M/2)+1,u=(N/2)+19一个简单的二维函数的中心谱

10

二维傅立叶变换的基本性质：分配性和比例变换性

傅立叶变换对加法具有分配性，对乘法没有：对于比例因子a和b，有：4.2傅立叶变换性质11

二维傅立叶变换的基本性质：周期性和对称性

周期性：共轭对称4.2傅立叶变换性质1314

二维傅立叶变换的基本性质：可分性其中对于每个x值，当v＝0,1,2,…,N-1时,该等式是完整的一维傅立叶变换。即F(x,v)是沿着f(x,y)的一行所进行的傅立叶变换。当y由0变为N-1时，沿着f(x,y)的所有行计算傅立叶变换。然而频率变量u仍然保持不变。为完成二维变换，将u值从0变到M-1.这涉及沿F(x,v)的每一列计算一维变换。可以通过先沿输入图像的每一行计算一维变换，然后沿中间结果的每一列再计算一维变换的方法来求二维变换。15图像频域滤波基本步骤174.2傅立叶变换性质18陷波滤波器滤波器可以置F(0,0)为零，因为滤波器除原点处有凹陷外，其它处均为常量函数

输出图像整体平均灰度降低，边缘突出

实际上显示的图像平均值并不是零，因为这样的话，一定存在灰度上的负值。图中是按照所有负值显示成黑色，而其它值按比例增加。19高通滤波的结果21空域滤波与频域滤波的对比22图像的平滑除了在空间域中进行外，也可以在频率域中进行。由于噪声主要集中在高频部分，为去除噪声改善图像质量，滤波器采用低通滤波器H(u,v)来抑制高频成分，通过低频成分，然后再进行逆傅立叶变换获得滤波图像，就可达到平滑图像的目的。4.3平滑的频域滤波器23理想低通滤波器25叠加的圆环具有5，15，30，80和230像素的半径，包含了92.0%，94.6%,96.4%,98%和99.5%的图像功率。理想低通滤波器26Butterworth低通滤波器

n阶Butterworth滤波器的传递函数为：它的特性是连续性衰减，而不象理想滤波器那样陡峭变化，即明显的不连续性。因此采用该滤波器滤波在抑制噪声的同时，图像边缘的模糊程度大大减小，没有振铃效应产生。29Butterworth低通滤波器3031巴特沃斯滤波器的空间解释32高斯低通滤波器Gauss滤波器的传递函数为：高斯低通滤波器的傅立叶反变换也是高斯的，这意味着反变换后高斯滤波器将没有振铃现象产生。

33高斯低通滤波器3435低通滤波实例36低通滤波实例374.4频域锐化滤波理想高通滤波器巴特沃斯高通滤波器高斯型高通滤波器频率域的拉普拉斯算子高频提升与高频加强滤波38理想高通滤波器巴特沃斯高通滤波器高斯高通滤波器39三种高通滤波器的空间解释40理想高通滤波结果41巴特沃斯高通滤波结果42高斯高通滤波结果434.4频率域锐化滤波器频率域的拉普拉斯算子:二维傅立叶变换后的拉普拉斯算子:一维傅立叶变换后的拉普拉斯算子:即:因此,频率域的拉普拉斯算子可以由如下滤波器实现:444.4频率域锐化滤波器钝化模板、高频提升滤波器和高频加强滤波钝化模板简单地由从一幅图像减去其自身模糊图像而生成的锐化图像构成。采用频率技术，这意味这从图像自身减去低通滤波后的图像而得到高通滤波的图像.高频提升过滤通过将f(x,y)乘上一个大于1的常熟A产生：可改写为：钝化模板在频率域中可由混合滤波器直接执行：类似的,高频提升过滤也可由混合滤波器直接执行：其中，A≥1，当A=1时，高频提升过滤为常规的高通过滤，A>1时，图象自身的贡献更加显著。45高频提升滤波(a)输入图像(b)(a)的拉普拉斯图像(c)A=2时高频提升滤波得到的图像(d)A=2.7时的图像与图3.43相同,但用频率域滤波464.4频率域锐化滤波器利用图像的高频成分强调增强的作用高频增强:在高通滤波器前简单地乘以一个常数,再增加一个偏移以便使零频率不被滤波器除掉.474.4频率域锐化滤波器(a)胸部X光图像(b)巴特沃思高通滤波的结果(c)高频增强滤波的结果(d)执行直方图均衡的结果48作业题1、近似一个离散导数的基本方法是对

f(x+1,y)-f(x,y)取差分。试找到空域一阶微分滤波器传递函数在频域中进行等价的操作H(u,v)；49答案答：已知空域滤波器在x方向的差分描述，则可写出y方向的差分描述f(x,y+1)-f(x,y) 则空域一阶微分滤波器可写为： g(x,y)=f(x+1,y)-f(x,y)+f(x,y+1)-f(x,y) 则对应频域表达式为： 则：50作业题2、给出一幅MⅹN大小的图像，要求做一个实验，它由截止频率为D0

的高斯型低通滤波器重复进行低通滤波。可以忽略计算的舍入误差