From the Bayesian regularization theory we derive a quadratic cost-function for integrating wave-front gradient fields. In the proposed cost-function, the term of conditional distribution uses a central-differences model to make the estimated function well consistent with the observed gradient field. As will be shown, the results obtained with the central-differences model are superior to the results obtained with the backward-differences model, commonly used in other integration techniques. As a regularization term we use an isotropic first-order differences Markov Random-Field model, which acts as a low-pass filter reducing the errors caused by the noise. We present simulated and real experiments of the proposal applied in the Foucault test, obtaining good results.