Interpolation multivariée

En analyse numérique, l'interpolation multivariée ou l'interpolation spatiale désigne l'interpolation numérique de fonctions de plus d'une variable.

Le problème est similaire à celui de l'interpolation polynomiale sur un intervalle réel : on connait les valeurs d'une fonction à interpoler aux points $(x_{i},y_{i},z_{i},\dots )$ et l'objectif consiste à évaluer la valeur de la fonction en des points $(x,y,z,\dots )$ .

L'interpolation multivariée est notamment utilisée en géostatistique, où elle est utilisée pour reconstruire les valeurs d'une variable régionalisée sur un domaine à partir d'échantillons connus en un nombre limité de points. Par exemple en météorologie, il s'agit de l'estimation de valeurs intermédiaires inconnues à partir de valeurs discrètes connues d'une variable dépendante, comme la température, sur une carte météorologique^[1].

Grille régulière[modifier | modifier le code]

Pour des fonctions connues sur une grille régulière (avec des intervalles prédéterminés, non nécessairement équidistants), les méthodes suivantes sont applicables.

Toute dimension[modifier | modifier le code]

Interpolation au plus proche voisin

2 dimensions[modifier | modifier le code]

Le redimensionnement d'image est l'application de l'interpolation dans le traitement d'images.

Trois méthodes sont ici appliquées sur un ensemble de 4x4 points.

Par proche voisin
Bilinéaire
Bicubique

Voir aussi les points de Padua pour l'interpolation polynomiale de deux variables.

3 dimensions[modifier | modifier le code]

Produit tensoriel en dimension N[modifier | modifier le code]

Les splines de Catmull-Rom peuvent être facilement généralisées en dimension quelconque. Les splines cubiques d'Hermite donnent $\mathrm {CINT} _{x}(f_{-1},f_{0},f_{1},f_{2})=\mathbf {b} (x)\cdot \left(f_{-1}f_{0}f_{1}f_{2}\right)$ pour un 4-vecteur $\mathbf {b} (x)$ donné, qui est donc une fonction de x, où $f_{j}$ est la valeur en $j$ de la fonction à interpoler.

En réécrivant cette approximation sous la forme

\mathrm {CR} (x)=\sum _{i=-1}^{2}f_{i}b_{i}(x)

cette formule peut être généralisée en dimension N^[2]

\mathrm {CR} (x_{1},\dots ,x_{N})=\sum _{i_{1},\dots ,i_{N}=-1}^{2}f_{i_{1}\dots i_{N}}\prod _{j=1}^{N}b_{i_{j}}(x_{j})

On remarque que des généralisations similaires peuvent être faites pour d'autres types d'interpolation par splines, dont les splines d'Hermite. En termes d'efficacité, la formule générale peut en effet être calculée comme une composition successive d'opérations de type CINT pour tout type de produit tensoriel de splines, comme dans le cas de l'interpolation tricubique. Cependant, il demeure que s'il y a n termes dans le terme en CR de la somme en dimension 1, il y aura alors n^N termes dans la somme en dimension N.

Grille irrégulière (données éparses)[modifier | modifier le code]

Les méthodes définies pour des données éparses sur une grille irrégulière peuvent être appliquées sur une grille régulière, ce qui permet de revenir à un cas connu.

Interpolation au plus proche voisin
Interpolation par voisins naturels sur un réseau irrégulier triangulé
Interpolation linéaire sur réseau irrégulier triangulé (une forme de fonctions linéaires par morceaux)
Pondération inverse à la distance
Krigeage
Interpolation par fonction de base radiale
Spline en plaque mince
Spline polyharmonique (la spline en plaque mince en est un cas particulier)
Spline par moindres carrés

Utilisations[modifier | modifier le code]

Interpolation et un lissage à deux dimensions

En météorologie, l'analyse d'une carte météorologique ;
En cartographie, une carte du niveau du terrain.

Notes et références[modifier | modifier le code]

(en) Cet article est partiellement ou en totalité issu de l’article de Wikipédia en anglais intitulé « Multivariate interpolation » (voir la liste des auteurs).

↑ Organisation météorologique mondiale, « Interpolation », sur Eumetcal (consulté le 16 novembre 2013)
↑ Two hierarchies of spline interpolations. Practical algorithms for multivariate higher order splines

Liens externes[modifier | modifier le code]

(en) Example C++ code for several 1D, 2D and 3D spline interpolations (including Catmull-Rom splines).
(en) Multi-dimensional Hermite Interpolation and Approximation, Prof. Chandrajit Bajaja, Université Purdue

Portail des mathématiques

[OMM-1] Organisation météorologique mondiale, « Interpolation », sur Eumetcal (consulté le 16 novembre 2013)

[2] Two hierarchies of spline interpolations. Practical algorithms for multivariate higher order splines

[1]

[2]