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');
matlab-digital-image.rar （45个子文件） 
matlab digital image matlab源代码 
road.m 1KB quzao.m 2KB lena测试图片 
本节算法去噪10%噪声.bmp 231KB lena5%.bmp 141KB x4.bmp 65KB lena40%噪声.bmp 65KB lena10%.bmp 141KB 40%噪声去除.bmp 65KB lena20%.bmp 141KB 3 3 10%.bmp 65KB p3.bmp 65KB x3.bmp 65KB lena10%噪声.bmp 65KB x5.bmp 65KB 30%噪声去除.bmp 65KB
lena1.jpg 63KB lena20%噪声.bmp 65KB 20%噪声去除.bmp 65KB lena40%.bmp 141KB untitled.bmp 231KB p6.bmp 65KB p1.bmp 65KB 5%噪声去除.bmp 65KB 经中值滤波后.bmp 65KB x6.bmp 65KB p4.bmp 65KB 7 7 10%.bmp 65KB lena50%.bmp 141KB lena2.bmp 65KB 9 9 10%.bmp 65KB lena5%噪声.bmp 65KB lena50%噪声.bmp 65KB 原始lena图像（含椒盐噪声）.bmp 65KB p5.bmp 65KB p2.bmp 65KB 5 5 10%.bmp 65KB lena30%噪声.bmp 65KB 10%噪声去除.bmp 65KB x2.bmp 65KB x1.bmp 65KB lena30%.bmp 141KB lena3.bmp 65KB lena.jpg 90KB 50%噪声去除.bmp 65KB a.bmp 65KB
评论0