中国行政区划数据下载

GADM 数据库

GADM，是一个全球行政区划数据库。包括了几乎全部国家和地区的行政区划界线数据。

下载数据

GADM 提供的数据格式有 6 种：ESRI file geodatabase、Shapefile、Geopackage(SpatiaLite)、R(SpatialPolygonsDataFrame)、Google Earth .kmz 和 ESRI personal geodatabase。这里，数据格式应选择 shapefile。

在 GADM 中，country 的定义为 “any entity with an ISO country code”。关于 ISO country code，可以参考维基百科。因而想要下载完整的中国数据，实际上需要四个文件：China、Hong Kong、Macao、Taiwan。

也可以在 http://www.gadm.org/version2 下载全球的行政区划数据，但非常不推荐。一方面是数据量偏大，另一方面是全球区划数据保存到一个文件中，难以整理。不如麻烦一点，需要哪个国家下哪个国家。

解压数据

以 china 数据为例

unzip CHN_adm_shp.zip -d CHN_adm

解压后得到一堆文件，其中 CHN_adm0.shp、CHN_adm1.shp、CHN_adm2.shp、 CHN_adm3.shp 为实际需要的 shapefile 数据。0、1、2、3 为第零、一、二、三级行政区划界线。对于 china 数据，基本相当于中国大陆地区的边界、省界、市界和区县界。如果是 Hong Kong 的数据，0 依然是边界，但是 1 理解为省界显然不合适，总之就是更低一级的行政区划界线。

格式转换

GMT 不能识别 shapefile 格式。需要使用 gdal 的 ogr2ogr 命令将 shapefile 格式转换为 GMT 可识别的格式。如果系统里找不到 ogr2ogr 命令，则说明 gdal 没有安装。Ubuntu 下的包名是 gdal-bin ，CentOS 下的是 gdal。Windows 用户则需要到 gdal 官网下载安装。另外，Windows 用户也可以用 ArcGIS 等软件实现格式的转换。

具体的转换命令是从网上找到的，CHN_adm0 为要生成的数据的文件名前缀，但是为什么要出现两次，表示很不解。

ogr2ogr -f GMT -nln CHN_adm0 CHN_adm0 CHN_adm0.shp
ogr2ogr -f GMT -nln CHN_adm1 CHN_adm1 CHN_adm1.shp
ogr2ogr -f GMT -nln CHN_adm2 CHN_adm2 CHN_adm2.shp
ogr2ogr -f GMT -nln CHN_adm3 CHN_adm3 CHN_adm3.shp

对于 Hong Kong、Macao、Taiwan 的数据做类似操作，最终生成了一堆以 gmt 结尾的文件。共计 12 个，如下：

CHN_adm0.gmt  CHN_adm2.gmt  HKG_adm0.gmt  MAC_adm0.gmt  MAC_adm2.gmt  TWN_adm1.gmt
CHN_adm1.gmt  CHN_adm3.gmt  HKG_adm1.gmt  MAC_adm1.gmt  TWN_adm0.gmt  TWN_adm2.gmt

绘图测试

PS1：数据为多段数据，在 GMT4 中需要使用 -m 选项，而 GMT5 已经可以自动处理多段数据，所以不需要使用 -m 选项。

PS2：数据中有领土纠纷的部分不符合中国政府的立场，笔者也不赞同。虽然有问题的部分不能在正规期刊上发表，但是没有问题的部分还是可以使用的。欢迎读者告诉我符合中国政府政策的数据。

中国全图

绘制中国全图需要全部的 0 级数据。

#!/bin/bash
R=72/136/15/54
J=M15c
PS=china.ps

gmt psxy -J$J -R$R -T -K -U > $PS
gmt psxy -R$R -J$J CHN_adm0.gmt -K -O >> $PS
gmt psxy -R$R -J$J HKG_adm0.gmt -K -O >> $PS
gmt psxy -R$R -J$J MAC_adm0.gmt -K -O >> $PS
gmt psxy -R$R -J$J TWN_adm0.gmt -K -O >> $PS
gmt psxy -R$R -J$J -T -O >> $PS

代码运行过程中 GMT 会出现如下的警告 (或错误？)

psxy: Bad OGR/GMT: @D record has more items than declared by @N

猜测是 ogr2ogr 转换的问题。在我的系统环境该警告不影响绘图效果，但评论中 @vv 指出在他的系统环境下会导致图形无法绘制。

效果图（缺了南海的九段线数据）：

1 级行政区划

与上面的代码几乎一样，1 级数据中包含了 0 级数据，所以直接绘制 1 级数据即可。

R=72/136/15/54
J=M15c
PS=china.ps

gmt psxy -J$J -R$R -T -K -U > $PS
gmt psxy -R$R -J$J CHN_adm1.gmt -K -O >> $PS
gmt psxy -R$R -J$J HKG_adm1.gmt -K -O >> $PS
gmt psxy -R$R -J$J MAC_adm1.gmt -K -O >> $PS
gmt psxy -R$R -J$J TWN_adm1.gmt -K -O >> $PS
gmt psxy -R$R -J$J -T -O >> $PS

效果图：

在上一篇博文《GMT 中 DCW 数据的使用》中利用 GMT 自带的 DCW 数据也生成了类似的图，如下图。查看全图，对比一下会发现，两张图的细节方面还是有些区别的，本文的数据绘制的似乎包含了更多的细节（主要是小的岛屿）。这个就得根据需求去选择了，当然也有可能两个都是有问题的。

2 级行政区划

转换出来的 2 级数据中包含了全国所有的市级边界，用编辑器打开查看内容就会发现，每条线段都有完整的注释，很容易从众多线段中提取出自己想要的部分。以安徽省为例，将与安徽有关的线段数据保存到文件 Anhui_adm2.gmt 中：

R=114/120/29/35
J=M10c
PS=anhui.ps

gmt psxy -J$J -R$R -T -K -U > $PS
gmt psxy -R$R -J$J Anhui_adm2.gmt -K -O >> $PS
gmt psxy -R$R -J$J -T -O >> $PS

上面的脚本有一个很不方便的地方：想要画一个省的 2 级数据，每次都要从 CHN_adm2.gmt 中手动提取该省的数据信息。下面的例子可以避免这种手动提取的过程，主要通过 DCW 数据和 psclip 命令，使用全国 2 级数据（CHN_adm2.gmt），但是只绘制安徽省的 2 级数据。

R=114/120/29/35
J=M10c
PS=anhui.ps

gmt psxy -J$J -R$R -T -K -U > $PS
gmt pscoast -FCN.34 -M > Anhui_bnd.gmt
gmt psclip -J$J -R$R Anhui_bnd.gmt -K -O >> $PS
gmt psxy -R$R -J$J CHN_adm2.gmt -K -O >> $PS
gmt psclip -C -K -O >> $PS
gmt psxy -R$R -J$J -T -O >> $PS

脚本利用 pscoast 命令，将安徽省（代码为 34）的省界数据导出到文件 Anhui_bnd.gmt 中，然后利用该文件进行 clip，psxy 绘图时虽然使用的是全国的 2 级数据 CHN_adm2.gmt，但是只有安徽省内的部分会被绘制出来，最后还需要再次调用 psclip 以结束 clip。

上面的这个脚本生成了一个中间文件 Anhui_bnd.gmt，有强迫症的人是无法忍受这个的，因而上面的代码利用管道可以进一步简化为：

R=114/120/29/35
J=M10c
PS=anhui.ps

gmt psxy -J$J -R$R -T -K -U > $PS
gmt pscoast -FCN.34 -M | gmt psclip -J$J -R$R -K -O >> $PS
gmt psxy -R$R -J$J CHN_adm2.gmt -K -O >> $PS
gmt psclip -C -K -O >> $PS
gmt psxy -R$R -J$J -T -O >> $PS

上面三个脚本的最终结果基本是一致的，效果图如下：

注：三个脚本的成图效果是有差的，但是目前没有体现出来，在下一段 “3 级行政区划” 中，可以更明显地看出区别。

3 级行政区划

这里还是以安徽省为例，实际上只用 3 级数据即可，这里同时用了 2 级数据和 3 级数据，并且用不同的粗细和颜色来区分。需要注意，由于 3 级数据中包含了 2 级数据，所以下面的例子先画了 3 级数据，再用 2 级数据覆盖。如果画的顺序反了，效果就会差很多。

R=114.8/120/29.3/36
J=M14c
PS=anhui.ps

gmt psxy -J$J -R$R -T -K -U > $PS
gmt psxy -R$R -J$J CHN_adm3.gmt -W0.5p,gray -K -O >> $PS
gmt psxy -R$R -J$J CHN_adm2.gmt -W1p -K -O >> $PS
gmt psxy -R$R -J$J -T -O >> $PS

效果图如下：

下面的脚本利用了前面提到的 psclip 的方法：

R=114.8/120/29.3/36
J=M14c
PS=anhui.ps

gmt psxy -J$J -R$R -T -K -U > $PS
gmt pscoast -FCN.34 -M | gmt psclip -J$J -R$R -K -O >> $PS
gmt psxy -R$R -J$J CHN_adm3.gmt -W0.5p,gray -K -O >> $PS
gmt psxy -R$R -J$J CHN_adm2.gmt -W1p -K -O >> $PS
gmt psclip -C -K -O >> $PS
gmt psxy -R$R -J$J -T -O >> $PS