close all;
clearvars;
clc;

image=imread('ccl1.png');
figure(1);
subplot(1,2,1);
imshow(image);
[XX,YY]=size(image);
L=1;
N=100;
id = zeros(N);
lut=zeros(N);
for i = 1:N
    id(i)=i;
end

for x = 2:XX-1
    for y = 2:YY-1
        X=image(x,y);
        if X~=0
            sasiedzi=[image(x-1,y-1) image(x-1,y) image(x-1, y+1) image(x,y-1)];
            suma = sum(sasiedzi);
            if suma==0
                image(x,y)=L;
                L=L+1;
            elseif suma>0
                sasiedziBezZer=nonzeros(sasiedzi);
                minimum=min(sasiedziBezZer);
                maksimum=max(sasiedziBezZer);
                id = union(minimum, maksimum, id);
                image(x,y)=minimum;
            end
        end
    end
end

for i = 1:N
    lut(i) = root(i,id);
end

for x = 2:XX-1
    for y = 2:YY-1
        X=image(x,y);
        if X~=0
            image(x,y) = root(image(x,y),id);
        end
    end
end

subplot(1,2,2);
% imshow(image);
imshow(label2rgb(image));