如何限制st_buffer经度范围，避免跨180°经线生成无效多边形？

针对端点位于180°经线的线要素使用st_buffer生成缓冲区时跨越日界线的问题，这里提供两种可靠的解决方案，避免直接修改坐标值导致的拓扑错误：

方案1：日界线包裹+东半球裁剪 #

核心思路是先确保线要素在日界线处正确渲染，再用东半球边界多边形截断缓冲区：

library(sf)
library(dplyr)

# 原始线要素（你的代码）
LineA<-st_linestring(matrix(c(179,179.2,180,66.5,66.2,66),3))%>%
  st_sfc()%>%
  st_as_sf()%>%
  st_set_crs(4326)

LineB<-st_linestring(matrix(c(179,179.8,180,66.9,66.4,66.5),3))%>%
  st_sfc()%>%
  st_as_sf()%>%
  st_set_crs(4326)

Lines<-bind_rows(LineA,LineB)

# 1. 处理日界线包裹，避免线要素跨日界线错误计算
Lines_wrapped <- Lines %>%
  st_wrap_dateline(options = c("WRAPDATELINE=YES", "DATELINEOFFSET=180"))

# 2. 创建东半球边界多边形（经度0-180，纬度-90-90）
east_hemisphere <- st_polygon(list(matrix(c(180, 90,
                                             180, -90,
                                             0, -90,
                                             0, 90,
                                             180, 90),
                                           ncol=2, byrow=TRUE))) %>%
  st_sfc(crs = 4326)

# 3. 生成缓冲区并裁剪到东半球
Buffer_fixed <- Lines_wrapped %>%
  st_buffer(500, endCapStyle='FLAT') %>%
  st_intersection(east_hemisphere)

原理说明 #

• st_wrap_dateline会调整跨180°经线的要素，确保几何图形在正确的半球范围内呈现，避免st_buffer计算时绕到地球另一侧。
• 用东半球多边形做交集操作，是基于拓扑的裁剪，不会破坏多边形的有效性，完美截断缓冲区在180°经线处。

方案2：高纬度投影转换处理 #

由于你的线要素位于高纬度（66°左右），转换到极地投影后再做缓冲，能从根源避免日界线问题，同时提升缓冲精度：

# 转换到北极极地立体投影（EPSG:3995）处理
Buffer_proj <- Lines %>%
  st_transform(3995) %>%  # 适合北极区域的平面投影
  st_buffer(500, endCapStyle='FLAT') %>%
  st_transform(4326) %>%  # 转换回WGS84
  st_intersection(east_hemisphere)  # 可选：确保缓冲区不跨180°

原理说明 #

• 极地投影下，180°经线是一条直线，缓冲计算不会跨越到地球另一侧；同时高纬度区域的平面缓冲精度比球面缓冲更高。
• 转换回WGS84后再裁剪，确保最终结果符合后续栅格统计需求。

关于你之前的尝试 #

直接修改坐标值（比如把-179.xxx改为179.9999999）容易破坏多边形的拓扑结构（比如环的闭合性、自相交），导致后续栅格统计出错。上述两种方案都是基于sf的拓扑操作，能保证几何有效性，更适合后续分析。

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