s110.m

简单的数字图像处理程序

%   产生高通滤波器
hd=zeros(20,20);
hd(1:10,1:10)=1;
[f1,f2]=freqspace(20,'meshgrid');
figure(2),mesh(f1,f2,hd);       %图示 高通滤波器
h=fwind1(hd,hanning(20));
figure(3),freqz2(h);    %如图示 高通滤波器
% 利用所设计滤波器进行图像滤波
I=imread('3.jpg');
I1=rgb2gray(I);
ff=filter2(h,I1);
figure(5),imshow(ff)  %如图示