Displaying 1 result from an estimated 1 matches for "cprat".
Did you mean:
carat
2012 Jul 02
0
Fit circle with R
...Now add some random noise
mult = 0.5;
xs = xs+mult*rnorm(rnorm(xs));
ys = ys+mult*rnorm(rnorm(ys));
plot(xs,ys,pch=19,cex=0.5,col="red",xlim=c(-10,20),ylim=c(-10,20),asp=1)
# real circle
draw.circle(x_c,y_c,radius=10,border="black")
points(x_c,y_c,,pch=4,col="black")
CPrat <- CircleFitByPratt(xs,ys)
draw.circle(CPrat[1],CPrat[2],radius=CPrat[3],border="blue")
points(CPrat[1],CPrat[2],pch=4,col="blue")
MyC <- fitCircle(xs,ys)
draw.circle(MyC[1],MyC[2],radius=MyC[3],border="green")
points(MyC[1],MyC[2],pch=4,col="green")...