如何避免ggplot2绘制sf对象时将坐标轴转换为经纬度角度

首先，问题的核心是coord_sf()默认会自动将坐标轴转换回WGS84经纬度（角度格式），哪怕你的sf对象已经是投影后的坐标系。要保留北东坐标（northing-easting），只需要在coord_sf中添加关键参数，再按需优化刻度显示即可。

修改后的可运行代码 #

library(rnaturalearth)
library(ggplot2)
library(sf)

# 获取世界多边形数据
world <- ne_countries()
# 重投影至兰伯特等面积投影
proj <- " +proj=laea +lat_0=0 +lon_0=-70 +x_0=0 +y_0=0 +datum=WGS84 +units=km +no_defs"
world_proj <- spTransform(world, CRS(proj))
# 转换为sf对象
world_proj <- st_as_sf(world_proj)

# 使用ggplot绘图——现在正确显示北东坐标（千米单位）
ggplot()+
  geom_sf(data = world_proj)+
  coord_sf(
    xlim = c(-3200, 4500), 
    ylim = c(-4900, 4900),
    datum = NA  # 核心参数：禁用自动转换到WGS84经纬度
  )+
  # 可选：给坐标轴加上单位标注，提升可读性
  scale_x_continuous(labels = function(x) paste0(x, " km"))+
  scale_y_continuous(labels = function(x) paste0(x, " km"))+
  labs(x = "Easting (km)", y = "Northing (km)")

关键细节解释 #

• datum = NA：这个参数直接告诉ggplot不要做坐标转换，完全沿用sf对象当前投影坐标系下的数值（也就是你设置的兰伯特投影的北东坐标，单位为千米）。而基础plot()函数之所以能正确显示，就是因为它默认不会自动转换坐标系统。
• 额外的scale_x/y_continuous：可以自定义刻度的显示格式，比如加上单位标注，让图表的信息传递更清晰。如果需要调整刻度的间隔，还可以用breaks参数进一步定制，例如：

  scale_x_continuous(
    breaks = seq(-3000, 4000, 1000),
    labels = function(x) paste0(x, " km")
  )
  

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