问题标签 [automap]

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

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

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

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

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

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

解决方案

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

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

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

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

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

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

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

我的数据如下所示：

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

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

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

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

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

我正在使用包中的 R 函数autokrige，automap但出现错误，我不知道如何解决。你有什么提示吗？

谢谢！

我正在尝试在 R 中编写使用 gstat 库的代码来创建插值。我已经阅读了 gstat 手册，并根据互联网上的一些示例，我设法编写了这段代码（这只是一部分）：

我的问题是，当我运行最后一个命令（预测）而不是使用普通克里金插值得到结果时，我得到的是反距离加权（IDW）。我在 gstat 手册中读到：“当没有指定变异函数时，反距离加权插值是默认操作。指定变异函数时，默认预测方法是普通克里金法。”

但是，正如您在我的代码中看到的那样，我指定了经验和理论变异函数。你知道为什么我一直使用 IDW 而不是普通的克里金法吗？它可以与我拥有的坐标类型有关吗？例如，如果我的坐标彼此接近，或者感兴趣的区域太大？任何帮助都会非常有用。

提前致谢 迪米特里斯

我是 R 的新手。我正在使用autofitVariogram来处理 50 个站点的每日降雨数据。示例数据如下所示。一些站点的缺失值由“NaN”值表示。

我的问题是关于 variogramfit。变异函数仅覆盖 60,000m 的距离。为什么没有绘制超过 60Km 的 bin 中的点。我从空间相关图中看到，与经纬度信息的最大距离> 200Km。

下面提供了纬度和经度信息的摘要。摘要（lonlat） lon lat
Min。:74.78 分钟 :15.77
1st Qu.:75.14 1st Qu.:16.04
Median :75.56 Median :16.33
Mean :75.54 Mean :16.37
第三Qu.:75.94 第三Qu.:16.66
Max. :76.31 最大。:17.23

Q2）如何将变异函数拟合 png 文件保存在一个循环中。我知道每次保存图形后都应该使用 dev.off()，但我无法保存图形。任何帮助，将不胜感激。

谢谢，

任何建议，将不胜感激？

我想对 geoR 包中的 ca20-Dataset 进行交叉验证。例如，使用 meuse 数据集，这工作正常，但对于这个数据集，我遇到了 SpatialPointsDataFrame 尺寸的奇怪问题。也许您可以自己尝试一下并解释为什么 autoKrige.cv 函数不起作用（我尝试了几个 nfold-values 但这只会更改错误消息的位置值...）：

我希望有人能重现这个问题，我的 R 版本是 2.15，所有包都是最新的（至少不超过一个月左右......）。

谢谢你的帮助！！

我正在使用 automap 库，我想在其中使用 autokrige 插值降雨数据。生成的插值图与我的预测点图具有不同的维度（不同的行数和列数）。这怎么可能？请在下面查看我的代码：

我的 autoKrig 函数有问题，我尝试在这里制作一个可重现的示例：

这将导致：

我的原始数据出现同样的错误。我很感激任何帮助。

我在 Windows 7 下将 sp 和 automap 库安装到我的 R 3.0.1 64 位（通过 install.packages 命令）。它们的安装没有显示任何错误，并且 library(sp) 工作正常但是当我尝试执行 library(automap) 时，我收到以下错误：

我从路径中查看，该路径中确实没有DESCRIPTION文件（或文件夹）。但是只有 libs 文件夹在 x64 文件夹下，其中包含文件 sp.dll

知道是什么原因造成的吗？