使用 Stata 绘制长三角 41 省市的区县地图

今天给大家分享一份使用 Stata 绘制长三角 41 省市的区县地图的数据和方法。

首先需要准备我特别设计的长三角 41 省市地图数据（包含指北针和比例尺）：

• jzmap 文件夹：地理矢量数据；
• labeldf.dta 数据：文本标签位置数据；

这份数据是使用 R 语言设计的，感兴趣的小伙伴可以参考附件中的 main.R 文件。

下面我们看一下如何在 Stata 中绘制长三角 41 省市的区县地图。

首先把 shp 文件转换成 dta 文件：

*- 把 shp 文件转换成 dta 文件
local name = "jzmap"
shp2dta using `name'/`name', database(`name'_db) coordinates(`name'_coord) genid(ID) gencentroids(centroids) replace
shp2dta using `name'/`name'_line, database(`name'_line_db) coordinates(`name'_line_coord) genid(ID) gencentroids(centroids) replace

这样就会得到四个文件：

• jzmap_db.dta
• jzmap_coord.dta
• jzmap_line_coord.dta
• jzmap_line_db.dta

其中前两个是底图的属性数据和坐标系数据，后面两个是线条的。

下面还需要对数据进行下调整。

首先是删除底图上的比例尺和指北针数据，因为这两个元素计划使用线条和多边形选项绘制：

use jzmap_db.dta, clear
* 删除比例尺和指北针
drop if index(县, "_")
save jzmap_db.dta, replace

再创建比例尺和指北针的多边形数据：

use jzmap_coord.dta, clear
*- 保留指北针和比例尺对应的观测值
keep if _ID > 306
gen value = 1
save jzmap_polyon, replace

再对线条进行分组：

*- 线条
use jzmap_line_db.dta, clear
*- 分组
gen group = 1
replace group = 2 if index(市, "边界")
replace group = 3 if index(省, "边界")
replace group = 4 if index(省, "_")
keep ID group
ren ID _ID
save temp, replace

use jzmap_line_coord.dta, clear
merge m:1 _ID using temp
drop _m
save jzmap_line_coord, replace

这里我把线条分为四类：

1. 区县边界；
2. 城市边界；
3. 省份边界；
4. 比例尺和指北针的线条。

后面的代码中也会根据不同的线条类型使用不同的颜色和粗细。

统计 2021 年每个区县的年均 PM2.5 浓度：

import excel using "1998~2021年中国各区县PM2_5面板数据.xlsx", clear first
keep if inlist(省, "长三角 41 省市省") & 年份 == 2021
drop 年份 sum min max sd
ren mean PM2_5
keep 县代码 PM2_5
save "PM2_5", replace

然后就可以绘制地图了：

此处代码需下载讲义材料查看～

点击这里跳转到 RStata 短书平台获取附件：使用 Stata 绘制长三角 41 省市的区县地图