2021版1:100万全国标准基础地理数据！你值得拥有

本文介绍国家标准地理数据的获取以及处理为 GMT 支持格式的方法。

具体包括：

• 国家标准1:100万（非涉密）基础地理数据的申请、下载
• 数据集内容简介
• 数据集要素合并、提取及格式转换

数据概况

国家基础地理信息中心于2017年11月份开始免费向公众提供全国1:100万非涉密的全国基础地理数据库。目前数据最新为2021版。该数据库：

• 共77幅DLG数据集，包括全国陆地范围以及台湾岛、海南岛、钓鱼岛、南海诸岛在内的主要岛屿及其临近海域。
• 包括水系、居民地及设施、公路铁路、行政边界、地名等5个数据集12个数据层。

数据更多介绍可访问全国地理信息资源目录服务系统。

数据获取

• 进入全国地理信息资源目录服务系统，点右上角“注册”。已注册用户直接登陆。
• 在“数据下载”栏目中，点击进入“1：100万全国基础地理信息数据库”。
• 该数据集下有77个图幅产品。可选择一两个图幅，或全部选中，加入成果车。
• 结算，填写订单中的个人信息，验证手机号，下载数据。 **注：**若未提示验证手机，换个firefox浏览器试试。

数据集内容简介

该数据集包括5类12个要素集，如下图所示：

每个要素集内又包括不同等级的自然地理要素，以BOUL要素集（图层）为例：

BOUL 中国及邻区各级行政界限（线）图层

该数据集包括各级行政区界线，以GB码进行区分：

其他各图层内不同要素的分类编码可参考基础地理信息要素分类与代码（GB/T 13923-2006）。

此外，常用的还有 AGNP 各级行政地名和城乡居民地坐标（点）图层

行政地名点以政府驻地来定位，无GB码，以CLASS类别区分：

数据处理

从网站上直接下载的数据是分图幅的 geodatabase(gdb) 数据库文件，如果需要用 GMT画“全国一张图”，需要：

• 合并77个gdb文件中的同类要素
• 提取某类要素中的某个等级数据
• 转为gmt识别的文本格式

gdb数据库文件合并

采用本社区一贯推荐的GDAL/ORG库将分幅的gdb数据进行批量合并。 首先解压所有下载的压缩包，得到以.gdb结尾的目录。目录名类似H51.gdb、I49.gdb等。 然后在这些目录的上级目录执行以下脚本：

Linux系统：

#!/bin/bash
file="./final"
for i in $(ls -d *.gdb)
do
 if [ ! -d "$file" ]
      then
           echo "creating final filefolder" 
           ogr2ogr -f 'ESRI Shapefile' $file $i -lco encoding=UTF-8
      else
           echo "merging……"
      ogr2ogr -f 'ESRI Shapefile' $file $i -update -append
fi
done

Windows系统。注意必须保存为 bat 文件执行，直接复制粘贴到 cmd 中会出现报错：

REM 将A49分幅的gdb数据转为一组shp格式的文件，文件夹命名为Merged
REM 编码为UTF-8以正确显示属性字段中的中文
ogr2ogr -f "ESRI Shapefile" -lco encoding=UTF-8 Merged A49.gdb

REM 将当前目录下剩余所有的gdb数据追加Merged，同类shp文件合并。
for /f %%i in ('dir /b "*.gdb"') do (
ogr2ogr -f "ESRI Shapefile" -update -append Merged %%i
)

pause

若安装有python环境，还可选择ogrmerge.py工具进行批量的格式转换和合并，详见ogrmerge使用手册。

该节将分幅的交通、水系、边界等12类要素分别合并为全国一张图，最后转为12个shp文件。

查询GB代码与CLASS类别

使用ogrinfo可以查询shp文件内包含哪些要素。例如以下命令查询BOUL.shp中包含的GB代码：

ogrinfo -sql "SELECT distinct GB FROM BOUL" BOUL.shp

而以下命令则查询AGNP.shp中包含的CLASS类别：

ogrinfo -sql "SELECT distinct CLASS FROM AGNP" AGNP.shp

提取某类要素

由于每个数据集（一个shp文件）中存在不同等级的要素，而实际可能只使用其中一种，则仍可利用ogr2ogr工具提供的sql查询选项进行提取。

例如：

行政区划BOUL图层中GB代码620201、620202、250200分别为已定国界、未定国界和海岸线，那么提取数据并转化为GMT格式的命令为：

ogr2ogr -f GMT -where "GB=620201 or GB=620202 or GB=250200" BOUL_sub.gmt BOUL.shp

最后GMT绘图如下：

gmt begin ChinaMap
    gmt plot BOUL_sub.gmt -R70/145/3/55 -JM105/35/6.5i 
gmt end show

绘图结果：

如果想要提取行政地名点AGNP图层中的建制村信息（对应CLASS为AK），则使用：

ogr2ogr -f GMT -where "CLASS='AK'" AGNP.gmt AGNP.shp

注意因为CLASS类别为字符串，所以此处AK需要加上引号。GB代码为整型数因而不需要加引号。

最后转换成的GMT格式数据如下：

...
# @NCLASS|NAME|PINYIN|GNID|XZNAME
# @Tstring|string|string|string|string
# FEATURE_DATA
# @DAK|岙里王村|Aoliwangcun|330226114037|胡陈乡
121.68517089 29.388652295
# @DAK|坝头村|Batoucun|331022000000|亭旁镇
121.310746425 29.0109724900001
# @DAK|白岩村|Baiyancun|331082111028|尤溪镇
120.971417165 28.7226002200001
# @DAK|柏树下村|Baishuxiacun|331082116042|邵家渡街道
121.257111875 28.8522451100001
...

以 Windows 系统为例，首先将以上数据的编码转换为GB2312，再使用以下代码绘制：

gmt begin AGNP jpg
  gmt set PS_CHAR_ENCODING Standard+

  gmt basemap -R120/121/28/29 -JM10c -Baf
  gmt grdimage @earth_relief_15s -Cgeo
  gmt plot -Sc0.2c -Gred -W0p,black AGNP.gmt
  chcp 936
  gmt convert AGNP.gmt -a2="NAME" | gmt text -F+f10p,40+jTL -D0.15c/0.3c -Gwhite
gmt end

绘图结果：