awimg=rgb2gray(im1)

tic

[accum, circen, cirrad] = CircularHough_Grd(rawimg, [20 30],5,50)

circen

toc

figure(1) imagesc(accum) axis image

title( Accumulation Array from Circular Hough Transform )

figure(2) imagesc(rawimg) colormap( gray ) axis image

hold on

plot(circen(:,1), circen(:,2), r+ )

for k = 1 : size(circen, 1),

DrawCircle(circen(k,1), circen(k,2), cirrad(k), 32, b- )

end

hold off

title([ Raw Image with Circles Detected , ...

(center positions and radii marked) ])

figure(3) surf(accum, EdgeColor , none ) a