Local hybrid approximation for scattered data fitting with bivariate splines