9

这是对这个问题的跟进。我现在已经设法使用对该问题的答案绘制数据的变异函数。我正在使用的代码（其中包含https://gist.github.com/2780792中给出的数据im的变量是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 手册中找不到任何关于此的内容，所以我假设他们使用了基础数据的单位，但情况似乎并非如此。

有任何想法吗？

4

1 回答 1

13

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

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

于 2012-05-26T21:54:09.403 回答