基于线交点创建多边形：通过四条线a、b、c、d生成围合区域多边形

要把你提供的四条线a、b、c、d围成的区域转换成多边形，可以借助sf包的拓扑处理函数实现，以下是完整可运行代码和步骤说明：

完整代码 #

library(sf)
library(dplyr)
library(ggplot2)

# 定义原始坐标数据
a <- data.frame (x  = c(1, 5, 9),
                 y = c(1, 0, 1))
b <- data.frame (x  = c(2, 2),
                 y = c(0, 5))
c <- data.frame (x  = c(1, 9),
                 y = c(4, 4))
d <- data.frame (x  = c(8, 8),
                 y = c(0, 5))

# 将所有线转换为LINESTRING并合并为单个几何集合
all_lines <- bind_rows(
  a %>% st_as_sf(coords = c("x", "y")) %>% summarize(do_union=FALSE) %>% st_cast("LINESTRING"),
  b %>% st_as_sf(coords = c("x", "y")) %>% summarize(do_union=FALSE) %>% st_cast("LINESTRING"),
  c %>% st_as_sf(coords = c("x", "y")) %>% summarize(do_union=FALSE) %>% st_cast("LINESTRING"),
  d %>% st_as_sf(coords = c("x", "y")) %>% summarize(do_union=FALSE) %>% st_cast("LINESTRING")
) %>% st_combine()

# 从线集合生成多边形
polygon <- st_polygonize(all_lines) %>% 
  st_cast("POLYGON") %>% 
  st_as_sf()

# 可视化验证：多边形填充色+原始线虚线标注
ggplot() +
  geom_sf(data = polygon, fill = "#87CEEB", alpha = 0.5) +
  geom_sf(data = all_lines, color = "darkred", linetype = "dashed") +
  theme_minimal()

关键步骤说明 #

1. 合并线对象：把四条独立的LINESTRING合并成一个几何集合，st_combine将多个线打包为整体，方便后续拓扑处理。
2. 生成多边形：st_polygonize()是核心函数，会自动识别线围成的闭合区域，转换为多边形几何。
3. 格式转换：用st_cast("POLYGON")确保输出为标准多边形类型，再转为sf对象，便于后续空间分析或可视化。

如果线存在微小间隙导致无法生成多边形，可先执行st_snap(all_lines, all_lines, tolerance = 0.01)对齐端点，再运行st_polygonize。

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