从数据框列表创建空间多边形数据框列表答案

【问题标题】：Create a list of Spatial Polygon dataframe from a list of dataframe从数据框列表创建空间多边形数据框列表
【发布时间】：2014-10-22 09:50:07
【问题描述】：

我想从数据框列表中创建空间多边形数据框列表。空间多边形列表称为list_sp_Tanzania，数据框列表称为list_df_Tanzania_Modis500。每个列表包含 61 个对象，每个对象包含几个多边形。

str(list_df_Tanzania_Modis500)
$ :'data.frame':    30 obs. of  11 variables:
  ..$ ID        : int [1:30] 296 298 321 323 324 330 331 361 419 453 ...
  ..$ LU_1990   : int [1:30] 11 11 11 11 11 11 11 11 11 11 ...
  ..$ LU_2000   : num [1:30] 12 12 12 12 12 12 12 12 12 12 ...
  ..$ CHLU_90_00: chr [1:30] "1112" "1112" "1112" "1112" ...
  ..$ LU_2005   : num [1:30] 12 12 12 12 12 12 12 12 12 15 ...
  ..$ CHLU_00_05: chr [1:30] "1212" "1212" "1212" "1212" ...
  ..$ Tile      : Factor w/ 1 level "S11_E039": 1 1 1 1 1 1 1 1 1 1 ...
  ..$ UNIQ_ID   : Factor w/ 30 levels "S11_E039_296",..: 1 2 3 4 5 6 7 8 9 10 ...
  ..$ AREA      : num [1:30] 219337 347133 393961 181875 105137 ...
  ..$ Sour_90_00: chr [1:30] "Modis500_2000" "Modis500_2000" "Modis500_2000" "Modis500_2000" ...
  ..$ Sour_00_05: chr [1:30] "Modis500_2005" "Modis500_2005" "Modis500_2005" "Modis500_2005" ...

str(list_sp_Tanzania)
[[61]]
class       : SpatialPolygons 
features    : 30 
extent      : 38.95413, 39.04577, -11.04522, -10.95469  (xmin, xmax, ymin, ymax)
coord. ref. : +proj=longlat +datum=WGS84 +ellps=WGS84 +towgs84=0,0,0

很遗憾，我无法通过dput 提供数据，因为它们是s4 objects，而且我找不到类似的示例。希望有人能以某种方式帮助我。

【问题讨论】：

我看不出你的两个列表是如何关联的，但基本上使用 30 行的 data.frame 和 30 个特征的 SpatialPolygons，你可以选择使用参数来控制 iD 匹配.如果“list_df_Tanzania_Modis500”是一个列表，您需要 df
[[61]] 部分是列表中的第 61 个空间多边形对象。它说 30 个特征，这意味着第 61 个空间多边形对象中有 30 个多边形。

标签： r dataframe spatial-data-frame

【解决方案1】：

这里是解决方案。

fun <- function(x, y) {
  SpatialPolygonsDataFrame(x, y, match.ID = F)
}

list_Spdf_Tanzania_Modis500 <- mapply(FUN = fun, 
                                      x = list_sp_Tanzania, 
                                      y = list_df_Tanzania_Modis500)

【讨论】：