MATLAN图像处理之高斯低通滤波

%图像平滑(低通滤波)

f=imread('D:\研究生课程\其他\MATLAB图像处理\亮度调节与空间滤波\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)));%高斯低通滤波器

g=dftfilt(f,H);
subplot(2,2,1)
imshow(fftshift(H),[])
title('图像频谱')

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

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

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