scipy.interpolate.LSQBivariateSpline.get_knots
scipy.interpolate.LSQBivariateSpline.integral
Return weighted sum of squared residuals of the spline approximation: sum ((w[i]*(z[i]-s(x[i],y[i])))**2,axis=0)