问题标签 [spatial-interpolation]

我正在分析来自鸡蛋调查的数据。数据来自北海的不同地点，一些站点在不同日期记录了两次。大海应该被 0.5 x 0.5 度的正方形覆盖。我有两个问题我还没有找到任何解决方案：

1. 如何将具有重复位置和不同日期的点替换为平均值？我知道如何删除重复项或如何用最大值或最小值替换它们，但找不到计算平均值的方法。

2. 如何根据相邻单元格计算缺失点的插值。只有当至少两个记录点相邻时，才应计算插值。

我尝试设置一个网格，但并没有走得很远，因为我找不到一种方法来告诉 R 什么时候插值，什么时候不插值。

样本数据：

非常感谢！！

我对ggplot比较陌生，所以如果我的一些问题真的很简单或根本无法解决，请原谅我。

我想做的是生成一个国家的“热图”，其中形状的填充是连续的。此外，我有国家的形状.RData。我使用hadley wickham 的脚本将我的 SpatialPolygon 数据转换为数据框。我的数据框的 long 和 lat 数据现在看起来像这样

这个长/纬度数据勾勒出德国的轮廓。数据框的其余部分在此省略，因为我认为不需要。对于某些长/纬度点，我还有第二个值数据框。这看起来像这样

我现在想做的是根据位于该点一定距离内的所有固定点的平均值为地图的每个点着色。这样我就可以（几乎）连续地为整个国家的地图着色。到目前为止，我所拥有的是用 ggplot2 绘制的国家地图

我的第一个想法是生成位于已绘制地图内的点，然后my_generated_point像这样计算每个生成点的值

我还没有找到产生这些点的方法。和整个问题一样，我什至不知道是否可以通过这种方式解决。我现在的问题是是否存在产生这些点的方法和/或是否有另一种方法来解决问题。

解决方案

多亏了保罗，我几乎得到了我想要的。以下是荷兰的示例数据示例。

我有在美国切萨皮克湾不同地点采集的物种丰富度调查数据，我想以图形方式将数据呈现为“热图”。

我有一个纬度/经度坐标和丰富度值的数据框，我将其转换为 aSpatialPointsDataFrame并使用autoKrige()automap 包中的函数生成插值。

首先，任何人都可以评论我是否正确实现了该autoKrige()功能？

其次，我无法绘制数据并覆盖该地区的地图。或者，我可以指定插值网格以反映海湾的边界（如此处建议）？关于我如何做到这一点以及我可以从哪里获得这些信息的任何想法？提供网格autoKrige()看起来很容易。

编辑：感谢保罗的超级有用的帖子！这就是我现在所拥有的。无法让 ggplot 同时接受插值数据和地图投影：

所有空间插值的实际例子我都能够通过采样额外的周围点来估计导数来找到工作。但是，如果导数已知，并且您只需要已知点中心的单个点的值（和导数），是否有更简单的方法？

举例说明：假设对于每个点(1, 1)、(-1, 1)、(-1, -1)，并且(1, -1)您知道f(x, y)、f'(x)、f''(x)、f'(y)和- 并且您想要在for 、、、和处的f''(y)插值。(0, 0)f(x, y)f'(x)f''(x)f'(y)f''(y)

我正在尝试使用克里金函数在海洋横断面中创建化学参数的垂直图，但我很难开始。

我的数据如下所示：

包中的autoKrige函数automap看起来应该为我完成这项工作，但它需要一个 class 对象SpatialPointsDataFrame。据我所知，spTransform包中的函数rgdal创建了SpatialPointsDataFrame对象，但是有两个问题：

• CRAN不提供此 OSX 二进制文件，而且我在 OXS 10.7 上运行的 RStudio 副本似乎无法安装它，并且

• 此功能似乎适用于经纬度数据和地球曲率的正确距离值。由于我正在处理垂直平面（以及短距离，数百米的规模），我不想更正我的距离。

这里有一个关于 R 中克里金法的精彩讨论，但由于上面列出的问题，我不太明白如何将它应用于我的具体问题。

我想要一个矩阵或数据框来描述具有我的化学参数的插值的点网格，然后我可以绘制（理想情况下使用 ggplot2）。我怀疑我的问题的解决方案比我想象的要容易得多。

我在 data.txt 中有 2 列 xy 数据，如下所示：

现在我想定义一个函数 f(x)，其中 x 是第一列， f(x) 是第二列，然后能够像这样打印这个函数的值：

哪个应该给我4。

我如何实现这一目标？

我有许多具有每周分辨率的地理参考水文数据：

Station name, Lat, Long, Week 1 average, Week 2 average ... Week 52 average

不幸的是，我也有一些只有每月分辨率的数据：

Station name, Lat, Long, January average, February average ... December average

与其“重新发明轮子”，任何人都可以推荐一个最喜欢的模块、包或技术，以提供从月值中合理插值周值的方法吗？线性会很好，但如果我们可以使用坐标来改进基于附近站点的插值，那就太好了。

我用 python 标记了这篇文章，因为它是我最近一直在使用的语言（尽管不是它的统计函数）。如果答案是“使用类似的统计程序r”，那就这样吧，但我很好奇 python 有什么。谢谢！

我一直在尝试在我的世界地图上创建一个网格，如下所示 的网格地图示例

我希望能够根据文件中的 rate2 列为网格着色。我得到的大多数解决方案仅根据国家/地区为世界着色。请问我该怎么做。我附上了我的csv。

热图.csv

我是 R 的新手，并试图在 R 中进行插值（对不起，如果这似乎是一个蹩脚的问题）。我有一些点（Lat、Lan 和每个点的一个整数值）需要对它们进行插值。我是 GIS 用户，我是在 ARCGIS 中完成的，但据我所知，ArcGIS 没有给我样条曲线的 RMSE，所以我想在 R 中做并学习一些 R。我已经学习了基础知识，例如绘图、、、 ,, 你能告诉我我应该用 R 语言做什么来插值（全局级别）并计算 RMSE 吗？

我真的很感谢汉克的任何帮助

我有一个列表点及其覆盖区域的值，如下所示：

我想编写一个函数interpolateValueAtCoordinate(x, y)，它可以很好地插值特定坐标处的值。

我能做什么？