Méthodes Open Source pour le krigeage?


42

Je voudrais un jeu de données de points que je voudrais Krige , idéalement en utilisant un logiciel open-source. Si possible, j'aimerais également choisir le modèle de semi-variogramme pendant le processus pour améliorer l'estimation.

Réponses:


45

Selon le type de Kriging que vous souhaitez appliquer, vous pouvez choisir parmi différents packages:

Krigeage ordinaire

La version la plus courante est implémentée par exemple dans:

Krigeage simple

Le Kriging simple utilise la moyenne de l'ensemble des données alors que le Kriging ordinaire utilise une moyenne locale. Par conséquent, le krigeage simple peut être moins précis, mais il produit généralement des résultats plus "lisses". Il est implémenté dans:

Krigeage Universel

Le krigeage universel permet de prendre en compte la dérive des données. Les implémentations sont incluses dans:

Autres types de krigeage

GRASS v.krige prend également en charge Block Kriging.

HPGL implémente un grand nombre de méthodes de Kriging moins connues (consultez le manuel pour plus d'informations sur celles-ci):

  • Indicateur de krigeage (IK)
  • Krigeage moyen local variable (LVM Kriging)
  • CoKrigeage simple (modèles de Markov 1 et 2)
  • Simulation par indicateurs séquentiels (SIS)
  • Corellogram Local Variable SIS (CLVM SIS)
  • Moyenne variable locale SIS (LVM SIS)
  • Simulation gaussienne séquentielle (SGS)
  • Simulation gaussienne tronquée (GTSIM) [dans la collection de scripts Python]

SAGA propose différentes versions du krigeage ordinaire et universel.

Gstat krige prend également en charge le krigeage en blocs et en point.


1
Grande réponse Kriging!
Ragi Yaser Burhum

1
Vos liens faunalia ne fonctionnent pas ...
Alex Leith

17

Il semble y avoir quelques options avec GRASS GIS. Consultez la page Wiki de GRASS Kriging: http://grass.osgeo.org/wiki/Kriging

Un projet Google Summer of Code en 2009 a produit V.krige: http://grass.osgeo.org/wiki/V.krige_GSoC_2009

Le paquet GPL gstat devrait fonctionner seul ou interfacé avec GRASS GIS. http://www.gstat.org/

Dylan Beaudette a un bel exemple de krigeage avec GRASS. http://casoilresource.lawr.ucdavis.edu/drupal/node/438 (son blog regorge d'exemples intéressants et intéressants d'utilisation de SIG OpenSource et d'outils statistiques!)


14

Le projet R contient un nombre important de progiciels de statistiques spatiales , mais la courbe d’apprentissage de R est plutôt abrupte.


Les gens disent toujours ça, mais je me demande: raide par rapport à quoi?
Matt Parker

J'ai vu le commentaire "courbe d'apprentissage abrupte" jeté à plusieurs reprises sur R - cela n'a tout simplement aucun sens pour moi. J'avais un an dans ma relation avec MATLAB lorsque j'ai découvert R. J'ai trouvé R tellement facile à apprendre que j'ai salué MATLAB avec le salut à un doigt et que j'ai rapidement cessé de l'utiliser.
Sharpie

1
Je pense que c'est parce que les gens essaient rarement de comprendre les statistiques, et à cause de cela que la syntaxe de la courbe d'apprentissage est raide, il y a rarement des problèmes pour le relever
dodobas

2
Je pense que la syntaxe est l’une des langues les plus faciles à apprendre. Voici un exemple de langage statistique facile à apprendre à partir de la ligne de commande. Je pense que les gens se plaignent parce que ce n'est pas Excel.
TheSteve0

1
C'est une étape en comparaison avec un programme basé sur une interface graphique. Si vous avez utilisé Windows toute votre vie et des programmes basés sur une interface graphique, vous serez exécuté lorsque vous verrez l'aspect de la ligne de commande. La comparaison Excel les fait effectivement fonctionner. Mais R est très simple à utiliser si quelqu'un peut vous montrer les astuces de base. Vous devez être prêt à apprendre de nouveaux concepts tels que vecteurs, matrices, fonctions, boucles, qui n'existent pas dans un monde Excel / Windows. Si vous utilisiez auparavant Linux, ce ne serait pas une courbe en étapes.
Niculita Mihai

14

Si vous êtes heureux de lire votre raster dans un tableau numpy ( gdal peut le faire), vous pouvez utiliser l' implémentation de la bibliothèque géostatistique à hautes performances de Python ou C / C ++.

HPGL implémente les algorithmes suivants:

  1. Krigeage simple (SK)
  2. Krigeage ordinaire (OK)
  3. Indicateur de krigeage (IK)
  4. Krigeage moyen local variable (LVM Kriging)
  5. CoKrigeage simple (modèles de Markov 1 et 2)
  6. Simulation par indicateurs séquentiels (SIS)
  7. Corellogram Local Variable SIS (CLVM SIS)
  8. Moyenne variable locale SIS (LVM SIS)
  9. Simulation gaussienne séquentielle (SGS)
  10. Simulation gaussienne tronquée (GTSIM) [dans la collection de scripts Python]

Je ne l'ai pas utilisé moi-même, mais j'ai entendu de bonnes choses à ce sujet, notamment en ce qui concerne la vitesse.



4

Je me souviens d’avoir utilisé SAGA pour faire cela il ya quelques années pour certaines sorties de modélisation des crues. Open Source et vaut bien un coup d'oeil.


4

gvSIG (un autre SIG gratuit) permet le krigeage en utilisant Sextante. C’est fondamentalement la même chose que d’utiliser SAGA, mais gvSIG fournit une expérience gis plus «typique» (c’est-à-dire semblable à ESRI).


3

Vous pouvez essayer le modèle Kriging dans Surfpack version 1.1 (je l’ai écrit alors que j’étais encore dans l’équipe DAKOTA), ou la dernière et meilleure version fournie avec la version "stable" de DAKOTA (Surfpack est un sous-package de DAKOTA). , il fait du krigeage universel du point de vue des fonctions de corrélation plutôt que des semi-variogrammes.

Récemment, un utilisateur, Joel Guerrero, l'a comparé à une foule d'autres implémentations et a déclaré: "Toujours lié au surfpack, nous le comparons à d'autres implémentations (y compris une implémentation commerciale), et jusqu'à présent, il les surperforme toutes, au point qu'il semble parfois que cela fait de la magie noire "


3

GSLIB (Geostatistical Software Library) est un logiciel de pointe basé sur les fichiers et les commandes, développé par l’Université Stanford et publié dans les années 1990, avec quelques mises à jour au cours de la dernière décennie. Le code source peut être téléchargé et compilé librement sur Linux / Windows à l'aide d'un compilateur Fortran. Il existe des ressources en ligne et un livre disponible.

Le krigeage est l'une des forces du logiciel:

  • Krigeage sur grille 1, 2 ou 3-D, validation croisée, jackknifing
  • SK, OK, UK, krigeage avec dérive externe
  • cokriguer
  • indicateur de krigeage
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.