#-*-?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