其他
数据呈现 | 让文稿shinly起来！地图绘制
有时候，我们在写文稿时，可能会涉及到数据的“统计制图”这一个环节。比如我们会遇到这样一批数据集，绘制全国31个省人均GDP的空间分布特征，以探索不同地区的经济发展水平差异。
这时候，我们会不假思索地很快想到使用地图来解决问题。如果你是学过空间计量的学生，也许会考虑使用Arcgis专业地图绘制，但大部分非空间计量学生可能会望而却步。
那么，不熟悉Arcgis的学生，如何使用常见的计量软件进行地图绘制呢？以热力地图绘制为例——热力地图是以颜色来表现数据强弱大小及分布趋势的可视化类型，可应用于经济发展水平差异分析、人口密度分析、活跃度分析等——基于全国31个省份人均GDP数据，下面我们利用两种软件Stata和R分别进行介绍。
数据集
表1 全国人均GDP数据集2018（单位：元）
| procode | prov | y |
|---|---|---|
| 110000 | 北京 | 140211 |
| 120000 | 天津 | 120711 |
| 130000 | 河北 | 47772 |
| 140000 | 山西 | 45328 |
| 150000 | 内蒙古 | 68302 |
| ... | ... | ... |
| 610000 | 陕西 | 63477 |
| 620000 | 甘肃 | 31336 |
| 630000 | 青海 | 47689 |
| 640000 | 宁夏 | 54094 |
| 650000 | 新疆 | 49475 |
左右滑动查看更多
stata地图绘制——spmap
stata绘制热力图，我们主要用到spmap命令：
**导入GDP数据集import excel "C:\Stata 14\map\prov.xlsx", sheet("Sheet1") firstrow save prov.dta,replace
**打开地图经纬度信息数据，并融合GDP数据集use "province_data.dta",clearmerge m:1 prov using prov.dta,nogen
**利用spmap命令绘图，请安装ssc install spmapspmap y using `"province_co"',id(_ID) title ("中国地图") label(label(prov) xcoord(x_coord) ycoord(y_coord) size(*.66)) fcolor(Purples)R地图绘制——REmap
R绘制热力图，我们主要用到REmap包：
#加载包library(openxlsx)library(REmap)
#导入数据pdata<-read.xlsx('prov.xlsx',1)
##绘制地图pdata<-pdata[,c(2,3)]
ph1<-remapC(pdata,maptype = 'china',color = c('gray','white'),##颜色设置theme = get_theme('none', lineColor='black', backgroundColor='white', titleColor='black', borderColor = "black"),title = "",subtitle = "",markPointTheme = markPointControl(),mindata = 30000,maxdata = 150000)
plot(ph1)左右滑动查看更多
获取更多R地图绘制（含动态地图），可查看：
https://blog.csdn.net/u013524655/article/details/72812181
►一周热文
数据Seminar
这里是大数据、分析技术与学术研究的三叉路口
欢迎扫描👇二维码添加关注