scipy.interpolate.LSQSphereBivariateSpline
scipy.interpolate.LSQSphereBivariateSpline.ev
Evaluate the spline at the grid ponts defined by the coordinate arrays theta, phi.