高频强调滤波增强图像

文章目录

• 前言
• 一、基础思路
• 二、例题
• 三、代码
• 四、结果展示
• 总结

 

---

前言

学习完高频滤波器的设计之后，紧接着就会接触到高频强调滤波的内容，这一块内容相对简单，只需要在高频滤波器的基础上稍加补充即可。

---

一、基础思路

高频强调就是在高频滤波的基础上对传递函数稍加修改，引入k1、k2两个变量。k1≥0给出了控制距原点的偏移量，k2≥0控制高频的贡献。

二、例题

给出一张胸部X射线图像，试用高频强调滤波增强该图像

三、代码

%由于在之前的文章中已经详细介绍了高斯滤波器的设计思路，在这里就不过多介绍，只介绍高频强调滤波增强部分
%本例体重采用2阶滤波器，D0=40，并令k1=0.5,k2=0.75
I1=imread('chest.tif');
[len,wed]=size(I1);
g=fft2(I1);
g=fftshift(g);
[M,N]=size(g);
m=fix(M/2);
n=fix(N/2);
D0=40;
n1=2
for i=1:M
    for j=1:N
        D=sqrt((i-m)^2+(j-n)^2);
        H1=1-exp((-1)*(D^2/(2*D0^2)));
        H2=0.5+0.75*H1;%根据给定的k1\k2两个系数对传递函数进行修改
        s1(i,j)=H1*g(i,j);
        s2(i,j)=H2*g(i,j);%高频强调滤波是利用修改后的函数进行频域相乘
    end
end
I2=im2uint8(real(ifft2(ifftshift(s1)/255)));
I3=im2uint8(real(ifft2(ifftshift(s2)/255)));
I4=histeq(I3);%高频强调滤波增强和直方图均衡化是常见的搭配处理方式，可以得到更好的图像效果
subplot(221),imshow(I1),title('原图');
subplot(222),imshow(I2),title('高斯高频滤波图');
subplot(223),imshow(I3),title('高斯高频强调滤波图');
subplot(224),imshow(I4),title('直方图均衡化');

四、结果展示

---

总结

可以看到原本普通的高频滤波之后图像的效果在引入k1、k2两个系数之后得到了明显的增强。

高频强调滤波和直方图均衡化联合使用得到的图像有奇效！