Stata spmap问题：内嵌缺失数据多边形渲染错误

核心原因 #

Stata的spmap默认按数据集内多边形的顺序渲染，若内嵌的小多边形排在外层大多边形之前，会被后者完全遮盖，导致缺失数据的灰色样式无法显示。

具体解决办法 #

方法1：调整多边形绘制顺序 #

1. 查看地图数据集中多边形的面积或标识，确认内嵌小多边形的位置：

   list _ID _AREA  // _ID为多边形标识，_AREA为面积字段，按需替换
   

2. 按面积从小到大排序，让小多边形（面积更小）排在大多边形之后：

   sort _AREA
   

   若有明确嵌套关系，也可手动调整数据集顺序，确保内嵌多边形处于更靠后的位置。
3. 重新执行spmap命令，此时小多边形会在大多边形之上渲染，缺失数据的灰色样式可正常显示。

方法2：分层叠加绘制 #

将缺失数据的内嵌多边形单独作为顶层图层绘制：

1. 先绘制除目标内嵌县外的所有区域：

   spmap your_var using map_data, if _ID != "内嵌县ID" ///
       fcolor(...) ocolor(...) ///
       ndfcolor(gs14) ndocolor(gs2) ndlabel("No data")
   

2. 叠加绘制目标内嵌县，指定缺失数据的样式：

   spmap using map_data, if _ID == "内嵌县ID" ///
       fcolor(gs14) ocolor(gs2) label("No data") ///
       addplot
   

   addplot选项会将此图层叠加在已有地图上方，确保内嵌县的灰色样式不被遮盖。

方法3：修复地图拓扑结构 #

若地图数据存在拓扑错误（如大多边形未正确挖空内嵌区域），可重新导入shapefile格式的地图文件：

spshape2dta "your_shapefile.shp", replace data(map_data) coord(map_coord)

重新导入后再执行spmap命令，可能解决遮盖问题。

内容的提问来源于stack exchange，提问作者Kay Ru