同态滤波去除光照不均匀的matlab代码

同态滤波是一种用于去除图像中光照不均匀的有效方法。在Matlab中，可以通过以下代码来实现同态滤波去除光照不均匀：

% 读取原始图像
img = imread('input.jpg');

% 将原始图像转换为灰度图像
gray_img = rgb2gray(img);

% 对灰度图像进行同态滤波处理
filtered_img = homomorphic_filter(gray_img);

% 显示处理后的图像
imshow(filtered_img);
title('同态滤波去除光照不均匀');

% 定义同态滤波函数
function filtered_img = homomorphic_filter(img)
    % 将灰度图像转换为双精度数据类型
    img_double = im2double(img);
    
    % 对图像进行对数变换
    log_img = log(1 + img_double);
    
    % 设定同态滤波参数
    gamma_l = 0.3;
    gamma_h = 2.0;
    c = 1;
    D0 = 10;
    
    % 设定图像频域变换的尺寸
    [M, N] = size(img);
    
    % 生成频域网格
    [U, V] = meshgrid(1:N, 1:M);
    
    % 计算频率域滤波器H
    H = (gamma_h - gamma_l) * (1 - exp(-c * (D0^2 ./ (U.^2 + V.^2)))) + gamma_l;
    
    % 对频域图像进行滤波
    filtered_log_img = H .* log_img;
    
    % 将滤波后的对数图像进行反变换
    filtered_img = exp(filtered_log_img) - 1;
    
    % 将处理后的图像转换为uint8数据类型
    filtered_img = uint8(filtered_img * 255);
end

以上代码中，首先读取原始图像，并将其转换为灰度图像。然后定义了一个名为`homomorphic_filter`的函数，利用同态滤波算法对灰度图像进行处理。在此算法中，将图像进行对数变换后，计算频域滤波器H并对频域图像进行滤波，最后对滤波后的对数图像进行反变换，得到处理后的图像。最后通过imshow函数显示处理后的图像，以可视化处理效果。