我有一个值为0或1的矩阵，我想获得一组相邻的1的列表。
例如矩阵
mat = rbind(c(1,0,0,0,0),
c(1,0,0,1,0),
c(0,0,1,0,0),
c(0,0,0,0,0),
c(1,1,1,1,1))
> mat
[,1] [,2] [,3] [,4] [,5]
[1,] 1 0 0 0 0
[2,] 1 0 0 1 0
[3,] 0 0 1 0 0
[4,] 0 0 0 0 0
[5,] 1 1 1 1 1
应该返回以下4个连接的组件：
C1 = {（1,1）;（2,1）}
C2 = {（2,4）}
C3 = {（3,3）}
C4 = {（5,1）;（5,2）;（5,3）;（5,4）;（5,5）}
有人知道如何在R中快速实现吗？我的真实矩阵确实相当大，例如2000x2000（但我希望连接的组件数量会相当少，即200）。
通过更新，您可以将二进制矩阵转换为栅格对象并使用丛集函数。然后，只是数据管理即可返回所需的确切格式。下面的例子：
library(igraph)
library(raster)
mat = rbind(c(1,0,0,0,0),
c(1,0,0,1,0),
c(0,0,1,0,0),
c(0,0,0,0,0),
c(1,1,1,1,1))
Rmat <- raster(mat)
Clumps <- as.matrix(clump(Rmat, directions=4))
#turn the clumps into a list
tot <- max(Clumps, na.rm=TRUE)
res <- vector("list",tot)
for (i in 1:tot){
res[i] <- list(which(Clumps == i, arr.ind = TRUE))
}
然后res在控制台输出：
> res
[[1]]
row col
[1,] 1 1
[2,] 2 1
[[2]]
row col
[1,] 2 4
[[3]]
row col
[1,] 3 3
[[4]]
row col
[1,] 5 1
[2,] 5 2
[3,] 5 3
[4,] 5 4
[5,] 5 5
如果有更好的方法从栅格对象到最终目标，我也不会感到惊讶。再次，2000年到2000年的矩阵对此并不重要。
旧的（错误的答案），但对于希望连接图的组件的人应该有用。
您可以使用igraph软件包将邻接矩阵转换为网络并返回组件。您的示例图形是一个组成部分，因此我删除了一条边进行说明。
library(igraph)
mat = rbind(c(1,0,0,0,0),
c(1,0,0,1,0),
c(0,0,1,0,0),
c(0,0,0,0,0),
c(1,1,1,1,1))
g <- graph.adjacency(mat) %>% delete_edges("5|3")
plot(g)
clu <- components(g)
groups(clu)
然后，最后一行在提示符处返回：
> groups(clu)
$`1`
[1] 1 2 4 5
$`2`
[1] 3
我使用此算法的经验非常快-因此我认为2,000 x 2,000不会有问题。
本文收集自互联网，转载请注明来源。
如有侵权，请联系[email protected] 删除。
我来说两句