Optimal 3-D field interpolation for spheroidal sources