matlab-digital-image.rar_Impulse noise image_去除图片噪声_去除脉冲噪声_脉冲噪声

clc clear; I=imread('f:\lena\lena2.bmp'); [m,n]=size(I); figure(1); subplot(2,2,1),imshow(I),title('源图像'); I1=imnoise(I,'salt & pepper',0.3); %添加椒盐脉冲噪声 subplot(2,2,2),imshow(I1);title('加入强度为0.3的脉冲噪声'); % %常规中值滤波 tic; Ic=medfilt2(I1); toc; subplot(2,2,3),imshow(Ic);title('常规中值滤波效果图'); imwrite(Ic,'F:\lena\x6.bmp','bmp'); tic; Temp=zeros(3,3); Ij=I1; Iz=I1; w=[1 2 3 4 5 4 3 2 1 ]; for i=2:m-1 for j=2:n-1 % clear Temp ; % clear T; k=0; Temp = Iz(i-1:i+1,j-1:j+1); T=Temp; while 1 [Tv,ind] = sort(double(T(:))); zmin=Tv(1); zmax=Tv(9); zmed=Tv(5); zxy=double(Iz(i,j)); a1=zmed-zmin; a2=zmed-zmax; if a1>0 & a2<0 b1=zxy-zmin; b2=zxy-zmax; if b1>0 & b2<0 Iz(i,j)=zxy; else Iz(i,j)=zmed; end break; else k=k+2; if i-(1+k)<1 || j-(1+k)<1 || i+(1+k)>m || j+(1+k)>n Iz(i,j)=zxy; break; else T=Iz(i-(1+k):i+(1+k),j-(1+k):j+(1+k)); if k>5 Iz(i,j)=zxy; break; end end end end end end toc; subplot(2,2,4),imshow(Iz);title('自适应开关中值滤波效果图'); imwrite(Iz,'F:\lena\x5.bmp','bmp');