我正在尝试使用 automap 包中的 autoKrige() 函数来简单地应用通用克里金法。我有一个不规则间隔的测量网格，我想在它们之间以精细的空间尺度进行插值。示例代码：

这会导致错误：

我尝试让 new_data 与原始数据具有相同的行数，甚至尝试使 new_data 中的坐标与原始数据完全相同，但我仍然收到此错误。我是地质统计学技术的新手，如果我犯了一个基本错误，我深表歉意。谁能告诉我哪里出错了？谢谢。

我的数据

我从一个虚拟字段中取样，现在想要预测相邻值。它是通过下面描述的两种方法完成的……都为每个点预测一个常数值。我不明白为什么。

请记住，我是新的 R 用户并且喜欢它。我一直在尝试按照我在网上找到的示例来执行我​​的任务。使用示例数据集，我可以毫无问题地实现他们的结果。但是，我无法将其应用于我自己的数据。请参阅此处的示例：http: //www.stat.ucla.edu/~nchristo/statistics_c173_c273/c173c273_lec11_w11.pdf

所有 35 个结果都返回 44.0833 #我就是不明白！！！

当有这么多零时，我不明白如何在字段中返回单个值。我希望看到单个区域的值递减。

想法？

我一直在使用包中的krigeST函数gstat，并尝试使用全套基于 SpatialPolygons 的数据 (STFDF) 进行本地克里金法。

当我尝试以下

问题似乎是在 krigeST.local 中，some_STFDF_with_SpatialPolygons_in_sp被强制转换为SFIDF. 在那里，有一个班级不允许的子集步骤SpatialPolygons，因为 中有重复的术语from@index[, 1]。

我想尝试不可分离的时空变异函数，因此不使用局部克里金法并不是一个真正的选择。有解决此问题的方法吗？

谢谢。

我正在尝试使用 r gstat 包使用回归克里金法拟合模型。我已经加载了一些称为流行率的疾病数据和一些协变量降雨量。降雨栅格数据作为 .grd 文件加载。我用来拟合此模型的代码是

我得到的错误结果是

有人可以帮忙吗？

我有一系列坐标表示二维标量字段中的数据：x[i]、y[i]、F[i]。但是，x[i] 和 y[i] 不在直线网格上，它们是“随机”放置的。我想对 f[i] 进行插值，这样我就可以在 x、y 中位置的直线网格上估计它的值。

这个问题应该是地形的相同高度插值（参见http://www.cs.uu.nl/geobook/interpolation.pdf）。例如，如果有一个地形的高度图，它不一定是在直线网格上测量的。但是，许多地图更喜欢直线网格上的高度。因此，必须将数据内插到直线网格中。

我使用 Matlab，它有一个griddata命令，可以对直线网格执行插值。该方法使用 Delaunay 三角剖分和一些未描述的插值。我可以猜到它是如何工作的，但在我看来，这种类型的代码应该已经由其他人编写了。

请参阅： http ://www.mathworks.com/help/matlab/ref/griddata.html ?refresh=true 以获取有关 griddata 如何从“分散”数据产生直线插值的参考。

我想将代码集成到 C# 应用程序中，因此我正在寻找要合并的库或源代码。ALGLIB 将执行从均匀间隔到均匀间隔的三次样条插值。我找不到任何适用于非均匀间隔数据的东西。

我有一组规则间隔的天气数据，需要插入到空间表面中。我想知道可用于此插值过程的最佳空间预测方法和软件。谢谢

我目前正在编写一个多重梯度下降算法，我使用克里格函数。我的问题是我找不到如何获得克里格函数的梯度（我尝试使用线性化但我不知道如何使其工作）。

我正在尝试使用普通克里金法使用 R 中的 gstat 或 automap 包根据预测变量在空间上预测动物将出现的数据。我有许多（超过 100 个）重复的坐标点，因为这些站点是采样的，所以我无法丢弃它们多年来多次。每次我为普通克里金运行下面的代码时，我都会收到一个 LDL 错误，这是由于重复点造成的。有谁知道如何在不丢弃数据的情况下解决这个问题？我已经尝试了 automap 包中应该纠正重复项的代码，但我无法让它工作。感谢您的帮助！

我是通过 Matlab 进行图像处理的新用户。我的第一个目标是应用这篇文章并比较我的结果和作者的结果。
这篇文章可以在这里找到：http: //arxiv.org/ftp/arxiv/papers/1306/1306.0139.pdf

第一个问题，图像质量：在图 7 中，定义了掩码，但我无法访问掩码数据集，并且我使用屏幕截图，因此图像质量很低。在我看来，它可以影响结果。有什么建议吗？

第二个问题，合并图像：我想在 Lena 上应用蒙版 1。但我不想使用paint =)另一方面，是否可以合并图像并保留lena？

我有一个名为“seoul032823”的 81 次观测的每小时 PM10 数据集。您可以从这里下载。我在这个数据集上执行了普通的克里金法，还得到了用于克里金法预测的空间图。我还可以在国家地图上显示观察数据点。但我不能在国家地图上重叠克里金空间预测图。

我想做什么：我想在韩国地图（不是整个韩国）上重叠我的空间预测地图。我感兴趣的区域是纬度 37.2N 到 37.7N 和经度 126.6E 到 127.2E。这意味着我需要从韩国地图上裁剪这个区域并将预测地图重叠在上面。我还需要根据浓度值显示原始观察数据点，这些数据点将遵循空间图的颜色。例如，我想要这种类型的地图：

我的克里金 R 代码，并在韩国地图上显示数据点：

我已经使用automap包进行克里金法和ggplot2绘制韩国地图。