r - 从拉长的缓冲区创建 SpatialLine

Question

我有一个SpatialPolygonsDataFrame是道路的缓冲区。我想通过找到SpatialLine适合该缓冲区的路径来重新创建这条道路的大致路径。

关于如何做到这一点的任何想法？

可重现的例子

library(rgeos)
library(sp)

line <- readWKT("LINESTRING(0 0, 0.2 0.1, 0.3 0.6, 0.4 0.1, .3 .3, 1 1)")

# buffer
  line_buff <- gBuffer(line, width = 0.01)
  plot(line_buff, col="red")

Answer 1

一个非常长的镜头，但你可以通过三角测量到正确的大小，然后找到穿过那些质心的路径，但我不喜欢路径实际上像这样自相交的机会.. .

tri <- sfdct::ct_triangulate(sf::st_as_sf(line_buff), a = 0.0005)
plot(st_coordinates(st_centroid(st_cast(tri$geometry))))