数字图像处理 高频提升加强

一.实验题目

参照下图实现图像的高频增强。

高频提升增强：a>=0,b>=a    H_hfe(u,v)=a+bH_hp(u,v)     其中a=0.5,b=2.0

二.实验分析

本实验主要是频域滤波，其算法过程如下：
①用(-1)^x+y乘以输入图像进行中心变换
②计算①中的DFT F(u,v)  ③用滤波器函数H(u,v)乘以F(u,v)
④计算③中结果的反DFT  ⑤得到④中结果的实部
⑥用(-1)^x+y乘以⑤中的结果，取消输入图像的乘数
我们可以通过fft2,fftshift,ifftshift,iftt2来实现相应操作，然后通过与巴特沃斯

滤波函数相结合便可以完成实验。

三.实验程序

clear;clc;figure;
I=imread('images\Fig4.30(a).jpg');
subplot(2,2,1);imshow(I);title('a图 原图');
f=double(I);[r,c]=size(f);
F=fft2(f);G=fftshift(F);
d0=15; %半径范围
n=2;%巴特沃斯阶次
a=0.5;b=2.0; %高频强调滤波传递函数系数
mu=floor(r/2);mv=floor(c/2);
for u=1:r
    for v=1:c
        d=sqrt((u-mu)^2+(v-mv)^2);
        Hlpbtw=1/(1+0.414*(d/d0)^(2*n));
        Hhpbtw=1-Hlpbtw;    Ghpbtw(u,v)=Hhpbtw*G(u,v);
        Hhfebtw=a+b*Hhpbtw; Ghfebtw(u,v)=Hhfebtw*G(u,v);
    end
end
ghpbtw=ifftshift(Ghpbtw);
fhpbtw=uint8(real(ifft2(ghpbtw)));
subplot(2,2,2);imshow(fhpbtw);title('b图 巴特沃斯高通');
ghfebtw=ifftshift(Ghfebtw);
fhfebtw=uint8(real(ifft2(ghfebtw)));
subplot(2,2,3);imshow(fhfebtw);title('c图 巴特沃斯高频强调滤波');
histeq_fhfebtw=histeq(fhfebtw,256);
subplot(2,2,4);imshow(histeq_fhfebtw);title('c图均衡化结果');
        

程序说明:

1.btw  butterworth,巴特沃斯滤波器

2.f,g为空间域,G,H为频域