gstat變異函數中的距離單位是多少？

這是this question的後續行動。我現在已經設法使用給出的答案來繪製數據的變異圖。我使用的代碼（其中im是包含在https://gist.github.com/2780792給出的dput數據的變量是gstat變異函數中的距離單位是多少？

library(gstat) 
point_data <- as(im, 'SpatialPointsDataFrame') 
v <- variogram(band1 ~ 1, data = point_data) 
plot(v)

這給了我下面的情節： enter image description here

正如你可以看到，在距離軸上的數字範圍從0到150左右。然而，數據的單位是度：

> head(coordinates(point_data)) 
      x  y 
[1,] -1.849353 52.06165 
[2,] -1.759728 52.06165 
[3,] -1.401227 52.06165 
[4,] -1.311602 52.06165 
[5,] -1.221977 52.06165 
[6,] -1.132352 52.06165

考慮到，什麼是對的距離度量單位X軸？我希望他們與座標處於同一單位。他們沒有意義是米，我猜他們可能是km。

我在gstat手冊中找不到任何關於此的內容，所以我假設他們使用了底層數據的單位，但似乎並非如此。

任何想法？

來源

2012-05-26 robintw

我不確定，但它們可能是以公里爲單位的大圓距離（例如前兩個觀測值之間的距離爲〜10km）。 –

如果你想手動檢查距離，這個博客可以幫助你設置功能來運行這些計算：http://www.r-bloggers.com/great-circle-distance-calculations-in-r/ PS：有是與GIS相關的問題http://gis.stackexchange.com。這是一個比SO小得多的社區，但肯定有一些R用戶！ –

我發現有些人在郵件列表中提出類似的問題，並重新檢查了幫助文件的細節。您的數據未投影（proj = + longlat）。由於the documentation和我在郵件列表中發現的信息表明，它應該在公里的基礎上great circle distances。退房的投影（邏輯）參數變差：

對於投影數據，歐幾里得距離計算，對未投影大圓距離（公里）。

來源

2012-05-26 21:54:09

gstat變異函數中的距離單位是多少？

回答

相關問題