如何在R中将大型ggtree/ggplot绘图拆分至多页PDF？

核心思路 #

你的问题在于生成了一个超大的单张树图（带多个facet），直接输出到单页PDF会导致压缩模糊。由于不能拆分树结构，我们可以将生成的完整gtable（绘图对象）按页面尺寸拆分为多个子对象，分别输出到PDF的不同页面。

修正原代码的错误 #

首先，你原代码中的assign(p11, facet_plot(...))是错误用法，应该直接赋值给p11，否则无法正确生成最终绘图对象：

# 替换原assign代码为直接赋值
p11 <- facet_plot(p11, panel=colnames(info_df[11]), geom=geom_text,
                  aes(x=0, label=unlist(info_df[11])), data=info_df)

完整实现代码 #

下面是调整后的代码，将生成的树图拆分为多页横向A4 PDF：

library(ggtree)
library(grid)

make_tree_plot <- function(phy_tree, tip_label){
   p <- ggtree(phy_tree) + geom_nodepoint(color='purple', size=2, alpha=0.2) + geom_tippoint() 

  # 保留你原来的facet_plot代码，最终生成p11
  ...

  # 修正赋值：直接将facet_plot结果赋值给p11
  p11 <- facet_plot(p11, panel=colnames(info_df[11]), geom=geom_text,
                    aes(x=0, label=unlist(info_df[11])), data=info_df)

  # 生成gtable并返回，而不是直接绘制
  gt <- ggplot_gtable(ggplot_build(p11))
  return(gt)
}

# 生成完整的树图gtable
gt <- make_tree_plot(tree, tips_info)

# 初始化多页PDF（横向A4，尺寸为11.69x8.27英寸，留边距）
pdf("multi_page_tree.pdf", width=11.69, height=8.27, paper='a4r')

# 计算总宽度和单页可容纳宽度（单位：英寸）
total_width <- sum(gt$widths)
page_width <- unit(11, "in")

# 计算拆分的列索引：找到累计宽度超过单页宽度的位置
cum_widths <- cumsum(gt$widths)
split_indices <- which(cum_widths > page_width)

# 根据拆分情况输出每页内容
if(length(split_indices) == 0){
  # 不需要拆分，直接输出整图
  grid.draw(gt)
} else {
  # 输出第一页：从第1列到第一个拆分点前一列
  gt_page1 <- gt[, 1:(split_indices[1]-1)]
  grid.draw(gt_page1)
  
  # 输出中间页
  prev_col <- split_indices[1]
  for(i in 2:length(split_indices)){
    gt_page <- gt[, prev_col:(split_indices[i]-1)]
    grid.draw(gt_page)
    prev_col <- split_indices[i]
  }
  
  # 输出最后一页：从最后一个拆分点到最后一列
  gt_last <- gt[, prev_col:ncol(gt)]
  grid.draw(gt_last)
}

dev.off()

关键调整说明 #

1. 返回gtable对象：修改函数使其返回完整的绘图gtable，而不是直接调用grid.draw，这样方便后续拆分。
2. 按页面尺寸拆分：由于你使用横向A4（a4r），树图是横向展开的，所以按列拆分gtable；如果是纵向树图，则改为按行拆分（调整为计算total_height和page_height，拆分gt的行）。
3. 边距适配：设置page_width为11英寸（略小于横向A4的11.69英寸），避免内容超出页面边距。你可以根据实际情况微调这个数值。
4. 兼容无拆分场景：如果树图宽度不超过单页尺寸，直接输出整图。

内容的提问来源于stack exchange，提问作者SixtyCorrupted