Stata绘制中国行政色彩分级地图

[唐基老爹]

在使用stata进行空间计量分析的时候，发现说它也可以用来绘制色彩分级地图（区域经济学和地理学中经常需要用到，可以将某一些区域性的数量指标在地图上变现出来，尤其是分行政区的地图），于是我就摸索使用stata绘制了中国广告产业积聚度的全国色彩分级地图。
主要使用到的命令有：spmap、shp2dta、mif2dta，如果没有这些命令就需要先安装。主要命令如下：
clear
*安装这几个必要的命令
*DEA-CCR原模型matlab模拟（非对偶形式）
*@经济地理学笔记，微信公众号分享空间计量经济学相关的模型与方法
*交流邮箱：[e-mail]
*参考:Stata之禅
ssc install spmap
ssc install shp2dta
ssc install mif2dta
*下载中国行政区域地图
findit china_map
*解压压缩文件，这一步我出现过命令不能别的的情况，后将路劲全部设置成英文就解决了
unzipfile china_map,replace
*转码成可识别的GB18030，这样方便后文合并数据
unicode encoding set gb18030
unicode translate china_label.dta
use china_label.dta,clear
*新生成一列省份的变量，删除行政区名中多余的字符，方便匹配
gen province= ustrregexra(name,"省|市|维吾尔|回族|壮族|自治区|特别行政区","")
save china_province.dta, replace
use advertising.dta, clear
*/按照province进行匹配
merge 1:m province using china_province.dta
*/列出未匹配的
list if _merge==2
spmap agglo using "china_map.dta", id(id)
最后的结果是这样的：


此外，还可以根据需要进行标题、颜色、边界等设定——