1

我希望在 2 维和 3 维域中进行克里格（分别）。我有很多观点，已在下面复制。

笔记：

我希望克里格的属性是 Y。这些是故意为负的，因为它们是指数值。
该域是 3D 的，因此从 xy 角度来看，一些点似乎是并置的。

# define kriging domain 
x <- sample(seq(0,300,by=5),10,replace=FALSE)
y <- sample(seq(0,300,by=5),10,replace=FALSE)
z <- sample(seq(0,30,by=0.5),100,replace=TRUE)
Y <- sample(seq(-3.0,-0.01,by=0.001),100*10*10,replace=TRUE)
(
mydata = data.frame(x=rep(x,each=10),y=rep(y,times=10),z=z,Y=Y)

grid = data.frame(x=rep(x,times=length(x)),y=rep(y,each=length(y)))

# 2_D Kriging
ok_2D <- krige(formula=Y~1, locations=~x+y, data=mydata, newdata=grid, model=fit.exp)

# 3_D Kriging
ok_3D <- krige(formula=Y~1, locations=~x+y+z, data=mydata, newdata=grid, model=fit.exp)

错误：

> # 2_D Kriging
> ok_2D <- krige(formula=Y~1, locations=~x+y, data=mydata, newdata=grid, model=fit.exp)
[using ordinary kriging]

"chfactor.c", line 131: singular matrix in function LDLfactor()
Error in predict.gstat(g, newdata = newdata, block = block, nsim = nsim,  : 
  LDLfactor
> 
> # 3_D Kriging
> ok_3D <- krige(formula=Y~1, locations=~x+y+z, data=mydata, newdata=grid, model=fit.exp)
Error in .subset(x, j) : only 0's may be mixed with negative subscripts
>

我的问题是：为什么奇异矩阵有问题？我的理解是，并置点在克里金法中不应该成为问题，因为权重将为其中一个并置点分配零权重。

或者，我还缺少另一个问题吗？谢谢。

4

0 回答 0