Kriging

Can someone kindly help me about kriging 3D interpolation in MATLAB? Is there an efficient code for such an interpolation? It should be much faster than TriScatterInterp. I am wondering why Mathworks does not include it in its program.

[Morteza Kiani]

Do you know how the prediction can be write out as matlab function in dace or kriging? Thanks

Answer 1

[Wolfgang Schwanghart]

To my knowlegde there is no comprehensive geostatistical software written in Matlab. I would always use packages such gstat or gslib. There are, however, some toolboxes that interface with gstat, e.g. mgstat http://mgstat.sourceforge.net/

You can also adapt my code http://www.mathworks.com/matlabcentral/fileexchange/29025-ordinary-kriging to work in 3d. If you work on large datasets, however, it will be very slow.

So, don't stick to Matlab if you want to work with geostatistics. There are much better solutions around.

HTH, Wolfgang

Answer 2

[Sean de Wolski]

You may wish to look at some of the results for kriging on the FEX. (I've never used any, but a few appear to be well rated)

http://www.mathworks.com/matlabcentral/fileexchange/?term=kriging

[AP]

Are you aware of any 3D interpolation by Kriging? Unfortunately, all those are 2D :(

Answer 3

[Thomas Clark]

This works in ND... http://www2.imm.dtu.dk/~hbni/dace/dace.pdf

[AP]

Thanks Thomas. Unfortunately, I had to write my own code. Hopefully, other users will benefit from the package you cited.