同态滤波增强算法

同态滤波增强算法是一种基于频域的图像增强方法，其原理是将图像在频域进行滤波，使得图像的低频分量得到增强，高频分量得到抑制。该方法的核心思想是在频域中将图像转换为可分离的两个部分：亮度和反射率。
亮度部分表示图像的整体亮度变化，而反射率则表示图像的细节信息。通过将亮度和反射率分别进行滤波处理，可以实现对图像的增强。

Image = imread('1.jpg');
% 对RGB三个通道分别进行同态滤波增强
for i = 1:3
    % 将当前通道的图像取对数
    logI = log(double(Image(:,:,i)) + 1);
    % 设定同态滤波参数
    sigma = 1.414;
    filtersize = [7 7];
    lowfilter = fspecial('gaussian', filtersize, sigma);
    highfilter = zeros(filtersize);
    highpara = 1;
    lowpara = 0.4;
    highfilter(ceil(filtersize(1,1)/2),ceil(filtersize(1,2)/2)) = 1;
    highfilter = highpara*highfilter - (highpara-lowpara)*lowfilter;
    % 进行同态滤波
    highpart = imfilter(logI, highfilter, 'replicate', 'conv');
    NewImage(:,:,i) = exp(highpart);
    % 对增强后的图像进行归一化处理
    top = max(NewImage(:));
    bottom = min(NewImage(:));
    NewImage(:,:,i) = (NewImage(:,:,i)-bottom)/(top-bottom);
    NewImage(:,:,i) = 1.5 * NewImage(:,:,i);
end
% 显示增强后的图像
imshow(NewImage);
title('彩色图像的同态滤波增强');
subplot(1,2,1),imshow(Image),title('原始图像')
subplot(1,2,2),imshow(NewImage),title('增强后的图像')