r - 最接近估计路径的点

Question

类似于我之前的问题:Closest point to a path

我希望能够找到最接近路径的所有中心。然而，路径中缺少一些数据，我想做线性段以在点之间进行插值以“估计”可能的路径，并且仍然找到可能接近该“估计路径”的中心。

set.seed(1)
n <- 10000
x <- 100*cumprod(1 + rnorm(n, 0.0001, 0.002))
y <- 50*cumprod(1 + rnorm(n, 0.0001, 0.002))

# original path
path <- data.frame(cbind(x=x, y=y))

# path with missing points/points every hundred
path.w.missing <- path[seq(1,n,by=100),]


centers <- expand.grid(x=seq(0, 500,by=0.5) + rnorm(1001), 
                       y=seq(0, 500, by=0.2) + rnorm(2501))

centers$id <- seq(nrow(centers))

没有模拟路径中给定点之间的数百万个线性点......我不确定人们会如何去做......

对我来说，这有点像找到一条线和一个细胞矩阵的交点......各种各样......但也许我离你很远......

任何帮助将不胜感激。

Answer 1

这是否有帮助:

## create data
set.seed(1)
n <- 10000
x <- 100 * cumprod(1 + rnorm(n, 0.0001, 0.002))
y <- 50 * cumprod(1 + rnorm(n, 0.0001, 0.002))

# original path
path <- data.frame(cbind(x=x, y=y))

# path with missing points/points every hundred
path.w.missing <- path[seq(1,n,by=100),]


## proposed solution
library(raster)

l <- SpatialLines(list(Lines(list(Line(path.w.missing)), "1")))
r <- raster(extent(bbox(l)), res = c(0.5, 0.2))
r <- rasterize(l, r, field = 1)

## xy will be a matrix
xy <- rasterToPoints(r)

plot(r)
points(xy)
points(path.w.missing, col = "red")