close all
clear all
I = imread('rice.png');

BW = I;
[X,Y]=size(I);
W2 = 15;
for i = 1:X
    for j = 1:Y
        if I(i,j) < meanLT(i,j,W2,I,X,Y)
            BW(i,j) = 0;
        else
            BW(i,j) = 255;
        end
    end
end
figure;
subplot(2,1,1);
imshow(I);
title('oryginalny')
subplot(2,1,2);
imshow(BW);
title('zbinaryzowany')
%%W przypadku obiektów kontrastowych które są blisko siebie może dochodzić
%%do połączeń.

BW = I;
BW2 = I;
[X,Y]=size(I);
W2 = 15;
for i = 1:X
    for j = 1:Y
        R = 128;
        k = 0.35;
        m = meanLT(i,j,W2,I,X,Y);
        if I(i,j) < m *( 1 - k * (stddevLT(i,j,W2,I,m,X,Y) /R -1 ))
            BW(i,j) = 0;
        else
            BW(i,j) = 255;
        end

        if I(i,j) < m *( 1 + k * (stddevLT(i,j,W2,I,m,X,Y) /R -1 ))
            BW2(i,j) = 0;
        else
            BW2(i,j) = 255;
        end
    end
end
figure;
subplot(2,2,1);
imshow(I);
title('oryginalny')
subplot(2,2,2);
imshow(BW);
title('zbinaryzowany Sauvoli +')
subplot(2,2,3);
imshow(BW2);
title('zbinaryzowany Sauvoli -')
%% - Zwieksza tolerancje progu, a + zmniejsza uzwgledniane są tylko sąsiednie