Hough transformation
Use the typical (?, ?) parametrisation for the space of lines (? is the angle of the normal to the line, ? is the distance along this normal from the line to the origin of the image)
All points on a line will fulfill ? = x cos(?) + y sin(?)
Record the results in a 2D histogram