<abstract><p>In this paper, a kind of bivariate multiquadric quasi-interpolant with the derivatives of a approximated function is studied by combining the known multiquadric quasi-interpolant with the generalized Taylor polynomials that act as the bivariate Lidstone interpolation polynomials. For practical purposes, a kind of improved approximation operator without any derivative of the approximated function is given by using bivariate divided differences to approximate the derivatives. It has the property of high-degree polynomial reproducing. In addition, the improved bivariate quasi-interpolation operators only demand information of the location points rather than the derivatives of the function approximated. Some error bounds in terms of the modulus of continuity of high order and Peano representations for the error are given. Several numerical comparisons with other existing methods are carried out to verify a higher degree of accuracy based on the obtained scheme. Furthermore, the advantage of our method is that the algorithm is very simple and easy to implement.</p></abstract>