R语言计算moran’I

最新推荐文章于 2022-11-20 07:15:26 发布

原创

最新推荐文章于 2022-11-20 07:15:26 发布 · 9.9k 阅读

47 ·

CC 4.0 BY-SA版权

文章标签：

#R语言 #moran‘I #邻接矩阵 #画地图

R语言计算moran’I

install.packages("maptools")#画地图的包
install.packages("spdep")#空间统计，moran'I
install.packages("tripack")
install.packages("RANN")
library("maptools")
library("spdep")
library("tripack")
library("RANN")

#读空间数据
rdata <- readShapePoly("D:/data/cairo.shp")#读取shp面数据
names(rdata)#显示数据字段名字
head(rdata@data)#显示属性表

#按照字段两幅画图
spplot(rdata[c("TFR96_03","TFR86_03") ],
       main = "spatial distribute of TFR96_03",#图名
       xlab = "X Coords",#横坐标名字
       ylab = "Y Coords",#纵坐标名字
       cut = 30#分段
       )

#按照字段一幅画图
spplot(rdata["TFR96_03"],
       main = "spatial distribute of TFR96_03",#图名
       xlab = "X Coords",#横坐标名字
       ylab = "Y Coords",#纵坐标名字
       cut = 30#分段
       )

#按边邻接角邻接生成邻居（方式一）
queen_nb <- poly2nb(rdata, queen = TRUE)#有8个
rook_nb <- poly2nb(rdata, queen = FALSE)#有4个

#获取中心点坐标编号
coords <- coordinates(rdata)

#地图数据转数据框-->生成ID
IDs <- row.names(as.data.frame(rdata))

#设置画参数
oopar <- par(mfrow = c(1,2