clc;
clear;
[obraz,map]=imread('cancer','bmp');
obraz=double(obraz);
[row col deep]=size(obraz);
figure(1);
imshow(uint8(obraz));
title('Obraz oryginalny');

[obrazcell,map]=imread('cell','bmp');
obrazcell=double(obrazcell);
[row1 col1 deep1]=size(obrazcell);
figure(2);
imshow(uint8(obrazcell));
title('Obraz komorki');

gray=obraz;
obrazR=obraz(:,:,1);
obrazG=obraz(:,:,2);
obrazB=obraz(:,:,3);
gray(:,:,1)=(obrazR+obrazG+obrazB)/3;
gray(:,:,2)=gray(:,:,1);
gray(:,:,3)=gray(:,:,2);
figure(3);
imshow(uint8(gray));
title('Obraz gray');

bin=gray;
%binaryzacja
for i=1:row
    for j=1:col
        for k=1:deep
            if gray(i,j,k)>80
                bin(i,j,k)=255;
            else
                bin(i,j,k)=0;
            end
        end
    end
end

figure(4);
imshow(uint8(bin));
title('Obraz po binaryzacji');


graycell=obrazcell;
obrazRcell=obrazcell(:,:,1);
obrazGcell=obrazcell(:,:,2);
obrazBcell=obrazcell(:,:,3);
graycell(:,:,1)=(obrazRcell+obrazGcell+obrazBcell)/3;
graycell(:,:,2)=graycell(:,:,1);
graycell(:,:,3)=graycell(:,:,2);

bincell=graycell;
%binaryzacja
for i=1:row1
    for j=1:col1
        for k=1:deep1
            if graycell(i,j,k)>80
                bincell(i,j,k)=255;
            else
                bincell(i,j,k)=0;
            end
        end
    end
end

max_r=row-row1+1;
max_c=col-col1+1;
A=zeros(max_r,max_c);

for r=1:max_r
    for c=1:max_c
        p_sum=0;
        for p_r=1:row1
            for p_c=1:col1
                p_sum=p_sum+255-abs(gray(r+p_r-1,c+p_c-1)-graycell(p_r,p_c));
            end
        end
        A(r,c)=p_sum;
    end
end

%normalizacja
A_min = min(min(A(:,:,1)));
A_max = max(max(A(:,:,1)));
A(:,:,1)=255*(A(:,:,1)-A_min)/(A_max-A_min);


figure(5);
imshow(uint8(A));
title('Akumulator');