Thin plate energy functional explained

The exact thin plate energy functional (TPEF) for a function

f(x,y) is \int ^{y _{1}} _{y _{0}} \int ^{x _{1}} _{x _{0}} (\kappa ^{2} _{1} + \kappa ^{2 )} _{2} \sqrt{g} dx dy

where

\kappa _{1} and \kappa _{2} are the principal curvatures of the surface mapping f at the point (x,y). ^{[1] [2]} This is the surface integral of \kappa ^{2} _{1} + \kappa ^{2 ,} _{2} hence the \sqrt{g} in the integrand.

Minimizing the exact thin plate energy functional would result in a system of non-linear equations. So in practice, an approximation that results in linear systems of equations is often used.^{[3] [4]} The approximation is derived by assuming that the gradient of

f is 0. At any point where f _{x} = f _{y} =0 , the first fundamental form g _{ij} of the surface mapping f is the identity matrix and the second fundamental form b _{ij} is\begin{pmatrix} f _{xx} & f _{xy} \ f _{xy} & f _{yy} \end{pmatrix} .H = b _{ij} g ^{ij} / 2 ^[5] to determine that H = (f _{xx} + f _{yy} )/ 2 and the formula for Gaussian curvature K = b/g (where b and g are the determinants of the second and first fundamental forms, respectively) to determine that K = f _{xx} f _{yy} - (f _{xy} ) ^{2.} Since H = (k _{1+k} _{2)/2} and K = k _{1k} _{2,} the integrand of the exact TPEF equals 4 H ^{2} - 2 K. The expressions we just computed for the mean curvature and Gaussian curvature as functions of partial derivatives of f show that the integrand of the exact TPEF is4 H ^{2} - 2 K = (f _{xx} + f _{yy} ) ^{2} - 2 (f _{xx} f _{yy} - f ^{2 )} _{xy} = f ^{2} _{xx} + 2 f ^{2} _{xy} + f ^{2 .} _{yy}

So the approximate thin plate energy functional is

J[f] = \int ^{y _{1}} _{y _{0}} \int ^{x _{1}} _{x _{0}} f ^{2} _{xx} + 2 f ^{2} _{xy} + f ^{2} _{yy} dx dy.

Rotational invariance

The TPEF is rotationally invariant. This means that if all the points of the surface

z(x,y) are rotated by an angle \theta about the z -axis, the TPEF at each point (x,y) of the surface equals the TPEF of the rotated surface at the rotated (x,y). The formula for a rotation by an angle \theta about the z -axis is

The fact that the

z value of the surface at (x,y) equals the z value of the rotated surface at the rotated (x,y) is expressed mathematically by the equation Z(X,Y) = z(x,y) = (z\circ xy)(X,Y)

where

xy is the inverse rotation, that is, xy(X,Y) = R ^{-1} (X, Y) ^{T} = R ^{T} (X,Y) ^{T} . So Z = z\circ xy and the chain rule implies

In equation,

Z _{0} means Z _{X,} Z _{1} means Z _{Y,} z _{0} means z _{x,} and z _{1} means z _{y.} Equation and all subsequent equations in this section use non-tensor summation convention, that is, sums are taken over repeated indices in a term even if both indices are subscripts. The chain rule is also needed to differentiate equation since z _{j} is actually the composition z _{j} \circ xy: Z _{ik} = z _{jl} R _{kl} R _{ij} .

Swapping the index names

j and k yields

Expanding the sum for each pair

i,j yields \begin{array}{lcl} Z _{XX} & = & R ^{2} _{00} z _{xx} + 2 R _{00} R _{01} z _{xy} + R ^{2} _{01} z _{yy} , \ Z _{XY} & = & R _{00} R _{10} z _{xx} + (R _{00} R _{11} + R _{01} R _{10} )z _{xy} + R _{01} R _{11} z _{yy} , \ Z _{YY} & = & R ^{2} _{10} z _{xx} + 2 R _{10} R _{11} z _{xy} + R ^{2} _{11} z _{yy} . \end{array}

Computing the TPEF for the rotated surface yields

Inserting the coefficients of the rotation matrix

R from equation into the right-hand side of equation simplifies it to z ^{2} _{xx} + 2 z ^{2} _{xy} + z ^{2 .} _{yy}

Data fitting

The approximate thin plate energy functional can be used to fit B-spline surfaces to scattered 1D data on a 2D grid (for example, digital terrain model data).^[6] Call the grid points

(x _{i,y} _{i)} for i =1... N (with x _{i} \in [a,b] and y _{i} \in [c,d] ) and the data values z _{i.} In order to fit a uniform B-spline f(x,y) to the data, the equation

(where

λ is the "smoothing parameter") is minimized. Larger values of λ result in a smoother surface and smaller values result in a more accurate fit to the data. The following images illustrate the results of fitting a B-spline surface to some terrain data using this method.

The thin plate smoothing spline also minimizes equation, but it is much more expensive to compute than a B-spline and not as smooth (it is only

C ^{1} at the "centers" and has unbounded second derivatives there).

Notes and References

1. Variational Design and Fairing of Spline Surfaces. Greiner. Günther. 1994. Eurographics '94. January 3, 2016.
2. Functional Optimization for Fair Surface Design. Moreton. Henry P.. 1992. Computer Graphics. January 4, 2016.
3. Automatic reconstruction of B-splines surfaces of arbitrary topological type. Eck. Matthias. 1996. Proceedings of SIGGRAPH 96, Computer Graphics Proceedings, Annual Conference Series. January 3, 2016.
4. Efficient, Fair Interpolation using Catmull-Clark Surfaces. Halstead. Mark. 1993. Proceedings of the 20th annual conference on Computer graphics and interactive techniques. January 4, 2016.
5. Book: Kreyszig, Erwin. Differential Geometry. limited. Dover. 1991. 0-486-66721-9. Mineola, New York. 131.
6. Multilevel Least Squares Approximation of Scattered Data over Binary Triangulations. Hjelle. Oyvind. 2005. Computing and Visualization in Science. January 14, 2016.