In this paper, several modifications of the quasi-interpolation operator of Scott and Zhang [30] are discussed. The modified operators are defined for non-smooth functions and are suited for application on anisotropic meshes. The anisotropy of the elements is reflected in the local stability and approximation error estimates. As an application, an example is considered where anisotropic finite element meshes are appropriate, namely the Poisson problem in domains with edges.