有关R中的stplanr软件包的问题

Question

在第89页上的stplanr软件包（https://cran.r-project.org/web/packages/stplanr/stplanr.pdf）的route_local函数方面，我需要您的帮助。

您可能会意识到，地图是从示例函数生成的，显示了两点之间的路径（我留下了代码和下面生成的图像）​​。我想做同样的事情。就我而言，这是考虑到我的道路的两点之间的路径。两者都是shapefile文件。我设法生成了要显示的道路（下面的代码），但是我想显示这些道路上任意两点之间的路线。有人能帮我吗？？ 我将其留在以下站点https://github.com/JovaniSouza/JovaniSouza5/blob/master/Example.zip上以下载shapefile。

library(geosphere)
library(sf)
library(stplanr)

roads<-st_read("C:/Users/Jose/Downloads/Example/Roads/Roads.shp")
p <- SpatialLinesNetwork(roads, uselonglat = FALSE, tolerance = 0)
plot(p)

代码生成的地图

示例

from <- c(-1.535181, 53.82534)
to <- c(-1.52446, 53.80949)
sln <- SpatialLinesNetwork(route_network_sf)
r <- route_local(sln, from, to)
plot(sln)
plot(r$geometry, add = TRUE, col = "red", lwd = 5)
plot(cents[c(3, 4), ], add = TRUE)
r2 <- route_local(sln = sln, cents_sf[3, ], cents_sf[4, ])
plot(r2$geometry, add = TRUE, col = "blue", lwd = 3)

Answer 1

尝试一下。为了使示例适合您的情况，您必须将axios的坐标系转换为roads shapefile（或反之）：

points