Raman spectroscopy provides information on the residual strain state of host-inclusion systems that, coupled with the elastic geobarometry theory, can be used to retrieve the P-T conditions of inclusion entrapment. In situ Raman measurements of zircon and coesite inclusions in garnet from the ultrahigh-pressure Dora Maira Massif show that rounded inclusions exhibit constant Raman shifts throughout their entire volume. In contrast, we demonstrate that Raman shifts can vary from the center to the edges and corners of faceted inclusions. Step-by-step polishing of the garnet host shows that the strain in both rounded and prismatic inclusions is gradually released as the inclusion approaches the free surface of the host. More importantly, our experimental results coupled with selected numerical simulations demonstrate that the magnitude and the rate of the strain release also depend on the contrast in elastic properties between the host and the inclusion and on the inclusion crystallographic orientation with respect to the external surface. These results allowed us to give new methodological guidelines for determining the residual strain in host inclusion systems.