запрос x, y совмещает внутренний / внешний ограничивающий эллипсоид

Question

У меня есть координаты x, y для ограничивающего эллипсоида в data.frame.И тогда у меня есть несколько query x, y координат в data.fame.Координаты 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, чтобы сделать это?Спасибо

Answer 1

Пакет 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)