points = [l.points for l in new_lines1]
for trial in xrange(3):
- line1, cons = ransac.estimate(points, 2, 800, Diagonal_model)
- points2 = map(lambda l: [(p if not p in cons else None) for p in l], points)
- line2, cons2 = ransac.estimate(points2, 2, 800, Diagonal_model)
+ model = Diagonal_model(points)
+ diag_lines = ransac.ransac_multi(2, points, 2, 800, model=model)
+ [(line1, cons), (line2, cons2)] = diag_lines
center = intersection(line1, line2)
data = sum(points, [])
diag1 = Line(line1)