数据呈现 | 让文稿shinly起来！地图绘制

有时候，我们在写文稿时，可能会涉及到数据的“统计制图”这一个环节。比如我们会遇到这样一批数据集，绘制全国31个省人均GDP的空间分布特征，以探索不同地区的经济发展水平差异。

这时候，我们会不假思索地很快想到使用地图来解决问题。如果你是学过空间计量的学生，也许会考虑使用Arcgis专业地图绘制，但大部分非空间计量学生可能会望而却步。

那么，不熟悉Arcgis的学生，如何使用常见的计量软件进行地图绘制呢？以热力地图绘制为例——热力地图是以颜色来表现数据强弱大小及分布趋势的可视化类型，可应用于经济发展水平差异分析、人口密度分析、活跃度分析等——基于全国31个省份人均GDP数据，下面我们利用两种软件Stata和R分别进行介绍。

数据集

下表1是全国人均GDP数据集，其中procode是省份代码，prov是省份名称，y是人均GDP。

表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",clear
merge m:1 prov using prov.dta,nogen

**利用spmap命令绘图，请安装ssc install spmap
spmap y using `"province_co"',id(_ID) title ("中国地图") label(label(prov) xcoord(x_coord) ycoord(y_coord) size(*.66))  fcolor(Purples)

左右滑动查看更多

stata热力图形展示：

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热力图形展示：

获取更多R地图绘制（含动态地图），可查看：

https://blog.csdn.net/u013524655/article/details/72812181

值得说明的是，stata和R软件的两种统计制图方法所达到的预期效果相差不大。其中，地图的底层坐标数据，stata来源于外部（含九段线），而R依赖于eCharts。

注：stata绘制地图详细方法，请点击阅读原文并下载说明，附命令和步骤（电脑客户端可直接打开下载，手机客户端需用浏览器打开喔）。

►一周热文

数据Seminar

这里是大数据、分析技术与学术研究的三叉路口

作者：简华（何年华）审阅：杨奇明编辑：青酱

    欢迎扫描👇二维码添加关注