如何在" krige"中完成矩阵求逆R工具的gstat包的功能

时间:2015-03-05 10:27:07

标签: r matrix-inverse gstat

我正在了解R工具中的gstat包如何实现kriging方法。我已经理解了经验半变异函数和拟合半变异函数模型的计算。但我不明白它如何实现矩阵求逆来计算克里金估计量的权重。我有一个包含50000 lat-long-precipitaion三胞胎的大型数据集。必须进行大小为50000x50000的矩阵的理论反演才能获得权重。虽然这个大矩阵需要几GB的人类记忆,这是特别不切实际的。

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

此致

CHANDAN

1 个答案:

答案 0 :(得分:1)

你没有告诉你的计算环境是什么,但我相信可以说它在一秒内没有解决50000点的克里金问题。为了理解它的作用,请提供更多信息,例如:您使用的命令和输出gstat给出的。