如何在ggplot2中使用pdf / ps作为背景图像？

Question

我正在尝试在我使用WMS（Geoserver）创建的地图（矢量图形）上叠加半透明的ggplot地图。

我可以使用grImport包导入R中的pdf（命令PostScriptTrace，readPicture和pictureGrob）。但现在我迷路了。我想在pdf背景图像上叠加ggplot图像。我可以指定ggplot将grob用作背景图像吗？

感谢baptiste让我指向“annotation_custom（）”的方向。 在他的帮助下，我设法为我的问题创建了一个更小的例子：

library("cluster")
library("lattice")
library("ggplot2")
library("grImport")

sink("test.ps")
cat("%!PS\n")
cat("/petal {\n")
cat("newpath\n")
cat("0 0 moveto\n")
cat("-5 10 lineto\n")
cat("-10 20 10 20 5 10 curveto\n")
cat("5 10 lineto\n")
cat("closepath\n")
cat("0 setgray\n")
cat("fill\n")
cat("} def\n")
cat("20 20 translate\n")
cat("5 {\n")
cat("petal 72 rotate\n")
cat("} repeat\n")
cat("showpage")
sink()
PostScriptTrace("test.ps")
test.ps.xml <- readPicture("test.ps.xml")
test.grob <- pictureGrob(test.ps.xml)

# following the example from 'Importing Vector Graphics: The grImport Package for R' by Paul Murrell I can test the graphic and add the .ps to a lattice plot
xyplot(V8 ~ V7, data = flower,
  xlab = "Height", ylab = "Distance Apart",
  panel = function(x, y, ...) {
    grid.symbols(test.ps.xml, x, y, units = "native", size = unit(5, "mm"))
  }
)

# ... but I would like to add the .ps as a background image to a plot with ggplot2
polygon.df <- data.frame(
  id = 1,
  x = c(-60,60,60,-60,-60),
  y = c(0,0,175,175,0)
)
ggplot() + annotation_custom(test.grob)+
  geom_polygon(data=polygon.df, aes(x=x, y=y), alpha=.2, fill="blue") +
  scale_x_continuous(limits=c(-70, 70), expand = c(0,0)) + 
  scale_y_continuous(limits=c(-10, 185), expand = c(0,0))

然而，ps-graphic没有绘制到ggplot。我确信我错过了一些非常微不足道的事情。

编辑：以下内容显示了grob：

qplot(x=polygon.df$x, y=polygon.df$y)+ 
annotation_custom(test.grob)

流程不显示grob：

ggplot() + geom_point(data=polygon.df, aes(x=x, y=y)) +
annotation_custom(test.grob)

我不太明白为什么dows用qplot（）显示grob，但不用ggplot（）。

Answer 1

annotation_custom非常错误，也许它无法处理pictureGrob。您可以尝试以下技巧将grob添加到gtable，

library(ggplot2)
library(gtable)

p <- qplot(1,1) + theme(panel.background=element_rect(fill=NA))
g <- ggplotGrob(p)

pos <- gtable_filter(g, "panel", trim=FALSE)$layout
g <- with(pos, gtable_add_grob(g, test.grob, t, l, b, r, z-1))

grid.newpage()
grid.draw(g)