- 1、本文档共20页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
数字图像处理三级项目
—咼通、低通、带通滤波命
摘要
在图像处理的过程中,消除图像的噪声干扰是一个非常重要的问题。利用 matlab软件,采用频域滤波的方式,对图像进行低通和高通滤波处理。低通滤波 是要保留图像中的低频分量而除去高频分量, 由于图像中的边缘和噪声都对应图
像傅里叶频谱中的高频部分,所以低通滤波可以除去或消弱噪声的影响并模糊边 缘轮廓;高通滤波是要保留图像中的高频分量而除去低频分量, 所以高通滤波可
以保留较多的边缘轮廓信息。低通滤波器有巴特沃斯滤波器和高斯滤波器等等, 本次设计使用的低通滤波器为****。高通滤波器有巴特沃斯滤波器、高斯滤波器、 Laplacian高通滤波器以及Unmask高通滤波器等等,本次设计使用巴特沃斯高 通滤波器。
1、频域低通滤波器:设计低通滤波器包括 butterworth and Gaussian ( 选择
合适的半径,计算功率谱比),平滑测试图像test1和2。
实验原理分析
根据卷积定理,两个空间函数的卷积可以通过计算两个傅立叶变换函数的 乘积的逆变换得到,如果f(x, y)和h(x, y)分别代表图像与空间滤波器, F(u, v)和H(u, v)分别为响应的傅立叶变换(H(u, v)又称为传递函数), 那么我们可以利用卷积定理来进行频域滤波。
在频域空间,图像的信息表现为不同频率分量的组合。如果能让某个范围 内的分量或某些频率的分量受到抑制,而让其他分量不受影响,就可以改变 输出图的频率分布,达到不同的增强目的。
频域空间的增强方法的步骤:
将图像从图像空间转换到频域空间;
在频域空间对图像进行增强;
将增强后的图像再从频域空间转换到图像空间
低通滤波是要保留图像中的低频分量而除去高频分量。图像中的边缘和
噪声都对应图像傅里叶频谱中的高频部分,所以低通滤波可以除去或消弱噪
声的影响并模糊边缘轮廓。理想低通滤波器具有传递函数:
其中DO为制定的非负数,D(u,v)为点(u,v)到滤波器中心的距离
功率谱比的定义:
功率谱比的定义:
Pg(u,v)
Pf(u,v)
其中,Pf (u,v)为滤波前图像的功率谱,Pg(u, V)为滤波后图像的功率谱。
频率计算公式为:Pf(u,v) F(u,v)2, Pg (u,v) G(u, v)2。
①Butterworth 滤波器设计:
理想低通滤波器在数学上定义得很清楚, 在计算机模拟中也可实现,但在截断频率处直上直下的理想低通滤波器是不能用实际的电子器件实现的。
n阶Butterworth 低通滤波器(BLPF的传递函数
n阶Butterworth 低通滤波器(BLPF的传递函数
离为D。)定义如下:
H(u,v)
D(u,v)/ Do
2n
(截止频率距原点的距
(1)
其中,D(u,v)= .(u-M/2)2+ (v-N/2)2。 ( 2)
不同于ILPF,BLPF变换函数在通带与被滤除的频率之间没有明显的截断。 对 于有平滑传递函数的滤波器,定义一个截止频率的位置并使 H(u,v)幅度降到其 最大值的一部分。在式(1)中,当D(u,v)=D0时,H(u,v)=0.5 (从最大 值降到它的50%。
一阶的巴特沃斯滤波器没有振铃,在二阶中振铃通常很微小,这是因为与理 想低通滤波器相比,它的通带与阻带之间没有明显的跳跃, 在高低频率间的过渡
比较光滑。巴特沃斯低通滤波器的处理结果比理想滤波器的要好, 但阶数增高时
振铃便成为一个重要因素。本次实验中设计实现了二阶巴特沃斯滤波器。
根据以上原理设计Butterworth低通滤波器,其处理结果如下图示:
源图像test1.pgmButterworth 低通 滤波图
源图像test1.pgm
Butterworth 低通 滤波图 test1.pgm
源图像
源图像test2.tif Butterworth 低通滤波图 test2.tif
源图像
源图像test2.tif Butterworth 低通滤波图 test2.tif
理想低通滤波器有明显的振铃现象,而巴特沃斯滤波器的效果较好。
计算得testl的功率谱比L=0.9939。test2的功率谱比为0.9902
②Gaussian滤波器设计:
二维高斯低通滤波器,其传递函数的形式为:
2 2
TOC \o "1-5" \h \z H(u,v)=1-e D (u,v)/2 ( 3)
其中,D(u,v)= ,( u-M/2)2 +( v-N/2)2。 表示高斯曲线扩展的程度。使 =D0,
可以将滤波器表示为:
H(u,v)=e D2(u,v)/2D02 (4)
其中,D0是截止频率。当D (u,v)=D0时,滤波器下降到它最大值的0.607倍 处。
由于高斯低通滤波器的傅里叶反变换也是高斯的,这就是说通过式( 3)或
式(4)的傅里叶反变换
文档评论(0)