我有这样的形象:
代表同心椭圆。这个图像包含在一个数字矩阵中,在不同椭圆周长的索引中有一个,在其他地方为零。
因为我想找出每一个椭圆的中心,我找到了这两个函数:
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 =