如何在同一x-y网格中精准叠加两幅Base R绘图？

问题核心在于image()和plot()默认的坐标轴扩展模式不同：

• image()默认使用xaxs="i"和yaxs="i"，坐标轴严格贴合数据范围，无额外空白
• plot()默认使用xaxs="r"和yaxs="r"，会在设置的坐标范围外自动添加5%的空白区域

即使设置了相同的xlim、ylim，这种扩展模式的差异也会导致两幅图错位。

修正后的代码 #

# 绘制初始热力图
image(t(volcano)[ncol(volcano):1,])
# 获取当前绘图的坐标参数
usr <- par("usr")
asp_val <- par("asp")

# 开启新绘图
par(new = TRUE)
# 绘制空图，关键添加xaxs和yaxs参数
plot(1, 
     type = "n",
     xlim = c(usr[1], usr[2]),
     ylim = c(usr[3], usr[4]),
     asp = asp_val,
     col.axis = "red",
     xaxs = "i",  # 与image保持一致的轴扩展模式
     yaxs = "i",
     xlab = "",   # 避免重复显示坐标轴标签
     ylab = ""
)

另一种更直接的写法 #

直接指定image()的坐标范围，再同步给plot()：

# 定义热力图的坐标范围
x_range <- c(1, ncol(volcano))
y_range <- c(1, nrow(volcano))

# 绘制热力图并指定坐标范围
image(x = x_range, y = y_range, z = t(volcano)[ncol(volcano):1,])
par(new = TRUE)

# 绘制匹配的空图
plot(1, 
     type = "n",
     xlim = x_range,
     ylim = y_range,
     asp = 1,
     col.axis = "red",
     xaxs = "i",
     yaxs = "i",
     xlab = "",
     ylab = ""
)

关键注意点 #

• 必须同步xaxs和yaxs参数为"i"，消除默认的空白扩展
• 若需要保留坐标轴标签，注意不要重复显示（设置xlab=""和ylab=""）
• 确保asp参数与初始绘图一致，维持宽高比

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