stata绘制地理分布图的spmap命令

前言

首先感谢连玉军老师的连享会，如果读者对spmap命令一点基础都没有的话,请首先阅读连享会的推文：
Stata：空间计量之用 spmap 绘制地图

本文主要是分享下我个人使用spmap绘制地图的一点点经验,主要是如何处理中国地图和当绘制城市层面的数据时如何让读者更清晰地识别城市。

如何同时清晰地展示省界和市界

反面例子

先上一个反面例子，下图中使用的数据是2005年城市级别的人口密度，白色区域表示数据缺失。可以看到这样的图片无法直观的区分大部分城市，对于读者是非常不友好的。我们可以通过加入省界的方式来区分省级区域，从而让读者比较好的识别具体城市位置。
这里使用的就是最基本的命令：

spmap 人口密度 using "city_map2",id(id)

其中city_map2就是地级市级别的coord文件。

进一步改善：加入省界

如下图所示，改善后的图片中地级市界限以灰色显示，省级界限以黑色显示（当然那些缺失数据的地区界限还需要进一步改善）。
该图使用的命令如下：

spmap 人口密度 using "