rtc_util.py 文件源码

def mk_tau(jmax,imax,f_list,eref,cref):
    taux=np.zeros(imax*jmax).reshape(jmax,imax)
    x=t_set 
    #t=cv2.getTickCount() 
    for j in range(jmax):
        for i in range(imax):
            fref=f_list[j][i]
            res=np.polyfit(x,fref(x,eref[j,i]),3) # 3rd order
        taux[j,i]=iestimate(res,cref[j,i])
        #if j % 100 == 0: print i,(cv2.getTickCount()-t)/cv2.getTickFrequency()
    temp=np.where(np.isnan(taux)==True)
    print np.nanmean(taux),np.nanstd(taux),len(temp[0])
    return taux

#5/25/2016