You can calculate the IR and Raman spectra using the Quantum Espresso package. If you are familiar with QE, it will be very simple for you. The procedure is as follows:
Optimize the wavefunction by performing a self-consistent field (scf) calculation with pw.x code.
Calculate the vibrational frequencies (normal modes/phonons) with ph.x code.
Extract the phonon information from ph.x output using dynmat.x code.
Parse the dynmat.x output section that contains the spectra data (frequencies and intensities) and plot it with gnuplot, to produce these two spectra.