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

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

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

数据集

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

表1 全国人均GDP数据集2018（单位：元）

stata地图绘制——spmap

stata绘制热力图，我们主要用到spmap命令：

**导入GDP数据集import excel "C:Stata 14mapprov.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地图绘制（含动态地图），可查看：

CSDN-专业IT技术社区-登录​blog.csdn.net

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

文章首发于公众号 数据Seminar

数据呈现 | 让文稿shinly起来！地图绘制​mp.weixin.qq.com

更多可视化内容请关注公众号 数据Seminar