高通滤波代码

%高通滤波
I=imread('lv.jpg');
I=rgb2gray(I);
[f1,f2]=freqspace(size(I),'meshgrid');  %分割频率空间
Hd=ones(size(I));
r=sqrt(f1.^2+f2.^2);   %计算到原点的r
Hd(r<0.1)=0;  %频率低于0.1就置变换矩阵为0
Y=fft2(I);
Y=fftshift(Y);
Ya=Y.*Hd;
Ya=ifftshift(Ya);
Ya=ifft2(Ya);

imshow(uint8(Ya));