如何更换坐标系且不丢失新转换绘图中的数据点？

核心原因 #

你的coords是普通DataFrame，没有绑定坐标系信息；直接使用coord_sf指定新CRS时，ggplot无法完成正确的坐标系转换，导致数据点要么超出绘图范围，要么坐标映射失败。

具体解决步骤 #

1. 将普通数据框转为带坐标系的sf对象 #

首先必须给你的坐标数据绑定原始坐标系（比如常见的WGS84经纬度对应EPSG:4326），转换为sf类：

library(sf)
# 假设你的x、y是WGS84经纬度，指定CRS为4326
df_sf <- st_as_sf(df, coords = c("x", "y"), crs = 4326)

2. 确认海岸线数据的坐标系 #

先检查coastlines_sf的CRS，确保后续转换逻辑一致：

st_crs(coastlines_sf)

3. 两种方式实现坐标系统一绘图 #

方式一：提前转换所有数据到目标CRS

先把数据和海岸线都转成目标坐标系，再绘图：

# 示例目标CRS为墨卡托投影（EPSG:3857）
target_crs <- 3857
df_sf_transformed <- st_transform(df_sf, crs = target_crs)
coastlines_sf_transformed <- st_transform(coastlines_sf, crs = target_crs)

ggplot() +
  geom_sf(data = coastlines_sf_transformed) +
  geom_sf(data = df_sf_transformed, color = "red", size = 2)

方式二：用coord_sf自动转换

确保数据已绑定原始CRS后，直接在coord_sf中指定目标CRS，ggplot会自动完成转换：

ggplot() +
  geom_sf(data = coastlines_sf) +
  geom_sf(data = df_sf, color = "red", size = 2) +
  coord_sf(crs = 3857) # 替换为你需要的目标CRS代码

4. 排查数据范围问题 #

如果仍不显示，检查转换后的数据点范围是否与海岸线匹配：

# 查看数据点和海岸线的边界范围
st_bbox(df_sf_transformed)
st_bbox(coastlines_sf_transformed)

若数据点范围超出海岸线显示范围，手动调整绘图边界：

ggplot() +
  geom_sf(data = coastlines_sf_transformed) +
  geom_sf(data = df_sf_transformed, color = "red", size = 2) +
  xlim(st_bbox(df_sf_transformed)[c(1,3)]) +
  ylim(st_bbox(df_sf_transformed)[c(2,4)])

关键注意事项 #

• 普通DataFrame无CRS属性，直接用coord_sf会导致转换失败，必须先转为sf对象并指定原始CRS。
• 原始CRS必须与坐标类型匹配（比如经纬度对应EPSG:4326，平面投影坐标对应对应EPSG代码）。

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