2015-04-24 67 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)' –