python椭圆拟合_将椭圆拟合到python中的点

实际上有一个帖子与此非常相似，但提出的问题与5年前不同。在

我有一组不是椭圆的点，我想用最小二乘法拟合椭圆。这些是我发现的用于计算的函数。我仔细检查了一下，但没有发现任何错误。在import numpy as np

import numpy.linalg as linalg

import matplotlib.pyplot as plt

def fitEllipse(x,y):

x = x[:,np.newaxis]

y = y[:,np.newaxis]

D = np.hstack((x*x, x*y, y*y, x, y, np.ones_like(x)))

S = np.dot(D.T,D)

C = np.zeros([6,6])

C[0,2] = C[2,0] = 2; C[1,1] = -1

E, V = linalg.eig(np.dot(linalg.inv(S), C))

#print E

n = np.argmax(np.abs(E))

a = V[:,n]

return a

def ellipse_center(a):

b,c,d,f,g,a = a[1]/2, a[2], a[3]/2, a[4]/2, a[5], a[0]

num = b*b-a*c

x0=(c*d-b*f)/num

y0=(a*f-b*d)/num

return np.array([x0,y0])

def ellipse_angle_of_rotation( a ):

b,c,d,f,g,a = a[1]/2, a[2], a[3]