MATLAN图像处理之高通滤波

结果如下

 

 

%锐化频域滤波器
%若有低通的传递函数H（u，v）  高通滤波器可由下式得出
%H（u,v）=1-Hlp(u,v)
%这副图像的边缘和其他亮度得到了变化  由于图像的平均值由 0，0给处
%而且至今为止，高通滤波器偏离了傅里叶变换的原点  图像失去了大部分原图像所呈现的背景色彩
%这个问题在高频强调滤波器中讲明（下一个博客）
f=imread('D：\亮度调节与空间滤波\5.jpg');
f=f(:,:,1);
PQ=[size(f,1)*2,size(f,2)*2];
[U,V]=dftuv(PQ(1),PQ(2));
D0=0.05*PQ(2);%使用D0的值等于填充后图像宽度的5% 当D=D0时 截至频率
F=fft2(f,PQ(1),PQ(2));
H=exp(-(U.^2+V.^2)/(2*(D0^2)));%高斯低通滤波器
H=1-H;
g=dftfilt(f,H);
subplot(2,2,1)
imshow(fftshift(H),[])
title('滤波器频谱')

subplot(2,2,2)
imshow(log(1+abs(fftshift(H))),[])
title('滤波器lg拉伸后的频谱')

subplot(2,2,3)
imshow(g,[])
title('频域高通滤波后的图像')

subplot(2,2,4)
imshow(f,[])
title('原图')