R中st_buffer的flat端帽样式无法生效问题求助

使用经纬度坐标创建线串并绘制缓冲区时，设置endCapStyle="FLAT"后，缓冲区仍超出线段端点，参数未按预期生效。

问题原因 #

你当前使用的是地理坐标系（EPSG:4326，WGS84），该坐标系下st_buffer会采用球面缓冲区算法，此时endCapStyle这类针对平面缓冲区的参数会被忽略。只有在平面投影坐标系（如UTM，单位为米）下，endCapStyle等参数才能正常发挥作用。

解决方案 #

先将线串转换为合适的平面投影坐标系（示例中坐标位于美国科罗拉多州，对应UTM 13N，EPSG:32613，需根据你的坐标区域调整），执行缓冲区操作后，再转回地理坐标系。

修改后的代码如下：

library(sf)
library(mapview)
library(tidyverse)

# 示例坐标
df_tbl <- tibble::tribble(
  ~lat,   ~lon,   ~lat2,   ~lon2,
  39.977917887303626, -105.28350982984472, 40.026678120413961, -105.20312377997453,
  39.934422269952776, -105.25986687400055, 39.902061443935025, -105.28605599432022)

# 创建起点和终点
beg_t <- df_tbl %>% select(lon, lat) %>% st_as_sf(coords = c("lon", "lat"), crs = 4326)  
end_t <- df_tbl %>% select(lon2, lat2) %>% st_as_sf(coords = c("lon2", "lat2"), crs = 4326)  

# 合并点、创建线串、转换坐标系并生成FLAT端点缓冲区
line_segments <- cbind(beg_t, end_t) %>% 
  mutate(
    path = st_sfc(purrr::map2(
      .x = geometry, 
      .y = geometry.1, 
      .f = ~ st_union(c(.x, .y)) %>% st_cast("LINESTRING")
    ), crs = 4326)) %>%
  # 转换为UTM平面坐标系（根据你的区域调整EPSG代码）
  st_transform(crs = 32613) %>%
  # 在平面坐标系下生成FLAT端点的缓冲区
  mutate(buffer = st_buffer(path, dist = 2000, endCapStyle = "FLAT")) %>%
  # 转回地理坐标系以便可视化
  st_transform(crs = 4326)

# 查看结果
mapview(line_segments$path) + mapview(line_segments$buffer)

关键说明 #

• 坐标系转换：必须切换到平面投影坐标系，地理坐标系的球面缓冲区逻辑不支持endCapStyle参数。
• UTM分带选择：UTM分带基于经度，需根据你的坐标所在区域选择对应EPSG代码（可通过st_transform自动匹配或手动查询）。
• 单位匹配：dist参数单位需与投影坐标系一致，UTM单位为米，示例中dist=2000即表示2000米的缓冲区宽度。

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