对于类似于GIS细化算法（例如r.thin，GRASS GIS）的栅格对象，R中是否存在细化或骨架化算法？

Question

我知道我可以在rgrass7包中执行GRASS GIS中的r.thin，但我不想为此任务建立一个草gis环境。那么R中还有其他解决方案吗？

小例子，期望的输出将是黑线：

library(raster)
#create synthetic raster
values<-c(1:10000)
r <- raster(ncol=100, nrow=100, xmn=0, xmx=100, ymn=0, ymx=100)
r[]<-0
r[]<-ifelse(values>4500 & values<5501,1,0)
plot(r)
lines(c(0,100),c(50,50),col="black")