Cubature formulae for two dimensional integrals based on multivariate quasi-interpolating splines