巴特沃斯高通滤波器教案资料

巴特沃斯高通滤波器

数字图像处理实验报告（二）

频域滤波增强之高通滤波

一、实验目的

通过高通滤波实现图像的频域滤波，加深对频域图像增强的理解；使用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);

for i=1:M

for j=1:N

d=sqrt((i-m)^2+(j-n)^2);

if (d==0)

h=0;

else

h=1/(1+0.414*(d0/d)^(2*nn));% 计算传递函数 end

result(i,j)=h*g(i,j);

end

end

result=ifftshift(result);

J2=ifft2(result);

J3=uint8(real(J2));

figure(3),imshow(J3); % 滤波后图像显示

四、效果图

原图像

灰度图像

当d0=5时所得图像