Implicit surface

implicit surface torus (R=40, a=15)

implicit surface of genus 2

implicit non algebraic surface (wineglas)

In mathematics an implicit surface is a surface in Euclidean space defined by an equation

$F(x,y,z)=0$

An implicit surface is the set of zeros of a function of 3 variables. Implicit means, that the equation is not solved for x or y or z.

The graph of a function is usually described by an equation $z=f(x,y)$ and is called an explicit representation. The third essential description of a surface is the parametric one: $(x(s,t),y(s,t),z(s,t))$, where the x-, y- and z-coordinates of surface points are represented by three functions $x(s,t)\,,y(s,t)\,,z(s,t)$ depending on common parameters $s,t$. The change of representations is unsually simple only, when the explicit representation $z=f(x,y)$ is given: $z-f(x,y)=0$ (implicit), $(s,t,f(s,t))$ (parametric).

Examples:

1. a plane $x+2y-3z+1=0$,
2. a sphere $x^{2}+y^{2}+z^{2}-4=0$,
3. a torus $(x^{2}+y^{2}+z^{2}+R^{2}-a^{2})^{2}-4R^{2}(x^{2}+y^{2})=0$,
4. surface of genus 2: $2y(y^{2}-3x^{2})(1-z^{2})+(x^{2}+y^{2})^{2}-(9z^{2}-1)(1-z^{2})=0$ (s. picture),
5. surface of revolution $x^{2}+y^{2}-(\ln(z+3.2))^{2}-0.02=0$ (s. picture wine glas).

For a plane, a sphere and a torus there exist simple parametric representations. This is not true for the 4. example.

The implicit function theorem describes conditions, under which an equation $F(x,y,z)=0$ can be solved (theoretically) for x , y or z. But in general the solution may not be feasible. This theorem is the key for the computation of essential geometric features of the surface: tangent planes, surface normals, curvatures (s. below). But they have an essential drawback: their visualization is difficult.

If $F(x,y,z)$ is polynomial in x,y and z, the surface is called algebraic.

Example 5. is non algebraic.

Despite the visualization of an implicit surface is difficult, they provide rather simple techniques to generate theoretically (e.g. Steiner surface) and practically (s. below) interesting surfaces.

Formulas

Throughout the following considerations the implicit surface is represented by an equation $F(x,y,z)=0$ where function $F$ meets the necessary conditions of differentiability. The partial derivatives of $F$ are $F_{x},...,F_{xx},...$.

Tangent plane and normal vector

A surface point $(x_{0},y_{0},z_{0})$ is called regular, if

• $(F_{x}(x_{0},y_{0},z_{0}),F_{y}(x_{0},y_{0},z_{0}),F_{z}(x_{0},y_{0},z_{0}))\neq (0,0,0)$ ,

otherwise the point is singular.

The equation of the tangent plane at a regular point $(x_{0},y_{0},z_{0})$ is

• $F_{x}(x_{0},y_{0},z_{0})(x-x_{0})+F_{y}(x_{0},y_{0},z_{0})(y-y_{0})+F_{z}(x_{0},y_{0},z_{0})(z-z_{0})=0$, and

$\mathbf {n} (x_{0},y_{0},z_{0})=(F_{x}(x_{0},y_{0},z_{0}),F_{y}(x_{0},y_{0},z_{0}),F_{z}(x_{0},y_{0},z_{0}))^{T}$ is a normal vector.

Normal curvature

In order to keep the formula simple the arguments $(x_{0},y_{0},z_{0})$ are omitted:

• $\kappa _{n}={\frac {\mathbf {v} ^{\top }H_{F}\mathbf {v} }{\|\operatorname {grad} F\|}}$ is the

normal curvature of the surface at a regular point for the unit tangent direction $\mathbf {v}$. $H_{F}$ is the Hessian matrix of $F$ (matrix of the second derivatives).

The proof of this formula relies (as in the case of an implicit curve) on the implict function theorem and the formula for the normal curvature of a parametric surface.

Applications of implicit surfaces

As in the case of implicit curves it is an easy task to generate implicit surfaces with desired shapes by applying algebraic operations (addition, multiplication) on simple premitives.

equipotential surface of 4 point charges

Equipotential surface of point charges

The electrical potential of a point charge $q_{i}$ at point $\mathbf {p} _{i}=(x_{i},y_{i},z_{i})$ generates at point $\mathbf {p} =(x,y,z)$ the potential (omitting physical constants)

$F_{i}(x,y,z)={\frac {q_{i}}{\|\mathbf {p} -\mathbf {p} _{i}\|}}$.

The equipotential surface for the potential value $c$ is the implicit surface $F_{i}(x,y,z)-c=0$ which is a sphere with center at point $\mathbf {p} _{i}$.

The potential of (e.g.) $4$ point charges is represented by

$F(x,y,z)={\frac {q_{1}}{\|\mathbf {p} -\mathbf {p} _{1}\|}}+{\frac {q_{2}}{\|\mathbf {p} -\mathbf {p} _{2}\|}}+{\frac {q_{3}}{\|\mathbf {p} -\mathbf {p} _{3}\|}}+{\frac {q_{4}}{\|\mathbf {p} -\mathbf {p} _{4}\|}}$.

For the picture the four charges are equal 1 and are located at the points $(\pm 1,\pm 1,0)$. The displayed surface is the equipotential surface (implicit surface) $F(x,y,z)-2.8=0$.

Constant distance product surface

A Cassini oval can be defined as the point set for which the product of the distances to two given points is constant (for an ellipse the sum is constant !). In a similar way implicit surfaces can be defined by a constant distance product to several fixed points.

In picture metamorphoses the upper left surface is generated by this rule: With

$F(x,y,z)={\sqrt {(x-1)^{2}+y^{2}+z^{2}}}\cdot {\sqrt {(x+1)^{2}+y^{2}+z^{2}}}\cdot {\sqrt {x^{2}+(y-1)^{2}+z^{2}}}\cdot {\sqrt {x^{2}+(y+1)^{2}+z^{2}}}$

the constant distance product surface $F(x,y,z)-1.1=0$ is displayd.

metamorphoses between two implicit surfaces (torus and a constant distance product surface)

Metamorphoses of implicit surfaces

A further simple method to generate new implicit surfaces is called metamorphoses of implicit surfaces:

For two implicit surfaces $F_{1}(x,y,z)=0,F_{2}(x,y,z)=0$ (in picture: a constant distance product surface and a torus) one defines new surfaces using the design parameter $\mu \in [0,1]$:

$F(x,y,z)=\mu F_{1}(x,y,z)+(1-\mu )F_{2}(x,y,z)=0$

For the picture the design parameter is: $\mu =0,\,0.33,\,0.66,\,1$ .

approximation of three tori (parallel projection)

PovRay-image (central projection) of an approximation of three tori

Smooth approximationen of several implicit surfaces

Analogously to the smooth approximation with implicit curves the equation

$F(x,y,z)=F_{1}(x,y,z)\cdot F_{2}(x,y,z)\cdot F_{3}(x,y,z)-c=0$

represents for suitable parameters $c$ smooth approximations of three intersecting tori with equations

$F_{1}=(x^{2}+y^{2}+z^{2}+R^{2}-a^{2})^{2}-4R^{2}(x^{2}+y^{2})=0$,

$F_{2}=(x^{2}+y^{2}+z^{2}+R^{2}-a^{2})^{2}-4R^{2}(x^{2}+z^{2})=0$,

$F_{3}=(x^{2}+y^{2}+z^{2}+R^{2}-a^{2})^{2}-4R^{2}(y^{2}+z^{2})=0$.

(For the picture the parameters are: $R=1,\,a=0.2,\,c=0.01$)

PovRay-Bild: metamorphoses between a sphere and a constant distance product surface (6 points)

Visualization of implicit surfaces

The visualization of implicit surfaces needs great efforts. Essentially there are two ideas visualizing an implicit surface. One generates a net of polygons which is visualized (see surface triangulation) and the second relies on ray tracing which determines intersection points of rays with the surface.

Weblinks

• Sultanow: Implizite Flächen
• Hartmann: Geometry and Algorithms for COMPUTER AIDED DESIGN
• GEOMVIEW
• K3Dsurf: 3d surface generator
• SURF: Visualisierung algebraischer Flächen

References

• Thorpe: Elementary Topics in Differential Geometry, Springer-Verlag, New York, 1979, ISBN 0-387-90357-7