栅格图裁切图RGB图像

Question

我有如下的R代码。 第一个plot()和第一个plotRGB()命令可以正常工作。 我想分割图像

在第二个代码块中plot()可以工作，但是plotRGB()失败。 在第二种情况下，如何使plotRGB工作？ 我的最终目标是获取原始颜色的裁剪图像

 r <- raster(system.file("external/rlogo.grd", package="raster"))
 r1 <- brick(system.file("external/rlogo.grd", package="raster"))
 plot(r)
 plotRGB(r1, interpolate = TRUE)

 plot(r,xlim=c(0,50),ylim=c(0,50),legend=FALSE,axes=FALSE,frame.plot=F)
 plotRGB(r1,xlim=c(0,50),ylim=c(0,50),legend=FALSE,axes=FALSE,frame.plot=F)

Answer 1

那是一个错误，但是可以通过“ crop”功能轻松解决

r1 <- brick(system.file("external/rlogo.grd", package="raster"))
x <- crop(r1, extent(0,50,0,50))
plotRGB(x)