使用 gCentroid 在 R 中查找多個多邊形的質心時保留 ID

Question

這似乎應該很簡單。 我在一個 shapefile 中包含多個多邊形，我正在使用 rgeos 的gCentroid() function 創建一堆質心點。 function 有一個id參數，如果未指定，它應該返回父幾何 ID，但這要么不起作用，要么我找錯了地方，或者我誤解了這個參數。

簡單的例子：

library(terra)
library(rgeos)
library(sp)

v <- vect(system.file("ex/lux.shp", package="terra"))
v <- as(v, "Spatial")

#Clearly there are IDs here (albeit not unique) 
v@data[["ID_1"]] 

所以當我在 go 上創建質心時

cents <- gCentroid(v, byid = TRUE)

我沒有看到任何關聯的“ID_1”插槽。 問題是我最終將使用這些質心從柵格中獲取值，並且需要 ID 來區分值來自哪個多邊形。

Answer 1

簡單的方法是

library(terra)
v <- vect(system.file("ex/lux.shp", package="terra"))
x <- centroids(v)

與sf

library(sf)
y <- st_centroid(st_as_sf(v))

使用 rgeos，您需要按照這些方式做一些事情

library(rgeos)
s <- as(v, "Spatial")
cents <- gCentroid(s, byid = TRUE)
atts <- data.frame(v)[as.integer(row.names(cents)), ]
s < SpatialPointsDataFrame(cents, atts)