> I have copied and assimilated the houghpeaks function into octave.
> %HOUGHPEAKS Detect peaks in Hough transform.
> % with peak neighborhood suppressed.
> % 0.5*max(H(:)). If NUMPEAKS is omitted, it defaults to 1.
> % Make sure the neighborhood size is odd.
> % Suppress this maximum and its close neighbors.
> % the rho direction.
> % along the rho axis for theta = +/- 90 degrees.
> % Convert to linear indices to zero out all the values.