ggplot2使用coord_sf时如何将x轴负坐标标签调整为正值

问题根源 #

使用coord_sf时容易出现自定义x轴标签不生效的情况，核心原因有两点：

• coord_sf默认会基于投影后的坐标空间自动生成西经为负值的刻度标签，如果只把投影基准传给datum参数，没有明确指定投影本身，自定义的刻度逻辑很容易被sf内置的默认刻度覆盖
• 直接修改原始数据的x值为正值会导致投影位置完全偏移，破坏EPSG:3978的投影对应关系

可直接复用的修正代码 #

全程保留要求的EPSG:3978投影配置，仅修改刻度显示逻辑，让西经以正值形式展示，匹配参考图的标注样式：

# 提前加载依赖包 ggplot2、sf
library(ggplot2)
library(sf)

plt <- ggplot(dts, aes(x = x, y = y)) +
  geom_tile(aes(fill = val), alpha = 1) +
  # 明确传入EPSG:3978投影，替换之前仅传datum的写法
  coord_sf(
    crs = st_crs(3978), 
    datum = st_crs(3978),
    expand = FALSE
  ) +
  # 自定义x轴刻度：断点位置保留原始负经纬度，仅将显示标签取绝对值去除负号
  scale_x_continuous(
    breaks = seq(-120, -100, by = 2),
    # R4.1及以上可直接用简写函数，低版本替换为 labels = function(x) abs(x)
    labels = \(x) abs(x),
    expand = c(0, 0)
  ) +
  xlab("Longitude (°W)") +
  ylab("Latitude (°N)")

plt

关键注意事项 #

• 关于coord_map使用EPSG:3978的可行性：不推荐操作，coord_map属于已经逐步淘汰的地图绘图接口，要适配EPSG:3978需要手动拼接冗长的proj4参数，容易出现投影变形误差，稳定性远不如原生支持EPSG编码的coord_sf
• 不要改动原始数据的x坐标值为正值：所有投影计算都依赖原始的经纬度坐标匹配EPSG:3978的转换规则，只需要修改刻度的显示文本，完全不影响数据在投影空间的位置

避坑提示：如果修改后刻度标签还是被覆盖，检查coord_sf是否传入了正确的crs参数，不要仅将EPSG:3978传给datum——datum只是坐标基准面，不是投影本身，缺少crs参数时sf会默认使用WGS84地理坐标系，导致自定义刻度失效。

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