我正在使用光栅功能,如下面的行所示。我的最后一行产生了一些输出。该输出有一行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)
答案 0 :(得分:5)
您可以通过设置xlim
和ylim
:
plot(r,xlim=c(0,299),ylim=c(0,199))
[<强>更新强>]
要摆脱白色背景,您可以尝试useRaster=F
参数:
plot(r,xlim=c(0,100),ylim=c(0,100),legend=FALSE,axes=FALSE,frame.plot=F,useRaster=F)
答案 1 :(得分:4)
如果理解你的问题,你有一个带有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)
另一种方法(对于较小尺寸的栅格)将使用索引和drop = 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,如mai
png('test.png', width=450, height=275)
plot(rc)
dev.off()
处理此问题的其他方法包括在rasterVis包中使用spplot或levelplot
spplot(rc)
library(rasterVis)
levelplot(rc)
获取原始图像上的红色矩形
plot(r)
plot(e, add=TRUE, col='red', lwd=2)