The algebraic distance of a point to a circle
is simply
.
Thus, we have a simple procedure to estimate the parameters:
find the parameters
that minimizes
Define the design matrix:
The least square solution is the eigenvector of the smallest eigenvalue of
the matrix .
Some matlab code for the fit is here
and for conversion to standard form is here.