Question

如何在SpatialPolygonsDataFrame中缓冲SpatialPolygons，并从其属性中获取每个SpatialPolygon的特定宽度。

下面的示例代码取自SpatialPolygons - Creating a set of polygons in R from coordinates，并稍加修改：

square <- t(replicate(50, {
  o <- runif(2)
  c(o, o + c(0, 0.1), o + 0.1, o + c(0.1, 0), o)
}))
ID <- paste0('sq', seq_len(nrow(square)))

# Create SP
polys <- SpatialPolygons(mapply(function(poly, id) {
  xy <- matrix(poly, ncol=2, byrow=TRUE)
  Polygons(list(Polygon(xy)), ID=id)
}, split(square, row(square)), ID))

# Create SPDF
polys.df <- (polys, data.frame(id=ID, row.names=ID))
polys.df$radius <- runif(50, 2.0, 8)

> head(polys.df)
      radius
sq1 7.397579
sq2 5.123130
sq3 4.607719
sq4 4.052619
sq5 6.571260
sq6 3.912436

每个SpatialPolygon应该使用其列“ radius”的值进行缓冲。

尝试以下循环将仅产生一个SpatialPolygon作为输出

for (i in length(polys.df)) {
  polys.df.buffer <- gBuffer(polys.df[i,], width=polys.df$radius[i])
}

以下循环返回所有缓冲区具有相同缓冲区宽度的多边形，但没有为“半径”列中指定的每个多边形使用定义的缓冲区宽度的缓冲区：

for (i in length(CrownsOut)) {
  polys.df.buffer <- gBuffer(polys.df[ ], width=polys.df$radius[i], byid=TRUE)
}

Answer 1

该代码几乎是正确的。您需要在for循环之前初始化一个列表，并在每次迭代中分配缓冲的Polygon。您还需要在循环条件中添加1:length(polys.df)，否则它将仅获取最后一个多边形，并且需要为每个多边形分配一个 ID 。

library(rgeos)

spList = vector("list", length(polys.df))
for (i in 1:length(polys.df)) {
  a <- gBuffer(polys.df[i,], width = polys.df$radius[i])
  a$id = i
  spList[[i]] <- a
}

spListAll <- do.call("rbind", spList)

plot(spListAll)

R-在循环中按属性缓冲SpatialPolygons

1 个答案: