Estimation (géostatistique)
# Estimation (géostatistique)

En géostatistique, l'estimation est la prédiction à partir d'une variable régionalisée pour pallier une lacune d'information.

## Estimation globale

Une estimation globale consiste à proposer une formule a priori de l'estimateur (généralement la moyenne des mesures) et de sa variance.

La variance d'estimation s'exprime par :

${\begin{aligned}{\sigma _{\mathrm {E} }}^{2}&={\frac {1}{[v]^{2}}}\int _{v}\int _{v}C\left(x-y\right)\mathrm {d} x\mathrm {d} y+{\frac {1}{N^{2}}}\sum _{i}\sum _{j}C\left(x_{i}-x_{j}\right)-{\frac {2}{N[v]}}\int _{v}\sum _{i}C\left(x_{i}-y\right)\mathrm {d} y\\\ &={\bar {C}}\left(v,v\right)+{\bar {C}}\left(v',v'\right)-2{\bar {C}}\left(v,v'\right)\end{aligned}}$

Dans les cas suivants, on suppose la géométrie connue (V connu). Quand cela n'est pas assuré, peuvent apparaître des effets de bord. Il peut être nécessaire de travailler alors en géostatistique transitive.

### Échantillonnage aléatoire pur

Si les échantillons sont implantés au hasard, indépendamment entre elles et uniformément dans le champ V à estimer, le problème revient à estimer $\scriptstyle Z_{V}={\frac {1}{V}}\int _{V}Z(x)\mathrm {d} x$ par la moyenne $\scriptstyle {\frac {1}{N}}\sum _{i}Z(x_{i})$.

La variance d'estimation s'écrit à l'aide d'erreurs partielles $Z(X_i)−Z_V$ sous la forme : ${\sigma _{\mathrm {E} }}^{2}={\frac {1}{N^{2}}}\mathbf {Var} \left[\sum _{i}\left(Z\left(X_{i}\right)-Z_{V}\right)\right]$

Sous l'hypothèse stationnaire ou sous l'hypothèse intrinsèque sans dérive, la variance d'estimation s'écrit :

${\sigma _{\mathrm {E} }}^{2}={\frac {1}{N}}\sigma ^{2}\left(o|V\right)$

Sous l'hypothèse stationnaire ou sous l'hypothèse intrinsèque sans dérive, la variance d'estimation s'écrit :

${\begin{aligned}{\sigma _{\mathrm {E} }}^{2}&=\mathbf {Var} \left[{\frac {1}{N}}\sum _{i}\left(Z\left(X_{i}\right)-Z_{V}\right)\right]&(1)\\\ &={\frac {1}{N^{2}}}\mathbf {Var} \left[\sum _{i}\left(Z\left(X_{i}\right)-Z_{V}\right)\right]&(2)\\\ &={\frac {1}{N^{2}}}\mathbf {E} \left[\left(\sum _{i}Z\left(X_{i}\right)-Z_{V}\right)^{2}\right]&(3)\\\ &={\frac {1}{N^{2}}}\mathbf {E} \left[\mathbf {E} \left[\left(\sum _{i}Z\left(X_{i}\right)-Z_{V}\right)^{2}\right]|Z\right]&(4)\\\ {\sigma _{\mathrm {E} }}^{2}\left(\bullet |Z=z\right)&={\frac {1}{N^{2}}}\mathbf {E} \left[\left(\sum _{i}z\left(X_{i}\right)-z_{V}\right)^{2}\right]&(5)\\\ &={\frac {1}{N^{2}}}\sum _{i}\mathbf {E} \left[\left(z\left(X_{i}\right)-z_{V}\right)^{2}\right]&(6)\\\ &={\frac {1}{N^{2}}}\sum _{i}{\frac {1}{V}}\int _{v}\left(z\left(x\right)-z_{V}\right)^{2}\mathrm {d} x&(7)\\\ &={\frac {1}{N^{2}}}\sum _{i}s^{2}(o|V)&(8)\\\ &={\frac {1}{N}}s^{2}\left(o|V\right)\\\ {\sigma _{\mathrm {E} }}^{2}&={\frac {1}{N}}\sigma ^{2}\left(o|V\right)&(9)\end{aligned}}$

* par définition de la variance d'estimation.
* les $Z(X_i)−Z_V$ sont les erreurs partielles.
* hypothèse stationnaire ou hypothèse intrinsèque sans dérive : les erreurs partielles sont d'espérance nulle.
* par espérance conditionnelle.
* à réalisation fixée de la fonction aléatoire (on travaille sur une variance conditionnelle).
* Les $X_i$ sont indépendantes ; les termes croisés sont des covariances de variables aléatoires indépendantes, donc nuls.
* La loi de $X_i$ dans V est uniforme.
* par définition de la variance statistique.
* en déconditionnant l'expression par rapport à Z.

### Échantillonnage aléatoire stratifié

Soit une partition $v_i,$ de volumes identiques v, du domaine à estimer V. Pour chaque sous-domaine est prélevé, indépendamment, un unique échantillon. La variance d'estimation est alors : ${\sigma _{\mathrm {E} }}^{2}={\frac {1}{N}}\sigma ^{2}\left(0|v\right)$

Cette variance d'estimation est inférieure à celle du cas précédent.

### Maille régulière à implantation préférentielle

Soit une partition $v_i,$ de volumes identiques v, du domaine à estimer V. Pour chaque sous-domaine est prélevé, en son centre, un échantillon. La variance d'estimation apparaît comme la somme de trois composantes:

* terme de ligne : variance de l'erreur commise en estimant un volume élémentaire par son échantillon central;
* terme de section : variance de l'erreur commise en estimant un plan par la moyenne pondérée des lignes qu'il contient;
* terme de tranche : variance de l'erreur commise en estimant le champ par la moyenne pondérée de ses sections.

La validité de ce principe de composition n'est pas forcée.

Une règle empirique est qu'un estimateur sera d'autant meilleur, si la fonction aléatoire très structurée, que les mesures seront régulièrement placées, et si la fonction aléatoire est peu structurée, qu'elles seront nombreuses.

## Estimation locale

Une estimation locale construit localement un estimateur à partir des données disponibles. En géostatistique linéaire, la quantité à estimer sera une fonctionnelle linéaire de la variable régionalisée ; de même, l'estimateur sera une combinaison linéaire des données, et l'erreur d'estimation une fonctionnelle linéaire sur la variable régionalisée. Les poids de la combinaison linéaire formant l'estimateur sont donnés par minimisation de la variance d'erreur. Cette estimation locale est dénommée krigeage.

