用matlab实现一种MSRCR图像增强算法，

MSRCR（Multi-Scale Retinex with Color Restoration）是一种图像增强算法，它在不同的尺度上对图像进行Retinex分解，并且通过颜色恢复来得到最终的增强图像。

以下是用MATLAB实现MSRCR算法的步骤：

1. 加载图像并将其转换为YCbCr颜色空间。

img = imread('image.jpg');
img_ycbcr = rgb2ycbcr(img);

2. 对亮度分量（Y通道）进行多尺度Retinex分解。

Y = double(img_ycbcr(:,:,1))/255;
sigma_list = [15,80,250];
G = fspecial('gaussian',[15 15],7);
Y_log = log(Y+1);
for i = 1:length(sigma_list)
    sigma = sigma_list(i);
    Gi = imfilter(Y,G*sigma,'replicate');
    Li = log(Y+1) - log(Gi+1);
    M(:,:,i) = imfilter(Li,G*sigma,'replicate');
end

3. 计算MSRCR权重。

epsilon = 0.01;
W = zeros(size(Y));
for i = 1:length(sigma_list)
    Wi = exp(-abs(M(:,:,i)-log(Y+1))/((sigma_list(i)^2)/epsilon));
    W = W + Wi;
end

4. 对MSRCR权重进行归一化处理。

W = W/length(sigma_list);

5. 进行颜色恢复。

Cb = imresize(img_ycbcr(:,:,2),size(Y));
Cr = imresize(img_ycbcr(:,:,3),size(Y));
R = Y.^(1-W).*Cb.^W.*Cr.^W;

6. 将增强后的图像转换回RGB颜色空间。

R = R.*255;
img_msrcr = cat(3,R,Cb,Cr);
img_msrcr = ycbcr2rgb(uint8(img_msrcr));

完整代码：

img = imread('image.jpg');
img_ycbcr = rgb2ycbcr(img);

Y = double(img_ycbcr(:,:,1))/255;
sigma_list = [15,80,250];
G = fspecial('gaussian',[15 15],7);
Y_log = log(Y+1);
for i = 1:length(sigma_list)
    sigma = sigma_list(i);
    Gi = imfilter(Y,G*sigma,'replicate');
    Li = log(Y+1) - log(Gi+1);
    M(:,:,i) = imfilter(Li,G*sigma,'replicate');
end

epsilon = 0.01;
W = zeros(size(Y));
for i = 1:length(sigma_list)
    Wi = exp(-abs(M(:,:,i)-log(Y+1))/((sigma_list(i)^2)/epsilon));
    W = W + Wi;
end

W = W/length(sigma_list);

Cb = imresize(img_ycbcr(:,:,2),size(Y));
Cr = imresize(img_ycbcr(:,:,3),size(Y));
R = Y.^(1-W).*Cb.^W.*Cr.^W;

R = R.*255;
img_msrcr = cat(3,R,Cb,Cr);
img_msrcr = ycbcr2rgb(uint8(img_msrcr));

imshow(img_msrcr);

注意：在实现过程中，可能需要根据图像的特点调整参数，如sigma_list、epsilon等。