Autokriging空間數據

我正嘗試使用克立格函數在海洋橫斷面上創建化學參數的垂直地圖，而且我很難開始使用。Autokriging空間數據

我的數據是這樣的：

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對象，但有兩個問題：

有克里格在R中here一個很好的討論，但由於我上面列出的問題不太明白如何將其應用到我的具體問題。

我想要一個矩陣或數據框來描述一個帶有內插值的網格點，然後我可以繪製我的化學參數（理想情況下使用ggplot2）。我懷疑我的問題的解決方案比我想要的要容易得多。

2012-06-13 Drew Steen

您應該指定您的MacOS版本。如果它小於10.6.x（因爲我的意思是和你錯誤地提出「沒有可用的二進制文件」），那麼就'rgdal'而言，你可能運氣不好。 –

編輯的問題反映我正在使用OSX 10.7。我無法在CRAN上找到rgdal的二進制文件（每個鏈接） - 他們是否可以在其他地方使用？ - 哦，@Paul Hiemstra提到他們在哪裏。 –

看到我的答案爲鏈接到rgdal二進制文件的mac，但我懷疑你需要他們的這個特定的問題。儘管rgdal是一件好事。 –

所以有一個你想幾個問題回答：

的spTransform功能不會產生SPDF的，但投影之間變換。要創建SPDF，您可以使用簡單的data.frame作爲開始。要變換df到SPDF：
```
coordinates(df) = c("horiz", "depth") 
```
的rgdal OS X的二進制文件可以在http://www.kyngchaos.com找到。但我懷疑你是否需要rgdal。
spTransform可以對latlong數據進行操作，也可以對預計數據進行操作。但我認爲你不需要rgdal或spTransform，另請參見第1點。
使用點1創建SPDF後，可以使用the info at the post you mentioned繼續。

2012-06-13 15:47:33

你能在你的答案中包含包含'coordinates（）'的包嗎？不知道這是如此簡單！ –

它來自SP包，它也提供了SPDF。 –

回答