test.m

来自「在一维信号中」· M 代码 · 共 38 行

% 读取原始图像并显示
   f0=imread('peppers.bmp');
   figure(1)
   imshow(f0); %如图4-1所示
   
  %加入高斯噪声并显示
  f1=imnoise (f0, 'gaussian', 0, 0.05);
  figure(2)
  imshow(f1); %如图4-2所示

  % 给定四种不同的滤波窗口
  h1=[-1 -1 -1; -1 9 -1; -1 -1 -1 ];
  h2=[1 1 1; 1 2 1; 1 1 1];
  h3=[1/9 1/9 1/9; 1/9 1/9 1/9; 1/9 1/9 1/9 ];
  h4=[-1 -1 -1; -1 8 -1; -1 -1 -1];

 %采用h1对图像f2进行卷积滤波
 f2=conv2(f1, h1,same);
 figure(3)
 imshow(f2); %如图4-3所示

%采用h2对图像f2进行卷积滤波
 f3=conv2(f1, h2,same);
 figure(4)
 imshow(f3); %如图4-4所示
 
%采用h3对图像f2进行卷积滤波
 f4=conv2(f1, h3,same);
 figure(5)
 imshow(f4); %如图4-5所示
   
%采用h4对图像f2进行卷积滤波
 f5=conv2(f1, h4,same);
 figure(6)
 imshow(f5); %如图4-6所示