jiaoyan yuantuzhugechuli.m

来自「图像椒盐噪声去除研究程序」· M 代码 · 共 53 行

A=imread('lena.bmp');
subplot(2,3,1);

imshow(A);

title('原图');
B = imnoise( A,'salt & pepper',0.56);
subplot(2,3,2);
imshow(B);

title('用噪声函数产生的椒盐噪声');
imwrite( B,'lenaNoiseFun.tif');
[M,N]=size(B);
outImg1=ones(M,N);
for i=2:M-1
    for j=2:N-1
        Model=B(i-1:i+1,j-1:j+1);
        ord=sort(Model(:),1);
        outImg1(i,j)=ord(5,1);
    end
end
subplot(2,3,3);
imshow(outImg1,[0,255]);
title('传统滤波处理后的图象');
[m,n]=size(B);
outImg2=ones(m,n);
for i=3:m-2
    for j=3:n-2
        Model=B(i-2:i+2,j-2:j+2);
        ord=sort(Model(:),1);
        outImg2(i,j)=ord(13,1);
    end
end
subplot(2,3,4);
imshow(outImg2,[0,255]);
title('5*5滤波处理后的图象');

[m1,n1]=size(B);
outImg3=ones(m,n);
for i=4:m-3
    for j=4:n-3
        Model=B(i-3:i+3,j-3:j+3);
        ord=sort(Model(:),1);
        outImg3(i,j)=ord(25,1);
    end
end
subplot(2,3,5);
imshow(outImg3,[0,255]);
title('7*7滤波处理后的图象');