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