3

自動マッピングを実行しようとしています。不確実性がアイソラインで表示される不確実性マップが必要です。これらのアイソラインとそれに関連付けられた値をエクスポートし、Web マップ サービスを介して表示したいと考えています。出来ますか？次のコードにより、等高線を生成できます

require(gstat)
require(lattice)
data(meuse)
data(meuse.grid)
coordinates(meuse) = ~x+y
gridded(meuse.grid) = ~x+y
m <- vgm(.59, "Sph", 874, .04)
# ordinary kriging:
x <- krige(log(zinc)~1, meuse, meuse.grid, model = m)
l=as.data.frame(x)
contourplot(var1.var ~ x+y,data=l,cuts=10)
4

1 に答える 1

1

これはあなたの質問に対する直接的な答えではありませんが、正しい方向に導くことができるかもしれません:

等高線図からオブジェクトを生成できます

CP<- contourplot(var1.var ~ x+y,data=l,cuts=10)

そして、その構造を調べてみてください

fix(CP)

アイソラインに関するデータは次の場所に保存されます。

x <- CP$panel.args.common[1]$x
y <- CP$panel.args.common[2]$y
z <- CP$panel.args.common[3]$z

しかし、私はそれを抽出する方法を今は理解できません。多分これはあなたを少し助けることができます。

于 2013-01-21T10:35:05.347 に答える