Examples
initial_dir = PWD; chdir (SIPDIR + 'images'); // ======= Example 1 im = imread('star.bmp'); im = bwborder(im); clf imshow(im,2); h = hough(im); clf imshow(h,[]); // theta varies horizontally from 0 to 180 ht = 1*(h>= 40); // threshold the hough space lim = ihough(ht,size(im)); // draw the detected lines clf imshow(lim + 2*im +1, hotcolormap(4)) // detected lines shown in yellow // ======= Example 2: how to obtain the parameters // // creating a empty picture with a line at y = -90 e = zeros(200,200); e(10,:) = 1; // (remember that the Y axis points downwards and is centered in the // middle of the image) // getting its hough transform, and finding the points // corresponding to y=10 [h, rrange] = hough(e); [r,c] = find(h == max(h)) // Gets the parameters of the line theta = c - 1 // 90 degrees rho = rrange(r) // -90 rho (upper half of image) // thx to Herve Lombaert for inspiring example #2 ! chdir(initial_dir); |  |  |