自动克里金插值空间数据

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对象，但存在两个问题：

• 这个在CRAN上没有OSX二进制文件, 我的运行在OXS 10.7上的RStudio似乎无法安装它，以及

• 这个函数似乎仅适用于经度/纬度数据，并且能够为地球曲率校正距离值。由于我处理的是一个垂直平面（和短距离，数百米级别），我不想校正距离。

这里有一个关于在R中进行Kriging的很好的讨论 here, 但由于上述问题，我不太明白如何将其应用到我的具体问题上。

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

- Drew Steen

---

你应该指定你的 MacOS 版本。如果它低于 10.6.x（就像我的版本一样，也像你错误地声称“没有可用的二进制文件”），那么对于 'rgdal' 来说，你可能会很不幸。 - IRTFM

修改问题以反映我正在使用OSX 10.7。我找不到CRAN上rgdal的二进制文件（根据链接） - 它们还在其他地方可用吗？- 噢，@Paul Hiemstra提到了它们的位置。 - Drew Steen

请查看我的答案，其中包含Mac平台上的rgdal二进制文件链接，但我怀疑你在这个特定问题上不需要它们。虽然无论如何，rgdal仍然是一个很棒的工具。 - Paul Hiemstra

Answer 1

所以你想要回答几个问题：

1. The spTransform function does not create SPDF's, but transforms between projections. To create a SPDF you can use a simple data.frame as a start. To transform df to a SPDF:

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

2. OS X binaries of rgdal can be found at http://www.kyngchaos.com. But I doubt if you need rgdal.

3. spTransform can operate on latlong data, but also on projected data. But I do not think you need rgdal, or spTransform, see also point 1.

4. After you create the SPDF using point 1, you can use the info at the post you mentioned to go on.

- Paul Hiemstra

---

дҪ иғҪеңЁеӣһзӯ”дёӯеҢ…еҗ«еҢ…еҗ«coordinates()еҮҪж•°зҡ„еә“еҗ—пјҹжҲ‘д»ҺжқҘдёҚзҹҘйҒ“е®ғиҝҷд№Ҳз®ҖеҚ•пјҒ - tim riffe

它来自于sp软件包，该软件包还提供了SPDF。 - Paul Hiemstra