使用行/列索引对栅格进行子集 [英] Subset a raster using row/column index

问题描述

当子设置矩阵或df时，可以引用行列，例如df1[1:5, 3:10]或df3[2:4, ]。

有什么方法可以使用栅格实现这一点吗？也就是说，我是否可以只裁剪500：700行，例如从栅格对象中裁剪？

我试过使用rasterFromCells()，但没有得到我想要的结果(似乎应该有一个更简单的解决方案，因为R的其他巧妙的子集方法)。

示例：

r <- raster(ncols = 50, nrow = 50)
r[] <- 1:ncell(r)

# I would like to subset the bottom 50 rows of cells, but keep it as a raster.
# However, this returns a numeric object.
rSub <- r[30:50, 1:50]  

谢谢！

推荐答案

我觉得这个问题不太清楚。

但是，这就是您要找的吗？

subR <- crop(r, extent(r, 30, 50, 1, 50))

plot(subR)

rasterPackage的函数crop()可以做到这一点，因为它允许您使用行和列设置栅格对象的子集。

这篇关于使用行/列索引对栅格进行子集的文章就介绍到这了，希望我们推荐的答案对大家有所帮助，也希望大家多多支持IT屋！