Interpolation bilinéaire

L'interpolation bilinéaire est une méthode d'interpolation pour les fonctions de deux variables sur une grille régulière. Elle permet de calculer la valeur d'une fonction en un point quelconque, à partir de ses deux plus proches voisins dans chaque direction. C'est une méthode très utilisée en imagerie numérique pour le redimensionnement d'image, qui permet d'obtenir de meilleurs résultats que l'interpolation par plus proche voisin, tout en restant de complexité raisonnable.

Principe général

La grille régulière à partir de laquelle on interpole. Les quatre points rouges sont les points existants, et le point vert est le point dont on cherche la valeur par interpolation.

Contrairement à ce que son nom suggère, la fonction d'interpolation n'est pas bilinéaire mais quadratique, qui peut se mettre sous la forme :

$f(x,y)=ax+by+cxy+d.$

La valeur Modèle:Math est la valeur interpolée au point de coordonnées Modèle:Math, et Modèle:Math sont des constantes déterminées à partir des quatre voisins Modèle:Math du point Modèle:Math dont on cherche la valeur. Connaissant les valeurs en ces points, on peut écrire un système de 4 équations à 4 inconnues :

$\{\begin{array}{c}f(x_1,y_1)=a{x}_1+b{y}_1+c{x}_1{y}_1+d\\ f(x_2,y_1)=a{x}_2+b{y}_1+c{x}_2{y}_1+d\\ f(x_1,y_2)=a{x}_1+b{y}_2+c{x}_1{y}_2+d\\ f(x_2,y_2)=a{x}_2+b{y}_2+c{x}_2{y}_2+d\end{array}$

L'interpolation bilinéaire peut s'interpréter comme une succession de deux interpolations linéaires, une dans chaque direction.

Solution du système

Exemple d'interpolation bilinéaire dans le carré formé par les points (0,0), (0,1), (1,1) et (1,0) et respectivement de valeur 0, 1, 0.5 et 1. Les valeurs de f(x,y) sont représentées par des variations de couleur.

Un changement de variable permet de simplifier considérablement le système à résoudre. Considérons les nouvelles variables suivantes :

$dx=x-x_1,dy=y-y_1,$

où Modèle:Math sont les coordonnées du coin inférieur gauche. La nouvelle fonction d'interpolation bilinéaire s'écrit alors :

$f(dx,dy)=adx+bdy+cdxdy+d.$

En introduisant les notations $\mathrm{\Delta }x=x_2-x_1$ et $\mathrm{\Delta }y=y_2-y_1$, la matrice à inverser devient :

$A=\left(\begin{array}{cccc}0& 0& 0& 1\\ \mathrm{\Delta }x& 0& 0& 1\\ 0& \mathrm{\Delta }y& 0& 1\\ \mathrm{\Delta }x& \mathrm{\Delta }y& \mathrm{\Delta }x\mathrm{\Delta }y& 1\end{array}\right)$ .

Il reste à introduire les notations suivantes :

$\mathrm{\Delta }{f}_x=f(x_2,y_1)-f(x_1,y_1),\mathrm{\Delta }{f}_y=f(x_1,y_2)-f(x_1,y_1),$

$\mathrm{\Delta }{f}_{xy}=f(x_1,y_1)+f(x_2,y_2)-f(x_2,y_1)-f(x_1,y_2).$

La fonction d'interpolation bilinéaire solution du problème vient alors directement :

$f(x,y)=\mathrm{\Delta }{f}_x\frac{dx}{\mathrm{\Delta }x}+\mathrm{\Delta }{f}_y\frac{dy}{\mathrm{\Delta }y}+\mathrm{\Delta }{f}_{xy}\frac{dx}{\mathrm{\Delta }x}\frac{dy}{\mathrm{\Delta }y}+f(x_1,y_1).$

Articles connexes

• Interpolation bicubique
• Filtrage bilinéaire
• Redimensionnement d'image