Rank: 178213 based on 0 downloads (last 30 days) and 0 file submitted
photo

Sébastien

E-mail

Personal Profile:

 

Watch this Author's files

 

Comments and Ratings by Sébastien
Updated File Comments Rating
03 Apr 2013 Ordinary Kriging 2D-interpolation using geostatistics Author: Wolfgang Schwanghart

Very nice submission. I was amazed by its speed.
However, two things to say:

1) there is a mistake in the code when expanding b with ones. It should be the following:
% expand b with ones
b = vstruct.func([vstruct.range vstruct.sill],b);
if ~isempty(vstruct.nugget)
b = b+vstruct.nugget;
end
b=[b;ones(1,chunksize)];
Indeed, the ones should be added AFTER adding the nugget. Otherwise, the nugget is also applied to the ones: an offset is then unfortunately observed on the results...

2) I have the same problem than Kirsten: it's working well with spherical variograms. However, when I apply a gaussian one, huge values are found. I'm actually investigating the problem. If I find something, I'll let you know. But I would appreciate if you find the solution first ;)

Contact us