Question

这是><这个问题。我现在已经设法用对该问题的答复绘制了数据的方略图。我使用的代码(这里的 im 是一个变量,包含dput 数据,在上提供的数据, https://gist.github.com/2780792。

library(gstat)
point_data <- as(im,  SpatialPointsDataFrame )
v <- variogram(band1 ~ 1, data = point_data)
plot(v)

This gives me the following plot: enter image description here

可以看到,距离轴上的数值介于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轴上的距离度量单位是什么?我期望它们与坐标站的单位相同。它们与坐标站的单位是多少。 3米是没有道理的, 我猜它们可能是千米。

在Gstat手册里我找不到任何关于这个的信息所以我估计他们用的是基础数据单位但情况似乎并非如此

有什么想法吗?

Answer 1

我发现一些人在邮件列表中提出了类似的问题,并重新核对了帮助文件的细节。您的数据没有预测(projlonglat ) 。 “ a href=” http://stat.ethz.ch/CRAN/web/packages/gstat/gstat.pdf ” rel= “noreferr” 文档。我在邮件列表中发现的信息显示,它应该以 http://en.wikipedia.org/wiki/Greater_terrate' rel=“noreferrer” > grandcour 距离 < /a> 。请查看 varigraph 中的预测(逻辑)参数:

For projected data, Euclidian distances are computed, for unprojected great circle distances (km).

友情链接