单尺度retinex增强的MATLAB实现代码

以下是单尺度Retinex增强的MATLAB实现代码：

function enhanced_img = single_scale_retinex(img, sigma)
% img: 待增强的图像
% sigma: 高斯核的标准差

% 将图像转换为双精度型
img = im2double(img);

% 计算高斯核
kernel = fspecial('gaussian', max(1, fix(6 * sigma)), sigma);

% 对原图像进行高斯滤波
log_img = log(img);
filtered_img = imfilter(log_img, kernel, 'replicate');

% 计算中心环绕函数
center_weight = center_weight_func(size(img));

% 对图像进行增强
enhanced_img = center_weight .* (log_img - filtered_img) + filtered_img;

% 将图像转换为灰度图像
if size(img, 3) == 3
    enhanced_img = rgb2gray(enhanced_img);
end

% 对图像进行灰度值缩放
enhanced_img = enhanced_img * 255;
enhanced_img = uint8(enhanced_img);

end

function weight = center_weight_func(size)
% size: 图像的尺寸

% 计算中心点的位置
center_x = fix(size(1) / 2) + 1;
center_y = fix(size(2) / 2) + 1;

% 计算中心环绕函数
[X, Y] = meshgrid(1:size(2), 1:size(1));
distance = sqrt((X - center_y).^2 + (Y - center_x).^2);
weight = 1 - exp(-(distance.^2) / (2 * (0.1 * size(1))^2));

end

使用方法：

img = imread('example.jpg');
enhanced_img = single_scale_retinex(img, 50); % 使用标准差为50的高斯核进行增强
imshow(enhanced_img);

以上代码只是单尺度Retinex的一种实现方式，具体的实现方式可能因人而异，需要根据实际情况进行调整和优化。