s108.m

来自「简单的数字图像处理程序」· M 代码 · 共 15 行

%矩形窗设计二维带阻滤波器
%产生带阻滤波器
hd=ones(11,11);
hd(5:8,5:8)=0;
[f1,f2]=freqspace(11,'meshgrid');
figure(2),mesh(f1,f2,hd);         % 带阻滤波器如图所示  
h=fwind1(hd,boxcar(11));
figure(3),freqz2(h);      % 带阻滤波器如图所示  
%利用所设计滤波器进行图像滤波
I=imread('3.jpg');
I1=rgb2gray(I);
figure(4),imshow(I1);       %原图像如图所示                                   
ff=filter2(h,I1);
figure(5),imshow(ff);     %    滤波后图像如图所示