R中能否实现ggplotGrob的稳定定位？（双异尺度Y轴场景）

你遇到的定位反常和不稳定，核心原因是对annotation_custom()的坐标逻辑、ggplot转grob后的结构理解不够清晰，下面拆解关键技术细节：

1. annotation_custom()的坐标陷阱 #

annotation_custom()默认使用主图的数据坐标系来定位嵌入的grob，而不是绘图面板的绝对/相对坐标。也就是说你设置的xmin=-0.4、ymin=-0.01这些值，是对应主图panel_D中x（log10的width）和y（velocity的value）的数据范围，不是面板的百分比位置：

• 因为主图x轴是log10尺度，数据坐标的数值和视觉位置是非线性的，这就导致你设置的x范围对应的视觉位置完全不符合预期；
• 一旦主图的轴范围调整，嵌入grob的位置会跟着偏移，这就是“不稳定”的根源。

2. ggplotGrob转换后的冗余元素 #

当你用ggplotGrob(panel_D1)把ggplot对象转成grob时，得到的是包含完整绘图元素的grob树：包括原图的x轴、标题、边距、面板等所有部分。而你只需要左侧的Y轴和散点，多余的元素（比如panel_D1中被隐藏的x轴）会占据额外空间，导致嵌入后出现偏移。

3. 可行的修正方案 #

步骤1：提取嵌入图的核心面板grob #

先从panel_D1的grob中只提取面板部分（包含散点、Y轴、刻度），去掉多余的边距和x轴元素：

library(grid)
# 转换为grob并提取面板部分
grob_D1 <- ggplotGrob(panel_D1)
# 找到面板的位置（通常是layout中的"panel"）
panel_pos <- grob_D1$layout$name == "panel"
# 提取面板grob，同时保留左侧Y轴相关元素
core_grob <- grob_D1$grobs[[which(panel_pos)]]

步骤2：用绝对视口（viewport）定位叠加 #

放弃annotation_custom()的数据坐标，改用grid包的视口（viewport）来指定嵌入grob的绝对位置，这样定位更稳定：

# 先绘制主图
print(panel_D)
# 创建左侧视口：宽度占主面板的20%，高度和主面板一致，左对齐
vp <- viewport(x = 0.1, y = 0.5, width = 0.2, height = 0.8, just = c("left", "center"))
# 在视口中绘制提取的核心grob
pushViewport(vp)
grid.draw(core_grob)
popViewport()

优化后的完整代码 #

library(ggplot2)
library(grid)
library(scales)

data_D <- structure(list(width = c(20, 6, 2, 1, 4, 0.6, 0.7484218, 0.9335587, 1.1644929, 1.4525532, 1.8118708, 2.2600728, 2.8191464, 3.5165178, 4.3863979, 5.4714599, 6.8249334, 8.5132152, 10.619127, 13.245977, 16.522631, 20.60983, 25.708077, 32.067477, 40), unit = structure(c(1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("shear_stress", "velocity" ), class = "factor"), value = c(0.00174, 0.01622, 0.16065, 0.65696, 0.02312, 0.01726544, 0.01507592, 0.02141866, 0.02538495, 0.04250598, 0.03747275, 0.05833458, 0.07434515, 0.09323235, 0.12565372, 0.11249062, 0.19224849, 0.21529261, 0.2861165, 0.29305465, 0.42021676, 0.55698482, 0.76004248, 0.783903, 0.83415913)), class = "data.frame", row.names = c(NA, -25L))

# 绘制左侧Y轴的子图
panel_D1 <- ggplot(data=subset(data_D,data_D$unit=="shear_stress"), aes(x=width,y=value)) + 
  geom_point(color="red") + 
  scale_x_log10(expand=c(0,0), limits=c(0.5,50)) + 
  scale_y_log10( expand = c(0, 0), labels = scales::trans_format('log10', math_format(10^.x)), breaks=c(0.001,0.01,0.1,1), limits = c(0.001,1) ) + 
  annotation_logticks(sides = "l") + 
  theme_classic() + 
  theme(axis.title.x=element_blank(), 
        axis.text.x=element_blank(), 
        axis.ticks.x=element_blank(), 
        axis.line.x = element_blank(),
        axis.title.y = element_text(vjust = -5),
        plot.margin = unit(c(0,0,0,0), "cm")) # 去掉子图的边距，避免叠加时出现空隙

# 绘制主图
panel_D <- ggplot(data=subset(data_D,data_D$unit=="velocity"), aes(x=width,y=value))+
  geom_point(color="blue") + 
  scale_x_log10(expand=c(0,0), breaks=c(0.5,1,2,5,10,20,50), labels=c(0.5,1,2,5,10,20,50), limits=c(0.5,50)) + 
  annotation_logticks(sides = "b") + 
  scale_y_continuous(position = "right", expand = c(0,0), breaks = seq(0,1,0.1), limits = c(0,1) ) + 
  theme_classic() + 
  theme(axis.text.y.right = element_text(margin = margin(10,10,10,10)), 
        axis.title.y.right = element_text(vjust = 5,angle = 90), 
        axis.title.x = element_text(vjust = 3), 
        plot.margin = unit(c(0.2,0,0,1), "cm"))

# 提取子图的核心面板grob
grob_D1 <- ggplotGrob(panel_D1)
panel_idx <- which(grob_D1$layout$name == "panel")
core_grob <- grob_D1$grobs[[panel_idx]]

# 叠加绘图
print(panel_D)
# 调整视口参数来精确控制位置，x是水平位置（0最左，1最右），width是宽度占比
vp <- viewport(x = 0.12, y = 0.5, width = 0.22, height = 0.8, just = c("left", "center"))
pushViewport(vp)
grid.draw(core_grob)
popViewport()

补充说明 #

• 视口的x、width参数需要根据你的绘图尺寸微调，确保左侧Y轴和主图的X轴对齐；
• 如果你想更便捷地实现双Y轴，也可以考虑ggside包（专门用于侧边轴），但既然你要手动叠加grob，上述方法更可控。

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