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和我在邮件列表中找到的信息所示，它应该以km为基础great circle distances。查看变异函数中的投影（逻辑）参数：

  

对于投影数据，计算欧几里德距离，用于未投影   大圆距离（km）。