使用r-raster向量化栅格砖块对象，这样我就可以对它们进行计数

Question

我有一个红色和蓝色边框圆圈列的图像，如下所示：

​

​

其中列交替为红色和蓝色(在本例中，第一列为红色)

我已经能够创建一块光栅砖，并在RGB层中绘制图像，但我想将这些列计数为一个向量，如上例所示。值1(红色)和2(蓝色)

1,1,1,1,2,2,2,1,1,2,1,1,1 ...

有没有可能清理出我不需要计数的砖块区域，并将砖块折叠成我可以转换成我想要的数字或标签的值？或者有没有更简单的方法让我无法定位？另外，从长远来看，我希望能够将程序指向几个图像，而不是自己打开它们。

编辑:为了弄清楚一些事情，我想从上到下，从左到右计算圆圈的数量。例如，一旦计算完第一列，我想从右边下一列的顶部重新开始。此外，我不确定我是否在朝着正确的方向前进，但我能够从图像中删除所有背景单元。给我留下了一个只有圆圈所在的值的图。

编辑2：

这是我为上面的图片编写的代码。

color.image <- brick("image")
color.image = dropLayer(color.image,4) #gets rid of a channel
plot(color.image)

e <- extent(-10, 240, 45, 84.8) #xmin,xmax, ymin,ymax
ccolor.image <- crop(color.image, e)

plot(ccolor.image)
#thresholding to simplify what I was dealing with
mini=ccolor.image[ccolor.image > 97] = NA
mini=ccolor.image[ccolor.image < 15] = NA
mini=ccolor.image[ccolor.image > 20] = 80
plot(ccolor.image)

mcolor = as.matrix(ccolor.image)
colSums(ccolor.image)
rowSums(ccolor.image)

编辑3：

我想通了！或者至少找到了一种可以做到这一点的方法，一旦我清理了一些代码，我会在稍后发布代码。然而，我仍然希望根据我的简化光栅砖块矩阵的值矩阵来创建一个向量。代码马上就要出来了！

Answer 1

对光栅中的值进行计数的最快方法是freq(x, merge=T)。这将为您提供一列中的值和与行一样多的列中的频率。通过这种方式，我们需要轮询感兴趣的值并对所有其他列(计数)求和。希望这能有所帮助！

freq_vals <- freq ( rasterbrick , merge = T )
sum( freq_vals [ which ( freq_vals$value == 1 ) , 2 : ncol ( freq_vals ) ] )