-def combine(image1, image2):
- im_l1 = image1.load()
- im_l2 = image2.load()
-
- on_both = []
-
- for x in xrange(image1.size[0]):
- for y in xrange(image1.size[1]):
- if im_l1[x, y] and im_l2[x, y]:
- on_both.append((x, y))
- return on_both
-
-def intersections_from_angl_dist(lines, size):
- intersections = []
- for (angl1, dist1) in sorted(lines[1], key=itemgetter(1)):
- line = []
- for (angl2, dist2) in sorted(lines[0], key=itemgetter(1)):
- if abs(angl1 - angl2) > 0.4:
- x = - ((dist2 / math.cos(angl2))-(dist1 / math.cos(angl1))) / (math.tan(angl1) - math.tan(angl2))
- y = (math.tan(angl1) * x) - (dist1 / math.cos(angl1))
- if (-size[0] / 2 < x < size[0] / 2 and
- -size[1] / 2 < y < size[1] / 2):
- line.append((int(x + size[0] / 2), int(y + size[1] / 2)))
- intersections.append(line)
- return intersections
-