2

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

我可以使用 grImport 包（命令 PostScriptTrace、readPicture 和 pictureGrob）在 R 中导入 pdf。但现在我迷路了。我想将 ggplot 图像覆盖在 pdf 背景图像上。我可以指定 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 图形并未绘制到 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)

我不太明白为什么它用 qplot() 显示 grob，但没有用 ggplot()。

4

1 回答 1

0

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)

在此处输入图像描述

于 2013-11-08T16:11:17.743 回答