Question

Trilinear interpolation approximates the value of a point (x, y, z) inside a cube using the values at the cube vertices. I´m trying to do an “inverse” trilinear interpolation. Knowing the values at the cube vertices and the value attached to a point how can I find (x, y, z)? Any help would be highly appreciated. Thank you!

Answer 1

You are solving for 3 unknowns given 1 piece of data, and as you are using a linear interpolation your answer will typically be a plane (2 free variables). Depending on the cube there may be no solutions or a 3D solution space.

I would do the following. Let v be the initial value. For each “edge” of the 12 edges (pair of adjacent vertices) of the cube look to see if 1 vertex is >=v and the other <=v – call this an edge that crosses v.

If no edges cross v, then there are no possible solutions.

Otherwise, for each edge that crosses v, if both vertices for the edge equal v, then the whole edge is a solution. Otherwise, linearly interpolate on the edge to find the point that has a value of v. So suppose the edge is (x1, y1, z1)->v1 <= v <= (x2, y2, z2)->v2.

s = (v-v1)/(v2-v1)
(x,y,z) = (s*(x2-x1)+x1, (s*(y2-y1)+y1, s*(z2-z1)+z1)

This will give you all edge points that are equal to v. This is a solution, but possibly you want an internal solution – be aware that if there is an internal solution there will always be an edge solution.

If you want an internal solution then just take any point linearly between the edge solutions – as you are linearly interpolating then the result will also be v.