2015-04-24 49 views
1

我有R代码如下。第一plot()和第一plotRGB()命令工作正常。我想分割图片R栅格修剪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) 

回答

2

这是一个错误裁剪后的图像,但它很容易通过“裁剪”功能

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

一个问题解决它。我想用红色边框突出(0,50,0,50)我的图像区域(没有填充)...我怎么能这样做? – user2543622

+1

@ user2543622继续罗伯特的例子,你可以做一些事情:'plot(extent(c(0,20,0,20)),lwd = 2,col =“red”,add = TRUE)' –