function HC = hough_circle_create(xcNb,xcMin,xcMax,ycNb,ycMin,ycMax,rNb,rMin,rMax)

HC.xcNb=xcNb;
HC.xcMin=xcMin;
HC.xcMax=xcMax;
HC.xcDelta=(xcMax-xcMin)/(xcNb-1);
HC.xc=linspace(xcMin,xcMax,xcNb);
HC.ycNb=ycNb;
HC.ycMin=ycMin;
HC.ycMax=ycMax;
HC.ycDelta=(ycMax-ycMin)/(ycNb-1);
HC.yc=linspace(ycMin,ycMax,ycNb);
HC.rNb=rNb;
HC.rMin=rMin;
HC.rMax=rMax;
HC.rDelta=(rMax-rMin)/rNb;
HC.r=linspace(rMin,rMax,rNb);
HC.stack=zeros(ycNb,xcNb,rNb);
for i=1:rNb
    nb=ceil(2*pi*HC.r(i)/sqrt(HC.xcDelta^2+HC.ycDelta^2)*2);
    [xc,yc]=pts_circle(HC.r(i),nb);
    HC.xcCircle{i}=xc;
    HC.ycCircle{i}=yc;
end
end


function [xc,yc]=pts_circle(r,nb)

theta=(0:nb-1)/nb*2*pi;
xc=r*cos(theta);
yc=r*sin(theta);
end




