我正在尝试与 建立邻里关系spdep::poly2nb
。我正在使用 77 个纽约市警察局,它们都分配了特定的号码:
当我创建邻域时,我得到了区域之间的邻域关系,但它们的编号是 1:77,而不是它们的实际编号。
我试图将警察辖区 ID 分配为 sf 或 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)
非常感谢任何帮助!欢呼 <3