r - R工具gstat包的“krige”函数中如何进行矩阵求逆

Question

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

我的问题是，克里格函数如何在一秒钟内完成所有这些工作？

问候，

钱丹

Answer 1

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