Génération de données pour suivre un variogramme donné

8

Il s'agit d'une approche simple ayant un ensemble de coordonnées (par exemple, en 2D comme {x,y}) et au moins une variable associée (par exemple v) pour calculer un variogramme comme descripteur de la dépendance spatiale de la variable à vtravers le champ étudié.

entrez la description de l'image ici entrez la description de l'image ici

La question qui m'est apparue est:
comment générer une réalisation d'un ensemble de données ayant un variogramme? (mouvement inverse!)
C'est-à-dire qu'il y a au moins un variogramme disponible mais ni jeu de données ni autre description n'est disponible et le but est de générer une réalisation de jeu de données original (inconnu) qui pourrait avoir un tel variogramme.
Quelle est la probabilité d'avoir une telle réalisation?

Mises à jour / commentaires: Dans variogramle contexte ci-dessus, je veux dire variogramme empirique. Je suppose que l'ajustement d'un modèle de variogramme n'est pas un problème au moins pour cette question. Le variogramme est également disponible sous forme de paires (h, gamma).

Développeur
la source
Par variogramme, vous voulez dire modèle de variogramme ?
Paul Hiemstra

Réponses:

7

Vous pouvez utiliser la simulation séquentielle pour générer des réalisations d'un champ aléatoire qui a la structure de covariance donnée dans le modèle de variogramme. Dans R, cela peut être fait en utilisant gstat. Voir démo (ugsim) et démo (uisim) des exemples de code R de gstat.

Paul Hiemstra
la source
1
+1 Tout d'abord, bien sûr, vous devez adapter un modèle de variogramme au variogramme empirique indiqué. Notez que geoRpermet également d'effectuer facilement des simulations.
whuber
Mais je suppose que l'OP a déjà un variogramme. Ou bien ce serait un problème de poulet et d'oeuf :).
Paul Hiemstra
Le graphique de la question, Paul, montre un variogramme empirique, pas un modèle de variogramme. Vous ne pouvez pas utiliser directement un variogramme empirique pour piloter une simulation stochastique (il ne serait pas défini positif, donc la simulation rencontrerait probablement des problèmes numériques). Mais il n'y a pas de problème poulet / œuf: on adapte le modèle de variogramme à ce variogramme empirique puis on procède. Certaines informations sont perdues - une bonne analyse variographique considère bien plus qu'un seul variogramme empirique in vacuo - mais cela ne peut pas être aidé ici.
whuber
Par poulet et œuf, je veux dire que sans données pour s'adapter au modèle de variogramme, il n'y a pas de modèle de variogramme pour générer de nouvelles données. Et par variogramme, je veux dire modèle de variogramme.
Paul Hiemstra
2
Une partie de ma thèse a été consacrée à l'ajustement automatique des variogrammes, donc je pense que c'est possible;)
Paul Hiemstra