matlab練習程序（圖像加噪，各種濾波）

cl;

img=imread('3.17.tif');
imgn=imnoise(img,'salt & pepper',0.02);
[m n]=size(img);

%h=fspecial('average',[3 3]);    %   3*3的均值濾波
%imgn1=imfilter(imgn,h,'replicate');

%imgn=medfilt2(imgn,[3 3]);     %    3*3的中值濾波

%%%%自己編寫3*3均值濾波
%{
tem=0;
for i=2:1:m-1
    for j=2:1:n-1     
   
        for p=i-1:1:i+1
            for q=j-1:1:j+1
                tem=tem+double(imgn(p,q));
            end
        end

        tem=tem/9;
        imgn(i,j)=tem;
        tem=0;
    end
end
%}

%%%%自己編寫3*3中值濾波
%{
s=zeros(1,9);
for i=2:1:m-1
    for j=2:1:n-1     
        
        he=1;
        for p=i-1:1:i+1
            for q=j-1:1:j+1            
                s(he)=imgn(p,q);
                he=he+1;               
            end
        end

        s=sort(s);
        imgn(i,j)=s(5);
        
    end
end
%}
imshow(imgn)