巴特沃斯高通滤波器

数字图像处理实验报告（二）频域滤波增强之高通滤波一、实验目的通过高通滤波实现图像的频域滤波，加深对频域图像增强的理解；使用MATLAB语言进行编程实现。二、实验原理1.计算需要增强的图像的傅里叶变换2.一个理想高通滤波器的转移函数满足下列条件H(u，v)=0当D(u，v)=d0时；H(u，v)=1当D(u，v)d0时3.巴特沃斯高通滤波器传递函数：H(u,v)=1/[1+(d0/D(u,v))^2n]4．利用巴特沃斯高通滤波器增强图像三、源程序I=imread('wll.jpg');figure(1),imshow(I)I=rgb2gray(I);figure(2),imshow(I);f=double(I);%数据类型转换，MATLAB不支持图像的无符号整型的计算g=fft2(f);%傅立叶变换g=fftshift(g);%转换数据矩阵[M,N]=size(g);nn=2;%二阶巴特沃斯(Butterworth)高通滤波器d0=5;m=fix(M/2);n=fix(N/2);fori=1:Mforj=1:Nd=sqrt((i-m)^2+(j-n)^2);if(d==0)h=0;elseh=1/(1+0.414*(d0/d)^(2*nn));%计算传递函数endresult(i,j)=h*g(i,j);endendresult=ifftshift(result);J2=ifft2(result);J3=uint8(real(J2));figure(3),imshow(J3);%滤波后图像显示四、效果图原图像灰度图像当d0=5时所得图像