hough变换检测多个圆 matlab

hough变换检测多个圆的matlab程序。就是对内存要求大点

function [Accumulator xdim ydim] = houghcircle2（Imbinaryr）
% Fast implementation of Circular HOUGH-Transform
% 
% parameters:
% 	Imbinary    ... a binary image （after applying an edge-detection algorithm）
% 	r           ... the radius of circle to detect within the image
% 	Accumulator ... result of the HOUGH-Transform - an image where all circles are accumulated 
% 	                the image is 2*r larger than the input image
% 	xdim ydim  ... dimensions of the Accumulator image
xdim = size（Imbinary 2） + ceil（2*r）;
ydim = size（Imbinary 1） + ceil（2*r）;
[circle_x circle_y] = circle（r r r）;
[edge_y edge_x] = find（Imbinary）;
clear Imbinary;
circlepoints = length（circle_x）;
edgepoints = length（edge_x）;
points = zeros（（circlepoints*edgepoints）1）;
pos = 1;
points = （ones（size（edge_x））*round（circle_x） + round（edge_x）*ones（size（circle_x））） * ydim + ...
		 （ones（size（edge_y））*round（circle_y） + round（edge_y）*ones（size（circle_y）））;
clear circle_x;
clear circle_y;
clear edge_x;
clear edge_y;
points = reshape（points 1 （size（points1） * size（points2）））;
Accumulator = histc（points （0:1:（max（points（:）））））;
clear points;
Accumulator（1） = 0; % delete counts from zeros
if （length（Accumulator） < （xdim*ydim））
	Accumulator = [Accumulator zeros（1 （xdim*ydim）-length（Accumulator））];
end;
Accumulator = reshape（Accumulator（1:（xdim*ydim）） ydim xdim）;
	

function [xy] = circle（center_x center_y radius）
% returns a vector with all circle-corner points containing the correct number of points
numpoints = ceil（2*radius*pi）;
theta = linspace（02*pi numpoints）;
rho = ones（1numpoints） * radius;
[xy] = pol2cart（theta rho）;
x = x+center_x;
y = y+center_y;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----
     文件        1711  2004-11-18 23:22  houghcircle2.m