R工具的gstat包的"krige"函数是如何进行矩阵求逆的

Question

我正在了解R工具中的gstat包如何实现克里金法。我已经了解经验半变异函数和拟合半变异函数模型的计算。但是我不明白它是如何实现矩阵求逆来计算克里金估计量的权重的。我有一个包含 50000 个经纬度降水三元组的大型数据集。理论上必须对大小为 50000x50000 的矩阵进行求逆才能获得权重。而这个大矩阵需要几个GB的人内存，这特别不切实际。

我的问题是 krige 函数如何在一秒钟内完成所有这些操作？

此致，

钱丹

Answer 1

你没有说出你的计算环境是什么，但我相信可以肯定地说它没有在一秒钟内解决 50000 点的克里金问题。为了了解它做了什么，请提供更多信息，例如您使用的命令以及 gstat 给出的输出。