0

我正在尝试与 建立邻里关系spdep::poly2nb。我正在使用 77 个纽约市警察局,它们都分配了特定的号码:

警察局的身份证

当我创建邻域时,我得到了区域之间的邻域关系,但它们的编号是 1:77,而不是它们的实际编号。

我试图将警察辖区 ID 分配为 sf 或 poly2nb 的行名,但我总是得到相同的结果:

poly2nb 输出

我的代码两次分配行名,但仍然得到与上面相同的输出:

admin_data <- read_sf("*/Data/NY_PP/nypp.shp")
pp <- admin_data[,c(1,3,4)]
pp_sf <- pp  %>% st_as_sf(wkt="geometry", crs=2263, remove=FALSE)
row.names(pp_sf) <- as.character(pp_sf$Precinct)
bnb <- poly2nb(pp_sf, row.names = pp_sf$Precinct, queen = TRUE)

pp_sf

非常感谢任何帮助!欢呼 <3

4

0 回答 0