图像处理(第五次).docxVIP

西安交通大学 图像与视频处理 第五次作业 姓名： 学号： 班级： 2011-4-24 摘要：本次作业主要针对图像的频率域滤波进行了研究。首先，介绍了频率域滤波的基本原理，其次介绍了典型的频率域滤波器，如平滑滤波器中的gaussian低通滤波器、巴特沃兹低通滤波器，锐化滤波器中的gaussian高通滤波器、巴特沃兹高通滤波器，unsharp mask, Laplace滤波，并对给定的图像进行测试，分析了图像的变化。最后，比较了空域与频域滤波的异同联系之处。  实验任务 频域低通滤波器：设计低通滤波器包括butterworth and Gaussian(选择合适的半径，计算功率谱比),平滑测试图像test1和2 频域高通滤波器：设计高通滤波器包括butterworth and Gaussian，在频域增强边缘。选择半径和计算功率谱比，测试图像test3,4： 其他高通滤波器：拉普拉斯和Unmask，对测试图像test3,4滤波； 比较并讨论空域低通高通滤波（Project3）与频域低通和高通的关系； 实验原理 二维离散傅里叶变换 令f(x,y)表示一幅大小为M*N的图像，其中处f的二维离散傅里叶变换可表示为F(u,v)，如下式所示： F 其中u=0,1,2…M-1和v=0,1,2…N-1。 在频率域滤波的步骤： 用(-1)x+y L 由（1）计算图像的DFT，即F(u,v)。 用滤波器函数H(u,v)乘以F(u,v)。 计算（3）中结果的反DFT。 得到（4）中结果的实部。 用-1x+y乘以（5 如下给出输出图像的傅里叶变换: G 被滤波的图像 如果两个滤波器是相同尺寸，那么通常在频率域进行滤波计算更为有效。但是，在空间域更适用于更小的滤波器。这正是我们所感兴趣的联系。滤波在频率域中更为直观，但在空间域适用更小的滤波器模板更为明智。我们可以在频率域制定滤波器，做反变换，然后再空间域使用结果滤波器作为在空间域构建更小空间滤波板得知道。 Butterworth低通滤波器 n阶巴特沃斯低通滤波器的传递函数定义如下： H 其中D(u,v)由式Du,v Gaussian低通滤波器 高斯低通滤波器的二维形式由下式给出： H 其中，如Du,v=[ 使σ=D H 其中，D0 Butterworth高通滤波器 n阶且截止频率距原点为D0 H 其中D(u,v)由式Du,v Gaussian高通滤波器 截止频率距原点为D0 H 其中D(u,v)由式Du,v 拉普拉斯高通滤波器 频域的拉普拉斯算子可以由如下滤波器实现： H unsharp mask和highboost filtering 思路为从原始图像中减去模糊地图像从而得到边缘锐化的图像。 g 其中Hlp 实验结果 设计butterworth低通滤波器，选择合适的半径，计算功率谱比,平滑测试图像test1和2。 截止频率为r=100的4阶巴特沃兹低通滤波器： 测试图像test1： 取半径为r=5时，功率谱比为p =0.9015 取半径为r=15时，功率谱比为p =0.9226 取半径为r=30，功率谱比为p = 0.9339 取半径为r=100，功率谱比为p =0.9537 测试图像test2： 取半径为r=15，功率谱比为p =0.9299 取半径为r=30，功率谱比为p =0.9349 取半径为r=100，功率谱比为p =0.9426 设计Gaussian低通滤波器，选择合适的半径，计算功率谱比,平滑测试图像test1和2 截止频率为r=100的高斯低通滤波器： 测试图像test1： 取半径为r=10，功率谱比为p =0.8615 取半径为r=30，功率谱比为p =0.9007 取半径为r=100，功率谱比为p =0.9327 取半径为r=15，功率谱比为p =0.8964 取半径为r=50，功率谱比为p =0.9126 取半径为r=100，功率谱比为p =0.9137 设计butterworth高通滤波器，在频域增强边缘。选择半径和计算功率谱比，测试图像test3,4： 巴特沃兹高通滤波器三维曲面： 测试test3： 取半径为r=5，功率谱比为p = 0.0154 取半径为r=15，功率谱比为p =0.0315 取半径为r=100，功率谱比为p= 0.1839 测试test4： 取半径为r=15，功率谱比为p =0.0122 取半径为r=30，功率谱比为p =0.0226 取半径为r=100，功率谱比为p =0.0635 设计Gaussian高通滤波器，在频域增强边缘。选择半径和计算功率谱比，测试图像test3,4： 高斯高通滤波器： 测试test3： 取半径为r=15，功率谱比为p =0.1250 取半径为r=30，功率谱比为p = 0.2207 取半径为r=100，功率谱比为p =