单尺度Retinex算法实现

function VV=SSR（S）
    if length（size（S））>2
        S=rgb2gray（S）;
    end
    [mn]=size（S）;
    S=double（S）;
    size1=m*n;
    mean=sum（sum（S）2）/size1;
    R=zeros（mn）;
    R=repmat（meanmn）;
    %N（ij）高斯函数
    c=500;
    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;
    G=conv2（SN1‘same‘）; %对与原始图像高斯卷积操作