## import the data luisData <- dget("luisData") olsFn <- dget("olsFn") lociLuis1 <- dget("lociLuis1") lociLuis <- matrix(c(-5.58,61.29,-5.92,62.04,-7.08,61.63),byrow=T,ncol=2) ## kriging luisKc <- krige.conv(luisData , locations = lociLuis , krige = krige.control(obj.m = ols.Fn)) ## plot kriging image(luisKc, loc = lociLuis, col = grey(seq(1, 0.1, l = 30)))