Tengo un conjunto de datos de puntos que me gustaría Krige , idealmente usando un paquete de software de código abierto. Si es posible, también me gustaría elegir el modelo de semi-variograma durante el proceso para mejorar la estimación.
Dependiendo del tipo de Kriging que desee aplicar, hay diferentes paquetes para elegir:
La versión más común se implementa, por ejemplo, en:
Kriging simple usa el promedio de todo el conjunto de datos, mientras que Kriging Ordinario usa un promedio local. Por lo tanto, el Kriging simple puede ser menos preciso, pero generalmente produce resultados "más suaves". Se implementa en:
Kriging universal permite considerar la deriva en los datos. Las implementaciones están incluidas en:
Otros tipos de Kriging
GRASS v.krige también es compatible con Block Kriging.
HPGL implementa una gran cantidad de métodos de Kriging menos conocidos (consulte el manual para obtener más información al respecto):
SAGA ofrece diferentes versiones de Kriging tanto ordinario como universal.
Gstat krige también es compatible con Block y Point Kriging.
Parece que hay algunas opciones con GRASS GIS. Echa un vistazo a la página Wiki de GRASS Kriging: http://grass.osgeo.org/wiki/Kriging
Un proyecto de Google Summer of Code en 2009 produjo V.krige: http://grass.osgeo.org/wiki/V.krige_GSoC_2009
El paquete GPL gstat debería funcionar solo o en interfaz con GRASS GIS. http://www.gstat.org/
Dylan Beaudette tiene un buen ejemplo de hacer kriging con GRASS. http://casoilresource.lawr.ucdavis.edu/drupal/node/438 (¡Su blog está lleno de excelentes e interesantes ejemplos del uso de OpenSource GIS y herramientas estadísticas!)
fuente
El proyecto R tiene un número considerable de paquetes de software de estadísticas espaciales , pero R tiene una curva de aprendizaje bastante empinada.
fuente
Si está contento de leer su ráster en una matriz numpy ( gdal puede hacer esto), entonces podría usar la implementación de la Biblioteca de geoestadística de alto rendimiento de Python o C / C ++.
HPGL implementa los siguientes algoritmos:
No lo he usado yo mismo, pero he escuchado cosas buenas al respecto, especialmente con respecto a la velocidad.
fuente
Consulte este libro gratuito, se trata de hacer geoestadística en R, y también contiene información sobre cómo hacerlo en SAGA y GRASS. http://spatial-analyst.net/book/ http://spatial-analyst.net/book/sites/default/files/Hengl_2009_GEOSTATe2c1w.pdf
fuente
Recuerdo haber usado SAGA para hacer esto hace unos años para obtener un resultado de modelado de inundación. Código abierto y vale la pena echarle un vistazo.
fuente
gvSIG (otro SIG gratuito) permite kriging, utilizando Sextante. Esto es básicamente lo mismo que usar SAGA, pero gvSIG proporciona una experiencia gis más 'típica' (es decir, similar a ESRI).
fuente
Puede probar el modelo Kriging en Surfpack versión 1.1 (lo escribí mientras todavía estaba en el equipo de DAKOTA), o la última y mejor versión que viene con la versión "estable" de DAKOTA (Surfpack es un sub-paquete de DAKOTA) , realiza Kriging universal desde la perspectiva de las funciones de correlación en lugar de semi-variogramas.
Recientemente, un usuario, Joel Guerrero, lo comparó frente a un montón de otras implementaciones y declaró que "Siempre relacionado con surfpack, lo estamos comparando con otras implementaciones (incluida una comercial), y hasta ahora las ha superado a todas, hasta el punto de que a veces parece que está haciendo magia negra "
fuente
GSLIB (Geostatistical Software Library) es un software de primera clase basado en archivos / comandos desarrollado por la Universidad de Stanford y lanzado en la década de 1990, con algunos mantenimientos en la última década. El código fuente se puede descargar y compilar libremente en Linux / Windows utilizando un compilador Fortran. Hay recursos en línea y un libro disponible.
Kriging es uno de los puntos fuertes del software:
fuente