单/多尺度Retinex（matlab）

    clc;clear;
    I=imread('DR_ck_0001.tiff');
     
    % Ir=I(:,:,1);%提取红色分量
    % Ig=I(:,:,2);%提取绿色分量
    % Ib=I(:,:,3);%提取蓝色分量
    % I1=SSR(Ir); %自定义单尺度Retinex处理函数
    % I2=SSR(Ig);
    % I3=SSR(Ib);
    % In=cat(3,I1,I2,I3);          %cat用于构造多维数组
    In = SSR(I);
    t = exp(-10);
    enhanced =  log(double(I)+t)-log(In+t);
    enhanced = normal(enhanced);
    subplot(1,3,1);imshow(I);
    title('A). 原始图像');  
    subplot(1,3,2);imshow(In,[]);
    title('B). 单尺度Retinex');
    subplot(1,3,3);imshow(enhanced,[]);
    title('C). 单尺度Retinex');
     
    %——————————