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）。メーリングリストで見つけたドキュメントと情報が示すように、大圏距離に基づいてkm単位である必要があります。バリオグラムの射影（論理）パラメーターを確認してください：

投影されたデータの場合、投影されていない大円距離（km）の場合、ユークリディアン距離が計算されます。

于 2012-05-26T21:54:09.403 に答える