function [acumulator,outputimage] = irhough(Im,pattern,fnd)
[h,w] = size(Im);
[hp,wp] = size(pattern);
acumulator = zeros(h,w);
outputimage(:,:,1)=uint8(Im);
outputimage(:,:,2)=uint8(Im);
outputimage(:,:,3)=uint8(Im);

%czesc zasadnicza tutaj
for x = 1:w-wp+1
    for y = 1:h-hp+1
        p_sum=sum(abs(Im(y:y+hp-1,x:x+wp-1)-pattern),'all');
        acumulator(y,x)=p_sum;
    end
end

%to tylko rysowanie
acumulator=acumulator(1:h-hp+1,1:w-wp+1);
acumulatortmp = acumulator;
for i = 1:fnd
    [hmin,wmin] = find(acumulatortmp == min(acumulatortmp(:)));

    outputimage(hmin:hmin+hp-1,wmin,2)=0;
    outputimage(hmin:hmin+hp-1,wmin+wp-1,2)=0;
    outputimage(hmin,wmin:wmin+wp-1,2)=0;
    outputimage(hmin+hp-1,wmin:wmin+wp-1,2)=0;

    outputimage(hmin:hmin+hp-1,wmin,3)=0;
    outputimage(hmin:hmin+hp-1,wmin+wp-1,3)=0;
    outputimage(hmin,wmin:wmin+wp-1,3)=0;
    outputimage(hmin+hp-1,wmin:wmin+wp-1,3)=0;

    acumulatortmp(hmin-(hp/2):hmin+(hp/2),wmin-(wp/2):wmin+(wp/2))=max(acumulatortmp(:));
end
end