R栅格图裁剪图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

这是一个错误，但很容易通过“裁剪”功能解决它

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