实验 4 图像的滤波处理与图像增强.pdf

实验四 图像的滤波处理与图像增强 实验目的 1. 了解MATLAB工具箱中的滤波器。 2. 掌握空间滤波 3. 学会对图像的空间变换 实验内容 A. 用滤波器祛除图象噪声 在数字图像处理中，常常会遇到图像中混杂有许多的噪声。因此，在进行图像处 理中，有时要先进行祛除噪声的工作。最常用的祛除噪声的方法是用滤波器进行 滤波处理。MATLAB 的图像处理工具箱里也设计了许多的滤波器。如均值滤波器、 中值滤波器、维纳滤波器等。 （分别用均值滤波，中值滤波，及维纳滤波器祛除加入高斯噪声的图象） I=imread(D:\图像库\DIP3E_CH04\FigP0438(left).tif ); J=imnoise(I,gaussian,0,0.002); %进行均值滤波 h=fspecial(average,3); I2=uint8(round(filter2(h,I))); %进行中值滤波 I3=medfilt2(J,[3,3]); %进行维纳滤波 I4=wiener2(J,[3,3]); %进行一次维纳滤波 I5=wiener2(I4,[3,3]);%进行二次维纳滤波 subplot(2,3,1),imshow(I),title(原图象) subplot(2,3,2),imshow(J),title(加噪声图象) subplot(2,3,3),imshow(I2),title(均值滤波后图象) subplot(2,3,4),imshow(I3),title(中值滤波后图象) subplot(2,3,5),imshow(I4),title(维纳滤波后图象) subplot(2,3,6),imshow(I5),title(两次维纳滤波后图象) B. 空间噪声滤波器 %用函数imnoise2生成具有表5.1中的CDF的随机数 function R=imnoise2(type,M,N,a,b) if nargin 1 a=0;b=1; M=1;N=1; elseif nargin 3 a=0;b=1; end switch lower(type) case uniform R=a+(b-a)*rand(M,N); case gaussian R=a+b*randn(M,N); case salt pepper if nargin =3 a=0.05;b=0.05; end if (a+b)1; error(The sum Pa+Pb must not exceed 1.) end R(1:M, 1:N) = 0.5; X=rand(M,N); c=find(X=a); R(c)=0; u=a+b; c=find(Xa X=u); R(c)=1; case rayleigh R=a+(-b*log(1-rand(M,N))).^0.5; case exponential if nargin =3; a=1; end if a=0 error(Parameter a must be positive for exponential type.) end k=-1/a; R=k*log(1-rand(M,N)); case erlang if nargin=3 a=2;b=5; end if (b~=round(b)|b=0) error(Parameter b must be a positive integer for Erlang) end k=-1/a; R=zeros(M,N); for j=