多尺度Retinex

function VV=MSR（S）
    if length（size（S））>2
        S=rgb2gray（S）;
    end
    [mn]=size（S）;
    S=double（S）;

    sizes=m*n;
    mean=sum（sum（S）2）/sizes;
    R=zeros（mn）;
    R=repmat（meanmn）;

    %N（ij）高斯函数
    c=15;%尺度一
    for i=1:m
        for j=1:n
            N（ij）=exp（-（i^2+j^2）/c^2）;
        end
    end
    K=1/sum（N（:））;
    N1=K*N;
    G1=conv2（SN1‘same‘）;  

    c=80;%尺度二
    for i=1:m
        for j=1:n
            N（ij）=exp（-（i^2+j^2）/c^2）;
        end
    end
    K=1/sum（N（:））;
    N1=K*N;
    G2=conv2（SN1‘same‘）;  

    c=200;%尺度三
    for i=1:m
        for j=1:n
            N（ij）=exp（-（i^2+j^2）/c^2）;
        end
    end
    K=1/sum（N（:））;
    N1=K*N;
    G3=conv2（SN1‘sam