我有一些块克里金问题。首先我创建一个网格（5000 * 5000米）。比我排除彼此非常接近的电台。下一步是创建变差函数和fit.variogram。 我指定投影和坐标。之后，krige函数正在为ordninary kriging工作。
x_2 <-as.integer(range(data$x))
y_2 <-as.integer(range(data$y))
grid_1 <-expand.grid(x=seq(from=x_2[1],to=x_2[2], by=5000),y=seq(y_2[1],to=y_2[2],by=5000))
coordinates(grid_1) <- ~x+y #Koordinaten zuweisen, Variablennamen festlegen, auf die man sich bezieht)
gridded(grid_1) =TRUE
grid.stats <- data.frame(data)[,1:2]
coordinates(grid.stats) <- ~x+y
ind.zd <- zerodist(grid.stats)
data.new <- data[-ind.zd,]
p.sill <- var(log(data.new$variable+0.1)) #0.4005127
nugget <- var(log(data.new$variable+0.1))/2 #0.2002564
vgm_data <- variogram(log(variable+0.1)~1, ~x+y, data.new)
m.fit <- fit.variogram(vgm_data, vgm(0.4, "Sph", 1000000, 0.2))
coordinates(data.new) <- ~x+y
krig_1 <- krige(log(variable+0.1)~x+y, data.new, grid_1, nmax=7, nmin=4, m=m.fit)
krig_block <- krige(log(variable+0.1)~1, data.new, grid_1, m=m.fit, block =c(5000/5000))
我尝试为块克里金添加“块”，但我没有得到任何结果。有谁知道我的问题？任何帮助都会非常有用，提前谢谢。 亲切的问候