在一个PNG图中组合ggplot和plotRGB

时间:2015-05-29 07:27:19

标签: r plot ggplot2 raster

我正在尝试将ggplot结果和空间栅格的plotRGB组合成一个图并将其写入PNG。我已经找到了如何组合ggplots或进行多重绘图但是这个特定的设置我不知道如何开始工作。此代码也集成在for循环中,这意味着需要使用print语句来实际绘制ggplot(在PNG或Rstudio中)。

library(raster)
library(ggplot2)

plotRGB(brick(raster(matrix(rnorm(10000),ncol=100)),
              raster(matrix(rnorm(10000),ncol=100)),
              raster(matrix(rnorm(10000),ncol=100))),r=1,g=2,b=3,stretch="lin")
ggdata=data.frame("x"=seq(100),"y"=rnorm(100))
ggplot(ggdata,aes(x=x,y=y))+
  geom_line()

你们中的任何人都知道如何将这两者合并为一个PNG(在for循环中有利地显示,但不是必需的)?

1 个答案:

答案 0 :(得分:5)

由于plotRGB正在使用base图片且ggplot2正在使用grid图片,要合并它们,我们可以使用gridBase包。这段代码似乎有用

require("raster")
require("ggplot2")
require("gridBase")
require("grid")
require("gridExtra")

png("gg.png", width = 600, height = 1200)

grid.newpage()
pushViewport(viewport(width=1, height=0.5, y = .75))
par(omi=gridOMI())
par(mfrow=c(2, 2), mfg=c(1, 1), mar=c(3, 3, 1, 0))
for (i in 1:3) {
  plotRGB(brick(raster(matrix(rnorm(10000),ncol=100)),
                raster(matrix(rnorm(10000),ncol=100)),
                raster(matrix(rnorm(10000),ncol=100))),r=1,g=2,b=3,stretch="lin")
}
upViewport()

pushViewport(viewport(width=1, height=0.5, y = .25))
ggdata=data.frame("x"=seq(100),"y"=rnorm(100))
grid.draw(ggplotGrob(
  ggplot(ggdata,aes(x=x,y=y))+
    geom_line() ))
upViewport()
dev.off()

这会产生以下png: enter image description here

我建议您熟悉已使用软件包的文档,以便能够以您希望的方式绘制所需内容:)