6 votos

Realización de Kriging ordinario en SAGA mediante el complemento SEXTANTE de QGIS

Estoy intentando aprender a aplicar el análisis de kriging en QGIS utilizando el plug-in SEXTANTE y las herramientas SAGA. Parece que tengo mis carpetas SEXTANTE y SAGA configuradas correctamente según la guía del usuario. Nunca he utilizado una función de kriging antes y no puedo encontrar un buen recurso de SAGA que explique los diferentes parámetros de entrada. ¿Alguien tiene experiencia en el uso de la herramienta de kriging ordinaria en SAGA que pueda dar consejos o una referencia? Actualmente estoy recibiendo el siguiente mensaje de error en QGIS:

An error has occured while executing Python code: 
Traceback (most recent call last):   
File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\gui\ParametersDialog.py", line 137, in accept if self.setParamValues():   
File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\gui\ParametersDialog.py", line 91, in setParamValue if not self.setParamValue(param, self.paramTable.valueItems[param.name]):   
File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\gui\ParametersDialog.py", line 129, in setParamValue return param.setValue(widget.getValue())   
File "C:/Documents and Settings/PS/.qgis/python/plugins\sextante\parameters\ParameterExtent.py", line 11, in setValue if value is None: NameError: global name 'value' is not defined

6voto

Mohit Jain Puntos 145

Si nunca ha hecho kriging antes, debe entender lo que está haciendo. Si no es el caso, consiga un libro de texto sobre geoestadística. De todos modos, la mayoría de los parámetros están ahí para determinar el variograma:

johan@cdh7:~$ saga_cmd libgeostatistics_kriging 5
_____________________________________________
  #####   ##   #####    ##
 ###     ###  ##       ###
  ###   # ## ##  #### # ##
   ### ##### ##    # #####
##### #   ##  ##### #   ##
_____________________________________________

library path:   /usr/lib/saga
library name:   libgeostatistics_kriging
module name :   Ordinary Kriging
author      :   (c) 2008 by O.Conrad
_____________________________________________
go...
Usage: 5 [-GRID <str>] [-VARIANCE <str>] -SHAPES <str> [-FIELD <str>] [-BVARIANCE] [-TARGET <str>] [-MODEL <str>] [-BLOCK] [-DBLOCK <str>] [-BLOG] [-NUGGET <str>] [-SILL <str>] [-RANGE <str>] [-LIN_B <str>] [-EXP_B <str>] [-POW_A <str>] [-POW_B <str>] [-MAXRADIUS <str>] [-NPOINTS_MIN <str>] [-NPOINTS_MAX <str>] [-USER_CELL_SIZE <str>] [-USER_FIT_EXTENT] [-USER_X_EXTENT_MIN <str>] [-USER_X_EXTENT_MAX <str>] [-USER_Y_EXTENT_MIN <str>] [-USER_Y_EXTENT_MAX <str>] [-SYSTEM_SYSTEM_NX <num>] [-SYSTEM_SYSTEM_NY <num>] [-SYSTEM_SYSTEM_X <str>] [-SYSTEM_SYSTEM_Y <str>] [-SYSTEM_SYSTEM_D <str>] [-GRID_GRID <str>] [-GRID_VARIANCE <str>]
  -GRID:<str>               Grid
    Data Object (optional output)
  -VARIANCE:<str>           Variance
    Data Object (optional output)
  -SHAPES:<str>             Points
    Shapes (input)
  -FIELD:<str>              Attribute
    Table field
  -BVARIANCE                Create Variance Grid
    Boolean
  -TARGET:<str>             Target Grid
    Choice
    Available Choices:
    [0] user defined
    [1] grid system
    [2] grid
  -MODEL:<str>              Variogram Model
    Choice
    Available Choices:
    [0] Spherical Model
    [1] Exponential Model
    [2] Gaussian Model
    [3] Linear Regression
    [4] Exponential Regression
    [5] Power Function Regression
  -BLOCK                    Block Kriging
    Boolean
  -DBLOCK:<str>             Block Size
    Floating point
    Minimum: 0.000000
  -BLOG                     Logarithmic Transformation
    Boolean
  -NUGGET:<str>             Nugget
    Floating point
    Minimum: 0.000000
  -SILL:<str>               Sill
    Floating point
    Minimum: 0.000000
  -RANGE:<str>              Range
    Floating point
    Minimum: 0.000000
  -LIN_B:<str>              Linear Regression
    Floating point
  -EXP_B:<str>              Exponential Regression
    Floating point
  -POW_A:<str>              Power Function - A
    Floating point
  -POW_B:<str>              Power Function - B
    Floating point
  -MAXRADIUS:<str>          Maximum Search Radius (map units)
    Floating point
    Minimum: 0.000000
  -NPOINTS_MIN:<str>        Min./Max. Number of m_Points
    Value range
  -NPOINTS_MAX:<str>        Min./Max. Number of m_Points
    Value range
  -USER_CELL_SIZE:<str>     Grid Size
    Floating point
    Minimum: 0.000000
  -USER_FIT_EXTENT          Fit Extent
    Boolean
  -USER_X_EXTENT_MIN:<str>  X-Extent
    Value range
  -USER_X_EXTENT_MAX:<str>  X-Extent
    Value range
  -USER_Y_EXTENT_MIN:<str>  Y-Extent
    Value range
  -USER_Y_EXTENT_MAX:<str>  Y-Extent
    Value range
  -SYSTEM_SYSTEM_NX:<num>   Grid System
    Grid system
  -SYSTEM_SYSTEM_NY:<num>   Grid System
    Grid system
  -SYSTEM_SYSTEM_X:<str>    Grid System
    Grid system
  -SYSTEM_SYSTEM_Y:<str>    Grid System
    Grid system
  -SYSTEM_SYSTEM_D:<str>    Grid System
    Grid system
  -GRID_GRID:<str>          Grid
    Grid (input)
  -GRID_VARIANCE:<str>      Variance
    Grid (optional input)

6 votos

+1 A mucha gente le gusta Isaaks y Srivastava . Si está utilizando el software de kriging sin tener los conocimientos contenidos en este libro de texto (a mínimo ), entonces no estás haciendo kriging, sólo te estás metiendo en problemas.

0 votos

Ja - Sí es cierto que me estoy metiendo en líos, pero al menos soy consciente de ello :) Voy a tomar un curso de postgrado en geoestadística el próximo otoño y vamos a leer "Geoestadística para la evaluación de los recursos naturales" por Goovaerts. Actualmente estoy explorando la posibilidad de utilizar QGIS para evaluar las plumas de agua subterránea en California para mi tesis a finales de año. Sin saber mucho sobre kriging, en este momento sólo estoy intentando un ensayo para ver si QGIS y el plugin SEXTANTE-SAGA es una opción viable. Así que con esa advertencia, si alguien sabe de algún tutorial para tontos como yo, por favor hágamelo saber. Gracias

2 votos

El libro de Goovaerts es definitivamente bueno, los cursos que hice se basaron en él :-). Definitivamente, yo investigaría primero la propia saga. Sus análisis serán mucho más rápido en saga que cuando se ejecuta a través de qgis/sextante. Pero en realidad, si es para una tesis, yo recomendaría saltar a R y gstat. Guía práctica de cartografía geoestática sería sin duda un buen comienzo (puedes descargarlo si tienes un presupuesto). [descargo de responsabilidad: utilicé R+saga para mi doctorado].

i-Ciencias.com

I-Ciencias es una comunidad de estudiantes y amantes de la ciencia en la que puedes resolver tus problemas y dudas.
Puedes consultar las preguntas de otros usuarios, hacer tus propias preguntas o resolver las de los demás.

Powered by:

X