gstat variogram中的距离单位是什么？

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手册中，我找不到任何关于此的信息，所以我本来会假设它们使用底层数据的单位，但事实并非如此。

有任何想法吗？

- robintw

---

我不确定，但它们可能是以公里为单位的大圆距离（例如，前两个观测点之间的距离约为10公里）。 - Gavin Simpson

如果您想手动检查距离，此博客可能会帮助您设置函数来运行这些计算：http://www.r-bloggers.com/great-circle-distance-calculations-in-r/ PS：有一个http://gis.stackexchange.com专门回答GIS相关问题。虽然它的社区比SO小得多，但那里肯定有一些R用户！ - Bryan Goodrich

Answer 1

我在邮件列表上发现有些人提出了类似的问题，重新检查了帮助文件的细节。你的数据没有投影（proj=+longlat）。根据文档和我在邮件列表中找到的信息，应该基于大圆距离以公里为单位计算。请查看variogram中的投影（逻辑）参数：

对于投影数据，将计算欧几里得距离，对于未投影的数据将计算大圆距离（公里）。

- Bryan Goodrich