###Example 1: Ordinary Kriging
## creating random variables first
## here, a grid is chosen, but does not matter
step <- 0.25
x <- seq(0,7,step)
param <- c(0,1,0,1)
model <- "exponential"
RFparameters(PracticalRange=FALSE)
p <- 1:7
points <- as.matrix(expand.grid(p,p))
data <- GaussRF(points, grid=FALSE, model=model, param=param)
## visualise generated spatial data
zlim <- c(-2.6,2.6)
colour <- rainbow(100)
image(p, p, xlim=range(x), ylim=range(x),
matrix(data,ncol=length(p)),
col=colour,zlim=zlim)
## now: kriging
krige.method <- "O" ## ordinary kriging
z <- Kriging(krige.method=krige.method,
x=x, y=x, grid=TRUE,
model=model, param=param,
given=points, data=data)
image(x,x,z,col=colour,zlim=zlim)
Run the code above in your browser using DataLab