Elastic geobarometry for host-inclusion systems can provide new constraints to assess the pressure and temperature conditions attained during metamorphism. Current experimental approaches and theory are developed only for crystals immersed in a hydrostatic stress field, whereas inclusions experience deviatoric stress. We have developed a method to determine the strains in quartz inclusions from Raman spectroscopy using the concept of the phonon-mode Grüneisen tensor. We used ab initio Hartree-Fock/Density Functional Theory to calculate the wavenumbers of the Raman-active modes as a function of different strain conditions. Least-squares fits of the phonon-wavenumber shifts against strains have been used to obtain the components of the mode Grüneisen tensor of quartz (γm1 and γm3) that can be used to calculate the strains in inclusions directly from the measured Raman shifts. The concept is demonstrated with the example of a natural quartz inclusion in eclogitic garnet from Mir kimberlite and has been validated against direct X-ray diffraction measurement of the strains in the same inclusion.