椭圆拟合-EllipseDirectFit.m

椭圆拟合-EllipseDirectFit.m
这是一个快速和非迭代椭圆拟合算法
. 用法：
   A = EllipseDirectFit
输入: XY 数组是n个点的坐标
            x=XY, y=XY
输出: A = [a b c d e f]' 时椭圆拟合的系数向量其方程方程为：:
   ax^2 bxy cy^2 dx ey f = 0,其中A是被归一化的 ||A||=1
可以转换输出的几何参数，比如（半轴，中心等），具体的理论公式在http://www.mathworks.com/matlabcentral/fx_files/2

function A = EllipseDirectFit（XY）;
%
%  Direct ellipse fit proposed in article
%    A. W. Fitzgibbon M. Pilu R. B. Fisher
%     “Direct Least Squares Fitting of Ellipses“
%     IEEE Trans. PAMI Vol. 21 pages 476-480 （1999）
%
%  Our code is based on a numerically stable version
%  of this fit published by R. Halir and J. Flusser
%
%     Input:  XY（n2） is the array of coordinates of n points x（i）=XY（i1） y（i）=XY（i2）
%
%     Output: A = [a b c d e f]‘ is the vector of algebraic 
%             parameters of the fitting ellipse:
%             ax^2 + bxy + cy^2 +dx + ey + f = 0
%             the vector A is normed so that ||A||=1
%
%  This is a fast non-iterative ellipse fit.
%
%  It returns ellipses only even if points are
%  better approximated by a hyperbola.
%