《第3章3.5-频域图象增强技术分解优秀PPT.ppt》由会员分享，可在线阅读，更多相关《第3章3.5-频域图象增强技术分解优秀PPT.ppt（73页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第三章图象增加技术频域图象增加 3.5 频域滤波增加频域滤波增加 3.5 3.5.1.1 概述概述概述概述频域滤波处理的一般方法：频域滤波处理的一般方法：频域滤波处理的一般方法：频域滤波处理的一般方法：G(u,v)=H(u,v)F(u,v)G(u,v)=H(u,v)F(u,v)f(x,y)f(x,y)F(u,v)F(u,v)G(u,v)G(u,v)g(x,y)g(x,y)变换变换变换变换乘乘乘乘H(u,v)H(u,v)反变换反变换反变换反变换因此：频域滤波处理的关键是选取合适的滤波器函数因此：频域滤波处理的关键是选取合适的滤波器函数因此：频域滤波处理的关键是选取合适的滤波器函数因此：频域滤波处
2、理的关键是选取合适的滤波器函数 H(u,v)H(u,v)!空间域与频率域空间域与频率域空间域空间域与频率域与频率域 3.5.2 3.5.2 低通滤波低通滤波低通滤波低通滤波一、志向低通滤波器一、志向低通滤波器一、志向低通滤波器一、志向低通滤波器H(u,v)=H(u,v)=式中式中式中式中D0D0为截止频率，是一个非负整数，为截止频率，是一个非负整数，为截止频率，是一个非负整数，为截止频率，是一个非负整数，D(u,v)D(u,v)是从点是从点是从点是从点(u,v)(u,v)到频率平面原点的距离。到频率平面原点的距离。到频率平面原点的距离。到频率平面原点的距离。1 1当当当当D(u,v)D(u,v
3、)D D0 00 0当当当当D(u,v)D(u,v)D D0 0志向低通滤波器的截面图志向低通滤波器的截面图志向低通滤波器的截面图志向低通滤波器的截面图0D0D(u,v)H(u,v)1H(u,v)H(u,v)作为距离函数作为距离函数D(u,v)D(u,v)的函数的截面图的函数的截面图D(u,v)=(u2+v2)D(u,v)=(u2+v2)1/21/2理想低通滤波器的三维透视图理想低通滤波器的三维透视图理想低通滤波器的三维透视图理想低通滤波器的三维透视图vuH(u,v)H(u,v)H(u,v)作为作为u u、v v的函数的三维透视图的函数的三维透视图特点：特点：特点：特点：物理上不行实现物理上不
4、行实现物理上不行实现物理上不行实现有抖动现象有抖动现象有抖动现象有抖动现象滤除高频成分使图象变模糊滤除高频成分使图象变模糊滤除高频成分使图象变模糊滤除高频成分使图象变模糊 3.5.2 3.5.2 低通滤波低通滤波低通滤波低通滤波志向低通滤波器的振铃（志向低通滤波器的振铃（志向低通滤波器的振铃（志向低通滤波器的振铃（ringring）现象：）现象：）现象：）现象：n n理想低通滤波器的截止频率的设计理想低通滤波器的截止频率的设计理想低通滤波器的截止频率的设计理想低通滤波器的截止频率的设计 先求出总的信号能量先求出总的信号能量先求出总的信号能量先求出总的信号能量P PT T：其中：其中：其中：其中
5、：P(u,v)=|F(u,v)|P(u,v)=|F(u,v)|2 2=R=R2 2(u,v)+I(u,v)+I2 2(u,v)(u,v)是能量模是能量模是能量模是能量模 假假假假如如如如将将将将变变变变换换换换作作作作中中中中心心心心平平平平移移移移，则则则则一一一一个个个个以以以以频频频频域域域域中中中中心心心心为为为为原原原原点点点点，r r为半径的圆就包含了百分之为半径的圆就包含了百分之为半径的圆就包含了百分之为半径的圆就包含了百分之 的能量的能量的能量的能量 求出相应的求出相应的求出相应的求出相应的D D0 0 D D0 0=r=(u=r=(u2 2+v+v2 2)1/21/2图图图图
6、(a)(a)为一幅包含不同细微为一幅包含不同细微为一幅包含不同细微为一幅包含不同细微环节的原始图象环节的原始图象环节的原始图象环节的原始图象图图图图(b)(b)为它的傅里叶频谱，为它的傅里叶频谱，为它的傅里叶频谱，为它的傅里叶频谱，其上迭加半径分别为其上迭加半径分别为其上迭加半径分别为其上迭加半径分别为5 5、1111、4545和和和和6868的圆周。这些圆周内的圆周。这些圆周内的圆周。这些圆周内的圆周。这些圆周内分别包含了原始图象中分别包含了原始图象中分别包含了原始图象中分别包含了原始图象中9090、9595、9999和和和和99.599.5的能量。的能量。的能量。的能量。图图图图(c)(c
7、)到图到图到图到图(f)(f)分别为用截止分别为用截止分别为用截止分别为用截止频率由各圆周半径确定的志频率由各圆周半径确定的志频率由各圆周半径确定的志频率由各圆周半径确定的志向低通滤波器进行处理得到向低通滤波器进行处理得到向低通滤波器进行处理得到向低通滤波器进行处理得到的结果。的结果。的结果。的结果。由图由图由图由图(c)(c)可见尽管只有可见尽管只有可见尽管只有可见尽管只有1010的的的的(高频高频高频高频)能量被滤除，但图能量被滤除，但图能量被滤除，但图能量被滤除，但图象中绝大多敌细微环节信息象中绝大多敌细微环节信息象中绝大多敌细微环节信息象中绝大多敌细微环节信息都丢失了，事实上这幅图已都
8、丢失了，事实上这幅图已都丢失了，事实上这幅图已都丢失了，事实上这幅图已无多少实际用途。无多少实际用途。无多少实际用途。无多少实际用途。由图由图由图由图(d)(d)可见当仅有可见当仅有可见当仅有可见当仅有5 5的的的的(高频高频高频高频)能量被滤除后，图象中仍有明显的振铃效应。能量被滤除后，图象中仍有明显的振铃效应。能量被滤除后，图象中仍有明显的振铃效应。能量被滤除后，图象中仍有明显的振铃效应。由图由图由图由图(e)(e)可见，假如只滤除可见，假如只滤除可见，假如只滤除可见，假如只滤除1 1的的的的(高频高频高频高频)能量，图象虽有确定程度的模糊但视觉效果尚可。能量，图象虽有确定程度的模糊但视觉
9、效果尚可。能量，图象虽有确定程度的模糊但视觉效果尚可。能量，图象虽有确定程度的模糊但视觉效果尚可。由图由图由图由图(f)(f)可见，滤除可见，滤除可见，滤除可见，滤除0505的的的的(高频高频高频高频)能量后所得到的滤波结果与原图象几乎无差别。能量后所得到的滤波结果与原图象几乎无差别。能量后所得到的滤波结果与原图象几乎无差别。能量后所得到的滤波结果与原图象几乎无差别。n n志向低通滤波器的分析志向低通滤波器的分析志向低通滤波器的分析志向低通滤波器的分析n n整个能量的整个能量的整个能量的整个能量的90%90%被一个半径为被一个半径为被一个半径为被一个半径为8 8的小圆周包含的小圆周包含的小圆周
10、包含的小圆周包含,大大大大部分尖锐的细微环节信息都存在于被去掉的部分尖锐的细微环节信息都存在于被去掉的部分尖锐的细微环节信息都存在于被去掉的部分尖锐的细微环节信息都存在于被去掉的10%10%的能量中的能量中的能量中的能量中n n小的边界和其它尖锐细微环节信息被包含在频谱小的边界和其它尖锐细微环节信息被包含在频谱小的边界和其它尖锐细微环节信息被包含在频谱小的边界和其它尖锐细微环节信息被包含在频谱的至多的至多的至多的至多0.5%0.5%的能量中的能量中的能量中的能量中n n被钝化的图像被一种特别严峻的振铃效果所影响被钝化的图像被一种特别严峻的振铃效果所影响被钝化的图像被一种特别严峻的振铃效果所影响
11、被钝化的图像被一种特别严峻的振铃效果所影响志向低通滤波器的处理效果：志向低通滤波器的处理效果：志向低通滤波器的处理效果：志向低通滤波器的处理效果：志向低通滤波器的处理效果：志向低通滤波器的处理效果：志向低通滤波器的处理效果：志向低通滤波器的处理效果：在空间域将模糊作为卷积过程来理解的关键是在空间域将模糊作为卷积过程来理解的关键是h(x,y)的特性。的特性。我们看到滤我们看到滤波器波器h(x,y)有两个主要特性：在原点处的一个主要成分，及中心成分四周有两个主要特性：在原点处的一个主要成分，及中心成分四周集中、呈周期性的成分。中心成分主要确定模糊，集中的成分主要确定集中、呈周期性的成分。中心成分主
12、要确定模糊，集中的成分主要确定了志向滤波器振铃现象的特性。中心成分的半径和距原点每单位距离上了志向滤波器振铃现象的特性。中心成分的半径和距原点每单位距离上周期的数量都与志向滤波器的截止频率值成反比。在顶部的插图是通过周期的数量都与志向滤波器的截止频率值成反比。在顶部的插图是通过空间滤波器中心的水平扫描线的灰度级剖面线。所示的轴表明白零幅度，空间滤波器中心的水平扫描线的灰度级剖面线。所示的轴表明白零幅度，所以，我们看到空间域滤波器有负值。通常这不是严峻的问题，因为较所以，我们看到空间域滤波器有负值。通常这不是严峻的问题，因为较大的中心成分确定着卷积的结果。然而，被滤波的图像能够有负值，因大的中心
13、成分确定着卷积的结果。然而，被滤波的图像能够有负值，因此通常要求做标定。此通常要求做标定。假定假定f(x,y)是一幅简洁图像，它由在黑色背景下五个光明的象素组成，如是一幅简洁图像，它由在黑色背景下五个光明的象素组成，如图图(c)所示。这些光明的点可被近似地看做冲激，其强度确定于点的亮度。所示。这些光明的点可被近似地看做冲激，其强度确定于点的亮度。然后，然后，h(x,y)和和f(x,y)的卷积仅仅是在每个冲激处的卷积仅仅是在每个冲激处“复制复制”h(x,y)的过程，的过程，此操作曲结果，如图此操作曲结果，如图4.13(d)所示，说明白原始点如何通过所示，说明白原始点如何通过f(x,y)和模糊滤和
14、模糊滤波器函数波器函数h(x,y)的卷积而变模糊。事实上，振铃在此种状况下很严峻，以的卷积而变模糊。事实上，振铃在此种状况下很严峻，以至于由相互之间的干扰而产生畸变。这些概念被扩展到更困难的图像上，至于由相互之间的干扰而产生畸变。这些概念被扩展到更困难的图像上，考虑将每个像素作为一个脉冲，而且其强度与象素的灰度级成比例。图考虑将每个像素作为一个脉冲，而且其强度与象素的灰度级成比例。图4.13底部的插图显示了通过被滤波图像中心的对角扫描线的灰度级剖面线。底部的插图显示了通过被滤波图像中心的对角扫描线的灰度级剖面线。(a)半径为半径为5的频率域低通滤波器的频率域低通滤波器(b)相应的空间滤波器相应
15、的空间滤波器(留意振铃留意振铃)(c)空间域的空间域的5个脉冲模拟个脉冲模拟5个象素值个象素值(d)空间域空间域(b)和和(c)的卷积的卷积二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器n n 为滤波器的阶次，为滤波器的阶次，为滤波器的阶次，为滤波器的阶次，D D0 0为滤波器的截止频率为滤波器的截止频率为滤波器的截止频率为滤波器的截止频率一个截止频率在与原点距离为一个截止频率在与原点距离为D0的的n阶阶Butterworth低通滤波器低通滤波器（BLPF）的变换函数如下：）的变换函数如下：尽管志向低通滤波器在数学上定义得很清晰，在计算机模拟中也可实现
16、，尽管志向低通滤波器在数学上定义得很清晰，在计算机模拟中也可实现，但在截止频率处直上直下的志向低通滤波器是不能用实际的电子器件实但在截止频率处直上直下的志向低通滤波器是不能用实际的电子器件实现的。现的。一阶巴特沃思低通滤波器转移函数三维图一阶巴特沃思低通滤波器转移函数三维图三阶巴特沃思低通滤波器转移函数三维图三阶巴特沃思低通滤波器转移函数三维图二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器n nButterworthButterworth低通滤波器的截面图低通滤波器的截面图低通滤波器的截面图低通滤波器的截面图02D(u,v)/D0H(u,v)1H(u,
17、v)H(u,v)作为作为D(u,v)/DD(u,v)/D0 0的函数的截面图的函数的截面图130.5低通巴特沃斯滤波器在凹凸频率间的过渡比较光滑，所以用巴特沃斯滤低通巴特沃斯滤波器在凹凸频率间的过渡比较光滑，所以用巴特沃斯滤波器得到的输出图象振铃效应不明显。波器得到的输出图象振铃效应不明显。二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器n nButterworthButterworth滤波器截止频率的设计滤波器截止频率的设计滤波器截止频率的设计滤波器截止频率的设计n n变换函数中不存在一个不连续点作为一个通过的和被滤波掉的变换函数中不存在一个不连续点作
18、为一个通过的和被滤波掉的变换函数中不存在一个不连续点作为一个通过的和被滤波掉的变换函数中不存在一个不连续点作为一个通过的和被滤波掉的截止频率的明显划分截止频率的明显划分截止频率的明显划分截止频率的明显划分n n通常把通常把通常把通常把H(u,v)H(u,v)起先小于其最大值的确定比例的点当作其截止频率起先小于其最大值的确定比例的点当作其截止频率起先小于其最大值的确定比例的点当作其截止频率起先小于其最大值的确定比例的点当作其截止频率点点点点n n有两种选择：有两种选择：有两种选择：有两种选择：n n选择选择选择选择1 1：H(u,v)=0.5 H(u,v)=0.5 当当当当 D0=D(u,v)D
19、0=D(u,v)时时时时n n选择选择选择选择2 2：H(u,v)=H(u,v)=当当当当 D0=D(u,v)D0=D(u,v)时时时时二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器二、巴特沃斯低通滤波器ButterworthButterworth低通滤波器的分析低通滤波器的分析低通滤波器的分析低通滤波器的分析在任何经在任何经在任何经在任何经BLPFBLPF处理过的图像中都没有明显的振铃效果，处理过的图像中都没有明显的振铃效果，处理过的图像中都没有明显的振铃效果，处理过的图像中都没有明显的振铃效果，这是滤波器在低频和高频之间的平滑过渡的结果这是滤波器在低频和高频之间的平滑过
20、渡的结果这是滤波器在低频和高频之间的平滑过渡的结果这是滤波器在低频和高频之间的平滑过渡的结果低通滤波是一个以牺牲图像清晰度为代价来削减干扰低通滤波是一个以牺牲图像清晰度为代价来削减干扰低通滤波是一个以牺牲图像清晰度为代价来削减干扰低通滤波是一个以牺牲图像清晰度为代价来削减干扰效果的修饰过程效果的修饰过程效果的修饰过程效果的修饰过程当图象由于量化不足产生虚假轮廓时常可用低通滤波进行平滑以改进图象质量。当图象由于量化不足产生虚假轮廓时常可用低通滤波进行平滑以改进图象质量。图图(a)为一幅由为一幅由256级灰度量化为级灰度量化为12个灰度级的图象，帽子和肩膀等处均有不同程度的个灰度级的图象，帽子和肩
21、膀等处均有不同程度的虚假轮廓存在。虚假轮廓存在。图图(b)和图和图(c)分别为用志向低通滤波器和用阶数为分别为用志向低通滤波器和用阶数为1的巴特沃斯低通滤波器进行平滑处的巴特沃斯低通滤波器进行平滑处理所得到的结果。所用两个滤波器的截止频率所对应的半径均为理所得到的结果。所用两个滤波器的截止频率所对应的半径均为30。比较两幅滤波结果图象，志向低通滤波器的结果图象中有较明显的振铃现象，而巴特比较两幅滤波结果图象，志向低通滤波器的结果图象中有较明显的振铃现象，而巴特沃斯滤波器的效果较好。沃斯滤波器的效果较好。例：频域低通滤波消退虚假轮廓例：频域低通滤波消退虚假轮廓巴特沃斯低通滤波器的处理效果巴特沃斯
22、低通滤波器的处理效果巴特沃斯低通滤波器的处理效果巴特沃斯低通滤波器的处理效果(a)(d)阶数为阶数为1，2，5和和20的巴特沃斯低通滤波器的空间表示及相的巴特沃斯低通滤波器的空间表示及相应的通过滤波器中心的灰度级剖画图应的通过滤波器中心的灰度级剖画图(全部的滤波器都有半径为全部的滤波器都有半径为5的截的截止频率止频率)，振铃作为滤波器阶数的函数越独创显。，振铃作为滤波器阶数的函数越独创显。一阶巴特沃斯滤波器没有振铃，在二阶中振铃通常很微小，一阶巴特沃斯滤波器没有振铃，在二阶中振铃通常很微小，但阶数增高时振铃便成为一个重要因素。但阶数增高时振铃便成为一个重要因素。三、指数形低通滤波器三、指数形低
23、通滤波器D(u,v)=D0，H(u,v)降为最大值的降为最大值的 。n为阶数。为阶数。一阶指数形低通滤波器转移函数三维图一阶指数形低通滤波器转移函数三维图一阶指数形低通滤波器转移函数剖面图一阶指数形低通滤波器转移函数剖面图三阶指数形低通滤波器转移函数三维图三阶指数形低通滤波器转移函数三维图三阶指数形低通滤波器转移函数剖面图三阶指数形低通滤波器转移函数剖面图四、梯形低通滤波器四、梯形低通滤波器梯形低通滤波器转移函数三维图梯形低通滤波器转移函数三维图梯形低通滤波器转移函数剖面图梯形低通滤波器转移函数剖面图 3.5.3 3.5.3 高通滤波高通滤波高通滤波高通滤波频域高通滤波器的基本思想频域高通滤波
24、器的基本思想频域高通滤波器的基本思想频域高通滤波器的基本思想G(u,v)=F(u,v)H(u,v)G(u,v)=F(u,v)H(u,v)F(u,v)F(u,v)是须要锐化图像的傅立叶变换形式。是须要锐化图像的傅立叶变换形式。是须要锐化图像的傅立叶变换形式。是须要锐化图像的傅立叶变换形式。目标是选取一个滤波器变换函数目标是选取一个滤波器变换函数目标是选取一个滤波器变换函数目标是选取一个滤波器变换函数H(u,v)H(u,v)，通过它削，通过它削，通过它削，通过它削减减减减F(u,v)F(u,v)的低频部分来得到的低频部分来得到的低频部分来得到的低频部分来得到G(u,v)G(u,v)。运用傅立叶逆变
25、换得到锐化后的图像。运用傅立叶逆变换得到锐化后的图像。运用傅立叶逆变换得到锐化后的图像。运用傅立叶逆变换得到锐化后的图像。一、志向高通滤波器一、志向高通滤波器一、志向高通滤波器一、志向高通滤波器定义定义定义定义:一个二维的志向高通滤波器（一个二维的志向高通滤波器（一个二维的志向高通滤波器（一个二维的志向高通滤波器（ILPFILPF）的转换函数满足）的转换函数满足）的转换函数满足）的转换函数满足（是一个分段函数）（是一个分段函数）（是一个分段函数）（是一个分段函数）其中：其中：其中：其中：D0 D0 为截止频率，为截止频率，为截止频率，为截止频率，D(u,v)D(u,v)为距离函数为距离函数为距
26、离函数为距离函数 D(u,v)=(u2+v2)1/2 D(u,v)=(u2+v2)1/2H(u,v)=H(u,v)=0 0当当D(u,v)D(u,v)D D0 01 1当当D(u,v)DD(u,v)D0 0 3.5.3 3.5.3 高通滤波高通滤波高通滤波高通滤波n n志向低通滤波器的截面图志向低通滤波器的截面图志向低通滤波器的截面图志向低通滤波器的截面图0D0D(u,v)H(u,v)1H(u,v)H(u,v)作为距离函数作为距离函数D(u,v)D(u,v)的函数的截面图的函数的截面图 3.5.3 3.5.3 高通滤波高通滤波高通滤波高通滤波n n志向高通滤波器的三维透视图志向高通滤波器的三维
27、透视图志向高通滤波器的三维透视图志向高通滤波器的三维透视图vuH(u,v)H(u,v)H(u,v)作为作为u u、v v的函数的三维透视图的函数的三维透视图 3.5.3 3.5.3 高通滤波高通滤波高通滤波高通滤波二、巴特沃斯高通滤波器二、巴特沃斯高通滤波器二、巴特沃斯高通滤波器二、巴特沃斯高通滤波器n n 为滤波器的阶次，为滤波器的阶次，为滤波器的阶次，为滤波器的阶次，D D0 0为滤波器的截止频率为滤波器的截止频率为滤波器的截止频率为滤波器的截止频率定义：一个截止频率在与原点距离为定义：一个截止频率在与原点距离为定义：一个截止频率在与原点距离为定义：一个截止频率在与原点距离为D0D0的的的
28、的n n阶阶阶阶ButterworthButterworth高通滤波器高通滤波器高通滤波器高通滤波器（BHPFBHPF）的变换函数如下：）的变换函数如下：）的变换函数如下：）的变换函数如下：三阶巴特沃思高通滤波器转移函数三维图三阶巴特沃思高通滤波器转移函数三维图 3.5.3 3.5.3 高通滤波高通滤波高通滤波高通滤波ButterworthButterworth高通滤波器的截面图高通滤波器的截面图高通滤波器的截面图高通滤波器的截面图02D(u,v)/D0H(u,v)1H(u,v)H(u,v)作为作为D(u,v)/DD(u,v)/D0 0的函数的截面图的函数的截面图130.5高通巴特沃斯滤波器在
29、通过和滤掉的频率之间也没有不连续的分界，高通巴特沃斯滤波器在通过和滤掉的频率之间也没有不连续的分界，由于在凹凸频率间的过渡比较光滑，所以用巴持沃斯滤波器得到的由于在凹凸频率间的过渡比较光滑，所以用巴持沃斯滤波器得到的输出图其振铃效应不明显。输出图其振铃效应不明显。3.5.3 3.5.3 高通滤波高通滤波高通滤波高通滤波 Butterworth Butterworth高通滤波器截止频率设计高通滤波器截止频率设计高通滤波器截止频率设计高通滤波器截止频率设计变换函数中不存在一个不连续点作为一个通过的和被滤掉的截止频变换函数中不存在一个不连续点作为一个通过的和被滤掉的截止频变换函数中不存在一个不连续点
30、作为一个通过的和被滤掉的截止频变换函数中不存在一个不连续点作为一个通过的和被滤掉的截止频率的明显划分率的明显划分率的明显划分率的明显划分通常把通常把通常把通常把H(u,v)H(u,v)起先小于其最大值（起先小于其最大值（起先小于其最大值（起先小于其最大值（1 1）的确定比例的点当作其截止频）的确定比例的点当作其截止频）的确定比例的点当作其截止频）的确定比例的点当作其截止频率点率点率点率点有两种选择：有两种选择：有两种选择：有两种选择：选择选择选择选择1 1：H(u,v)=0.5 H(u,v)=0.5 当当当当 D0=D(u,v)D0=D(u,v)时时时时选择选择选择选择2 2：H(u,v)=H
31、(u,v)=当当当当 D0=D(u,v)D0=D(u,v)时时时时D0 /D(u,v)3.5.3 3.5.3 高通滤波高通滤波高通滤波高通滤波n nButterworthButterworth高通滤波器的分析高通滤波器的分析高通滤波器的分析高通滤波器的分析n n问题：低频成分被严峻地消弱了，使图像失去层问题：低频成分被严峻地消弱了，使图像失去层问题：低频成分被严峻地消弱了，使图像失去层问题：低频成分被严峻地消弱了，使图像失去层次次次次n n改进措施：改进措施：改进措施：改进措施：n n加一个常数到变换函数加一个常数到变换函数加一个常数到变换函数加一个常数到变换函数 H(u,v)+A H(u,v
32、)+A n n这种方法被称为高频强调这种方法被称为高频强调这种方法被称为高频强调这种方法被称为高频强调n n为了解决变暗的趋势，在变换结果图像上再进行为了解决变暗的趋势，在变换结果图像上再进行为了解决变暗的趋势，在变换结果图像上再进行为了解决变暗的趋势，在变换结果图像上再进行一次直方图均衡化。这种方法被称为后滤波处理一次直方图均衡化。这种方法被称为后滤波处理一次直方图均衡化。这种方法被称为后滤波处理一次直方图均衡化。这种方法被称为后滤波处理巴特沃斯高通滤波器的处理效果巴特沃斯高通滤波器的处理效果巴特沃斯高通滤波器的处理效果巴特沃斯高通滤波器的处理效果图图(a)为一幅比较模糊的图象为一幅比较模糊
33、的图象图图(b)给出用阶数为给出用阶数为1的巴特沃斯高通滤波器进行处理所得到的结果。因为高通处理的巴特沃斯高通滤波器进行处理所得到的结果。因为高通处理后低频重量大部分被滤除，所以虽然图中各区域的边界得到了较明显的增加，但图后低频重量大部分被滤除，所以虽然图中各区域的边界得到了较明显的增加，但图中原来比较平滑区域内部的灰度动态范围被压缩，整幅图比较昏暗。中原来比较平滑区域内部的灰度动态范围被压缩，整幅图比较昏暗。图图(c)给出变换函数给出变换函数 H(u,v)+A(A=0.5)的结果，不仅边缘得到了增加，整个图象层次的结果，不仅边缘得到了增加，整个图象层次也比较丰富。也比较丰富。三、指数形高通滤
34、波器三、指数形高通滤波器三阶指数形高通滤波器转移函数三维图三阶指数形高通滤波器转移函数三维图三阶指数形高通滤波器转移函数剖面图三阶指数形高通滤波器转移函数剖面图四、梯形高通滤波器四、梯形高通滤波器梯形高通滤波器转移函数三维图梯形高通滤波器转移函数三维图梯形高通滤波器转移函数剖面图梯形高通滤波器转移函数剖面图D1D0n n 为滤波器的阶次，为滤波器的阶次，为滤波器的阶次，为滤波器的阶次，D D0 0为滤波器的放射中心，为滤波器的放射中心，为滤波器的放射中心，为滤波器的放射中心，WW为阻带宽度为阻带宽度为阻带宽度为阻带宽度巴特沃斯带阻滤波器巴特沃斯带阻滤波器巴特沃斯带阻滤波器巴特沃斯带阻滤波器 3
35、.5.4 3.5.4 巴特沃斯巴特沃斯巴特沃斯巴特沃斯带通与带阻滤波带通与带阻滤波带通与带阻滤波带通与带阻滤波巴特沃斯带通滤波器巴特沃斯带通滤波器巴特沃斯带通滤波器巴特沃斯带通滤波器带通滤波器允许确定频率范围内的信号通过而阻挡其它频军范围内带通滤波器允许确定频率范围内的信号通过而阻挡其它频军范围内的信号通过，与此相对应，带阻滤波器阻挡确定频率范围内的信号的信号通过，与此相对应，带阻滤波器阻挡确定频率范围内的信号通过而允许其它频率范围内的信号通过。通过而允许其它频率范围内的信号通过。3.5.5 3.5.5 同态滤波同态滤波同态滤波同态滤波同态滤波器的基本思想同态滤波器的基本思想同态滤波器的基本思
36、想同态滤波器的基本思想一个图像一个图像一个图像一个图像f(x,y)f(x,y)可以依据它的照度和反射重量的乘积来表示可以依据它的照度和反射重量的乘积来表示可以依据它的照度和反射重量的乘积来表示可以依据它的照度和反射重量的乘积来表示 f(x,y)=i(x,y)r(x,y)f(x,y)=i(x,y)r(x,y)其中：其中：其中：其中：i(x,y)i(x,y)为照度函数，为照度函数，为照度函数，为照度函数，r(x,y)r(x,y)反射重量函数反射重量函数反射重量函数反射重量函数通过同时实现压缩亮度范围和增加对比度，来改进图像的表现通过同时实现压缩亮度范围和增加对比度，来改进图像的表现通过同时实现压缩
37、亮度范围和增加对比度，来改进图像的表现通过同时实现压缩亮度范围和增加对比度，来改进图像的表现定义：因为两个函数乘积的傅立叶变换不是可分别的，也即：定义：因为两个函数乘积的傅立叶变换不是可分别的，也即：定义：因为两个函数乘积的傅立叶变换不是可分别的，也即：定义：因为两个函数乘积的傅立叶变换不是可分别的，也即：Ff(x,y)Fi(x,y)Fr(x,y)Ff(x,y)Fi(x,y)Fr(x,y)假设我们定义假设我们定义假设我们定义假设我们定义z(x,y)=ln f(x,y)z(x,y)=ln f(x,y)=ln(i(x,y)r(x,y)=ln(i(x,y)r(x,y)=ln i(x,y)+ln r(
38、x,y)=ln i(x,y)+ln r(x,y)同态滤波是一种在频域中同时将图象亮度范围进行压缩和将图象对同态滤波是一种在频域中同时将图象亮度范围进行压缩和将图象对比度进行增加的方法。比度进行增加的方法。3.5.5 3.5.5 同态滤波同态滤波同态滤波同态滤波那么有：那么有：Fz(x,y)=Fln f(x,y)=Fln i(x,y)+Fln r(x,y)或或Z(u,v)=I(u,v)+R(u,v)其中其中I(u,v)和和R(u,v)分别是分别是ln i(x,y)和和ln r(x,y)的傅立叶变换的傅立叶变换用滤波器函数用滤波器函数H(u,v)的方法处理的方法处理Z(u,v)，有：，有：S(u,
39、v)=H(u,v)Z(u,v)=H(u,v)I(u,v)+H(u,v)R(u,v)其中其中S(u,v)是结果图像的傅立叶变换是结果图像的傅立叶变换在空域中：在空域中：s(x,y)=F-1S(u,v)=F-1H(u,v)I(u,v)+F-1H(u,v)R(u,v)通过设：通过设：i(x,y)=F-1H(u,v)I(u,v)r(x,y)=F-1H(u,v)R(u,v)可以表示为：可以表示为：s(x,y)=i(x,y)+r(x,y)最终，通过最终，通过i(x,y)和和 r(x,y)的逆操作（指数操作）产生增加后的的逆操作（指数操作）产生增加后的图像图像g(x,y)3.5.5 3.5.5 同态滤波同态
40、滤波同态滤波同态滤波也即：也即：g(x,y)=exps(x,y)=expi(x,y)expr(x,y)=i0(x,y)r0(x,y)其中其中i0(x,y)=expi(x,y)和和r0(x,y)=expr(x,y)是输出图像的明度和反射重量是输出图像的明度和反射重量g0(x,y)=i0(x,y)r0(x,y)利用前述概念进行增加的方法可以归纳为：利用前述概念进行增加的方法可以归纳为：在此特定应用中，问题的关键在于将照度和反射重量用进行分在此特定应用中，问题的关键在于将照度和反射重量用进行分别。同态滤波器函数别。同态滤波器函数H(u,v)能够分别对这两部分进行操作。能够分别对这两部分进行操作。ln
41、FFTH(u,v)(FFT)-1expf(x,y)g(x,y)3.5.5 3.5.5 同态滤波同态滤波同态滤波同态滤波同态滤波器的效果分析同态滤波器的效果分析同态滤波器的效果分析同态滤波器的效果分析图像的照度重量的特点是平缓的空域变更，而反射重量则图像的照度重量的特点是平缓的空域变更，而反射重量则图像的照度重量的特点是平缓的空域变更，而反射重量则图像的照度重量的特点是平缓的空域变更，而反射重量则近于陡峭的空域变更近于陡峭的空域变更近于陡峭的空域变更近于陡峭的空域变更这些特性使得将图像的对数的傅立叶变换的低频部分对应这些特性使得将图像的对数的傅立叶变换的低频部分对应这些特性使得将图像的对数的傅立
42、叶变换的低频部分对应这些特性使得将图像的对数的傅立叶变换的低频部分对应于照度重量，而高频部分对应于反射重量于照度重量，而高频部分对应于反射重量于照度重量，而高频部分对应于反射重量于照度重量，而高频部分对应于反射重量尽管这种对应关系只是一个粗略的近似，但它们可以用于尽管这种对应关系只是一个粗略的近似，但它们可以用于尽管这种对应关系只是一个粗略的近似，但它们可以用于尽管这种对应关系只是一个粗略的近似，但它们可以用于优化图像的增加操作优化图像的增加操作优化图像的增加操作优化图像的增加操作一个好的限制可以通过用同态滤波器比照度和反射重量分一个好的限制可以通过用同态滤波器比照度和反射重量分一个好的限制可
43、以通过用同态滤波器比照度和反射重量分一个好的限制可以通过用同态滤波器比照度和反射重量分别操作得到别操作得到别操作得到别操作得到这个限制要求指定一个滤波器函数这个限制要求指定一个滤波器函数这个限制要求指定一个滤波器函数这个限制要求指定一个滤波器函数H(u,v)H(u,v)，它对于傅立叶，它对于傅立叶，它对于傅立叶，它对于傅立叶变换的低频和高频部分的影响是不同的变换的低频和高频部分的影响是不同的变换的低频和高频部分的影响是不同的变换的低频和高频部分的影响是不同的 3.5.5 3.5.5 同态滤波同态滤波同态滤波同态滤波同态滤波器的截面图同态滤波器的截面图同态滤波器的截面图同态滤波器的截面图0D(u
44、,v)H(u,v)1H(u,v)H(u,v)作为作为D(u,v)D(u,v)的函数的截面图的函数的截面图H HL 3.5.5 3.5.5 同态滤波同态滤波同态滤波同态滤波同态滤波器的效果分析同态滤波器的效果分析假如参数假如参数L和和H的选取使得的选取使得L 1 则滤波器函数将削减低频部分、扩大高频部分，则滤波器函数将削减低频部分、扩大高频部分，最终的结果将是既压缩了有效范围，又扩大了对最终的结果将是既压缩了有效范围，又扩大了对比度。比度。3.5.5 3.5.5 同态滤波同态滤波同态滤波同态滤波同态滤波增加效果同态滤波增加效果左图为一幅人脸图象，由于单侧光照明的缘由使得人脸在图象的右侧产左图为一
45、幅人脸图象，由于单侧光照明的缘由使得人脸在图象的右侧产生阴影，头发的发际线很不清晰。右图为用生阴影，头发的发际线很不清晰。右图为用HL0.5，HH2.0进行同进行同态滤波得到的增加结果。图象增加后，人脸与头独创显分开，衣领也看态滤波得到的增加结果。图象增加后，人脸与头独创显分开，衣领也看出来了。同时使动态范围压缩出来了。同时使动态范围压缩(如眼睛处如眼睛处)并使对比度增加并使对比度增加(加入脸与头加入脸与头发交界处发交界处)。原始图象像中墙体内部的细微环节被墙外的光所遮挡，图原始图象像中墙体内部的细微环节被墙外的光所遮挡，图(b)显示显示了用同态滤波器处理这一图像的结果，其中图滤波函数中的了用
46、同态滤波器处理这一图像的结果，其中图滤波函数中的L0.5,H=2。减小亮度的动态范围，增加对比度可以使墙体内的细。减小亮度的动态范围，增加对比度可以使墙体内的细微环节显示出来，并且平衡墙外的灰度，增加后的图像同样更加锐微环节显示出来，并且平衡墙外的灰度，增加后的图像同样更加锐化。化。TermsHighpass filter:Highpass filter:高通滤波器高通滤波器高通滤波器高通滤波器Bandpass filter:Bandpass filter:带通滤波器带通滤波器带通滤波器带通滤波器Bandreject filterBandreject filter、Bandstop filte
47、r:Bandstop filter:带阻滤波器带阻滤波器带阻滤波器带阻滤波器Ideal filter:Ideal filter:志向滤波器志向滤波器志向滤波器志向滤波器Butterworth filter:Butterworth filter:巴特沃思滤波器巴特沃思滤波器巴特沃思滤波器巴特沃思滤波器Exponential filter:Exponential filter:指数滤波器指数滤波器指数滤波器指数滤波器Trapezoidal filter:Trapezoidal filter:梯形滤波器梯形滤波器梯形滤波器梯形滤波器Transfer function:Transfer functio
48、n:传递函数传递函数传递函数传递函数TermsFrequency response:Frequency response:频率响应频率响应频率响应频率响应Cut-off frequency:Cut-off frequency:截止频率截止频率截止频率截止频率Spectrum:Spectrum:频谱频谱频谱频谱Amplitude spectrum:Amplitude spectrum:幅值谱幅值谱幅值谱幅值谱Phase spectrum:Phase spectrum:相位谱相位谱相位谱相位谱Power spectrum:Power spectrum:功率谱功率谱功率谱功率谱Blur:Blur:模糊模糊模糊模糊High frequency enhancement filter:High frequency enhancement filter:高频加强滤波器高频加强滤波器高频加强滤波器高频加强滤波器作业：编程实现以下算法作业：编程实现以下算法1.1.志向低通滤波志向低通滤波志向低通滤波志向低通滤波2.2.巴特沃斯低通滤波巴特沃斯低通滤波巴特沃斯低通滤波巴特沃斯低通滤波3.3.志向高通滤波志向高通滤波志向高通滤波志向高通滤波4.4.巴特沃斯高通滤波巴特沃斯高通滤波巴特沃斯高通滤波巴特沃斯高通滤波5.5.同态滤波同态滤波同态滤波同态滤波