kruh.m

图像对齐（inage registration）工具箱

function circle=kruh(vel,radius);
% make circle of ones in the zero-matrix sized SxS with the radius RADIUS 
% circle=kruh(s,radius);
% B.Zitova, [e-mail]

pom=1:floor(vel/2);
if floor(vel/2)*2 == vel 
	pom1=[ -1*fliplr(pom) pom];
else 	
   pom1=[-1*fliplr(pom) 0 pom];
end		
	
pom2=ones(vel,1);
pom2=pom2*pom1;

pom2=pom2.*pom2;
pom1=pom2';
pom2=pom2+pom1;
pom2=pom2-radius*radius;
pom2=sign(pom2)-1;
circle=pom2&ones(vel,vel);