Evaluate the spline at the points (theta[i], phi[i]), i=0,...,len(theta)-1
scipy.interpolate.SmoothSphereBivariateSpline.__call__
scipy.interpolate.SmoothSphereBivariateSpline.get_coeffs