i
i
i
i
i
i
i
i
6 1. Surface Representations
Since our surface representations are supposed to support efficient pro-
cessing, a natural choice is to restrict functions to the class of polynomi-
als because those can be evaluated by elementary arithmetic operations.
Another justification for the restriction to polynomials is the well-known
Weierstrass theorem that guarantees that each smooth function can be
approximated by a polynomial up to any desired precision [Ross 80].
From calculus we know that a C
∞
function g with bounded derivatives
can be approximated over an interval of length h by a polynomial of degree
p such that the approximation error behaves like O(h
p+1
) (e.g., Taylor’s
theorem or generalized mean value theorem) [Rudin 02]. As a consequence
there are, in principle, two possibilities to improve the accuracy of an ap-
proximation with piecewise polynomials. We can either raise the degree of
the polynomial (p-refinement) or we can reduce the size of the individual
segments and use more segments for the approximation (h-refinement).
In geometry processing applications, h-refinement is usually preferred
over p-refinement since, for a discretely sampled input surface, we can-
not make reasonable assumptions about the boundedness of higher-order
derivatives. Moreover, for piecewise polynomials with higher degree, the
C
k
smoothness conditions between segments are sometimes quite difficult
to satisfy. Finally, with today’s computer architectures, processing a large
number of very simple objects is often much more efficient than processing a
smaller number of more complex ones. This is why the somewhat extremal
choice of C
0
piecewise linear surface representations, i.e., polygonal meshes,
have become the widely established standard in geometry processing.
While, for parametric surfaces, the O(h
p+1
) approximation error esti-
mate follows from the mean value theorem in a straightforward manner,
a more careful consideration is necessary for implicit representations. The
generalized mean value theorem states that if a sufficiently smooth function
g over an interval [a, a + h] is interpolated at the abscissae t
0
, . . . , t
p
by a
polynomial f of degree p, then the approximation error is bounded by
|f(t) − g(t)| ≤
1
(p + 1)!
max f
(p+1)
p
Y
i=0
(t
i
− t) = O(h
p+1
).
For an implicit representation G: IR
3
→ IR and the corresponding polyno-
mial approximant F , this theorem is still valid; however, here the actual
surface geometry is not defined by the function values G(x), for which this
theorem gives an error estimate, but by the zero level set of G, i.e., by
S = {x ∈ IR
3
|G(x) = 0}.
Consider a point x on the implicit surface defined by the approximating
polynomial F , i.e., F (x) = 0 within some voxel. We can find a correspond-
ing point x + d on the implicit surface defined by G, i.e., G(x + d) = 0
by shooting a ray in normal direction to F , i.e., d = d ∇F/k∇F k. For a