ossfim: Kriging standard errors as function of grid spacing and block...

ossfimR Documentation

Kriging standard errors as function of grid spacing and block size


Calculate, for a given variogram model, ordinary block kriging standard errors as a function of sampling spaces and block sizes


ossfim(spacings = 1:5, block.sizes = 1:5, model, nmax = 25, debug = 0)



range of grid (data) spacings to be used


range of block sizes to be used


variogram model, output of vgm


set the kriging neighbourhood size


debug level; set to 32 to see a lot of output


data frame with columns spacing (the grid spacing), block.size (the block size), and kriging.se (block kriging standard error)


The idea is old, simple, but still of value. If you want to map a variable with a given accuracy, you will have to sample it. Suppose the variogram of the variable is known. Given a regular sampling scheme, the kriging standard error decreases when either (i) the data spacing is smaller, or (ii) predictions are made for larger blocks. This function helps quantifying this relationship. Ossfim probably refers to “optimal sampling scheme for isarithmic mapping”.


Edzer Pebesma


Burrough, P.A., R.A. McDonnell (1999) Principles of Geographical Information Systems. Oxford University Press (e.g., figure 10.11 on page 261)

Burgess, T.M., R. Webster, A.B. McBratney (1981) Optimal interpolation and isarithmic mapping of soil properties. IV Sampling strategy. The journal of soil science 32(4), 643-660.

McBratney, A.B., R. Webster (1981) The design of optimal sampling schemes for local estimation and mapping of regionalized variables: 2 program and examples. Computers and Geosciences 7: 335-365.

See Also



## Not run: 
x <- ossfim(1:15,1:15, model = vgm(1,"Exp",15))
levelplot(kriging.se~spacing+block.size, x, 
  main = "Ossfim results, variogram 1 Exp(15)")

## End(Not run)
# if you wonder about the decrease in the upper left corner of the graph,
# try the above with nmax set to 100, or perhaps 200.

gstat documentation built on April 6, 2023, 5:21 p.m.