s114.m

简单的数字图像处理程序

[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)