如何高效修改sf中大规模LINESTRING对象的首尾坐标？

需求背景 #

我有一个包含约10万行LINESTRING几何类型的大型sf对象，需要修改每个LINESTRING的首尾坐标点。目前的实现方法是先将线串转为点、修改首尾点后再转回线串，但这种步骤冗余且计算量高，希望找到更直接的实现方式。

示例数据与目标效果 #

先定义示例数据和预期的修改效果：

library(sf)
library(tidyverse)

# 示例数据：生成一条LINESTRING
mypoints <- data.frame(id = 1, lon = 1:5, lat = 1:5)
line1 <- mypoints %>% 
  st_as_sf(coords = c("lon", "lat"), crs = 4326) %>%
  group_by(id) %>% 
  summarise() %>% 
  st_cast("LINESTRING") 

# 预期修改效果：将首坐标(1,1)改为(1,2)，尾坐标(5,5)改为(5,4)
mypoints2 <- data.frame(id = 1, lon = c(1:5), lat = c(2,2:4,4))
line2 <- mypoints2 %>% 
  st_as_sf(coords = c("lon", "lat"), crs = 4326) %>%
  group_by(id) %>% 
  summarise() %>% 
  st_cast("LINESTRING")

# 可视化对比：原线（蓝） vs 修改后线（红）
ggplot() + 
  geom_sf(data = line1, color = "blue", linewidth = 6) +
  geom_sf(data = line2, color = "red", linewidth = 3)

优化后的直接修改方法 #

无需转点再转回线串，直接操作LINESTRING的坐标矩阵：

# 提取原线串的坐标矩阵
coords <- st_coordinates(line1)

# 修改首尾坐标
coords[1, c("X", "Y")] <- c(1, 2)  # 修改首点
coords[nrow(coords), c("X", "Y")] <- c(5, 4)  # 修改尾点

# 重新构建LINESTRING对象
line_modified <- st_sf(
  id = line1$id,
  geometry = st_sfc(st_linestring(coords), crs = st_crs(line1))
)

# 查看结果
line_modified
# Simple feature collection with 1 feature and 1 field
# Geometry type: LINESTRING
# Dimension:     XY
# Bounding box:  xmin: 1 ymin: 2 xmax: 5 ymax: 4
# Geodetic CRS:  WGS 84
# # A tibble: 1 × 2
#      id                  geometry
#   <dbl>          <LINESTRING [°]>
# 1     1 (1 2, 2 2, 3 3, 4 4, 5 4)

批量处理多线串的方法 #

如果是包含多个LINESTRING的sf对象（比如10万行），可以用purrr批量处理：

# 生成多条LINESTRING的示例数据
multi_lines <- map_dfr(1:3, function(i) {
  data.frame(id = i, lon = 1:5 + i, lat = 1:5 + i) %>%
    st_as_sf(coords = c("lon", "lat"), crs = 4326) %>%
    group_by(id) %>%
    summarise() %>%
    st_cast("LINESTRING")
})

# 定义每条线的首尾修改规则（这里示例为每个id对应不同的首尾点）
new_points <- tibble(
  id = 1:3,
  start_x = 1:3 + 1,
  start_y = 1:3 + 2,
  end_x = 5:7,
  end_y = 5:7 - 1
)

# 批量修改
multi_lines_modified <- multi_lines %>%
  left_join(new_points, by = "id") %>%
  mutate(
    geometry = pmap(list(geometry, start_x, start_y, end_x, end_y), function(line, sx, sy, ex, ey) {
      coords <- st_coordinates(line)
      coords[1, c("X", "Y")] <- c(sx, sy)
      coords[nrow(coords), c("X", "Y")] <- c(ex, ey)
      st_linestring(coords)
    }) %>% st_sfc(crs = st_crs(multi_lines))
  ) %>%
  select(id, geometry)

# 查看结果
multi_lines_modified

方法优势 #

• 避免了st_cast转点再转回线串的冗余操作，减少了内存占用和计算时间
• 直接操作坐标矩阵，逻辑清晰，适合处理大型数据集

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