Он «привязывается» (проецируется) к ближайшему сегменту независимо от того, как далеко он находится . Таким образом, вам, вероятно, следует подготовить свои данные еще с некоторыми шагами перед вызовом lpp
. Ниже приведен простой пример, который вы можете попробовать имитировать, если все ваши точки находятся в ppp
, называемом X
, а основной компонент вашего linnet
называется L. Конечно, вам нужно адаптировать пороги расстояния и c , но, надеюсь, это поможет вам в правильном направлении.
library(spatstat)
set.seed(42)
X <- runifpoint(10)
L <- simplenet
proj <- project2segment(X, as.psp(L))
i <- which(proj$d<.05)
Xclose <- X[i]
Xfinal <- lpp(Xclose, simplenet)
plot(L)
plot(X, add = TRUE, col = "red")
plot(Xclose, add = TRUE, pch = 2, col = "blue")
plot(Xfinal)