We developed a fast algorithm to model a rough planetary surface by means of a statistical multi-facets approach. Consistently with Hapke theory (2012) we model the surface with a distribution of non-spatially resolved facets, being the distribution of their slopes completely described by roughness parameter