Отображение в RGB

close all;
clear all;
%function y = Retinex(inImage)
InputFileName = 'C:\\Images_tr\00003\00036_00029.jpg';
inImage = imread(InputFileName);
figure(20); imshow(uint8(inImage));
 
inImage = double(inImage);
 
[m,n] = size(inImage);
 
 
for j=1:n % create a log version of input image
    for i=1:m
        if inImage(i,j) == 0 % aproximate 0 with a small number to avoid NaN issues
            L(i,j) = 0.00001;
        else
        L(i,j)=log(inImage(i,j));  
        end
    end
end
 
imCopy = L;
Maximum = max(L(:));% maximum intensity value in the image
 
[sRow, sCol] = size(L);
shift = 2^(fix(log2(min(sRow, sCol)))-1); % initial place to begin
lastProduct = Maximum*ones(sRow, sCol); % initialize Last Product
while (abs(shift) >= 1) % iterate intil spiral is one unit large
    for i = 1:5  
        currRow = 0; % move horizontaly
        currCol = shift;
        interProduct = lastProduct;
        if (currRow + currCol > 0) % bottom right spiral part
            interProduct((currRow+1):end, (currCol+1):end) = lastProduct(1:(end-currRow), 1:(end-currCol)) + imCopy((currRow+1):end, (currCol+1):end) - imCopy(1:(end-currRow), 1:(end-currCol)); % propigate inward (adding is log mulitply)
        else
            interProduct(1:(end+currRow), 1:(end+currCol)) = lastProduct((1-currRow):end, (1-currCol):end) + imCopy(1:(end+currRow),1:(end+currCol)) - imCopy((1-currRow):end, (1-currCol):end); % propigate inward (adding is log mulitply)
        end
        interProduct(interProduct > Maximum) = Maximum; % % if values excede the max image value, rescale
        product = (interProduct + lastProduct)/2; % average with the previous Product
        lastProduct = product;
         
        currRow = shift; % move verticly
        currCol = 0;
        interProduct = lastProduct;
        if (currRow + currCol > 0) % bottom right spiral part
            interProduct((currRow+1):end, (currCol+1):end) = lastProduct(1:(end-currRow), 1:(end-currCol)) + imCopy((currRow+1):end, (currCol+1):end) - imCopy(1:(end-currRow), 1:(end-currCol)); % propigate inward (adding is log mulitply)
        else
            interProduct(1:(end+currRow), 1:(end+currCol)) = lastProduct((1-currRow):end, (1-currCol):end) + imCopy(1:(end+currRow),1:(end+currCol)) - imCopy((1-currRow):end, (1-currCol):end); % propigate inward (adding is log mulitply)
        end
         
        interProduct(interProduct > Maximum) = Maximum; % if values excede the max image value, rescale
 
        product = (interProduct + lastProduct)/2; % average with the previous Product
        lastProduct = product;
        
    end
    shift = -shift/2; % tighten spiral
end
y = product;
 
for j=1:n
    for i=1:m
        y(i,j)=exp(y(i,j));  % rescale output to be non logerithmic
    end
end
maxVal = max(max(y));
for j=1:n
    for i=1:m
        y(i,j)= y(i,j)/ maxVal;  
    end
end
 
y = uint8(round(y * 255));
cl=n/3;
ry=y(1:m,1:cl);gy=y(1:m,cl+1:2*cl);by=y(1:m,2*cl+1:n);
figure, imshow(uint8(round(inImage))); %This commented code shows the input
figure, imshow(uint8(y)); %and output images in the acgoritm