Accurately estimating all strain components in quasi-static ultrasound elastography is crucial for the full analysis of biological media. In this paper, 2D strain tensor imaging is investigated, using a partial differential equation (PDE)-based regularization method. More specifically, this method employs the tissue property of incompressibility to smooth the displacement fields and reduce the noise in the strain components. The performance of the method is assessed with phantoms and in vivo breast tissues. For all the media examined, the results showed a significant improvement in both lateral displacement and strain but also, to a lesser extent, in the shear strain. Moreover, axial displacement and strain were only slightly modified by the regularization, as expected. Finally, the easier detectability of the inclusion/lesion in the final lateral strain images is associated with higher elastographic contrast-to-noise ratios (CNRs), with values in the range [0.68 - 9.40] vs [0.09 - 0.38] before regularization.