A computational procedure for numerical integration of 2-D integrals based on local bivariate C^1 quasi-interpolating splines