昆明理工大学数字图像处理实验二重点分析.docxVIP

昆明理工大学信息工程与自动化学院学生实验报告 （ 2012—2013学年 第 二 学期 ） 课程名称：数字图像处理 开课实验室：444 2015年 4月 17 日 年级、专业、班学号姓名成绩实验项目名称频率域滤波指导教师王剑教师评语 教师签名： 年 月 日实验目的及内容 内容： 频率域滤波 获得图像的傅里叶图 1、平滑 采用低通平滑滤波器（理想、buttworth、Gaussian）完成图像的平滑功能，并用mesh函数绘制滤波器 2、锐化 采用高通平滑滤波器（理想、buttworth、Gaussian）完成图像的平滑功能，并用mesh函数绘制滤波器 实验原理 在傅里叶图像中，高频位于图像的四角，低频位于图像中心，图像的边缘和细节主要在高频部分显示，对应到四角，而图像的背景和基调对应为图像的低频部分，在图像的中心。 频率域滤波的基本步骤： F=fft2(f) F:傅里叶图像 Fc=fftshift(F) Fc:中心化后的傅里叶图像 Gc=H(u,v)*Fc G:频域的输出图像 G=ifftshift(Gc) g:空间域的输出图像 g=ifft2(G) H:频域滤波器 g=real(g) real:输出图像的实部 ilpf(理想低通滤波器)：由于高频成分包含大象的图像边缘信息，因此ilpf在去噪的同事会导致边缘信息的损失，同时造成图像边缘模糊。 BLPF（巴特沃斯低通滤波器）：连续性衰减没有ILPF那样陡峭不连续，因此图像边缘模糊大大减少，但大大增加了计算量。 GLPF（高斯低通滤波器）：图像边缘模糊成都较ILPF小，而比BLPF大 理想低通滤波器的传递函数：Hu,v=1 Du,v≤D00 Du,vD0 N阶巴特沃斯低通滤波器的传递函数：Hu,v=11+[D(u,v)/D0]2n 高斯低通滤波器的传递函数：Hu,v=e-D2(u,v)/2D02 理想高通滤波器的传递函数：Hu,v=0 Du,v≤D01 Du,vD0 巴特沃斯高通滤波器的传递函数：Hu,v=11+[D0/D(u,v)]2n 高斯高通滤波器的传递函数：Hu,v=1-e-D2(u,v)/2D02 实验源代码及结果 平滑： 理想低通平滑滤波器 ff=imread(zoro.png) %读取zoro.png图像,赋值给ff f=rgb2gray(ff); %将rgb图像转换成灰度图像，赋值给f subplot(2,2,1),imshow(f),title(原图); F=fft2(f); %对图像进行傅里叶变换 subplot(2,2,2),imshow(F),title(傅里叶图); Fc=fftshift(F); [M N]=size(Fc); %分别返回Fc的行数到M中，列数到N中 m=floor(M/2); %对M/2进行取整 n=floor(N/2); %对N/2进行取整 D0=5; %初始化D0 for i=1:M for j=1:N d=sqrt((i-m)^2+(j-n)^2);%点（i,j）到傅里叶变换中心的距离 if d=D0; %点（i,j）在通带内的情况 H(i,j)=1; %通带变换函数 else H(i,j)=0; %阻带变换函数 end end end Gc=H.*Fc; %理想低通滤波器滤波后的频域表示 G=ifftshift(Gc); %去中心化 g=ifft2(G); %反傅里叶变换 g=uint8(real(g)); %显示图像的实部 subplot(2,2,3),mesh(H(1:10:300,1:10:300)),title(滤波器);%绘制滤波器 subplot(2,2,4),imshow(g),title(理想平滑图); 巴特沃斯平滑滤波器 ff