We introduce a method for estimation of rotation invariant local shape descriptors for 3D models. This method follows a successful idea commonly used to obtain rotation invariant descriptors in 2D images, and improves it by tackling the difficulty of the 3 degrees of freedom that exists in 3D models. Our method is simple, yet it achieves high levels of invariance after rotation transformations, and it produces short descriptors that can be efficiently used in several tasks. Such is the case of automatic classification of 3D surfaces with archaeological value, in which the proposed method attains state-of-the-art results in shorter times when compared with previous methods.