Question

我目前正在尝试从列表列表中创建矩阵列表。列表列表的组成部分是多边形，即五个顶点中每个顶点的顶点坐标：

$HEW18
     [,1]    [,2]
[1,] 595165.9 5688133
[2,] 595065.9 5688133
[3,] 595065.9 5688233
[4,] 595165.9 5688233
[5,] 595165.9 5688133

$HEW19
     [,1]    [,2]
 [1,] 593512.7 5672780
 [2,] 593412.7 5672781
 [3,] 593414.0 5672881
 [4,] 593514.0 5672880
 [5,] 593512.7 5672780

此文件中总共有150个多边形

现在，为了创建空间多边形数据框，我必须从“内部”列表或矩阵列表中创建矩阵。

我设法对一个多边形进行了如下操作：

#turn list to matrix
matrixHEW39 = matrix(c(forest$HEW39[1:5,]), nc=2, byrow = FALSE)

我现在想做的是节省一些时间，而不是单独为所有150个多边形编写此代码，而是为其创建一个for循环。

以下代码代表了我未定义的“想法”：

forest_list <-list()

for (i in forest) {
   matrix[i] = matrix(c(forest$i[1:5], nc =2, byrow = FALSE)
   forest_list[i]<-matrix[i]
}

您看到的是这里的问题，或者其中至少一个是“ forest $ i [1：5]”

如何像单个多边形代码中那样在每行1 -5内选择一个多边形？另外，我需要维护每个多边形的名称吗？

Answer 1

我认为您在这里有些困惑（学习R！时很容易发生）。这是基于看到您使用forest_list[i]<-matrix[i]。 Lists可以使用单括号[i]，双括号[[i]]或您曾经使用的美元名称约定来访问：forest$HEW39。单括号包含一个“切片”，该切片返回数据，但它仍然是一个列表。双括号返回对象在列表位置'i'，与dollar-name返回具有该名称的列表项相同。注意，您也可以使用[["HEW39"]]。

我认为这是混乱的根源。我认为您没有列表列表，只有普通列表，但是尝试使用单括号访问此列表会返回矩阵的列表（长度为1）。因此，它看起来像一个列表列表。当然，您尝试从中创建一个矩阵。但是矩阵在那里，您只需要使用双括号将其访问。

为回答您的问题，我列出了与您相似的5个矩阵，并将它们命名为：

forest <- lapply(c(1,2,3,4,5), function(f){replicate(2, rnorm(5))})
names(forest) <- paste0("HEW", 1:5)

NB：

class(forest[1])
[1] "list"
class(forest[[1]])
[1] "matrix"
class(forest$HEW1)
[1] "matrix"

要为这些内容建立索引，我们可以使用for循环，并捕获到新列表中

forest_list <- as.list(names(forest))

for (i in 1:length(forest)) {
    matrixi <- forest[[i]]
    forest_list[[i]] <- matrixi
}

如果您确实有一个列表列表，例如：

forest_lol <- list(1)
forest_lol[[1]] <- forest_list

您可以遍历所需的主列表（例如[[1]]）

forest_list <- as.list(names(forest))

for (i in 1:length(forest_lol[[1]])) {
    matrixi <- forest_lol[[1]][[i]]
    forest_list[[i]] <- matrixi
}

在新列表中保留名称仅意味着在for循环期间将其捕获到向量中，或更简单地根据前一个列表进行命名：

name_vector <- c()
for (i in 1:length(forest)) {
    name_vector <- c(name_vector, names(forest)[i])
}

遍历列表中的列表

1 个答案: