後で R でのクリギングに使用する空間データ ポイントで構成される補間グリッドを作成したいと考えています。私は sp および gstat パッケージで作業してきました。私が作業しているコードは機能しますが、グリッドにすべてのポイントが含まれているわけではなく、それを拡張したいのですが、うまくいかないようです。
私が使用しているコードは次のとおりです (私の data.frame は kr.data です):
x.range <- as.integer(range(kr.data@coords[,1]))
y.range <- as.integer(range(kr.data@coords[,2]))
data.grd <- expand.grid(x=seq(from=x.range[1], to=x.range[2], by=0.5),
y=seq(from=y.range[1], to=y.range[2], by=0.5))
coordinates(data.grd) <- ~x+y
gridded(data.grd) <- TRUE
plot(data.grd, cex=0.5)
points(kr.data, pch=1, col='red', cex=0.7)
出力を追加して、グリッドから離れた右下隅に集中しているポイントを表示したかったのですが、私は初めてで、これを追加する方法がわかりません。空間データを使用してグリッドを作成するのに役立つ場所や、コードを変更する場所について誰か提案できますか? コードを変更してグリッドを拡張すると、データ ポイントが表示されなくなります。Rで空間データを扱うのは初めてです。