The surface properties of the icy bodies in the saturnian system have been investigated by means of the CASSINI-VIMS (Visual Infrared Mapping Spectrometer) hyperspectral imager which operates in the 0.35-5.1 m wavelength range. In particular, we have analyzed 111 full disk hyperspectral images of Rhea ranging in solar phase between 0.08° and 109.8°. These data have been previously analyzed by Filacchione et al. (2007Filacchione et al. ( , 2010 to study, adopting various "spectral indicators" (such as spectral slopes, band depth, continuum level, etc.), the relations among various saturnian satellites.As a further step we proceed in this paper to a quantitative evaluation of the physical parameters determining the spectrophotometric properties of Rhea's surface. To do this we have applied Hapke (1993) IMSA model (Isotropic Multiple Scattering Approximation) which allow us to model the phase function at VIS-IR (visible-infrared) wavelengths as well as the spectra taking into account various types of mixtures of surface materials. Thanks to this method we have been able to constrain the size of water ice particles covering the surface, the amount of organic contaminants, the large scale surface roughness and the opposition effect surge. From our analysis it appears that wavelength dependent parameters, e.g. opposition surge width (h) and single-particle phase function parameters (b,v), are strongly correlated to the estimated single-scattering albedo of particles. For Rhea the best fit solution is obtained by assuming:1) an intraparticle mixture of crystalline water ice and a small amount (0.4%) of Triton tholin; 2) a monodisperse grain size distribution having a particle diameter a m = 38 µm; and 3) a surface roughness parameter value of 33°. The study of phase function shows that both Shadow Hiding and Coherent Backscattering contribute to the opposition surge. This study represents the first attempt, in the case of Rhea, to join the spectral and the photometric analysis. The surface model we derived gives a good quantitative description of both spectrum and phase curve of the satellite. The same approach and model, with appropriate modifications, shall be applied to VIMS data of the other icy satellites of Saturn, in order to reveal similarities and differences in the surface characteristics to understand how these bodies interact with their environment.