如何在r中裁剪（切碎）；

Question

我想裁剪一个栅格，然后将每个和平与栅格的各个部分保存到一个列表（或其他方式）中。

想法：

• 保持所有棋子的覆盖面积相等。
• 用户可以设置要裁剪的范围区域的百分比。即总面积的10%、20%、30%应该是网格的范围。

下面可以看到一段可重现的代码：它是由 Chris 和 罗伯特·希曼斯在此问题。

library(terra)
r <- rast(volcano)
n <- 10

获取感兴趣的起始单元格

rows <- seq(1, nrow(r), by=n)
cols <- seq(1, ncol(r), by=n)    
cells <- cellFromRowColCombine(r, rows, cols)

获取坐标

起始单元格的左上角坐标

xy <- xyFromCell(r, cells)
rs <- res(r)
xy[,1] <- xy[,1] - rs[1]/2
xy[,2] <- xy[,2] + rs[2]/2

添加结束单元格的右下角坐标

xy <- cbind(xy[,1], xy[,1] + n*rs[1], xy[,2] - n*rs[2], xy[,2])

并循环

x <- lapply(1:nrow(xy), function(i) {
         crop(r, xy[i,])
       })

验证

e <- lapply(x, \(i) ext(i) |> as.polygons()) |> vect()
plot(r)
lines(e, col="blue", lwd=2)

注意裁剪后的栅格的右侧。对这些多边形进行统计似乎很糟糕。

还： 我想：

• 开头的 n 表示要在小栅格中裁剪的区域的百分比
• 用户可以输入 10、20 或 30 作为数字，例如，表示百分比决定要裁剪的“窗口”大小的命令
• 为了在窗口上保留相等的空间，为此，该命令应该在 10 到 20 之间稍微变化，例如，以保持对称性。

Answer 1

您确实应该尝试专注于手头的问题，并展示您已经尝试过的内容。因为仍然很难准确地知道你想要什么。无论如何，当 n=20% 和 40（列）乘 60（行）的矩阵时，

cols = 40
rows = 60
n = 20 

您可以找到这样的块大小：

h = sqrt(n/100) * rows
w = sqrt(n/100) * cols 

“proof”

p = h * w
100 * p / (cols * rows) 
# 20

当然 h 和 w不是整数

> h
#[1] 26.83282
> w
#[1] 17.88854

，你会得到

c(seq(1, cols, w) |> round(), cols)
#[1]  1 19 37 40
c(seq(1, rows, h) |> round(), rows)
[1]  1 28 55 60

另一种方法。首先计算块的数量。 20%就是五。在大多数情况下，这实际上应该是偶数，除非您采用整行或整列，因此您无法准确得到该数字。

# number of blocks
nb <- 100/n 
# block rows and cols
brow <- round(nb * rows/(rows+cols))
brow
#[1] 3
bcol <- round(nb / brow)
bcol
#[1] 2

现在你也许可以使用：

round(c(seq(1, rows, rows/brow), rows))
#[1]  1 21 41 60
round(c(seq(1, cols, cols/bcol), cols))
#[1]  1 21 40