%%
he = imread('land.jpg');
imshow(he), title('H&E image');
cform = makecform('srgb2lab');
lab_he = applycform(he,cform);
ab = double(lab_he(:,:,2:3));
nrows = size(ab,1);
ncols = size(ab,2);
ab = reshape(ab,nrows*ncols,2);
nColors = 2;
[cluster_idx cluster_center] = kmeans(ab,nColors,'distance','sqEuclidean','Replicates',3);
pixel_labels = reshape(cluster_idx,nrows,ncols);
figure, imshow(pixel_labels,[]),title('image labeled by clusterindex');
segmented_images = cell(1,3);
rgb_label = repmat(pixel_labels,[1 1 3]);
for k = 1:nColors
color = he;
color(rgb_label ~= k) = 0;
segmented_images{k} = color;
end
imwrite(pixel_labels,'obr.jpg');
painted = imread('land_porovnanie.jpg');
segmentated = imread('obr.jpg');
A = (painted == segmentated);
m1 = size(painted,1);
m2 = size(painted,2);

pocetPixelov = m1*m2;
sucet = sum(A(:) == 1);
print(sucet/pocetPixelov);