We present a full-wave inversion approach with total variation regularization for elastography. The proposed method is based on the minimization of an error in constitutive equations functional augmented with a least squares data misfit term referred to as MECE for “modified error in constitutive equations.” The main theme of this paper is to demonstrate several key strengths of the proposed method on experimental data. In addition, some illustrative examples are provided where the proposed method is compared with a common shear wave elastography (SWE) approach. To this end, ultrasonically tracked displacement data from an acoustic radiation force (ARF) pulse are used in different phantoms including phantom with layered inclusion and triangle inclusion. The results indicate that the MECE approach can produce accurate shear modulus reconstructions in comparison with SWE, especially around the sharp edges in the layered and triangle inclusions. We compare shear modulus reconstruction using MECE and SWE with original inclusion shapes using two-dimensional normalized zero mean cross correlation, edge preservation index and dice coefficient similarity index. [Work supported by NIH Grant R01 CA174723.]