R语言制作长三角城市群空间权重矩阵及作图显示

长三角一体化进入新的发展阶段，下面通过制作最新长三角城市群地图，并运用R语言形成空间权重矩阵，同时画出了空间连接权重矩阵图示。

library("raster")
library("spdep")
#p1<- shapefile("E:\new_delta_map\\new_delta_map_2.shp")
p=shapefile("E:\\new_delta_map\\deltacity_2.shp",encoding="UTF-8")
cityname=c('安庆市',    '池州市',    '滁州市',    '宣城市',    '常州市',    '南京市',    '南通市',    '苏州市',    '泰州市',    '无锡市',    '盐城市',    '扬州市',    '镇江市',    '杭州市',    '湖州市',    '嘉兴市',    '金华市',    '宁波市',    '绍兴市',    '台州市',    '舟山市',    '马鞍山市',   '芜湖市',    '合肥市',    '上海市',    '铜陵市')
p$count_name=citynamelatlon=read.csv("\\data\\zuobiao.csv")
latlon=as.data.frame(latlon)
xy[,1]=latlon[,2]#换成城市实际坐标
xy[,2]=latlon[,3]
###############
wr <- poly2nb(p, row.names=p$ID_1, queen=TRUE)
wr[[18]]=sort(as.integer(c(wr[[18]],21)))
wr[[21]]=sort(as.integer(c(18)))
wm <- nb2mat(wr, style='B', zero.policy = TRUE)
write.csv(wm,"\\new_delta_map\\deltacitycluster-1.csv") #相邻矩阵##############
#Distance based:
wd100 <- dnearneigh(xy, 0, 100, longlat=TRUE)
wmwd100<- nb2mat(wd100, style='B', zero.policy = TRUE)
write.csv(wmwd100,"\\new_delta_map\\deltacitycluster-wd100.csv") #相距离100km矩阵
#Nearest neighbors:
k3 <- knn2nb(knearneigh(xy, k=3, RANN=FALSE))
wmk3<- nb2mat(k3, style='B', zero.policy = TRUE)
write.csv(wmk3,"\\new_delta_map\\deltacitycluster-k3.csv") #相邻3个矩阵
plot(p, col='white', border='blue')
plot(wr, xy, col='red', lwd=2, add=TRUE)#作图，以相邻矩阵为连接
text(xy,p$count_name,cex=.7, halo=TRUE, hc="white", col='blue', hw=0.0)#halo是晕圈，设置hw=0则去除晕圈

以100km为门槛建立空间权重矩阵的图示如下：

以相邻与否（wij=1，if adj=1）建立空间权重矩阵的图示如下：

以三个相近邻居（k=3）建立空间权重矩阵的图示如下：

R语言制作长三角城市群空间权重矩阵及作图显示相关推荐

手把手教你使用R语言制作临床决策曲线
DCA(Decision Curve Analysis)临床决策曲线是一种用于评价诊断模型诊断准确性的方法,在2006年由AndrewVickers博士创建,我们通常判断一个疾病喜欢使用ROC曲线的A ...
r语言 html gif,推荐：6步教你用R语言制作动图
原标题:推荐 :6步教你用R语言制作动图引言数据可视化可能是数据科学领域最重要却通常最少被提及的部分. 我这么说是因为创造数据故事和可视化对你的顾客最终怎么看待你的工作有巨大的影响.数据科学不仅仅 ...
R语言----制作数据分布图（直方图+概率密度曲线）
R语言--制作group1和group2的分布图并输出tiff格式300dpi保存制作分布图我选择的是制作直方图+概率密度曲线~ library(ggplot2) data<-read.csv ...
R语言 REmap安装以及城市信息获取失败问题
R语言 REmap安装以及城市信息获取失败问题 (一)针对REmap安装问题首先进行如下四步 1.install.packages("devtools") 2.library(d ...
实例演示R语言制作限制性立方条图
限制性立方样条函数(RCS)在比较非线性关系中很常用.既往我们已经讲过R语言制作限制性立方条图,但是讲得比较简单,中间有些环节没写出来,我也不是很满意,今天重新来说一下.主要是要用到rms包的rcs函 ...
长三角城市数字经济发展水平评估白皮书
我国高度重视数字经济发展.十九届五中全会.<中华人民共和国国民经济和社会发展第十四个五年规划和2035年远景目标纲要>强调,要推动数字经济和实体经济深度融合,加快数字化发展,建设数字中国. ...
R语言ggplot2可视化：使用ggplot2按照热力图的方式显示全是分类变量的dataframe数据、并自定义因子（factor）的图例颜色legend
R语言ggplot2可视化:使用ggplot2按照热力图的方式显示全是分类变量的dataframe数据.并自定义因子(factor)的图例颜色legend(use heatmap to visuali ...
R语言ggplot2可视化将图像标题（title）居中显示实战
R语言ggplot2可视化将图像标题(title)居中显示实战目录 R语言ggplot2可视化将图像标题(title)居中显示实战
利用R语言制作好看的Meta分析文献偏倚风险图
在<Meta踩坑集>里,向大家分享了使用RevMan进行基本的文献偏移风险图是如何制作的.在很多文献里,我们看到的偏倚风险图都长这样: 但也有的文献偏移风险图长这样: 那这些花里胡哨,看起 ...
R语言制作Meta分析偏倚风险评估（ROB）图
我们进行偏倚风险评估一般都是用Cochrane的偏倚风险评估工具,在Revman里面会做出下面的图. 这个图在R语言里面完全可以重现,下面就给大家演示怎么用R来制作ROB图. 首先准备数据,在Exce ...

R语言制作长三角城市群空间权重矩阵及作图显示

R语言制作长三角城市群空间权重矩阵及作图显示相关推荐

最新文章

热门文章