testing scripts

[imago.git] / imago.py

diff --git a/imago.py b/imago.py
index ab44f4f..cf04f53 100755 (executable)
--- a/imago.py
+++ b/imago.py
@@ -100,7 +100,10 @@ def main():
     if verbose:
         print >> sys.stderr, "second hough transform"
 
-    hough2 = Hough(im_h2.size)
+    hough2 = Hough(im_h2.size) 
+    # im_hough might be used instead im_h2, but at the moment it brings a lot of
+    # noise to the second transform, which later confuses the center-finding
+    # mechanism (which is not very robust yet)
     im_hough2 = hough2.transform(im_h2)
     if show_all:
         do_something(im_hough2, "second hough transform")
@@ -109,26 +112,40 @@ def main():
     if show_all:
         do_something(im_h3, "third high pass filter")
      
-    im_h3 = filters.half_centers(im_h3)
+    im_h3 = filters.components(im_h3)
     if show_all:
         do_something(im_h3, "half centers")
 
     if verbose:
         print >> sys.stderr, "finding the grid"
 
-    lines_m = hough2.all_lines(im_h3)
+    lines_m = hough2.all_lines_h(im_h3)
     lines = []
+    im_c = im_h2.convert('RGB').convert('RGB', (1, 0.5, 0.5, 0))
+    draw_c = ImageDraw.Draw(im_c)
 
-    for line in lines_m:
+    for line_l in lines_m:
         im_line = Image.new('L', im_h2.size)
         draw = ImageDraw.Draw(im_line)
-        draw.line(line_from_angl_dist(line, im_h2.size), fill=255, width=5)
-        if show_all:
-            do_something(im_line, "line")
-        im_c = combine(im_h2, im_line)
-        if show_all:
-            do_something(im_c, "hough x lines")
-        lines.append(hough1.all_lines(im_c))
+        line_points = set()
+        for line in line_l:
+            draw.line(line_from_angl_dist(line, im_h2.size), fill=255, width=5)
+            draw_c.line(line_from_angl_dist(line, im_c.size), fill=(70, 70, 70), width=5)
+            for p in combine(im_h2, im_line):
+                line_points.add(p)
+        for point in line_points:
+            draw_c.point(point, fill=(120, 255, 120))
+        lines.append(hough1.lines_from_list(line_points))
+
+    if show_all:
+        do_something(im_c, "hough x lines")
+
+    image_g = image.copy()
+    draw = ImageDraw.Draw(image_g)
+    for line in [l for s in lines for l in s]:
+        draw.line(line_from_angl_dist(line, image.size), fill=(120, 255, 120))
+    if show_all:
+        do_something(image_g, "the grid")
 
     intersections = intersections_from_angl_dist(lines, image.size)
     image_g = image.copy()
@@ -142,7 +159,7 @@ def main():
                        line])
 
     if show_all:
-        do_something(image_g, "the grid")
+        do_something(image_g, "intersections")
 
     return 0
 
@@ -172,14 +189,13 @@ def combine(image1, image2):
     im_l1 = image1.load()
     im_l2 = image2.load()
 
-    im_n = Image.new('L', image1.size)
-    im_nl = im_n.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]:
-                im_nl[x, y] = 255
-    return im_n
+                on_both.append((x, y))
+    return on_both
 
 def line_from_angl_dist((angle, distance), size):
     x1 = - size[0] / 2
@@ -192,7 +208,7 @@ def intersections_from_angl_dist(lines, size):
     intersections = []
     for (angl1, dist1) in sorted(lines[1], key=itemgetter(1)):
         line = []
-        for (angl2, dist2) in lines[0]:
+        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))