Nothing
## area2point in space and time
library(sp)
library(gstat)
data("meuse")
coordinates(meuse) <- ~x+y
data("meuse.grid")
coordinates(meuse.grid) <- ~x+y
gridded(meuse.grid) <- TRUE
meuse.coarse.grid <- SpatialGrid(GridTopology(meuse.grid@grid@cellcentre.offset, c(600,600), c(5, 7)))
separableModel <- vgmST("separable",
space=vgm(0.85,"Exp", 831, 0.15),
time =vgm(0.9,"Exp", 3.25, 0.1),
sill=135000)
attr(separableModel,"temporal unit") <- "days"
library(spacetime)
stf <- STF(meuse[sample(155, 25),], Sys.time()-2:0*24*3600)
stf_grid <- STF(geometry(meuse.coarse.grid), stf@time)
krigedSim <- krigeSTSimTB(newdata = stf_grid, modelList = separableModel, nsim = 1, nLyrs = 100)
# area-to-point kriging:
a2pST = krigeST(sim1 ~ 1, krigedSim, stf, modelList = vgmAreaST, ndiscr = 9,
model = separableModel, # point variogram,
verbose = TRUE)
p1 <- stplot(krigedSim, color.key=F)
p2 <- stplot(a2pST, color.key=F)
print(p1, position=c(0,0.5,1,1), more=TRUE)
print(p2, position=c(0,0,1,0.5))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.