作业1（高斯模糊，椒盐噪声，中值滤波去噪）

%% imread 得到黑白原图im(图一)

>> im=rgb2gray(im);

>> im=im2double(im);

>> imshow(im);

 

%% 高斯核模糊 得到模糊后的图片B(图二)

方法一

>> sigma=11;

>> hsize=10;

>>h=fspecial('gaussian',hsize,sigma);

>> B=imfilter(im,h,'symmetric');

>> imshow(B);

方法二： 

    sigma=11;

>> hsize=10;

       x=floor((hsize+1)/2);  

     for i=1:sigma  

       forj=1:sigma  

       h(i,j)=exp(-((i-x)^2+(j-x)^2)/(sigma*sigma))/(2*pi*sigma*sigma);

       end  

     end  

B=zeros(1184,720);

tem=zeros(11,11);

tem2=zeros(11,11);

for m=6:1184-5

for n=6:720-5

       tem=im(m-5:m+5,n-5:n+5);

           tem2=tem.*h;

       B(m,n)=sum(tem2(:));

end

end

imshow(B,[]);

%% 原图以30%概率加上椒盐噪声  得到图C(图三)

方法一

>> C=imnoise(im,'salt &pepper',0.3);

>> imshow(C);

 

方法二

>> C=im;

>>

for i=1:1184;

for j=1:720;

if rand >=0.7

  ifrand>=0.5

 C(i,j)=255;

 else

 C(i,j)=0;

  end

end

end

end

>> imshow(C);

 

 

%% 将C中值滤波,得到图median_img(图四)

>> disp(size(C));

median_img=zeros(1184,720);

temp=zeros(5,5);

for i=3:1184-2

for j=3:720-2

       temp=C(i-2:i+2,j-2:j+2);

       median_img(i,j)=median(temp(:));

end

end

imshow(median_img);