t4_5

function[] = circleHoughFunc(image)
    e = edge(image, 'canny');
    radii = 20:1:50;
    h = circle_hough(e, radii, 'same', 'normalise');
    peaks = circle_houghpeaks(h, radii, 'nhoodxy', 41, 'nhoodr', 51, 'npeaks', 14);
    imshow(image);
    hold on;
    for peak = peaks
        [x, y] = circlepoints(peak(3));
        plot(x+peak(1), y+peak(2), 'g-');
    end
    hold off
end