Based on a linearization approximation coupled with path integral formalism, we propose a method derived from the propagation of quasi-classical trajectories to simulate resonance Raman spectra. This method is based on a ground state sampling followed by an ensemble of trajectories on the mean surface between the ground and excited states. The method was tested on three models and compared to quantum mechanics solution based on a sum-overstates approach: harmonic and anharmonic oscillators and the HOCl molecule (hypochlorous acid). The method proposed is able to correctly characterize resonance Raman scattering and enhancement, including the description of overtones and combination bands. The absorption spectrum is obtained at the same time and the vibrational fine structure can be reproduced for long excited state relaxation times. The method can be applied also to dissociating excited states (as is the case for HOCl).