冈萨雷斯数字图像处理精校版

1、会计学1冈萨雷斯数字图像处理精校版冈萨雷斯数字图像处理精校版第1页/共59页变化最慢的频率成分变化最慢的频率成分(u=v=0)(u=v=0)对应一幅图像的平均灰度对应一幅图像的平均灰度级级( , )( , )( , )( .)M 1 N 1x 0 y 0F 0 0f x yMN f x y4 621 当从变换的原点移开时，对低频对应着图像的慢变化分当从变换的原点移开时，对低频对应着图像的慢变化分量，如图像的平滑部分量，如图像的平滑部分进一步离开原点时，较高的频率对应图像中变化越来进一步离开原点时，较高的频率对应图像中变化越来越快的灰度级，如边缘或噪声等尖锐部分越快的灰度级，如边缘或噪声等尖锐部

2、分4.7.1、频率域的其他特性：、频率域的其他特性：()( , )( , )ux vyM 1 N 1i2MNx 0 y 0F u vf x y e 第2页/共59页从幅度谱中我们可以看出明亮线和原始图像中对应的轮从幅度谱中我们可以看出明亮线和原始图像中对应的轮廓线是垂直的。如果原始图像中有圆形区域那么幅度谱廓线是垂直的。如果原始图像中有圆形区域那么幅度谱中也呈圆形分布。中也呈圆形分布。( , )( , )( , )i u vF u vF u v e 第3页/共59页傅里叶频谱显示了傅里叶频谱显示了45450 0的强边缘，在垂直轴偏左的的强边缘，在垂直轴偏左的部分有垂直成分（对应两个氧化物突起）

3、。部分有垂直成分（对应两个氧化物突起）。第4页/共59页频率域滤波的基本步骤频率域滤波的基本步骤思想：通过滤波器函数以某种方式来修改图像变换，思想：通过滤波器函数以某种方式来修改图像变换，然后通过取结果的反变换来获得处理后的输出图像然后通过取结果的反变换来获得处理后的输出图像前处前处理理f (x , y)g (x , y)F (u , v)H (u , v) F (u , v)后处理后处理DFT 滤波器滤波器H (u , v)IDFT第5页/共59页 1( , )( , ) ( , )g x yH u v F u v 4.7.2、频率域滤波基础：、频率域滤波基础：前处前处理理f (x , y)

4、g (x , y)F (u , v)H (u , v) F (u , v)后处理后处理DFT 滤波器滤波器H (u , v)IDFT第6页/共59页陷波滤波器（带阻滤波）陷波滤波器（带阻滤波）0( , )(,)( , )221M Nu vH u v 其其它它设置设置F(0,0)=0(F(0,0)=0(结果图像的平均值为零结果图像的平均值为零) )，而保留，而保留其它傅里叶变换的频率成分不变其它傅里叶变换的频率成分不变由于图像平均值为由于图像平均值为0 0而产生整体平均灰度级的降而产生整体平均灰度级的降低，低，因此几乎没有平滑的灰度级细节因此几乎没有平滑的灰度级细节第7页/共59页低通滤波器：低

5、通滤波器：被低通滤波的图像比原始图像少了尖锐的细节部被低通滤波的图像比原始图像少了尖锐的细节部分而突出了平滑过渡部分分而突出了平滑过渡部分使低频通过，高频衰减的滤波器使低频通过，高频衰减的滤波器第8页/共59页高通滤波器：高通滤波器：被高通滤波的图像比原始图像少了灰度级的平被高通滤波的图像比原始图像少了灰度级的平滑过渡而突出了边缘等细节部分滑过渡而突出了边缘等细节部分使高频通过，低频衰减的滤波器使高频通过，低频衰减的滤波器第9页/共59页高通滤波结果高通滤波结果高通滤波改进结果高通滤波改进结果陷波滤波器将原点陷波滤波器将原点设置为设置为0，平均灰，平均灰度为度为0，负灰度置，负灰度置为为0。陷

6、波滤波结果陷波滤波结果 在高通滤波器中加在高通滤波器中加入常量，以使入常量，以使F(0,0)F(0,0)不被完全消不被完全消除。（防止直流项除。（防止直流项消除，保持色调）消除，保持色调）该高通滤波器原点该高通滤波器原点为为0，因此几乎没有，因此几乎没有平滑的灰度级细节平滑的灰度级细节，且图像较暗。且图像较暗。第10页/共59页错误的填充图像会导致错误的结错误的填充图像会导致错误的结果果第11页/共59页4.7.3 频率域的滤波步骤：频率域的滤波步骤：1 1、对要滤波的图像、对要滤波的图像 进行填充得到进行填充得到 ，典型地：典型地：P=2M,Q=2NP=2M,Q=2N( , )M Nfx y

7、 ( , )P Qfx y 2 2、填充图像，用、填充图像，用 乘以输入图像进行中心变换乘以输入图像进行中心变换(1)xy ( , ) ( 1)(,)22x yP QPQfx yF uv 3 3、变换到频域、变换到频域( , )( , )( 1)x yP QF u vfx y 4 4、生成一个实的、中心对称的滤波器、生成一个实的、中心对称的滤波器 ，中心在，中心在 ( , )P QHu v (,)22P Q( , )( , )( , )G u vF u v H u v频域滤波频域滤波: : 1( , )( , )( , )P Qgx yF u v H u v 5 5、变换到空间域变换到空间域:

8、 : ( , )P Qreal gx y 6 6、取实部取实部: : ( , )( , )( 1)x ypP Qgx yreal gx y 7 7、取消输入图像的乘数取消输入图像的乘数: :8 8、提、提取取 区域区域: :MN ( , )( , )M NP Qgx ygx y 的的对对应应部部分分第12页/共59页( , )M Nfx y ( , )P Qfx y ( 1)( , )xyP Qfx y ( , )P QFu v ( , ):LPHu vP Q ( , )( , )LPP QHu vFu v ( 1)( , )xyP Qreal gx y ( , )M Ngx y 图图4.36

9、第13页/共59页1100( , )( , )(, ) (,)(4.623)MNmnf x yh x yf m n h xm yn对比空间域滤波：在对比空间域滤波：在M MN N的图像的图像f f上，用上，用m mn n的滤波器进的滤波器进行线性滤波行线性滤波( , )( , ) (,)( .)absa tbg x yw s t f xs yt 3 41(4.6-23)和和(3.4-1)本质上是相似的；相差之处只在于：常数、负号及求和的本质上是相似的；相差之处只在于：常数、负号及求和的上、下限；上、下限；在实践中，我们宁愿使用在实践中，我们宁愿使用(3.4-1)和较小的滤波器模板来实现滤波处理

10、；和较小的滤波器模板来实现滤波处理；滤波在频率域中更为直观，可以在频率域指定滤波器，做反变换，然后在滤波在频率域中更为直观，可以在频率域指定滤波器，做反变换，然后在空间域使用结果滤波器作为在空间域构建小滤波器模板的指导；空间域使用结果滤波器作为在空间域构建小滤波器模板的指导；大小为大小为M MN N的两个函数的两个函数f(x,y)f(x,y)和和h(x,y)h(x,y)的频率域滤波表示为：的频率域滤波表示为： 由卷积定理，该运算对应的空间域运算为：由卷积定理，该运算对应的空间域运算为：( , )( , )H u vF u v 第14页/共59页对应空间域高斯低通滤波器为对应空间域高斯低通滤波器

11、为A B , 1 2对应空间域高斯高通滤波器为对应空间域高斯高通滤波器为22/2( )(4.75)uH uAe 2 2 22( )2(4.7 6)xh xAe 222212/2/2( )(4.77)uuH uAeBe 222222122212( )22(4.78)xxh xAeBe 频率域高斯低通滤波器函数频率域高斯低通滤波器函数频率域高斯高通滤波器函数频率域高斯高通滤波器函数第15页/共59页图图4.374.37频域高斯低通滤波器频域高斯低通滤波器频域高斯高通滤波器频域高斯高通滤波器空域高斯低通滤波器及模板空域高斯低通滤波器及模板空域高斯高通滤波器及模板空域高斯高通滤波器及模板 第16页/共

12、59页图图4.38( , ):600 600f x y ( , ) :600 600F u v 图图4.393 3( , )hx y 3 3( , )Hu v ( , ):602 602H u v 1( , )( , ) ( , )g x yH u v F u v 3 3602 602( , )( , )hx yfx y 空域线性滤波空域线性滤波的结果的结果例例4.15第17页/共59页4.8.1理想低通滤波器理想低通滤波器频率域的中心在频率域的中心在 ，从点，从点(u,v)(u,v)到中心（原点）的到中心（原点）的距离如下距离如下1222( , )()()(4.82)22PQD u vuv

13、截断傅里叶变换中的所有处于指定距离截断傅里叶变换中的所有处于指定距离D D0 0之外的高频成分之外的高频成分0001( , )( , ),0(4.81)0( , )ILPFD u vDHu vDD u vD (,)2 2P Q第18页/共59页理想低通滤波器理想低通滤波器说明：在半径为说明：在半径为D D0 0的圆内，所有频率没有衰减地通过的圆内，所有频率没有衰减地通过滤波器，而在此半径的圆之外的所有频率完全被衰滤波器，而在此半径的圆之外的所有频率完全被衰减掉减掉第19页/共59页其中其中原点在频率域的中心，半径为原点在频率域的中心，半径为D D0 0的圆包含的圆包含 %的功率的功率100(

14、, )/(4.84)TuvP u vP 理想低通滤波器理想低通滤波器1100( , )(4.83)QPTuvPP u v 222( , )( , )( , )( , )P u vF u vR u vIu v 总图像功率值总图像功率值P PT T其中其中: :第20页/共59页理想低通滤波器举例理想低通滤波器举例87%87%以上的功率以上的功率( (能量能量) )集中在半径小于集中在半径小于1010的圆周内的圆周内; ;随滤波器半径的增加随滤波器半径的增加, ,越来越少的功率被滤出掉越来越少的功率被滤出掉, ,使模使模糊减弱糊减弱; ;第21页/共59页原始图D0=10的ILPF滤波 损失能量为

15、8%D0=30的ILPF滤波 损失能量为5.4%D0=460的ILPF滤波 损失能量为0.5%D0=160的ILPF滤波 损失能量为2%D0=60的ILPF滤波 损失能量为3.6%理想低通滤波器举例理想低通滤波器举例第22页/共59页a ) a ) 半径为半径为1010的频的频域域ILPFILPFb ) b ) 半径为半径为1010空域空域ILPFILPFc ) c ) 图像图像b b）的水平）的水平 扫描线灰度变化扫描线灰度变化理想低通滤波器举例理想低通滤波器举例具有振铃现象具有振铃现象第23页/共59页n n阶布特沃思低通滤波器阶布特沃思低通滤波器(BLPF)(BLPF)定义如定义如下下2

16、01( , )(4.58)1 ( , )BLPFnHu vD u vD D D0 0为截至频率距原点的距离，为截至频率距原点的距离，D(u,v)D(u,v)是点是点(u,v)(u,v)距原点的距离。距原点的距离。1222( , )()()(4.82)22PQD u vuv 第24页/共59页它的特性是连续性衰减，而不象理想滤波器那样陡峭变化它的特性是连续性衰减，而不象理想滤波器那样陡峭变化，即明显的不连续性。因此采用该滤波器滤波在抑制噪声的即明显的不连续性。因此采用该滤波器滤波在抑制噪声的同时，图像边缘的模糊程度大大减小，没有振铃效应产生同时，图像边缘的模糊程度大大减小，没有振铃效应产生第25

17、页/共59页原始图D0=10的BLPF滤波D0=30的BLPF滤波 D0=460的BLPF滤波D0=160的BLPF滤波D0=60的BLPF滤波布特沃斯低通滤波器举例布特沃斯低通滤波器举例第26页/共59页布特沃斯低通滤波器举例布特沃斯低通滤波器举例振铃现象振铃现象注：二阶注：二阶BLPFBLPF处于有效低通滤波和可接受的振铃特征之间处于有效低通滤波和可接受的振铃特征之间阶数阶数n=n=1 1无振铃和负值无振铃和负值阶数阶数n=n=2 2轻微轻微振铃和负值振铃和负值 阶数阶数n=n=5 5明显明显振铃和负值振铃和负值阶数阶数n=20n=20与与ILPFILPF相似相似第27页/共59页二维高斯

18、低通滤波器二维高斯低通滤波器(GLPF)(GLPF)定义如下定义如下220( , ) 2( , )(4.87)Du vDH u ve 当当D(u,v)=DD(u,v)=D0 0时，滤波器下降到它最大值的时，滤波器下降到它最大值的0.6070.607处处第28页/共59页原始图D0=10的GLPF滤波D0=30的GLPF滤波 D0=460的GLPF滤波D0=160的GLPF滤波D0=60的GLPF滤波高斯低通滤波器举例高斯低通滤波器举例第29页/共59页字符识别举例字符识别举例原始图像原始图像D0=80的高斯低通滤波器的高斯低通滤波器修复字符修复字符用于机器识别系统识别断裂字符的预处理用于机器识

19、别系统识别断裂字符的预处理D0=80GLPF第30页/共59页人脸图像处理人脸图像处理原图像原图像D0=100的的GLPF滤波滤波,细纹减少细纹减少D0=80的的GLPF滤波，滤波，细纹减少细纹减少第31页/共59页D0=30的GLPF滤波图D0=10的GLPF滤波图佛罗里达亮墨西哥湾暗奥基乔 比湖卫星、航拍图像处理卫星、航拍图像处理目的：尽可能模糊细节，而保留大的可识别特征目的：尽可能模糊细节，而保留大的可识别特征第32页/共59页频率域锐化滤波器主要有：频率域锐化滤波器主要有：理想高通滤波器理想高通滤波器布特沃思高通滤波器布特沃思高通滤波器高斯高通滤波器高斯高通滤波器频率域的拉普拉斯算子频

20、率域的拉普拉斯算子钝化模板、高频提升滤波和高频加强滤波钝化模板、高频提升滤波和高频加强滤波第33页/共59页几种高通滤波器图示：几种高通滤波器图示：IHPFBHPFGHPF第34页/共59页IHPFBHPFGHPF几种高通滤波器空域图示：几种高通滤波器空域图示：第35页/共59页4.9.1理想高通滤波器（理想高通滤波器（IHPF）截断傅里叶变换中所有处于指定距离截断傅里叶变换中所有处于指定距离D D0 0之内的低频成之内的低频成分分0000( , )( , ),0(4.92)1( , )IHPFD u vDHu vDD u vD 频率域的中心在频率域的中心在 ，从点，从点(u,v)(u,v)到

21、中心（原点）到中心（原点）的距离如下的距离如下1222( , )()()(4.82)22PQD u vuv (,)2 2P Q第36页/共59页D0=30D0=60D0=160结论：图结论：图a a和和b b的振铃问题十分明显的振铃问题十分明显理想高通滤波示例：理想高通滤波示例：第37页/共59页n n阶巴特沃思高通滤波器阶巴特沃思高通滤波器(BHPF)(BHPF)定义如定义如下下201( , )(4.93)1 ( , )BHPFnHu vDD u v 20202001( , )1( , )11( , )/( , )/11( , )/1/( , )hplpnnnHu vHu vD u vDD

22、u vDD u vDDD u v 第38页/共59页D0=30D0=60D0=160二阶巴特沃思高通滤波示例：二阶巴特沃思高通滤波示例：第39页/共59页截频距原点为截频距原点为D0的高斯高通滤波器的高斯高通滤波器(GHPF)定义定义为为220( , )/2( , )1(4.94)Du vDGHPFHu ve 第40页/共59页D0=30D0=60D0=160高斯高通滤波示例：高斯高通滤波示例：第41页/共59页1026 962 0962 5%50D 二值化的结果二值化的结果第42页/共59页( , )H u v第43页/共59页频率域的拉普拉斯算子定义频率域的拉普拉斯算子定义:( )( 2)

23、( )nnnd f xju F udx ( , )( , )()( , )()( , )() ( , )nnd f x yd f x yju F u vjv F u vdxdyuv F u v 22222222242222( , )4()( , )f x yuvF u v 222( , )4()(4.95)H u vuv 第44页/共59页原点从原点从(0,0)(0,0)移到移到(P/2,Q/2),(P/2,Q/2),所以，滤波函数平移为所以，滤波函数平移为22222( , )4(2)(2)4( , )(4.96)H u vuPvQD u v 21( , )( , )( , )(4.97)f

24、x yH u vF u v 从原始图像中减去拉普拉斯算子部分，形成从原始图像中减去拉普拉斯算子部分，形成g(x,y)g(x,y)的增强图的增强图像像2122( , )( , )( , )(4.98)14( , ) ( , )(4.99)g x yf x yf x yD u v F u v 空间域拉普拉斯算子过滤后的图像可由计算空间域拉普拉斯算子过滤后的图像可由计算 的的反傅里叶变换得到反傅里叶变换得到( , ) ( , )H u v F u v4.9.4 频率域的拉普拉斯算子频率域的拉普拉斯算子第45页/共59页拉普拉斯锐化举例说明拉普拉斯锐化举例说明.例例4.20( , )f x y2( ,

25、 )( , )f x yf x y 21( , )( , )( , )(4.97)f x yH u vF u v 其其中中：第46页/共59页钝化模板钝化模板(锐化或高通图像锐化或高通图像)：在频率域，即从图像本身减去低通滤波（模糊）后的在频率域，即从图像本身减去低通滤波（模糊）后的图像而得到高通滤波（锐化）的图像。图像而得到高通滤波（锐化）的图像。1( , )( , )( , )(4.9 10)( , )( , ) ( , )(4.9 11)maskLPLPLPgx yf x yfx yfx yHu v F u v 钝化模板和高频提升过滤：钝化模板和高频提升过滤：11( , )( , )(

26、, )(4.912)1*1( , )( , )(4.913)( , )1*( , )( , )(4.914)maskLPHPg x yf x ykgx ykHu vF u vg x ykHu vF u v 第47页/共59页更一般的高频提升加强：更一般的高频提升加强：用图像的高频成分进行增强用图像的高频成分进行增强增加增加k1的目的是使零频率不被滤波器过滤的目的是使零频率不被滤波器过滤当当k21，高频得到加强，高频得到加强 12112( , )*( , )( ,00)(4.915)HPg x ykkHu vF u vkk 其其中中：控控制制距距原原点点的的偏偏移移量量；控控制制高高频频的的贡贡

27、献献第48页/共59页高频提升过滤举例高频提升过滤举例频率域滤波频率域滤波图图a: 胸部胸部X光图像光图像图图b: Gaussian高通滤波的结果高通滤波的结果图图c: 高频增强滤波的结果高频增强滤波的结果图图d: 图图c直方图均衡化的结果直方图均衡化的结果高频加强高频加强 直方图均衡化的混合图像增强方法直方图均衡化的混合图像增强方法第49页/共59页4.9.6 同态滤波同态滤波( , )( , )()f x yi x yr x,y( , )ln ( , )ln ( , )ln ()z x yf x yi x yr x,y( , )( , )( , )irZ u vF u vF u v( ,

28、)ln ( , )rF u vr x y( , )( , ) ( , )S u vH u v Z u v( , )ln ( , )iF u vi x y( , )( , )Z u vz x y( , ) ( , )( , )( , )irH u v F u vH u v F u v一幅图像可以表示为照射量和反射量的乘积：一幅图像可以表示为照射量和反射量的乘积：令：令：则有：则有：对对Z(uZ(u，v)v)进行滤波：进行滤波：第50页/共59页( , )( , )S u vs x y( , )( , )( , )rr x yH u v F u v( , )( , ) ( , )ii x yH u

29、 v F u v( , )( , )( , )s x yr x yi x y同态滤波增强后的图像同态滤波增强后的图像 g g（x x，y y）为：）为：( , )( , )( , )( , )s x yr x yi x yg x yeee( , )( , )i x yr x y00在空间域，滤波后的图像为：在空间域，滤波后的图像为：第51页/共59页同态滤波中滤波器函数的选择：同态滤波中滤波器函数的选择：滤波器函数能够实现对照射分量和反射分量的不同操作滤波器函数能够实现对照射分量和反射分量的不同操作而加强图像，以不同的方式影响傅里叶变换的高频和低而加强图像，以不同的方式影响傅里叶变换的高频和低

30、频成分。滤波器函数如图示：频成分。滤波器函数如图示：1L减少低频减少低频 增加高频增加高频1H( , )/)( , )()c Du v DHLLH u ve2201一种同态滤波器函数表达式：一种同态滤波器函数表达式：第52页/共59页变化最慢的频率成分变化最慢的频率成分(u=v=0)(u=v=0)对应一幅图像的平均灰度对应一幅图像的平均灰度级级( , )( , )( , )( .)M 1 N 1x 0 y 0F 0 0f x yMN f x y4 621 当从变换的原点移开时，对低频对应着图像的慢变化分当从变换的原点移开时，对低频对应着图像的慢变化分量，如图像的平滑部分量，如图像的平滑部分进一步离开原点时，较高的频率对应图像中变化越来进一步离开原点时，较高的频率对应图像中变化越来越快的灰度级，如边缘或噪声等尖锐部分越快的灰度级，如边缘或噪声等尖锐部分4.7.1、频率域的其他特性：、频率域的其他特性：()( , )( , )ux vyM 1 N 1i2MNx 0 y 0F u vf x y e 第53页/共59页高通滤波器：高通滤波器：被高通滤波的图像比原始图像少了灰度级的平被高通滤波的图像比原始图像少了灰度级的平滑过渡而突出了边缘等细节部分滑过渡而突出了边缘等细节部分使高频通过，低频衰减