如何在R中可靠可视化跨日界线的缓冲区？

针对你遇到的跨日界线缓冲区可视化异常问题，最稳定的解决方案是先转换到适合的投影创建缓冲区，再转换回WGS84处理日界线，避免直接在地理坐标系（EPSG:4326）下操作导致的球面几何可视化冲突。

核心思路 #

地理坐标系（WGS84）下的球面缓冲区在平面可视化时，跨日界线的几何会被底图的平面投影切割，导致显示异常。通过转换到无日界线问题的专用投影创建缓冲区，能保证几何形状的完整性，再转换回WGS84并切割日界线，就能正常可视化。

针对高纬度点的实现代码 #

你的两个案例都位于南极高纬度区域，推荐使用南极极地立体投影（EPSG:3031），该投影下无日界线限制，缓冲区为标准圆形，转换回WGS84后适配可视化：

案例1：坐标(178.4, -80.1) #

library(sf)
library(mapview)

sf_use_s2(TRUE)

# 创建点对象
coord <- c(178.4, -80.1)
point <- st_sfc(st_point(coord), crs = 4326)

# 转换到南极极地投影创建缓冲区
point_proj <- st_transform(point, crs = 3031)
buffer_proj <- st_buffer(point_proj, 2000000)  # 2000000米即2000公里

# 转换回WGS84并处理日界线
buffer <- st_transform(buffer_proj, crs = 4326) %>%
  st_wrap_dateline(options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180"))

# 可视化
mapview(buffer) + mapview(point)

案例2：坐标(78.4, -80.1) #

仅需修改坐标，其余逻辑完全一致：

library(sf)
library(mapview)

sf_use_s2(TRUE)

coord <- c(78.4, -80.1)
point <- st_sfc(st_point(coord), crs = 4326)

point_proj <- st_transform(point, crs = 3031)
buffer_proj <- st_buffer(point_proj, 2000000)

buffer <- st_transform(buffer_proj, crs = 4326) %>%
  st_wrap_dateline(options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180"))

mapview(buffer) + mapview(point)

通用适配建议 #

• 高纬度区域：北极点附近用EPSG:3413（北极极地立体投影），南极点附近用EPSG:3031。
• 中低纬度跨日界线点：可选用EPSG:4087（等距圆柱投影，无日界线切割）或对应区域的UTM投影，创建缓冲区后再转换回WGS84处理日界线。
• 避免直接在EPSG:4326下创建缓冲区后可视化，球面几何到平面底图的转换容易引发显示异常，st_shift_longitude()仅能调整经度范围，无法解决高纬度几何的完整性问题。

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