图像处理第三次实验.docxVIP

实验名称 ：图像的空域滤波处理（2015-6-10 1:30-4:30pm 计算中心）一、实验目的：在图像的空域滤波实现去噪与边缘提取的原理。实验平台：Matlab2014a实验内容： 1、具体要求：对一幅灰度图像（lena.bmp）分别进行均值低通滤波和拉普拉斯高通滤波。 2、实验要点：理解图像像素的数学运算以及图像滤波的原理。 3、实验原理：频域低通滤波的实现：先对合适大小的图像进行二维离散傅里叶变换，然后利用滤波器进行滤波；均值：利用3*3 5*5 7*7模板，例：K1=filter2(fspecial(average,3),I)/255;对图像进行均值滤波；拉普拉斯，用拉普拉斯模板进行运算实验结果Matlab程序：频域低通滤波：clc;clear all;close all;f=imread(D:/lena.bmp);[row,col]=size(f);g=fft2(f);gls=fftshift(g);H=zeros(row,col);x0=floor(row/2);y0=floor(col/2);u0=x0;v0=y0;D0=50;n=0.2for u=1:row for v=1:col D=(u-u0)^2+(v-v0)^2; H(u,v)=1/(1+(D/D0)^(2*n)); endendrss=gls.*H;rs=ifftshift(rss);r=real(ifft2(rs));subplot(2,2,1);imshow(f);subplot(2,2,2);imshow(uint8(r));subplot(2,2,3);H1=H(x0,y0:col);plot(H1);axis([0 y0 min(H1) max(H1)]);均值：clear all;I=imread(D:/lena.bmp);%读入预处理图像imshow(I)%显示预处理图像K1=filter2(fspecial(average,3),I)/255;%进行3*3均值滤波K2=filter2(fspecial(average,5),I)/255;%进行5*5均值滤波K3=filter2(fspecial(average,7),I)/255;%进行7*7均值滤波subplot(2,2,1) imshow(I) title(原图像) subplot(2,2,2) imshow(K1)title(3*3) subplot(2,2,3) imshow(K2) title(5*5) subplot(2,2,4) imshow(K3)title(7*7)拉普拉斯：clear all;clcf=imread(D:/lena.bmp);w=fspecial(laplacian,0);g1=imfilter(f,w,replicate);f2=im2double(f);g2=imfilter(f2,w,replicate);figure,imshow(f,[])figure,imshow(g2,[])g=f2-g2;figure,imshow(g)w8=[1 1 1;1 -8 1;1 1 1];w8=[1 1 1 1 -8 1 1 1 1];g8=f-imfilter(f,w8,replicate);figure,imshow(g8)实验截图：频域低通滤波：均值：拉普拉斯：实验心得：对理解图像像素的数学运算以及图像滤波的原理有了更为清晰直观的认识，对于图像处理这门学科的学习有很大的帮助，实验中学习了这次实验所涉及的函数，能较为熟练的应用，对各个对图像进行的不同处理能起到的不同作用也有了直观的认识，能更为准确的判断在不同的噪声下应该选用什么进行处理。