This paper describes the application of a first order regularization technique to the reconstruction of visible surfaces. Our approach is a computationally efficient first order method that simultaneously achieves approximate invariance and preservation of discontinuities. It is also robust with respect to the smoothing parameter /spl lambda/. The robustness property to /spl lambda/ allows a free choice of /spl lambda/ without struggling to determine an optimal /spl lambda/ that provides the best reconstruction. A new approximately invariant first order stabilizing function for surface reconstruction is obtained by employing a first order Taylor expansion of a nonconvex invariant stabilizing function that is expanded at the estimated value of the squared gradient instead of at zero. The data compatibility measure is the squared perpendicular distance between the reconstructed surface and the constraint surface. This combination of stabilizing function and data compatibility measure is necessary to achieve invariance with respect to rotations and translations. Sharp preservation of discontinuities is achieved by a weighted sum of adjacent pixels. The results indicate that the proposed methods perform well on sparse noisy range data. In addition, the volume between two surfaces normalized by the surface area (interpreted as average distance between two surfaces) is proposed as an invariant measure for the comparison of reconstruction results.< <ETX xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">></ETX>