我有相关经度和纬度的数据.如何根据点之间的大圆距离获得此数据的变异函数？

这个简单的例子包含了赤道上的所有数据：

require(geoR)

long <- seq(-179, 180)
x <- sin(pi * long / 180) + rnorm(length(long))
V <- variog(data=x, coords=cbind(long, 0))
# variog: computing omnidirectional variogram
plot(V)

第一个和最后一个点实际上只相差1度,但是我天真的尝试导致变化认为他们被分开了359度.

最佳答案 你应该使用
Semi-variogram from nmle.它允许你指定一个距离矩阵,你可以自己轻松地工作.