python利用scipy的optimize实现非线性最小二乘算法进行球心拟合

#-*- coding: utf-8 -*-
import numpy as np
import scipy.optimize as opt

def sphere_fit（P）:
    best = np.zeros（P.shape[0]）
    fitafitb=opt.curve_fit（fmaxP.Tbestp0=[1.50.41.50.4]）
    x0 = fita[0]
    y0 = fita[1]
    z0 = fita[2]
    Rsphere = fita[3]
    Rsphere = abs（Rsphere）
    o=[x0y0z0]
    return o Rs