Home » Tags » Figures

Fig. 5.13: Four simulations of liming requirements (indicator variable) using ordinary kriging.

T.Hengl - Tue, 12/08/2009 - 17:57
Tags:
Fig. 5.13: Four simulations of liming requirements (indicator variable) using ordinary kriging.
# geostatistical simulations (lime)
lime.ovgm <- fit.variogram(variogram(I(lime == 1) ~ 1, meuse), vgm(1, "Sph", 800, 1))
plot(variogram(I(lime == 1) ~ 1, meuse), lime.ovgm)
lime.sim <- krige(I(lime == 1) ~ 1, meuse, meuse.grid, lime.ovgm, nsim=4, indicators = TRUE, nmax=40)
spplot(lime.sim, col.regions=grey(c(rev(seq(0,0.9,0.05)), 0)), main="OK simulations of lime")