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繼續。