如何生成一个 GeoBUGS 格式的中国地图

第一节 导言

之前有些对 GeoBUGS 感兴趣的同学发邮件询问我有没有 GeoBUGS 的中国地图，以用于分析中国国内的一些空间数据。我想有必要将如何生成 GeoBUGS 格式的地图的方法分享给大家。这样的话， GeoBUGS 就可以真正为我们所用，从而对于其他 GeoBUGS 没有自带的地图，我们也可以轻松生成了。本节不涉及统计分析，仅为 GeoBUGS 的研究使用者提供一个软件使用的技术参考。关于 GeoBUGS 的统计的书，国外实在是很多了，但关于这块地图定制的参考资料较少，故提出来供大家参阅。

目前分析用地图普遍采用的是 shp 格式，该格式可以用 ESRI 公司开发的 ArcInfo 工具进行编辑和分析。由于其通用性，故很容易在网上找到相应的资源。比如可以在国家基础地理信息系统的网站 http://nfgis.nsdi.gov.cn/nfgis/chinese/ 上下载到有用的 shp 文件。我们主要需要的是其中的国界和省界的 shp 文件。可以点击链接下载中国省级行政区域地图 shp 版。我们要用到的是其中的 bou2_4p.shp 文件。

第二节 用 ArcInfo 生成一个 GeoBUGS 格式的中国地图

第一步，打开 ArcInfo，选中 Layers（图层），点击右键，在打开的菜单中选择 Add Data…（添加数据…），并定位到 bou2_4p.shp，确定。

第二步，打开位于图层工具箱上方的 Editor 按钮菜单，在其中选中 Start Editing。

第三步，选中 bou2_4p 图层，右键选中 Open Attribute Table，

你会发现这个地图的精度过高，如浙江省的舟山群岛，多边形区域数目很多，而在一般的 GeoBUGS 分析中，这些多边形区域可以不予显示。所以我们可以通过删减，得到只含有主要部分的中国行政区域地图。一个简便的方法是对区域面积进行排序，小于一定阈值的区域予以删除。例如可以以 Area=.066 的香港特别行政区为界，删除所以面积小于 0.066 的区域。

第四步，打开位于图层工具箱上方的 Editor 按钮菜单，在其中选中 Stop Editing，并保存你的修改。

第五步，在 ArcToolBox 中选择 Conversion Tools  -> To Coverage -> Feature Class To Coverage。

将 Input Feature classes 定位到修改后的 bou2_4p.shp。在 XY Tolerance（optional）中设定精度，这里我将其设定为 0.005 Decimal degrees。其他参数取默认值。

第六步，在 ArcToolBox 中选择 Coverage Tools  -> Conversion -> From Coverage -> Ungenerate。将 Input Coverage 定位到 bou2_4p_feat。Feature Type 选 region. bou2_4p。Numeric Format 取 FIXED。其他选项取默认值。

第七步，编辑 GeoBUGS 支持的格式，在 bou2_4p_feat1.txt 中添加一段头，如下

map: 33
1 Heilongjiang
2 NeiMonggol
3 Xingjiang
4 Jilin
5 Liaoning
6 Gansu
7 Hebei
8 Beijing
9 Shanxi
10 Tianjin
11 Shaanxi
12 Ningxia
13 Qinghai
14 Shandong
15 Xizang
16 Henan
17 Jiangsu
18 Anhui
19 Sichuan
20 Hubei
21 Chongqing
22 Shanghai
23 Zhejiang
24 Hunan
25 Jiangxi
26 Yunnan
27 Guizhou
28 Fujian
29 Guangxi
30 Taiwan
31 Hainan
32 Guangdong
33 Hongkong

regions
1 Heilongjiang
2 NeiMonggol
3 Xingjiang
4 Jilin
5 Liaoning
6 Gansu
7 Hebei
8 Beijing
9 Shanxi
10 Tianjin
11 Hebei
12 Shaanxi
13 Ningxia
14 Qinghai
15 Shandong
16 Xizang   
17 Henan
18 Jiangsu
19 Anhui
20 Sichuan
21 Hubei
22 Chongqing
23 Shanghai
24 Shanghai
25 Zhejiang
26 Hunan
27 Jiangxi
28 Yunnan
29 Guizhou
30 Fujian
31 Guangxi
32 Taiwan
33 Hainan
34 Guangdong
35 Hongkong
END

map: 33 表示这个中国地图有 33 个行政区域，而 regions: 35 表示后面给出的通过经纬坐标描述的多边形区域隶属于那个行政区域。

第八步，在 GeoBUGS 中打开这个 txt 文件，然后选择 Map -> Import ArcInfo。

第九步，保存成 GeoBUGS 的. map 文件

第十步，重启 GeoBUGS，恭喜你可以在 GeoBUGS 中使用中国地图进行分析了。

生成的地图可以在这里下载中国行政区域地图 GeoBUGS 版。我们可以在许多政府官网得到许多可以分析的数据资源，如 http://www.moh.gov.cn/publicfiles//business/htmlfiles/zwgkzt/pwstj/index.htm

那么就下载一些数据用 GeoBUGS 进行分析吧

WinBUGS 在统计分析中的应用 第四部分完

敬告各位友媒，如需转载，请与统计之都小编联系（直接留言或发至邮箱：[email protected]），获准转载的请在显著位置注明作者和出处（转载自：统计之都），并在文章结尾处附上统计之都微信二维码。

发表 / 查看评论