R光栅识别黑色光栅图像

Question

下面的代码在我的图像上生成两个框。我正计划进一步分析这些框内的像素。

在下面的例子中，在红色方块的情况下，我不想继续下去，因为它的右上角有黑色像素。而我想继续在绿色方块的情况下，因为它没有一个黑色像素沿着它的边缘。

library(raster)
r1 <- brick(system.file("external/rlogo.grd", package="raster"))
x <- crop(r1, extent(0,50,0,50))
plotRGB(x)
plot(extent(c(0,20,0,20)), lwd=2, col="red", add=TRUE)
plot(extent(c(21,35,0,10)), lwd=2, col="Green", add=TRUE)

Answer 1

这并不是很好的定义，因为在本例中颜色是由RGB值组成的。但这里有一个您可以适应的通用解决方案。我通过取平均值将这些“平坦化”为单个通道，然后在边界处测试低于阈值的最小值（白色为RGB中的255、255、255，黑色为0、0、0

proceed <- function(f, e, threshold) {
    lns <- as(as(e, 'SpatialPolygons'), 'SpatialLines')
    v <- unlist(extract(f, lns))
    ifelse( min(v, na.rm=TRUE) < threshold, FALSE, TRUE)
}

# flat <- mean(x) # not sophisticated see
# http://stackoverflow.com/questions/687261/converting-rgb-to-grayscale-intensity
flat <- sum(x * c(0.2989, 0.5870, 0.1140))
proceed(flat, extent(c(0,20,0,20)), 100)
proceed(flat, extent(c(21,35,0,10)), 100)

（在看到JBaums的解决方案后有了很大改善；现在已经不存在了）