Question

我一直在尝试在各个点之间生成LCP网络。我的一些数据集取得了一些成功，但是每当我尝试运行该函数时，其他集都会导致错误。数据都遵循相同的格式并驻留在同一个CRS中，因此我不明白问题所在。我尝试了一些解决方法但到目前为止还没有采取任何措施。我已经包含了我使用的基本脚本。

library(rgdal)
library(gdistance)
TRI<-raster("pathway.tif")
points<-readOGR("pathway.shp")
cost<-transition(TRI,sum,8)
LCP<-shortestPath(cost,points[1,],points[-1,],output="SpatialLines")

这些参数适用于我的数据集的某些部分，但是在我输入上面的最后一行后，其中一些显示以下错误消息，我不确定修复是什么，因为我对R很新： / p>

Error in validObject(.Object) : 
invalid class “SpatialLines” object: bbox should never contain infinite values
In addition: Warning messages:
1: In .Call("R_igraph_get_shortest_paths", graph, as.igraph.vs(graph,  :
 At structural_properties.c:4517 :Couldn't reach some vertices
2: In min(x) : no non-missing arguments to min; returning Inf
3: In max(x) : no non-missing arguments to max; returning -Inf
4: In min(x) : no non-missing arguments to min; returning Inf
5: In max(x) : no non-missing arguments to max; returning -Inf

非常感谢任何帮助或建议。谢谢！

Answer 1

这不是一个可重复的例子。您似乎在TransitionLayer中创建了一些未连接的组件。检查gdistance vignette以了解避免和检测这种情况的方法。请参阅https://www.jstatsoft.org/article/view/v076i13

的第8页

LCP生成错误

1 个答案: