Trilinear interpolation

From Wikipedia, the free encyclopedia
Jump to: navigation, search

Trilinear interpolation is a method of multivariate interpolation on a 3-dimensional regular grid. It approximates the value of an intermediate point (x, y, z) within the local axial rectangular prism linearly, using data on the lattice points. For an arbitrary, unstructured mesh (as used in finite element analysis), other methods of interpolation must be used; if all the mesh elements are tetrahedra (3D simplices), then barycentric coordinates provide a straightforward procedure.

Trilinear interpolation is frequently used in numerical analysis, data analysis, and computer graphics.

Compared to linear and bilinear interpolation[edit]

Trilinear interpolation is the extension of linear interpolation, which operates in spaces with dimension D=1, and bilinear interpolation, which operates with dimension D=2, to dimension D=3. The order of accuracy is 1 for all these interpolation schemes, and it requires (1 + n)^D = 8 adjacent pre-defined values surrounding the interpolation point. There are several ways to arrive at trilinear interpolation, it is equivalent to 3-dimensional tensor B-spline interpolation of order 1, and the trilinear interpolation operator is also a tensor product of 3 linear interpolation operators.

Method[edit]

Eight corner points on a cube surrounding the interpolation point C
Depiction of 3D interpolation
A geometric visualisation of trilinear interpolation. The product of the value at the desired point and the entire volume is equal to the sum of the products of the value at each corner and the partial volume diagonaly opposite the corner.

On a periodic and cubic lattice, let x_d, y_d, and z_d be the differences between each of x, y, z and the smaller coordinate related, that is:

 \ x_d = (x - x_0)/(x_1 - x_0)
 \ y_d = (y - y_0)/(y_1 - y_0)
 \ z_d = (z - z_0)/(z_1 - z_0)

where  x_0 indicates the lattice point below  x , and   x_1 indicates the lattice point above  x and similarly for y_0, y_1, z_0 and z_1.

First we interpolate along x (imagine we are pushing the front face of the cube to the back), giving:

 \ c_{00} = V[x_0,y_0, z_0]  (1 - x_d) + V[x_1, y_0, z_0] x_d
 \ c_{10} = V[x_0,y_1, z_0]  (1 - x_d) + V[x_1, y_1, z_0] x_d
 \ c_{01} = V[x_0,y_0, z_1]  (1 - x_d) + V[x_1, y_0, z_1] x_d
 \ c_{11} = V[x_0,y_1, z_1]  (1 - x_d) + V[x_1, y_1, z_1] x_d

Where V[x_0,y_0, z_0] means the function value of  (x_0,y_0,z_0). Then we interpolate these values (along y, as we were pushing the top edge to the bottom), giving:

 \ c_0 = c_{00}(1 - y_d) + c_{10}y_d
 \ c_1 = c_{01}(1 - y_d) + c_{11}y_d

Finally we interpolate these values along z(walking through a line):

 \ c = c_0(1 - z_d) + c_1z_d .

This gives us a predicted value for the point.

The result of trilinear interpolation is independent of the order of the interpolation steps along the three axes: any other order, for instance along x, then along y, and finally along z, produces the same value.

The above operations can be visualized as follows: First we find the eight corners of a cube that surround our point of interest. These corners have the values C000, C100, C010, C110, C001, C101, C011, C111.

Next, we perform linear interpolation between C000 and C100 to find C00, C001 and C101 to find C01, C011 and C111 to find C11, C010 and C110 to find C10.

Now we do interpolation between C00 and C10 to find C0, C01 and C11 to find C1. Finally, we calculate the value C via linear interpolation of C0 and C1

In practice, a trilinear interpolation is identical to three successive linear interpolations, or a bilinear interpolation combined with a linear interpolation:

 C \approx\ l( b(C_{000}, C_{010}, C_{100}, C_{110}), b(C_{001}, C_{011}, C_{101}, C_{111}))

See also[edit]

External links[edit]

  • pseudo-code from NASA, describes an iterative inverse trilinear interpolation (given the vertices and the value of C find Xd, Yd and Zd).
  • Paul Bourke, Interpolation methods, 1999. Contains a very clever and simple method to find trilinear interpolation that is based on binary logic and can be extended to any dimension (Tetralinear, Pentalinear, ...).