clear,clc,close all
I=imread('lena.bmp');figure,imshow(I);
%對圖像進行掃描
[m,n]=size(I);%圖像的寬度和高度
W=ones(3,3);W=W/9;%濾波模板
J=zeros(m,n);
%對邊緣不處理
J(1:m,1)=I(1:m,1);
J(1:m,n)=I(1:m,n);
J(1,2:n-1)=I(1,2:n-1);
J(m,2:n-1)=I(m,2:n-1);
%對不是邊緣點處理
for i=2:m-1
    for j=2:n-1
        k=1;
        for a=-1:1
            for b=-1:1
                S(k)=I(i+a,j+b);
                k=k+1;
            end
        end
        J(i,j)=median(S);
    end
end
figure,imshow(uint8(J));%處理后的圖像