conicfit (version 1.0.4)

CircleFitBySpath: Geometric circle fit by Spath

Description

CircleFitBySpath applies the Geometric circle fit by Spath

Usage

CircleFitBySpath(XY, ParIni = NA, epsilon = 1e-06, IterMAX = 50)

Arguments

XY

array of sample data

ParIni

initial guess (a, b, R)

epsilon

tolerance (small threshold)

IterMAX

maximal number of iterations, with a bad initial guess it may take >100 iterations

Value

vector(a, b, R)

vector with the values for the circle: center (a,b) and radius R

Source

Nikolai Chernov, 2014 Fitting ellipses, circles, and lines by least squares http://people.cas.uab.edu/~mosya/cl/ Nikolai Chernov, 2010 Circular and linear regression: Fitting circles and lines by least squares Chapman & Hall/CRC, Monographs on Statistics and Applied Probability, Volume 117

References

Nikolai Chernov, 2014 Fitting ellipses, circles, and lines by least squares http://people.cas.uab.edu/~mosya/cl/

Nikolai Chernov, 2010 Circular and linear regression: Fitting circles and lines by least squares Chapman & Hall/CRC, Monographs on Statistics and Applied Probability, Volume 117

Examples

xy<-calculateCircle(0,0,200,50,randomDist=TRUE,noiseFun=function(x) (x+rnorm(1,mean=0,sd=50)))
plot(xy[,1],xy[,2],xlim=c(-250,250),ylim=c(-250,250));par(new=TRUE)
c5 <- CircleFitBySpath(xy)
xyc5<-calculateCircle(c5[1],c5[2],c5[3])
plot(xyc5[,1],xyc5[,2],xlim=c(-250,250),ylim=c(-250,250),col='magenta',type='l');par(new=TRUE)