怎么在matlab上对灰度图像进行Retinex处理

在 MATLAB 中，可以使用以下步骤对灰度图像进行 Retinex 处理：

1. 读取灰度图像：

img = imread('image.jpg');
gray_img = rgb2gray(img);

2. 对灰度图像进行 Retinex 处理：

% 参数设置
sigma = 50;
k = 0.1;
epsilon = 0.0001;

% Retinex 处理
retinex_img = zeros(size(gray_img));
for i=1:size(gray_img,1)
    for j=1:size(gray_img,2)
        % 高斯模糊
        blur_img = imgaussfilt(double(gray_img(i,j)), sigma);
        
        % 对数变换
        log_img = log(k + gray_img(i,j));
        
        % Retinex 处理
        retinex_img(i,j) = exp(log_img - log(blur_img + epsilon));
    end
end

3. 显示 Retinex 处理后的图像：

imshow(retinex_img, []);

其中，`sigma` 是高斯模糊的参数，可以根据实际情况进行调整；`k` 和 `epsilon` 是对数变换的参数，`k` 用于控制对数变换的幅度，`epsilon` 用于防止出现除 0 错误。