我在使用R中的spdep
库来绘制Queen Contiguity Neighbors时,会对以下结果感到困惑。
library(spdep)
使用世界形状文件wrld_simp获取质心点(空间坐标)
data(wrld_simpl)
plot(wrld_simpl)
coord.w <- coordinates(wrld_simpl)
ID <- wrld_simpl$IDb
为了让女王邻居,我做
queen_nb <- poly2nb(wrld_simpl, row.names=ID, queen=TRUE)
然后我绘制形状文件
plot(wrld_simpl, border="grey60",xlab="",ylab="",xlim=NULL, axes=TRUE)
title(main="Queen contiguity neighbors")
和女王邻居:
plot(queen_nb, coord.w, add=TRUE, col="red", pch=19, cex=0.6)
你绘制它(或者你只是运行summary(queen_nb)
),你可以看到许多邻居没有连接(例如智利是孤立的)。
你有任何提示吗?