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

【问题标题】：What are the units of distance in gstat variogram?gstat 变异函数中的距离单位是什么？
【发布时间】：2012-05-27 05:54:09
【问题描述】：

这是this question 的后续行动。我现在已经设法使用对该问题的答案绘制数据的变异函数。我正在使用的代码（其中im 是一个包含dput 数据的变量，https://gist.github.com/2780792 是

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 手册中找不到任何关于此的内容，所以我假设他们使用了基础数据的单位，但事实并非如此。

有什么想法吗？

【问题讨论】：

我不确定，但它们可能是以公里为单位的大圆距离（例如，前两次观测之间的距离约为 10 公里）。
如果您想手动检查距离，此博客可以帮助您设置函数来运行这些计算：r-bloggers.com/great-circle-distance-calculations-in-r PS：有一个gis.stackexchange.com 用于 GIS 相关问题。这是一个比 SO 小得多的社区，但那里肯定有一些 R 用户！

标签： r spatial

【解决方案1】：

我发现有些人在邮件列表中提出了类似的问题，并重新检查了帮助文件的详细信息。您的数据未投影（proj=+longlat）。正如the documentation 和我在邮件列表中找到的信息所示，它应该基于great circle distances 以公里为单位。查看变异函数中的投影（逻辑）参数：

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

【讨论】：