R 寻找图片边缘

发布于04月18日

我有一张图片，它由每个像素中的坐标(x，y)和一些强度参数组成.它的形状不规则.它还可能有一些洞.这是我拥有的矩阵的一个例子:

set.seed(20)
test <- expand.grid(1:25, 1:50)
test$Int <- rnorm(nrow(test), 7, 1)
edge_rows <- which(test$Var1 %in% c(1, 25) | test$Var2 %in% c(1, 50))
rows_to_remove <- sample(edge_rows, 50) 
test <- test[-rows_to_remove, ]
edge_rows <- which(test$Var1 %in% c(2, 24) | test$Var2 %in% c(2, 49))
rows_to_remove <- sample(edge_rows, 25) 
test <- test[-rows_to_remove, ]
rows_to_remove <- sample(edge_rows, 75) 
test <- test[-sample(nrow(test), 50), ]
ggplot(test) +
    geom_tile(aes(Var1, Var2, fill = Int))

没有(x，y，Int)注释且被其他像素包围的像素被认为是洞，我正在寻找仅针对边缘的索引.

提前感谢您

推荐答案

我有点不清楚您的"结果边缘指数"，但会try /建议:

library(terra)
test_rast = terra::rast(test) # your data as above

plot(boundaries(test_rast, classes = FALSE, inner = TRUE)) # or
# test_bound_inner <- boundaries(test_rast, classes = FALSE, inner = TRUE)
# plot is wrapping the boundaries function but doesn't result in an object, just the plot
# plot(test_bound_inner) 
# retains test_bound_inner as an object in your session

或者:

plot(boundaries(test_rast, classes = FALSE, inner = FALSE))
# test_bound_inF <- boundaries(test_rast, classes = FALSE, inner = FALSE)

关于x/y索引的一些 idea -这里使用上面不太喜欢的第二张图像.

test_bound_inF_mtx <- as.matrix(test_bound_inF, wide = TRUE)
> which(test_bound_inF_mtx == 0, arr.ind = TRUE)
        row col
   [1,]   1   1
   [2,]   3   1
   [3,]   5   1
   [4,]  10   1
   [5,]  12   1
   [6,]  13   1
   [7,]  14   1
   [8,]  15   1
   [9,]  16   1
  [10,]  17   1
<---snip--->

但是，如果总体轨迹进一步对齐连续切片，(例如如果构建3D精细化)，您可能需要?terra::?xyFromCell的某种变体，以给出像素重心的x/y.切片可以简单地与my_many_slices <- c(slice1, slice1, slice3 & etc)"堆叠".