Вам определенно нужно предоставить воспроизводимый пример , чтобы получить лучший ответ;тем не менее, я думаю, что приведенный ниже скрипт будет служить цели:
library(spatstat)
# setting seeds
set.seed(222)
# two different point patterns
X <- runifpoint(15)
Y <- runifpoint(20)
plot(X, pch=19, main="")
plot(Y, col="red", pch=19, add=T)
#you can get both which and dist from nncross
#N.which <- nncross(X,Y, k=1:20, what="which")
#N.dist <- nncross(X,Y, k=1:20, what="dist")
out <- subset(X, nncross(X,Y, k=1:20, what="dist") < 0.1) # you may change 0.1
plot(out, col="blue", pch=19, add=T)
![enter image description here](https://i.stack.imgur.com/Zf3WW.jpg)
Для приведенного выше графика черные точки представляют X
, а красные точки представляютY
.Синий - это точки пересечения, которые находятся на расстоянии 0,1 .Это расстояние может быть дополнительно изменено.Для получения дополнительной информации см. spatstat
, чтобы вычислить расстояния между двумя различными наборами данных, используя nncross
.