Precise measurement of the sound source directivity not only requires special equipment, but also is time-consuming. Alternatively, one can reduce the number of measurement points and apply spatial interpolation to retrieve a high-resolution approximation of directivity function. This paper discusses the interpolation error for different algorithms with emphasis on the one based on spherical harmonics. The analysis is performed on raw directivity data for two loudspeaker systems. The directivity was measured using sampling schemes of different densities and point distributions (equiangular and equiareal). Then, the results were interpolated and compared with these obtained on the standard 5° regular grid. The application of the spherical harmonic approximation to sparse measurement data yields a mean error of less than 1 dB with the number of measurement points being reduced by 89%. The impact of the sparse grid type on the retrieval error is also discussed. The presented results facilitate optimal sampling grid choice for low-resolution directivity measurements.