2011-03-15 3 views
0

나는 data.frame에서 경계 타원체에 대한 x, y 좌표를가집니다. 그리고 나서 나는 여러 개의 query x, y 좌표를 data.fame에 가지고 있습니다. 바운딩 타원체의 X, Y 좌표는 쿼리 x, y는 내부/외부 경계 타원체를 공동 좌표합니다.

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

날 ....

이와 같은 플롯을 제공

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

같은 것을 .... ... 함수를 사용하여 계산 하였다

enter image description here

이 같은 쿼리를 가지고 ....

 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 함수가 있습니까? 감사합니다

답변

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) 

enter image description here

+0

는 @Roman ... 감사합니다. 'sp '패키지는 내가하고 싶은 일에 완벽한 피트이다. – user645600