如何在ggplot2的geom_sf()中对坐标点添加经度方向抖动？

我有一个名为DRCac的坐标点数据框，结构如下：

该数据框包含近400个点，使用以下代码绘制时，叠加在名为System的shapefile上的点重叠情况严重：

DRCmap <- ggplot() +
  geom_sf(data = System) +
  geom_sf(data = DRCac, aes(color = Species)) +
  coord_sf(ylim = c(36.246, 36.253), xlim = c(-92.241,-92.245), expand = FALSE)

请问是否可以在geom_sf()函数内对这些点添加抖动效果？若不可行，是否需要使用其他函数？另外，图例位于东侧，抖动需沿经度方向进行。

1. 直接修改sf对象实现经度定向抖动 #

geom_sf本身没有内置抖动参数，但可以通过调整几何对象的经度坐标实现定向抖动，纬度保持不变：

library(sf)
library(ggplot2)

# 设置随机种子保证抖动效果可重复
set.seed(123)

# 生成带经度抖动的新sf对象
DRCac_jittered <- DRCac %>%
  mutate(
    # 提取原始经纬度
    lon = st_coordinates(.)[,1],
    lat = st_coordinates(.)[,2],
    # 仅对经度添加随机抖动，sd控制抖动幅度（根据坐标范围调整）
    lon_jitter = lon + rnorm(nrow(.), mean = 0, sd = 0.0001),
    # 重新生成几何对象
    geometry = st_sfc(st_point(cbind(lon_jitter, lat)), crs = st_crs(DRCac))
  ) %>%
  select(-lon, -lat, -lon_jitter) # 清理中间变量

# 绘制带抖动的图层
DRCmap <- ggplot() +
  geom_sf(data = System) +
  geom_sf(data = DRCac_jittered, aes(color = Species)) +
  coord_sf(ylim = c(36.246, 36.253), xlim = c(-92.241,-92.245), expand = FALSE) +
  theme(legend.position = "east")

2. 转换为普通数据框用geom_jitter实现 #

如果不想修改原始sf对象，可以先提取坐标转为普通数据框，再使用geom_jitter定向抖动：

library(sf)
library(ggplot2)

set.seed(123)

# 将sf对象转为普通数据框，提取经纬度
DRCac_df <- DRCac %>%
  mutate(
    lon = st_coordinates(.)[,1],
    lat = st_coordinates(.)[,2]
  ) %>%
  st_drop_geometry()

# 绘制时仅抖动经度（x轴），height=0表示不抖动纬度
DRCmap <- ggplot() +
  geom_sf(data = System) +
  geom_jitter(data = DRCac_df, aes(x = lon, y = lat, color = Species),
              width = 0.0001, height = 0) +
  coord_sf(ylim = c(36.246, 36.253), xlim = c(-92.241,-92.245), expand = FALSE) +
  theme(legend.position = "east")

关键注意点 #

• 抖动幅度（sd或width）需要根据你的坐标范围微调，避免点抖动后超出地图可视范围
• 两种方法都能实现仅沿经度方向的抖动，且不影响图例位置设置
• 设置set.seed()可以让抖动效果固定，方便后续复现

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