使用SpatVector创建独立线段时生成空对象的问题

问题原因 #

1. 语法错误导致数据生成异常：你的代码中dupz <- dupz[-c(1,nrow(dupz)]缺少右括号，这会导致代码执行失败，无法正确生成目标数据框dupz，后续创建的线对象自然不符合预期。
2. vect()参数未明确指定，列解析混乱：传入dupz时包含hole列，但未明确告知vect()哪些列是几何相关字段（id、part、x、y），哪些是属性字段，导致属性列未被正确识别，几何分组也出现错误。
3. 几何要素计数不符预期：按逻辑应生成15条独立线段（原16个点，相邻点组成15条线段），但输出显示16个几何，本质是语法错误导致的数据预处理异常。

解决步骤 #

步骤1：修正语法错误，生成正确的线段数据 #

先修正括号问题，重新生成符合要求的dupz：

# 原参考代码
x1 <- rbind(c(-180,-20), c(-140,55), c(10, 0), c(-140,-60))
x2 <- rbind(c(-10,0), c(140,60), c(160,0), c(140,-55))
x3 <- rbind(c(-125,0), c(0,60), c(40,5), c(15,-45))
hole <- rbind(c(80,0), c(105,13), c(120,2), c(105,-13))
z <- rbind(cbind(object=1, part=1, x1, hole=0), cbind(object=2, part=1, x3, hole=0),
           cbind(object=3, part=1, x2, hole=0), cbind(object=3, part=1, hole, hole=1))
colnames(z)[3:4] <- c('x', 'y')
z[z[, "hole"]==1, "object"] <- 4

# 修正后的线段数据生成代码
z[,1] <- 1:nrow(z)
duptimes <- rep(2,nrow(z))
idx <- rep(1:nrow(z), duptimes)
dupz <- z[idx,]
dupz <- dupz[-c(1,nrow(dupz))]  # 修复括号缺失问题
dupz[,1] <- rep(1:(nrow(z)-1), rep(2, nrow(z)-1))
dupz[,2] <- 1  # 每个线段仅一个部分，统一设为1

步骤2：明确参数创建SpatVector #

调用vect()时明确指定几何相关列，确保属性和几何被正确解析：

# 明确指定id、part、坐标列，剩余列作为属性
river <- vect(dupz, type = "lines", id = "object", part = "part", x = "x", y = "y")

# 查看结果
river
plot(river)

步骤3：验证结果 #

执行后river的输出应类似：

class       : SpatVector 
geometry    : lines 
dimensions  : 15, 1  (geometries, attributes)
extent      : -180, 160, -60, 60  (xmin, xmax, ymin, ymax)
coord. ref. :  
names       : hole
type        : num 
values      : 0, 0, 0, 0, 0, ...

此时plot(river)会正确绘制所有15条独立线段。

补充说明 #

• 每条独立线段对应唯一的object（id），仅包含两个相邻点，完全匹配“连续点间独立线段”的需求。
• 明确指定vect()参数可避免列解析歧义，确保hole这类属性列被正确保留。

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