filfs.m

filfs滤波器的程序

% 读取原图像
   g0=imread('eight.tif');
   figure(1)
   imshow(g0) ; %如图4-7所示

   % 加入椒盐噪声
   g1=imnoise(g0, 'salt and peppers', 0.02);
   g1=im2double(g1);
   figure(2)
   imshow(g1); %如图4-8所示
 
   % 进行高斯低通滤波
   h1=fspecial('gaussian', 4, 0.3);
   g2=filter2(h1,g1, 'same');
   figure(3)
   imshow(g2); %如图4-9所示

   % 进行sobel 滤波
   h2=fspecial('sobel');
   g3=filter2(h2,g1, 'same');
   figure(4)
   imshow(g3); %如图4-10所示
 
  % 进行 prewitt滤波
   h3=fspecial('prewitt');
   g4=filter2(h3,g1, 'same');
   figure(5)
   imshow(g4); %如图4-11所示

   % 进行拉普拉斯滤波
   h4=fspecial('laplacian', 0.5)
   g5=filter2(h4,g1, 'same');
   figure(6)
   imshow(g5); %如图4-12所示
 
   % 进行高斯拉普拉斯滤波
   h5=fspecial('log', 4, 0.3)
   g6=filter2(h5,g1, 'same');
   figure(7)
   imshow(g6); %如图4-13所示

   % 进行均值滤波
   h6=fspecial('average')
   g7=filter2(h6,g1, 'same');
   figure(8)
   imshow(g7); %如图4-14所示

  % 进行模糊增强滤波
   h7=fspecial('unsharp',0.3)
   g8=filter2(h7,g1, 'same');
   figure(9)
   imshow(g8); %如图4-15所示

  % 进行高斯高通滤波
   h8=[0 -1 0; -1 5 -1; 0 -1 0];
   g9=filter2(h8,g1, 'same');
   figure(10)
   imshow(g9); %如图4-16所示