中国行政区划数据下载
地球物理相关资源
文章目录
- GADM 数据库
- 绘图测试
- 修订历史
GADM 数据库
GADM,是一个全球行政区划数据库。包括了几乎全部国家和地区的行政区划界线数据。
主页: http://www.gadm.org/
下载数据
下载地址: http://www.gadm.org/country
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 会出现如下的警告 (或错误?)
猜测是 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
效果图如下:
将这两张图对比一下,容易发现,省界的部分线段明显变细了,这算是 clip 的一个缺点,使用省界数据进行 clip,同时又要绘制省界数据,如何判断省界数据点是否在 clip 区域内部是个问题。
修订历史
- 2013-11-23:初稿;
- 2013-11-28:绘制 2 级和 3 级边界时,利用 clip 的方法以减少人工操作。Thanks to @yangtze。
- 2013-12-05:删除了数据包中的冗余隐藏文件,重新打包,提供多种格式下载。
- 2014-01-19:不推荐下载全球行政区划数据;
- 2015-04-11:存在个别 IP 恶意下载数据,导致我流量异常高,因而删除其他格式的数据压缩包,仅留下 7z 格式的,若仍出现流量异常,则删除该数据;
- 2015-05-01:流量依然异常,猜测是搜索引擎会对齐定时抓取导致的;虽然已经在七牛限制了文件的非本站访问,但似乎对压缩文件无效。故删除该数据,有需要的可以自己去转换格式;
- 2016-04-13:王亮添加了 GADM 数据库新增的 Macou 1、2 级数据,并修正了部分理解和表述;