0

data.frame内の境界楕円体のx、y座標があります。そしてquery、data.fameにいくつかのx、y座標があります。境界楕円体のx、y座標は、関数...を使用して計算されました。

exy <- ellipsoidhull(X[,1:2])

そのような....

plot(predict(exy), xlim=c(-0.018, 0.015), ylim=c(-0.018,0.015), 
     cex=0.1, type="l")

私にこのようなプロットを与えます...

ここに画像の説明を入力してください

私はこのようなクエリを持っています...

     V2      V3
-0.0167 -0.0137
-0.0159 -0.0127
-0.0150 -0.0127
-0.0164 -0.0137
-0.0164 -0.0134
-0.0173 -0.0131

query境界楕円体の内側/外側にあるものをどのように見つけることができますか？これを行うためのR関数はありますか？ありがとう

4

1 に答える 1

2

パッケージmgcvはそのような機能を提供します (ただし、それだけではありません。空間オブジェクトについて知りたい場合は、たとえば、を参照してくださいsp::overlay)。これはin.out()関数の例です。

library(mgcv)
data(columb.polys)
bnd <- columb.polys[[2]]
plot(bnd,type="n")
polygon(bnd)
x <- seq(7.9,8.7,length=20)
y <- seq(13.7,14.3,length=20)
gr <- as.matrix(expand.grid(x,y))
inside <- in.out(bnd,gr)
points(gr,pch=as.numeric(inside)+1)

ここに画像の説明を入力

于 2011-03-15T08:32:50.137 に答える