Displaying 1 result from an estimated 1 matches for "rsa2_grd".
2011 Jan 05
1
Prediction error for Ordinary Kriging
...l +units=m +no_defs"))
x2 <- readShapeSpatial("ptna2", CRS("+proj=tmerc +lat_0=39.66666666666666
+lon_0=-8.131906111111112 +k=1 +x_0=0 +y_0=0 +ellps=intl +units=m +no_defs"))
bb <- bbox(rsa2)
cs <- c(1, 1)
cc <- bb[, 1] + (cs/2)
cd <- ceiling(diff(t(bb))/cs)
rsa2_grd <- GridTopology(cellcentre.offset = cc,cellsize = cs, cells.dim = cd)
getClass("SpatialGrid")
p4s <- CRS(proj4string(rsa2))
x2_SG <- SpatialGrid(rsa2_grd, proj4string = p4s)
x2_SP <- SpatialPoints(cbind(x2$X, x2$Y))
v <- variogram(log1p(tsport_ace) ~ 1, x2, cutoff=100, widt...