我想计算空间中任何点到水文网络的最短距离。结果应该是标量场，标量=到水文网络的最短距离。
我已经尝试了执行作业的光栅包的功能距离
r <- raster(ncol=200,nrow=200)
r[] <- NA
r[c(2000:2060, 15000:15070)]=1
dist <- distance(r)
plot(dist / 1000)
plot(r,add=TRUE,col="blue")
但是当光栅和水文网络非常大时（栅格光栅6000x6500），它非常慢。有没有人有解决方案来加速它？