halfcircle.m

图像经过二维傅立叶频谱分析后

function [xc,yc]=halfcircle(r,x0,y0)
theta=91:270;
theta=theta*pi/180;
[xc,yc]=pol2cart(theta,r);
xc=round(xc)'+x0;
yc=round(yc)'+y0;