Gauss map computation for free-form surfaces

The Gauss map of a smooth doubly-curved surface characterizes the range of variation of the surface normal as an area on the unit sphere. An algorithm to approximate the Gauss map boundary to any desired accuracy is presented , in the context of a tensor-product polynomial surface patch, r(u, v) for (u, v) is an element of [0, 1] x [0, 1]. Boundary segments of the Gauss map correspond to variations of the normal along the patch boundary or the par abolic lines (loci of vanishing Gaussian curvature) on the surface, To comp ute the latter, points of vanishing Gaussian curvature are identified with the zero-set of a bivariate polynomial, expressed in the numerically-stable Bernstein basis-the subdivision and variation-diminishing properties then govern an adaptive quadtree decomposition of the (it, v) parameter domain t hat captures the zero-set of this polynomial to any desired accuracy. Loci on the unit sphere corresponding to the patch boundaries and parabolic line s are trimmed at their mutual or self-intersection points (if any), and the resulting segments are arranged in a graph structure with the segment end- points as nodes. By appropriate traversal of this graph, the Gauss map boun dary segments may then be identified in proper order, and extraneous segmen ts (lying in the Gauss map interior) are discarded. The symmetrization of t he Gauss map (by identification of antipodal points) and its stereographic projection onto a plane are also discussed. (C) 2001 Elsevier Science B.V. All rights reserved.