这是巴特沃斯高通滤波代码 有点问题 显示后不像锐化了的图像 图像变匾了变长了 还出现了大量竖线 麻烦高手指点一下！~~谢谢！
另求理想高通滤波代码~~

I=imread('4.jpg');
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,imshow(J3); % 滤波后图像显示

高手指点一下啊～ 比较急～～ 另外我改了一下代码
I=imread('4.jpg');
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
g(i,j)=h*g(i,j);
end
end
g=ifftshift(g);
J2=ifft2(g);
J3=uint8(real(J2));

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


改后图片变成绿色 还有我更改截至频率 为什么图象看不出什么变化 只是绿色变淡了而已了那？？

[此贴子已经被作者于2007-6-5 13:56:19编辑过]