圖像光柵R包 - 沒有背景和邊框和圖例的光柵打印

Question

我正在使用柵格函數，如下面的行所示。我的最後一行產生一些輸出。該輸出有一行說dimensions : 240, 320, 76800 (nrow, ncol, ncell)。我想重印那張圖片，但只說前200行和前300列。我怎樣才能做到這一點？下面的第二個最後一行繪製整個圖像

f <- "pictures/image1-1421787394.jpeg" 
f 
r <- raster(f) 
plot(r); 
r 

============================= UPDATE1

我做了png(filename = '~/x.png');par(mar=rep(0, 4), xpd = TRUE, oma=rep(0, 4),bty='n') ; plot(r,xlim=c(0,200),ylim=c(0,200),legend=FALSE,axes=FALSE); dev.off()來保存裁剪後的圖像。我能夠擺脫傳說，斧頭和黑匣子。但問題是保存的圖像比剪切的部分多 - 例如圖像周圍的白色部分。我只想保存原稿的裁剪部分（保持圖像大小200 * 200像素）。請讓我知道該怎麼做？

此外，我怎麼能添加一個紅色正方形對應於上述裁剪部分原始圖像？我的意思是我想在原始圖像的頂部獲得一個紅色正方形（僅邊），然後將此圖像（原始圖像+正方形）保存爲新圖像。

我該怎麼做？

UPDATE2 ++++++++++++++++++++++++++++++++++++++++++++++++

添加重複的例子來說明我所說的白色背景

以下地塊的最後一行裁剪圖像的意思。我希望該圖像爲100 * 100，因爲我的xlim和ylim是100.但是我看到如下例所示的白色背景。 （你看不到的後臺，但如果你在你的機器上運行的代碼，打開圖像，你會看到它。）

library(raster) 
r <- raster(nrow=240, ncol=320) 
values(r) <- 1:ncell(r) 
plot(r) 
plot(r,xlim=c(0,100),ylim=c(0,100),legend=FALSE,axes=FALSE,frame.plot=F) 

Answer 1

您可以通過設置xlim和ylim做到這一點：

plot(r,xlim=c(0,299),ylim=c(0,199)) 

[UPDATE] 爲了擺脫白色背景的，你可以嘗試useRaster=F參數：

plot(r,xlim=c(0,100),ylim=c(0,100),legend=FALSE,axes=FALSE,frame.plot=F,useRaster=F) 

Answer 2

如果很好地理解你的問題，你有一個帶dim（r）c（240,320,1）的RasterLayer r，你想把它裁剪到前200行和300列，然後繪製它白色空間。

總是提供示例數據。在這種情況下，這很容易做到。

library(raster) 
r <- raster(nrow=240, ncol=320) 
values(r) <- 1:ncell(r) 

有不同的方式按行/列裁剪。例如，您可以創建一個擴展盤對象並使用它。

e <- extent(r, 1, 200, 1, 300) 
rc <- crop(r, e) 

另一種方式（對於更小尺寸的柵格）是使用索引和降= FALSE

rc <- r[1:200, 1:300, drop=FALSE] 

爲了使地圖，你可以使用「陰謀」或「圖像」。也許，圖像更自己的喜好（更少的空白區域，但沒有傳說）

image(rc) 

與陰謀，你可以繪製之前設置該設備的大小。

dev.new(height=nrow(r), width=ncol(r)) 
plot(rc, legend=FALSE) 

你也可以繪製到像png這樣的文件來避免空白;這取決於你如何設置你的「par'ameters如麥

png('test.png', width=450, height=275) 
plot(rc) 
dev.off() 

其他的方式來處理這包括使用spplot，或levelplot在rasterVis包

spplot(rc) 
library(rasterVis) 
levelplot(rc) 

要獲得關於紅色矩形原始圖像

plot(r) 
plot(e, add=TRUE, col='red', lwd=2)