使用ggplot2的sec_axis()实现波数转微米时轴刻度标签异常问题

问题核心：主坐标轴是从4000到400递减的波数，而转换后的波长是从2.5μm到25μm递增的，默认sec_axis()会直接沿用主轴的刻度位置计算标签，导致次轴标签顺序混乱或显示不符合预期。

方案1：反转主坐标轴（推荐，符合光谱图常规展示） #

将主坐标轴改为从400到4000递增，此时次轴的波长标签会自然呈现从25μm到2.5μm递减的逻辑，和光谱图的常规布局一致（短波对应高波数，位于左侧）：

library(ggplot2)

# 示例数据
data <- data.frame(
  wavenumber = seq(4000, 400, by = -100),
  intensity = runif(37)
)

# 绘图：反转主x轴，次轴转换逻辑不变
p <- ggplot(data, aes(x = wavenumber, y = intensity)) +
  geom_line() +
  scale_x_continuous(
    name = "Wavenumber (cm⁻¹)",
    limits = c(400, 4000),  # 反转主轴范围
    sec.axis = sec_axis(~ 10000 / ., name = "Wavelength (μm)")
  ) +
  theme_minimal()

print(p)

方案2：手动指定次轴刻度（保留主轴反向） #

如果必须保留主轴从4000到400的反向显示，需要手动计算次轴刻度对应的主轴位置和标签，确保次轴标签递增显示：

library(ggplot2)

data <- data.frame(
  wavenumber = seq(4000, 400, by = -100),
  intensity = runif(37)
)

# 定义次轴需要显示的波长刻度
wavelength_breaks <- c(2.5, 5, 10, 15, 20, 25)
# 转换为对应的主轴波数刻度位置
wavenumber_breaks <- 10000 / wavelength_breaks

p <- ggplot(data, aes(x = wavenumber, y = intensity)) +
  geom_line() +
  scale_x_continuous(
    name = "Wavenumber (cm⁻¹)",
    breaks = wavenumber_breaks,  # 主轴刻度对应次轴关键点
    sec.axis = sec_axis(
      ~ 10000 / .,
      name = "Wavelength (μm)",
      breaks = wavenumber_breaks,  # 次轴刻度位置与主轴对齐
      labels = wavelength_breaks   # 手动指定次轴标签为波长值
    )
  ) +
  theme_minimal()

print(p)

验证转换逻辑：

• 4000 cm⁻¹ → 10000/4000 = 2.5 μm
• 400 cm⁻¹ → 10000/400 = 25 μm
  两种方案均能保证转换计算正确，同时刻度标签显示符合预期。

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