[Gstat-info] conditional sequential simulation, and value honoring
giohappy at gmail.com
Tue Apr 8 13:10:56 CEST 2008
I've done some work with krige() to run conditional simulations on my
datas. It's a simulation grid with 6458 nodes on 50x50 meters cells,
and 198 observed values (sparsed all over the grid). I run the
simulation with nsim=500, nmax=100 and range=700, choosed on the base
my variogram model range.
At the end of the 500 simulations of done the mean for every cell, and
then I've overlayed the values in the 198 locations. If I compare the
observed values and the values retrieved from the cells, I obtain very
different values. I.e. on a scale of 0-70 dS/m, if the observed value
is 52, I may obtain a simulated value of 15 in the cell it falls in,
even if it's the only observed value ina range of 200-300 meters.
Shouldn't it weigh much more in the simulated value? I would expect
something that honor more the true value...
Am I wrong in expecting this from the conditional simulation algorithm?
More information about the Gstat-info