The reconstruction of surfaces from speckle interferometry data is a demanding data-analysis task that involves edge detection, edge completion, and image reconstruction from noisy data. We present an approach that makes optimal use of the experimental information to minimize the hampering influence of the noise. The experimental data are then analyzed with a combination of wavelet transform and Bayesian probability theory. Nontrivial examples are presented to illustrate the proposed technique.