Numerical and experimental methods are presented for multiple scattering of light in discrete random media of densely-packed particles. The methods are utilized to interpret the observations of asteroid (4) Vesta and comet 67P/Churyumov-Gerasimenko.