ggplot2 中面的旋转

Question

我想我有一个棘手的案例。我正在绘制进化植物疾病水平，使用 geom_raster : x 和 y 是任意场坐标，z 是在几个时间点测量的疾病水平，我希望将每个日期绘制在不同的方面。
到目前为止，没有问题。下面是一个模拟数据集和代码:

library(ggplot2)
data <- data_frame(month=factor(rep(c("march","april","may","june"), each=100), levels=c("march","april","may","june")),
              x=rep(rep(1:10, each=10), 4),
              y=rep(rep(1:10, 10), 4),
              z=c(rnorm(100, 0.5, 1), rnorm(100, 3, 1.5), rnorm(100, 6, 2), rnorm(100, 9, 1)))
ggplot(data, aes(x=x, y=y, fill=z)) +
  geom_raster(color="white") +
  scale_fill_gradient2(low="white", mid=mean(range(dat$z)), high="red") +
  scale_x_discrete(limit=1:10, expand = c(0, 0)) +
  scale_y_discrete(limit=1:10, expand = c(0, 0)) +
  coord_equal() +
  facet_wrap(~month)

但我真正想要的是让每个面以特定角度(例如 15°)旋转，以反射(reflect)我的场没有完全按照北方定向的事实(即，顶部不是北，底部不是不是南)。
ggplot2 或任何与网格相关的工具是否有可能自动执行此操作？即使是将各个方面保存到图像、旋转它们并在新页面上打印旋转图像的自动方法也足以满足我的需要。这是我想要获得的图像示例(在图像编辑器中将刻面旋转 15°):
http://imgur.com/RYJ3EaR

Answer 1

这是一种独立旋转刻面的方法。我们为 month 的每个级别创建一个包含单独旋转图的列表。 ，然后使用 grid.arrange将四个地块布置在一起。我还从各个图中删除了图例并分别绘制了图例。下面的代码包含一个帮助函数来提取图例。

我将图例对象提取到 lapply 中的全局环境中下面的函数(更不用说多次重复提取了)。可能有更好的方法，但这种方法很快。

library(gridExtra)

# Helper function to extract the legend from a ggplot
# Source: http://stackoverflow.com/questions/12539348/ggplot-separate-legend-and-plot
g_legend<-function(a.gplot){
  tmp <- ggplot_gtable(ggplot_build(a.gplot))
  leg <- which(sapply(tmp$grobs, function(x) x$name) == "guide-box")
  legend <- tmp$grobs[[leg]]
  legend
}

# Create a list containing a rotated plot for each level of month
pl = lapply(unique(data$month), function(m) {

  # Create a plot for the current level of month
  p1 = ggplot(data[data$month==m,], aes(x=x, y=y, fill=z)) +
    geom_raster(color="white") +
    scale_fill_gradient2(low="white", high="red", 
                         limits=c(floor(min(data$z)), ceiling(max(data$z)))) +
    scale_x_discrete(limit=1:10, expand = c(0, 0)) +
    scale_y_discrete(limit=1:10, expand = c(0, 0)) +
    coord_equal() +
    facet_wrap(~month) 

  # Extract legend into global environment
  leg <<- g_legend(p1)

  # Remove legend from plot
  p1 = p1 + guides(fill=FALSE)

  # Return rotated plot
  editGrob(ggplotGrob(p1), vp=viewport(angle=-20, width=unit(0.85,"npc"), 
                                       height=unit(0.85,"npc")))                    
})

# Lay out the rotated plots and the legend and save to a png file
png("rotated.png", 1100, 1000)
grid.arrange(do.call(arrangeGrob, c(pl, ncol=2)),
             leg, ncol=2, widths=c(0.9,0.1))
dev.off()