%Accumulator
max_ro=floor((sqrt(row*row+col*col)/2));
min_ro=-max_ro;
rozmiar_aku=max_ro-min_ro+1;
max_alfa=180;
A=zeros(rozmiar_aku,max_alfa);

%Hough transform
for x=1:col
    for y=1:row
        if IG0(row-y+1,x,1)==0
            for ro=min_ro:max_ro
                for alfa=1:max_alfa
                    alf=pi*alfa/180.0;
                    if abs(ro)==round(x.*cos(alf)+y.*sin(alf))
                        A(floor(rozmiar_aku/2)-ro+1,alfa)=A(floor(rozmiar_aku/2)-ro+1,alfa)+1;
                    end;
                end;
            end;
        end; %if
    end;
end;


IA=double(zeros(rozmiar_aku,max_alfa,3));
IA(:,:,1)=A;
minIA=min(min(IA(:,:,1)));
maxIA=max(max(IA(:,:,1)));
IA(:,:,1)=255*(IA(:,:,1)-minIA)/(maxIA-minIA);
IA(:,:,2)=IA(:,:,1);
IA(:,:,3)=IA(:,:,1);