Abstract. The multiple signal classification (MUSIC) imaging method is applied to determine the locations of a collection of small anisotropic spherical scatterers in the framework of the total internal reflection tomography. Multiple scattering between scatterers is considered and the inverse scattering problem is nonlinear, which, however, is solved by the proposed fast analytical approach where no associated forward problem is iteratively evaluated. The paper also discusses the role of the polarization of incidence waves, the incidence angle, the separation of scatterers from the surface of the substrate, and the level of noise on the resolution of imaging.