这是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 手册中找不到关于此的任何内容,因此我假设他们使用了基础数据的单位,但情况似乎并非如此。
有任何想法吗?
最佳答案
我发现有些人在邮件列表上提出了类似的问题,并重新检查了帮助文件的详细信息。您的数据未投影(proj=+longlat)。如 the documentation我在邮件列表中找到的信息表明,它应该以公里为单位,基于 great circle distances .查看变异函数中的投影(逻辑)参数:
For projected data, Euclidian distances are computed, for unprojected great circle distances (km).
关于r - gstat变异函数中的距离单位是多少?,我们在Stack Overflow上找到一个类似的问题: https://stackoverflow.com/questions/10769884/