图像的频域滤波-高通滤波

图像的频域滤波-高通滤波

介绍

• 这里的高通滤波主要是巴特沃兹高通和理想高通，其频率特性正好和对应的低通相反，这里就不再赘述。

效果图

• 可见，理想高通仍有振铃效应的产生，而巴特沃兹就好很多。
• 加强滤波是为了弥补一下低频信息，让图像看上去更有意思一点哈。

示例代码

clc;
clear all;
close all;
J=imread('lena.jpg');
if size(J, 3)==3
    J = rgb2gray(J);
end
subplot(2,3,1);imshow(uint8(J));
xlabel('(a) 原图像');

J=double(J);
f=fft2(J);      %采用傅里叶变换
g=fftshift(f);   %数据矩阵平衡
[M,N]=size(f);
n1=floor(M/2);
n2=floor(N/2);
d0=20;
for i=1:M    %进行理想高通滤波和理想高通加强滤波
    for j=1:N
        d=sqrt((i-n1)^2+(j-n2)^2);
        if d>=d0
            h1=1;
            h2=1+0.5;
        else
            h1=0;
            h2=0.5;
        end
        g1(i,j)=h1*g(i,j);
        g2(i,j)=h2*g(i,j);
    end
end
g1=ifftshift(g1);
g1=uint8(real(ifft2(g1)));   
subplot(2,3,2);imshow(g1);  %显示理想高通滤波结果
xlabel('(b) 显示理想高通滤波');
g2=ifftshift(g2);
g2=uint8(real(ifft2(g2)));
subplot(2,3,3);imshow(g2); %显示理想高通加强滤波结果
xlabel('(c) 显示理想高通加强滤波');
n=2;
d0=20;
for i=1:M        %进行Bufferworth高通滤波和Bufferworth高通加强滤波
    for j=1:N
        d=sqrt((i-n1)^2+(j-n2)^2);
        if d==0
            h1=0;
            h2=.5;
        else
            h1=1/(1+(d0/d)^(2*n));
            h2=1/(1+(d0/d)^(2*n))+0.5;
        end
        gg1(i,j)=h1*g(i,j);
        gg2(i,j)=h2*g(i,j);
    end
end
gg1=ifftshift(gg1);
gg1=uint8(real(ifft2(gg1))); 
subplot(2,3,4);imshow(gg1);   %显示Bufferworth高通滤波
xlabel('(d) Bufferworth高通滤波');
gg2=ifftshift(gg2);
gg2=uint8(real(ifft2(gg2))); 
subplot(2,3,6);imshow(gg2);   %显示Bufferworth高通加强滤波
xlabel('(e) Bufferworth高通加强滤波');