创建空空间线对象答案

【问题标题】：Create empty spatial lines object创建空空间线对象
【发布时间】：2018-03-25 06:55:29
【问题描述】：

我正在尝试创建一个空的 SpatialLines 对象。使用多边形很容易：

SpatialPolygons(list())

对于空间线，这不起作用：

SpatialLines(LinesList = list())
Error in bb[1, ] : incorrect number of dimensions
SpatialLines(LinesList = Lines(list(),ID = "a"))
Error in as.list.default(X) : 
  no method for coercing this S4 class to a vector
SpatialLines(LinesList = Lines(slinelist = Line(coords = cbind(x = c(), y = c())), ID = c()))
Error in (function (classes, fdef, mtable)  : 
  unable to find an inherited method for function ‘coordinates’ for signature ‘"NULL"’

有人知道我如何创建一个空的 SpatialLines 对象吗？

解决方法

我找到了一个可能不是最好的解决方法。我生成了一条没有长度的空间线：

SpatialLines(list(Lines(Line(coords = cbind(x = c(0,0), y = c(0,0))), ID = "A")))

【问题讨论】：

标签： r spatial lines sp empty-list

【解决方案1】：

有趣的问题！

我能够解决的唯一方法是创建一个虚拟线并像这样删除它：

sl <- SpatialLines(LinesList = list(Lines(Line(matrix(0, ncol = 2)), ID = NA)))
sl <- sl[0]
length(sl)
# [1] 0

添加虚拟线时，长度按预期返回为 1：

length(rbind.SpatialLines(sl, SpatialLines(list(Lines(Line(coords = cbind(x = c(0,0), 
                                                                          y = c(0,0))), 
                                                           ID = "A")))))
# [1] 1

【讨论】：