0

複雑なポリゴンによって制約された多数のポイント (緯度と経度) があります。ただし、一部のポイントはポリゴンの境界内にありません。これらのポイントを元のデータ フレーム (以下で説明する ppp オブジェクトではありません) からサブセット化したいと考えています。

#Simple example of a polygon and points.
ex.poly <- data.frame(x=c(0,5,5,2.5,0), y=c(0,0,5,2.5,5))
df <- data.frame(x=c(0.5, 2.5, 4.5, 2.5), y=c(4,1,4, 4))

bound <- owin(poly=data.frame(x=ex.poly$x, y=ex.poly$y))

test.ppp <- ppp(x=df$x, y=df$y, window=bound)

#plotting example, to show the one out of the bound owin object
plot(bound)
points(df$x, df$y)

予想通り、エラーメッセージ1 point was rejected as lying outside the specified windowが表示されます。元のデータフレームをサブセット化して、df拒否されたポイントを見つけるにはどうすればよいですか?

4

2 に答える 2

4

ポイントの境界とリストを生成する

ex.poly <- data.frame(x=c(0,5,5,2.5,0), y=c(0,0,5,2.5,5))
df <- data.frame(x=c(0.5, 2.5, 4.5, 2.5), y=c(4,1,4, 4))

パッケージspatstatを使用して、ポリゴンの内側と外側のポイントを特定してプロットします

library(spatstat)
bound <- owin(poly=data.frame(x=ex.poly$x, y=ex.poly$y))
isin<-inside.owin(x=df$x,y=df$y,w=bound)
point_in <- df[isin,]
point_out <- df[!isin,]
plot(bound)
points(df)
points(point_out,col="red",cex = 3 )
points(point_in,col="green",cex = 3 )

多角形外のspatstat

または、パッケージ splancs を使用する

library(splancs)
pts<-as.matrix(df)
poly<-as.matrix(ex.poly)
point_in<-pip(pts,poly,out=FALSE)
point_out<-pip(pts,poly,out=TRUE)
plot(pts, ylim=c(0,5), xlim=c(0,5))
polygon (poly)
points(point_out,col="red",cex = 3 )
points(point_in,col="green",cex = 3 )

PointsInOutOfPolygon

于 2014-10-12T15:34:21.867 に答える