使用行/列索引子集栅格

gva*_*van 1 gis r raster

子集矩阵或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]  
Run Code Online (Sandbox Code Playgroud)

谢谢!

Sey*_*our 5

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

但是,这是您要找的东西吗?

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

plot(subR)
Run Code Online (Sandbox Code Playgroud)

包功能之所以有用crop()raster是因为它允许您使用行和列对栅格对象进行子集化。