Tell me more ×
Geographic Information Systems Stack Exchange is a question and answer site for cartographers, geographers and GIS professionals. It's 100% free, no registration required.
up vote 26 down vote favorite
14

I have a point dataset which I'd like to Krige, ideally using an open-source software package. If possible, I'd also like to choose the semi-variogram model during the process to improve the estimation.

share|improve this question

7 Answers

up vote 18 down vote accepted

Depending on which Kriging type you want to apply, there are different packages to choose from:

Ordinary Kriging

The most common version is implemented for example in:

Simple Kriging

Simple Kriging uses the average of the entire data set while Ordinary Kriging uses a local average. Therefore, Simple Kriging can be less accurate, but it generally produces "smoother" results. It's implemented in:

Universal Kriging

Universal Kriging allows for consideration of drift in data. Implementations are included in:

Other Kriging Types

GRASS v.krige and QGIS Kriging as part of the SDA4PP plugin also support Block Kriging.

HPGL implements a big number of less known Kriging methods (check the manual for more information on those):

  • Indicator Kriging (IK)
  • Local Varying Mean Kriging (LVM Kriging)
  • Simple CoKriging (Markov Models 1 & 2)
  • Sequential Indicator Simulation (SIS)
  • Corellogram Local Varying Mean SIS (CLVM SIS)
  • Local Varying Mean SIS (LVM SIS)
  • Sequential Gaussian Simulation (SGS)
  • Truncated Gaussian Simulation (GTSIM) [in Python scripts collection]

SAGA offers different versions of both Ordinary and Universal Kriging.

Gstat krige additionally supports Block and Point Kriging.

share|improve this answer
Great Kriging answer! – Ragi Yaser Burhum Nov 22 '11 at 4:37
up vote 14 down vote

It looks like there are a few options with GRASS GIS. Check out the GRASS Kriging Wiki page: http://grass.osgeo.org/wiki/Kriging

A Google Summer of Code project in 2009 produced V.krige: http://grass.osgeo.org/wiki/V.krige_GSoC_2009

The GPL gstat package should work by itself or interfaced with GRASS GIS. http://www.gstat.org/

Dylan Beaudette has a nice example of doing kriging with GRASS. http://casoilresource.lawr.ucdavis.edu/drupal/node/438 (His blog is full of great and interesting examples of using OpenSource GIS and statistical tools!)

share|improve this answer
up vote 11 down vote

If you are happy to read your raster into a numpy array (gdal can do this), then you could use the High Performance Geostatistics Library implementation from Python or C/C++.

HPGL implements the following algorithms:

  1. Simple Kriging (SK)
  2. Ordinary Kriging (OK)
  3. Indicator Kriging (IK)
  4. Local Varying Mean Kriging (LVM Kriging)
  5. Simple CoKriging (Markov Models 1 & 2)
  6. Sequential Indicator Simulation (SIS)
  7. Corellogram Local Varying Mean SIS (CLVM SIS)
  8. Local Varying Mean SIS (LVM SIS)
  9. Sequential Gaussian Simulation (SGS)
  10. Truncated Gaussian Simulation (GTSIM) [in Python scripts collection]

I haven't used it myself but have heard good things about it, especially with respect to speed.

share|improve this answer
up vote 10 down vote

The R-project has substantial number of spatial statistics software packages, but R has rather steep learning curve.

share|improve this answer
People always say that, but I wonder: steep relative to what? – Matt Parker Aug 11 '10 at 19:20
I've seen the "steep learning curve" comment thrown at R a few times- it just doesn't make sense to me. I was a year into my relationship with MATLAB when I discovered R. I found R so easy to learn that I gave MATLAB the one finger salute and promptly quit using it heavily. – Sharpie Aug 12 '10 at 8:33
i think it's because people rarely try to understand statistics, and because of that it has a steep learning curve syntax wise, there are rarely problems picking it up – dodobas Aug 12 '10 at 8:51
1  
I think syntax wise it is one of the easier languages to learn. What is an example of a statistical language that is easy to learn from the command line. I think people complain because it is not Excel. – TheSteve0 Aug 14 '10 at 4:39
It is step in comparison with a GUI based program. If you used Windows all your life, and GUI based programs, you will run when you see the command line look. The Excel comparison indeed make them run. But R is very simple to use if someone can show you the basic tricks. You must be prepared to learn new concepts like vectors, matrices, functions, loops, which in an Excel/Windows world does not exist. If you previously used Linux, it wouldn't be a step curve. – Niculita Mihai Feb 1 at 12:27
up vote 5 down vote

Check this free book, it's about doing geostatistics in R, and contains some info on doing it in SAGA and GRASS as well. http://spatial-analyst.net/book/ http://spatial-analyst.net/book/sites/default/files/Hengl_2009_GEOSTATe2c1w.pdf

share|improve this answer
up vote 3 down vote

I remember using SAGA to do this a few years back for some flood modelling output. Open Source and well worth a look.

share|improve this answer
up vote 2 down vote

gvSIG (another free GIS) does allow kriging, using Sextante. This is basically the same as using SAGA, but gvSIG provides a more 'typical' (i.e. ESRI-like) gis experience.

share|improve this answer

Your Answer

 
discard

By posting your answer, you agree to the privacy policy and terms of service.

Not the answer you're looking for? Browse other questions tagged or ask your own question.