资源描述:
《数字图像处理高通滤波器.docx》由会员上传分享,免费在线阅读,更多相关内容在行业资料-天天文库。
1、数字图像处理高通滤波器姓名:***学号:**********高通滤波是常见的频域增强的方法之一。高通滤波与低通滤波相反,它是高频分量顺利通过,使低频分量受到削弱。这里考虑三种高通滤波器:理想高通滤波器、巴特沃斯高通滤波器和高斯高通滤波器。这三种滤波器涵盖了从非常尖锐(理想)到非常平坦(高斯)范围的滤渡器函数,其转移函数分别为:1、理想高通滤波器(IHPF)2、巴特沃斯高通滤波器(BHPF)3、高斯高通滤波器(GHPF)一、数字图像高通滤波器的实验过程:1、理想高通滤波器程序clearall;clc;image=imread
2、('test.jpg');gimage_15=func_ihpf(image,15);gimage_30=func_ihpf(image,30);gimage_80=func_ihpf(image,80);figuresubplot(221),imshow(image);title('Original');subplot(222),imshow(gimage_15);title('d0=15');subplot(223),imshow(gimage_30);title('d0=30');subplot(224),imsho
3、w(gimage_80);title('d0=80');%被调函数子函数G(u,v)=F(u,v)H(u,v)functiongimage=func_ihpf(image,d0)image=double(image);f=fftshift(fft2(image));[M,N]=size(f);a0=fix(M/2);b0=fix(N/2);fori=1:Mforj=1:ND=sqrt((i-a0)^2+(j-b0)^2);if(D>d0)h=1;elseh=0;endg(i,j)=h*f(i,j);endendgimage
4、=uint8(real(ifft2(ifftshift(g))));2、巴特沃斯高通滤波器程序clearall;clc;image=imread('test.jpg');gimage_15=func_bhpf(image,15);gimage_30=func_bhpf(image,30);gimage_80=func_bhpf(image,80);figuresubplot(221),imshow(image);title('Original');subplot(222),imshow(gimage_15);title('
5、d0=15');subplot(223),imshow(gimage_30);title('d0=30');subplot(224),imshow(gimage_80);title('d0=80');%被调函数子函数G(u,v)=F(u,v)H(u,v)functiongimage=func_bhpf(image,d0)image=double(image);f=fftshift(fft2(image));[M,N]=size(f);nn=1;%巴特沃斯阶数a0=fix(M/2);b0=fix(N/2);fori=1:Mf
6、orj=1:ND=sqrt((i-a0)^2+(j-b0)^2);if(D==0)h=0;elseh=1/(1+0.414*(d0/D)^(2*nn));endg(i,j)=h*f(i,j);endendgimage=uint8(real(ifft2(ifftshift(g))));3、高斯高通滤波器程序clearall;clc;image=imread('test.jpg');gimage_15=func_ghpf(image,15);gimage_30=func_ghpf(image,30);gimage_80=fun
7、c_ghpf(image,80);figuresubplot(221),imshow(image);title('Original');subplot(222),imshow(gimage_15);title('d0=15');subplot(223),imshow(gimage_30);title('d0=30');subplot(224),imshow(gimage_80);title('d0=80');%被调函数子函数G(u,v)=F(u,v)H(u,v)functiongimage=func_ghpf(image,
8、d0)image=double(image);f=fftshift(fft2(image));[M,N]=size(f);a0=fix(M/2);b0=fix(N/2);fori=1:Mforj=1:ND=sqrt((i-a0)^2+(j-b0)^2);h=1-exp(-(D.^2)./(2*(d0^2