说明：  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