同态滤波器matlab代码

版权声明：<–本博客所有内容均为自己在学习工作中的总结、摘录等-- --转载请注明出处-- --如有侵权请联系删除–> https://mp.csdn.net/mdeditor/89389530#
function [I6] =HomomorphicFilitering (rL,rH,c,image_input , D0)
clc
I = image_input; % tun.bmp
figure(1),subplot(121),imshow(I); title(‘原图’)
I=im2double(I); %转换数据类型为double型
[M,N]=size(I);
P = 2M; Q = 2N;
I2 = zeros(P,Q);
for i = 1:M
for j =1:N
I2(i,j) = I(i,j); %对图像进行填充
end
end
I2=log(I2+1); %取对数
FI=fft2(I2); %傅里叶变换
n1=floor(P/2);
n2=floor(Q/2);
for u=1:P
for v=1:Q
D(u,v)=sqrt(((u-n1).^2+(v-n2).2)); %频率域中点（u，v）与频率矩形中心的距离
H(u,v)=(rH-rL).(1-exp(-c.(D(u,v)^2./D02)))+rL; %高斯同态滤波
end
end
H=ifftshift(H); %对H做反中心化
I3=ifft2(H.*FI); %傅里叶逆变换
I4=real(I3);
I5 =I4(1:M, 1:N); %截取一部分
I6=exp(I5)-1; %取指数
figure(1),subplot(122),imshow(I6,[]);title(‘同态滤波增强后’);

结果图