La interpolación trilineal es un método de interpolación multivariable en un retículo regular tridimensional.^[1]​ Se aproxima linealmente al valor de una función en un punto intermedio ${\displaystyle (x,y,z)}$ dentro de un prisma rectangular axial local, utilizando los datos de la función en los puntos del retículo. Para un retículo no estructurado arbitrario (como se usa en el análisis del método de los elementos finitos), se deben usar otros métodos de interpolación. Si todos los elementos de la malla son tetraedros (símplices 3D), entonces las coordenadas baricéntricas proporciona un procedimiento de cálculo más sencillo.

La interpolación trilineal se utiliza con frecuencia en el análisis numérico, el análisis de datos y la computación gráfica.

La interpolación trilineal es la extensión de la interpolación lineal,^[1]​ que opera en espacios con dimensión ${\displaystyle D=1}$, y de la interpolación bilineal, que opera con la dimensión ${\displaystyle D=2}$, a la dimensión ${\displaystyle D=3}$. Todos estos sistemas de interpolación utilizan polinomios de orden 1, lo que proporciona una precisión de orden 2, y requieren ${\displaystyle 2^{D}=8}$ valores predefinidos adyacentes que rodeen el punto de interpolación. Hay varias formas de llegar a la interpolación trilineal, que es equivalente a la interpolación tridimensional tensorial de orden 1 con B-splines, y el operador de interpolación trilineal también es un producto tensorial de 3 operadores de interpolación lineal.^[2]​

En una red periódica y cúbica, sean ${\displaystyle x_{\text{d}}}$, ${\displaystyle y_{\text{d}}}$ y ${\displaystyle z_{\text{d}}}$ sean las diferencias entre cada uno de ${\displaystyle x}$, ${\displaystyle y}$, ${\displaystyle z}$ y la coordenada más pequeña relacionada, es decir:^[3]​

${\displaystyle {\begin{aligned}x_{\text{d}}={\frac {x-x_{0}}{x_{1}-x_{0}}}\\y_{\text{d}}={\frac {y-y_{0}}{y_{1}-y_{0}}}\\z_{\text{d}}={\frac {z-z_{0}}{z_{1}-z_{0}}}\end{aligned}}}$

donde ${\displaystyle x_{0}}$ indica el punto de red debajo de ${\displaystyle x}$, y ${\displaystyle x_{1}}$ indica el punto de red por encima de ${\displaystyle x}$ y de manera similar para ${\displaystyle y_{0},y_{1},z_{0}}$ y ${\displaystyle z_{1}}$.

Primero se interpola en ${\displaystyle x}$ (imagínese que se está "empujando" la cara del cubo definida por ${\displaystyle C_{0jk}}$ hacia la cara opuesta, definida por ${\displaystyle C_{1jk}}$), dando:

${\displaystyle {\begin{aligned}c_{00}&=c_{000}(1-x_{\text{d}})+c_{100}x_{\text{d}}\\c_{01}&=c_{001}(1-x_{\text{d}})+c_{101}x_{\text{d}}\\c_{10}&=c_{010}(1-x_{\text{d}})+c_{110}x_{\text{d}}\\c_{11}&=c_{011}(1-x_{\text{d}})+c_{111}x_{\text{d}}\end{aligned}}}$

donde ${\displaystyle c_{000}}$ significa el valor de la función de ${\displaystyle (x_{0},y_{0},z_{0}).}$ A continuación, se interpolan estos valores (en ${\displaystyle y}$, "empujando" de ${\displaystyle C_{i0k}}$ a ${\displaystyle C_{i1k}}$), dando:

${\displaystyle {\begin{aligned}c_{0}&=c_{00}(1-y_{\text{d}})+c_{10}y_{\text{d}}\\c_{1}&=c_{01}(1-y_{\text{d}})+c_{11}y_{\text{d}}\end{aligned}}}$

Finalmente, se interpolan estos valores en ${\displaystyle z}$ (recorriendo una recta):

${\displaystyle c=c_{0}(1-z_{\text{d}})+c_{1}z_{\text{d}}.}$

Esto da un valor previsto para el punto.

El resultado de la interpolación trilineal es independiente del orden de los pasos de interpolación en los tres ejes: cualquier otro orden, por ejemplo en ${\displaystyle x}$, luego en ${\displaystyle y}$ y finalmente en ${\displaystyle z}$, produce el mismo valor.

Las operaciones anteriores se pueden visualizar de la siguiente manera: Primero se localizan las ocho esquinas de un cubo que rodean al punto de interés. Estas esquinas tienen los valores ${\displaystyle c_{000}}$, ${\displaystyle c_{100}}$, ${\displaystyle c_{010}}$, ${\displaystyle c_{110}}$, ${\displaystyle c_{001}}$, ${\displaystyle c_{101}}$, ${\displaystyle c_{011}}$, ${\displaystyle c_{111}}$.

A continuación, se realiza una interpolación lineal entre ${\displaystyle c_{000}}$ y ${\displaystyle c_{100}}$ para encontrar ${\displaystyle c_{00}}$, ${\displaystyle c_{001}}$ y ${\displaystyle c_{101}}$ para encontrar ${\displaystyle c_{01}}$, ${\displaystyle c_{011}}$ y ${\displaystyle c_{111}}$ para encontrar ${\displaystyle c_{11}}$, ${\displaystyle c_{010}}$ y ${\displaystyle c_{110}}$ para encontrar ${\displaystyle c_{10}}$.

Ahora se interpola entre ${\displaystyle c_{00}}$ y ${\displaystyle c_{10}}$ para encontrar ${\displaystyle c_{0}}$, ${\displaystyle c_{01}}$ y ${\displaystyle c_{11}}$ para encontrar ${\displaystyle c_{1}}$. Finalmente, se calcula el valor ${\displaystyle c}$ mediante interpolación lineal de ${\displaystyle c_{0}}$ y ${\displaystyle c_{1}}$.

En la práctica, una interpolación trilineal es idéntica a dos interpolaciones bilineales combinadas con una interpolación lineal:^[3]​

${\displaystyle c\approx l\left(b(c_{000},c_{010},c_{100},c_{110}),\,b(c_{001},c_{011},c_{101},c_{111})\right)}$

Una forma alternativa de escribir la solución al problema de interpolación es^[4]​

${\displaystyle f(x,y,z)\approx a_{0}+a_{1}x+a_{2}y+a_{3}z+a_{4}xy+a_{5}xz+a_{6}yz+a_{7}xyz}$

donde los coeficientes se encuentran resolviendo el sistema lineal

${\displaystyle {\begin{aligned}{\begin{bmatrix}1&x_{0}&y_{0}&z_{0}&x_{0}y_{0}&x_{0}z_{0}&y_{0}z_{0}&x_{0}y_{0}z_{0}\\1&x_{1}&y_{0}&z_{0}&x_{1}y_{0}&x_{1}z_{0}&y_{0}z_{0}&x_{1}y_{0}z_{0}\\1&x_{0}&y_{1}&z_{0}&x_{0}y_{1}&x_{0}z_{0}&y_{1}z_{0}&x_{0}y_{1}z_{0}\\1&x_{1}&y_{1}&z_{0}&x_{1}y_{1}&x_{1}z_{0}&y_{1}z_{0}&x_{1}y_{1}z_{0}\\1&x_{0}&y_{0}&z_{1}&x_{0}y_{0}&x_{0}z_{1}&y_{0}z_{1}&x_{0}y_{0}z_{1}\\1&x_{1}&y_{0}&z_{1}&x_{1}y_{0}&x_{1}z_{1}&y_{0}z_{1}&x_{1}y_{0}z_{1}\\1&x_{0}&y_{1}&z_{1}&x_{0}y_{1}&x_{0}z_{1}&y_{1}z_{1}&x_{0}y_{1}z_{1}\\1&x_{1}&y_{1}&z_{1}&x_{1}y_{1}&x_{1}z_{1}&y_{1}z_{1}&x_{1}y_{1}z_{1}\end{bmatrix}}{\begin{bmatrix}a_{0}\\a_{1}\\a_{2}\\a_{3}\\a_{4}\\a_{5}\\a_{6}\\a_{7}\end{bmatrix}}={\begin{bmatrix}c_{000}\\c_{100}\\c_{010}\\c_{110}\\c_{001}\\c_{101}\\c_{011}\\c_{111}\end{bmatrix}},\end{aligned}}}$

dando el resultado

${\displaystyle {\begin{aligned}a_{0}={}&{\frac {-c_{000}x_{1}y_{1}z_{1}+c_{001}x_{1}y_{1}z_{0}+c_{010}x_{1}y_{0}z_{1}-c_{011}x_{1}y_{0}z_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}}+{}\\&{\frac {c_{100}x_{0}y_{1}z_{1}-c_{101}x_{0}y_{1}z_{0}-c_{110}x_{0}y_{0}z_{1}+c_{111}x_{0}y_{0}z_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}},\\[4pt]a_{1}={}&{\frac {c_{000}y_{1}z_{1}-c_{001}y_{1}z_{0}-c_{010}y_{0}z_{1}+c_{011}y_{0}z_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}}+{}\\&{\frac {-c_{100}y_{1}z_{1}+c_{101}y_{1}z_{0}+c_{110}y_{0}z_{1}-c_{111}y_{0}z_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}},\\[4pt]a_{2}={}&{\frac {c_{000}x_{1}z_{1}-c_{001}x_{1}z_{0}-c_{010}x_{1}z_{1}+c_{011}x_{1}z_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}}+{}\\&{\frac {-c_{100}x_{0}z_{1}+c_{101}x_{0}z_{0}+c_{110}x_{0}z_{1}-c_{111}x_{0}z_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}},\\[4pt]a_{3}={}&{\frac {c_{000}x_{1}y_{1}-c_{001}x_{1}y_{1}-c_{010}x_{1}y_{0}+c_{011}x_{1}y_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}}+{}\\&{\frac {-c_{100}x_{0}y_{1}+c_{101}x_{0}y_{1}+c_{110}x_{0}y_{0}-c_{111}x_{0}y_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}},\\[4pt]a_{4}={}&{\frac {-c_{000}z_{1}+c_{001}z_{0}+c_{010}z_{1}-c_{011}z_{0}+c_{100}z_{1}-c_{101}z_{0}-c_{110}z_{1}+c_{111}z_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}},\\[4pt]a_{5}=&{\frac {-c_{000}y_{1}+c_{001}y_{1}+c_{010}y_{0}-c_{011}y_{0}+c_{100}y_{1}-c_{101}y_{1}-c_{110}y_{0}+c_{111}y_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}},\\[4pt]a_{6}={}&{\frac {-c_{000}x_{1}+c_{001}x_{1}+c_{010}x_{1}-c_{011}x_{1}+c_{100}x_{0}-c_{101}x_{0}-c_{110}x_{0}+c_{111}x_{0}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}},\\[4pt]a_{7}={}&{\frac {c_{000}-c_{001}-c_{010}+c_{011}-c_{100}+c_{101}+c_{110}-c_{111}}{(x_{0}-x_{1})(y_{0}-y_{1})(z_{0}-z_{1})}}.\end{aligned}}}$

• Interpolación lineal
• Interpolación bilineal
• Interpolación tricúbica
• Interpolación lineal esférica

• pseudocódigo de la NASA, describe una interpolación trilineal inversa iterativa (dados los vértices y el valor de C, encuentre Xd, Yd y Zd).
• Paul Bourke, Métodos de interpolación, 1999. Contiene un método muy inteligente y simple para encontrar la interpolación trilineal que se basa en la lógica binaria y se puede extender a cualquier dimensión (tetralineal, Pentalineal, ...).
• Kenwright, Deformación del tetraedro de forma libre. Simposio Internacional sobre Computación Visual. Springer International Publishing, 2015 [1].