如何在R语言中使沉积物岩芯照片与绘图Y坐标精准对齐？

我完全懂你现在的困扰——好不容易画出了多面板数据图，也放了岩芯照片，但两者的Y轴刻度对不上，根本没法直观对应层位和元素峰值。别担心，核心思路就是强制统一两者的Y轴坐标系统，让照片的深度刻度和绘图的深度刻度完全匹配。下面给你一套可落地的步骤和代码：

步骤1：准备关键对齐信息 #

首先你需要明确两个核心参数，这是对齐的基础：

• 你的沉积物岩芯的实际深度范围（比如从0cm表层到100cm底层）；
• 用图像工具（比如Windows画图、GIMP）查看岩芯照片，找到卷尺刻度对应的像素位置：
  • 记录表层（0cm）对应的像素Y坐标（比如photo_y_top = 120）；
  • 记录底层（比如100cm）对应的像素Y坐标（比如photo_y_bottom = 2120）。

步骤2：绘制统一Y轴的多面板数据图 #

用ggplot2绘制多面板图时，必须固定Y轴范围并关闭自动扩展，避免额外空白干扰对齐：

# 加载所需工具包
library(ggplot2)
library(grid)
library(patchwork)
library(tidyr)
library(jpeg)

# 替换成你的实际岩芯数据
set.seed(123)
depth <- seq(0, 100, by = 1)
element1 <- rnorm(101, mean = 50, sd = 10) + depth*0.2
element2 <- rnorm(101, mean = 30, sd = 8) + (100-depth)*0.15
core_data <- data.frame(depth, element1, element2)

# 转成长格式方便分面绘图
core_long <- pivot_longer(core_data, cols = c(element1, element2), 
                          names_to = "element", values_to = "concentration")

# 绘制多面板图，核心是固定Y轴范围、关闭扩展
p <- ggplot(core_long, aes(x = concentration, y = depth)) +
  geom_line(linewidth = 0.8) +
  facet_wrap(~element, ncol = 1) +
  scale_y_continuous(limits = c(0, 100), expand = c(0, 0)) + # 强制Y轴范围匹配岩芯实际深度
  labs(y = "Depth (cm)", x = "Concentration") +
  theme_bw() +
  theme(
    strip.text = element_text(size = 12, face = "bold"),
    panel.spacing = unit(0.5, "lines")
  )

步骤3：处理照片并实现精准对齐 #

把照片转换成图形对象，根据之前记录的像素刻度，缩放并放置到绘图右侧，确保Y轴完全贴合：

# 读取岩芯照片（替换成你的照片路径）
core_photo <- readJPEG("core_photo.jpg")

# 计算照片缩放比例，匹配绘图Y轴范围
photo_y_top <- 120    # 你测量的表层像素Y坐标
photo_y_bottom <- 2120 # 你测量的底层像素Y坐标
total_depth <- 100    # 岩芯实际总深度
pixels_per_cm <- (photo_y_bottom - photo_y_top) / total_depth

# 缩放照片，让其高度刚好匹配绘图的Y轴范围
scaled_photo <- rasterGrob(core_photo, 
                           height = unit(total_depth, "native"),
                           interpolate = TRUE)

# 如果照片的深度方向和绘图相反（比如照片表层在底部），翻转照片
# scaled_photo <- flipVertically(scaled_photo)

# 用patchwork组合图，确保Y轴严格对齐
final_plot <- p + 
  inset_element(scaled_photo,
                left = 1.05, right = 1.3, # 调整照片在右侧的显示宽度
                bottom = 0, top = 1,       # 上下完全贴合绘图Y轴范围
                align_to = "panel")        # 对齐到面板坐标系统，避免偏移

print(final_plot)

关键注意事项 #

• 像素测量要精准：尽量放大照片查看刻度对应的像素点，这是对齐的核心；
• 方向一致性：如果照片的深度方向和绘图相反，记得用flipVertically()翻转；
• base绘图适配：如果你习惯用base绘图系统，可通过layout()设置布局，再用grid.raster()放置照片，原理同样是统一Y轴范围。

内容的提问来源于stack exchange，提问作者Jordi F. Pagès