s114.m

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

[f1,f2]=freqspace(21,'meshgrid');
Hd=ones(21);
r=sqrt(f1.^2+f2.^2);
Hd(r<0.5)=0;
win=fspecial('gaussian',21,2);
win=win./max(win(:));
h=fwind2(Hd,win);
figure(1),freqz2(h);
a=imread('3.jpg');
a1=rgb2gray(a);
ff=filter2(h,a1);
figure(2),imshow(ff)