Autokriging空间数据

Question

我正尝试使用克立格函数在海洋横断面上创建化学参数的垂直地图，而且我很难开始使用。

我的数据是这样的：

horiz=rep(1:5, 5) 
depth=runif(25) 
value = horiz+runif(25)/5 
df <- data.frame(horiz, depth, value) 

在automap包中的autoKrige功能看起来像它应该为我做的工作，但它需要SpatialPointsDataFrame类的一个对象。据我所知，在包rgdal功能spTransform创建SpatialPointsDataFrame对象，但有两个问题：

• 这aren't available从CRAN的OSX的二进制文件，并运行我的RStudio的副本上OXS 10.7不似乎可以安装它，并且

• 此函数似乎对纬度/长度数据和地球曲率的正确距离值有效。由于我正在处理一个垂直平面（和几百米的短距离，比例尺），我不想纠正我的距离。

有克里格在R中here一个很好的讨论，但由于我上面列出的问题不太明白如何将其应用到我的具体问题。

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

Answer 1

所以有一个你想几个问题回答：

1. 的spTransform功能不会产生SPDF的，但投影之间变换。要创建SPDF，您可以使用简单的data.frame作为开始。要变换df到SPDF：

   coordinates(df) = c("horiz", "depth") 
   

2. 的rgdal OS X的二进制文件可以在http://www.kyngchaos.com找到。但我怀疑你是否需要rgdal。

3. spTransform可以对latlong数据进行操作，也可以对预计数据进行操作。但我认为你不需要rgdal或spTransform，另请参见第1点。

4. 使用点1创建SPDF后，可以使用the info at the post you mentioned继续。