for i=1:10
    rom = max(A(:));

    alfm = pi*(max_alpha-90)/180.0;
    for x=1:col
        for y=1:row
            value_temp = round(x.*cos(alfm)+y.*sin(alfm));

            if rom==value_temp
                lin_image(row-y+1,x,1)=255;
                lin_image(row-y+1,x,2)=0;
                lin_image(row-y+1,x,3)=0;
            end
        end
    end

    % remove global maximum
    r = 5;
    [rho,alpha] = find(A == max(A(:)));
    A(rho-r:rho+r, alpha-r:alpha+r) = 0;

end