在下面的矩阵中,我试图识别由 1 组成的每个最大矩形,如图所示。只有超过3列的矩形才能只有一行。只有超过3行的矩形才可以只有一列。矩形应该少于maxcol列和maxrow行,多于mincol和minrow。图中mincol=minrow=2。最大行=最大列=6
尝试使用子集函数,但它提供了 1 的子矩阵,而不是每个叠加的矩形。
mdat <- matrix(c(0,1,1,1,0,0,0,1,0,1,1,1,1,1,1,1,1,1,0,1,0,1,1,0), nrow = 4, ncol = 6, byrow = TRUE)
dimnames(mdat) <- list( c("row1", "row2", "row3","row4"), c("col1", "col2", "col3", "col4","col5","col6"))