问题标签 [kriging]

同时试图从克里格苯值

我像这样生成一个网格（这些井所在的属性）

看看变异函数..（不是很漂亮，但正在努力）

使用 ksline 生成克里格值...

但后来它拿着手机，当我尝试对结果进行成像时出错！！！！

这似乎完全超出了左侧字段，因为我已经多次检查过......我用谷歌搜索了错误，似乎我正在尝试调用函数的子集，而 90% 的答案是我的语法某处是错误的，但我已经检查了所有内容，但我无法弄清楚......任何帮助将不胜感激。

谢谢

ZR

我正在尝试应用克里金法来插值空气污染浓度（目标变量）当我运行如下克里格函数时，R 返回错误。

RSPAVE 是目标变量；air 是包含 RSPAVE 的数据集；TPU 是 shapefile

这可能是因为我的网格数据是一个 shapefile。但我不想阻止克里金法，如何将多边形转为点，并应用普通克里金法？

非常感谢！

我是 R 的新手，我想测试不同的参数来执行反距离加权 (IDW) 插值。

我想多次重复以下功能，只需更改变量的n值（比如说 from1到20）nmax并分别存储结果以执行模型敏感性分析。

我想我需要某种基本的循环。有人能帮我吗？

非常感谢！

我遇到了一个问题，我的特定输入会产生错误。

老实说，我不知道出了什么问题。我的格式似乎与我看到的示例代码非常一致。

以下代码用于使用 gstat 包在 R 中使用克里金法预测三个位置的 v 值。

我想知道为什么为了拟合变异函数模型，我们需要预先使用 vgm() 方法提供窗台、块金和范围值。根据克里金理论，我认为我们必须通过最小化 WLS 目标函数来计算这些值。

〜问候，钱丹

以下 R 程序使用 gstat 包中的 walker Lake 数据使用 470 个数据点创建插值曲面。

该程序为 2688 个点中的每个点计算以下内容

gstat 包是否正在使用一些聪明的方法进行计算。我从之前的 stackoverflow 查询中知道它使用 cholesky 分解进行矩阵求逆。一台机器算这么快算正常速度吗。

我想对二进制存在-不存在和主机网格数据执行回归克里金法 (RK) 作为常数预测器。我已经使用逻辑函数来估计二元结果和预测变量之间的关系，但是我认为它没有通过 RK 假设？预测变量在模型中并不显着。有没有其他方法可以解决它？

代码数据：https ://drive.google.com/folderview?id=0B7-8DA0HVZqDYk1BcFFwSkZCcjQ&usp=sharing

然后，当涉及到实际的克里金法时，我已经从教程中构建了这段代码，但似乎 glm 的实际残差没有输入到krige函数中。可以在gstat中改进吗？

有没有一种（简单的）方法可以将 SpatialPixelsDataFrame（来自 krige）转换为例如 SpatialPolygonsDataFrame（矢量图形而不是像素图）。

设置值范围并将栅格插入多边形或使用另一种生成 SpatialPolygonsDataFrame 的 krige 方法会很好。我期待一个简单的例子。

krige 示例：即书中https://oscarperpinan.github.io/spacetime-vis/ 8.1.5 Spatial Interpolation，完整源码https://github.com/oscarperpinan/spacetime-vis/blob/master/bubble.R

我想知道：克里金法时是否有必要重新定义自己的列？下面的错误似乎表明了这一点：

有我没有看到的问题吗？或者，我需要定义自己的列吗？谢谢。

以下程序从这里开始是完全可重现和可运行的：

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

笔记：

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

错误：

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

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