glpf.rar_GLPF_gaussian low pass_matlab 低通滤波器_滤波器_高斯低通滤波

function f=glpf(filename,r) % I=imread('E:\数字图像处理\images_chapter_04\Fig4.11(a).jpg'); I=imread(filename); Im=double(I); F=fft2(Im); F_result=fftshift(F); [m,n]=size(F_result); M=fix(m/2); N=fix(n/2); for u=1:m for v=1:n D=sqrt((u-M)^2+(v-N)^2); % H=1/(1+(D/r)^4); H=exp(-D^2/(2*r^2)); F_result(u,v)=F_result(u,v)*H; end end G_result=ifftshift(F_result); g_result=ifft2(G_result); f=real(g_result); f=uint8(f); % imshow(uint8(f)); %高斯低通滤波器 I=imread('E:\数字图像处理\images_chapter_04\Fig4.11(a).jpg'); subplot(121);imshow(I);title('原图(a)'); f_1=glpf('E:\数字图像处理\images_chapter_04\Fig4.11(a).jpg',5); subplot(122);imshow(f_1);title('截止频率为5进行GLPF的结果(b)'); f_2=glpf('E:\数字图像处理\images_chapter_04\Fig4.11(a).jpg',15); figure;subplot(121);imshow(f_2);title('截止频率为15进行GLPF的结果(c)'); f_3=glpf('E:\数字图像处理\images_chapter_04\Fig4.11(a).jpg',30); subplot(122);imshow(f_3);title('截止频率为30进行GLPF的结果(d)'); f_4=glpf('E:\数字图像处理\images_chapter_04\Fig4.11(a).jpg',80); figure;subplot(121);imshow(f_4);title('截止频率为80进行GLPF的结果(e)'); f_5=glpf('E:\数字图像处理\images_chapter_04\Fig4.11(a).jpg',230); subplot(122);imshow(f_5);title('截止频率为230进行GLPF的结果(f)');