R 如何估计地球仪上数据的变异函数？

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

此简单示例包含赤道上的所有数据：

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)

require（geoR）
long从R-sig-geo（R中专门用于空间数据的邮件列表）上的帖子中，我似乎记得R中没有支持大圆周距离的现成功能：


我的建议是投影数据，然后对投影数据执行插值。
您应该使用。它允许您指定一个距离矩阵，您可以自己计算。
据我所知，没有保持距离的投影。您的目标是插值一个全局数据集？也许在R-sig-geo邮件列表上发布一个问题。不，我想绘制一个变异函数。与插值无关。投影数据不是一个选项吗？我想你可能会在投影过程中得到一些误差，但这会对变异函数的总体趋势产生多大的影响。这只是一个选择，如果我能找到一个合适的投影。