如何将bordering_polygons外边界扩展至eez线且避免重叠？

数据对象 #

• eez：RDS格式空间线数据
• bordering_polygons：RData格式面数据

需求目标 #

将每个bordering_polygons的外边界扩展至eez线，且扩展后的边界不重叠。

尝试的代码 #

v <- vect(bordering_polygons)
r <- rast(eez, res=.1)
# 计算每个栅格到各多边形的距离，找到最近的多边形
rr <- lapply(1:nrow(v), \(i) distance(r, v[i])) |> rast() |> which.min()

z <- erase(as.polygons(rr), v)
values(z) <- values(v)[unlist(z[[1]]),  ]
x <- rbind(v, z)
a <- aggregate(x, names(x))
plot(v)
plot(a, "ADM1_EN", type="classes")
lines(v, lwd=3, col="gray")
lines(eez, col="red", lwd=3)

结果可视化 #

当前问题 #

生成的扩展区域未与现有多边形合并。

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