gstat变异函数中的距离单位是什么？

Question

这是this question的后续版本。现在，我已经成功地使用这个问题的答案绘制了数据的方差函数。我使用的代码(其中im是包含https://gist.github.com/2780792中给定的dput数据的变量)是

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

这给我提供了以下图：

​

如您所见，距离轴上的数字范围从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轴上的距离度量的单位是什么？我希望他们与坐标在相同的单位。对于他们来说，米是没有意义的，我猜他们可能是公里。

我在gstat手册中找不到任何关于这方面的内容，所以我以为他们使用了底层数据的单位，但事实似乎并非如此。

有什么想法吗？

Answer 1

我发现一些人在邮件列表上提出了类似的问题，并重新检查了帮助文件的详细信息。您的数据未投影(proj=+longlat)。正如我在邮件列表中找到的the documentation和信息所表明的那样，它应该是基于great circle distances的km。检查变差函数中的投影(逻辑)参数

对于投影的数据，计算欧几里得距离，对于未投影的大圆距离(Km)。