geoR克里金法错误: solve.default

Question

我正在尝试在geoR中为相当大的区域(约100万千米^2)做克里金法。这是为了我的论文，所以不幸的是我不能分享数据。我已经检查了坐标和实际数据中的重复项，看看这是否是求解中出现错误的原因。我已经尝试了不同的方式来写趋势，例如手写或“第二”。我见过类似的错误，但它们都是回归建模的例子。我的"test.geo“数据与坐标和数据在一起是有序的，并且我的变异函数工作。所以我知道.geo数据的初始部分是起作用的，我模型的克里金法部分就是不想起作用。

我的代码是：

    my.bdr <- map.shp2@polygons[[1]]@Polygons[[1]]@coords
    x <- my.bdr[,1]
    y <- my.bdr[,2]
    krige.grid <- expand.grid(seq(min(x), max(x), l=1000), 
                      seq(min(y), max(y), l=1000))
    test.ply <- as.matrix(cbind(x, y))   
    test <- krige.conv(test.geo, 
                 krige = krige.control(
                   type.krige = "OK",           
                   trend.d = trend.spatial(~ 1 + test.geo$coords[,1] + test.geo$coords[,2] 
                                           + I(test.geo$coords[,1]* test.geo$coords[,1]) 
                                           + I(test.geo$coords[,1]*test.geo$coords[,2])
                                           + I(test.geo$coords[,2]*test.geo$coords[,2])),    
                   trend.l = trend.spatial(~ 1 + krige.grid[,1] + krige.grid[,2] 
                                           + I(krige.grid[,1]* krige.grid[,1]) 
                                           + I(krige.grid[,1]*krige.grid[,2])
                                           + I(krige.grid[,2]*krige.grid[,2])), 
                   obj.model = test.svm),
                 locations = krige.grid, borders=test.ply)

我一直收到这个错误：

    Error in solve.default(ttivtt, crossprod(ivtt, as.vector(data))) : 
    system is computationally singular: reciprocal condition number = 4.81502e-36

我会感谢任何建议或帮助，以解决这个问题！谢谢!

Answer 1

这是一个古老的话题，但也许有人会发现它很有用。您使用的是什么CRS？如果您使用的是UTM，请尝试将coordenates除以1000。这是包的创建者提出的解决方案。它适用于我的线性趋势(而不是二次多项式趋势)