Home » Tags » Figures
warning: Parameter 2 to gmap_gmap() expected to be a reference, value given in /home/spatiala/public_html/book/includes/module.inc on line 497.

Fig. 2.6: Comparison of ordinary kriging and regression-kriging using a simple example with 5 points

T.Hengl - Tue, 12/08/2009 - 14:49
Tags:
Fig. 2.6: Comparison of ordinary kriging and regression-kriging using a simple example with 5 points


library(maptools)
library(gstat)
points5 <- as.data.frame(matrix(c(2,3,9,6,5,2,7,9,5,3,3,4,2,4,6),nrow=5,ncol=3,dimnames=list(c("p1","p2","p3","p4","p5"),c("X","Y","Z"))))
coordinates(points5) <- ~X+Y
points5.plt <- bubble(points5, "Z", fill=F, scales=list(draw=T))
print(points5.plt)
grid10 = expand.grid(x=seq(0,10,10), y=seq(0,10,10))
gridded(grid10) <- ~x+y
newpoint <- as.data.frame(matrix(c(5,5),nrow=1,ncol=2,dimnames=list(c("x1"),c("X","Y"))))
coordinates(newpoint) <- ~X+Y
krige(Z~1, points5, newpoint, vgm(nugget=2.5,"Sph",psill=7.5,range=10))