Arcgis(一) 制作全国行政区shp文件(精确到县级)

        最近在忙毕业设计，需要用行政区的shp面文件来做空间相关分析。查了一些资料，发现并没有现成的某个省份或直辖市的shp文件，不少还需要RMB下载，唉，走了不少弯路，这里分享一下制作方法。

一、Shapefile说明

Shapefile是一种用于存储地理要素的几何位置和属性信息的非拓扑简单格式，它可以在Arcgis中进行使用和编辑。通常包括以下这些文件：

文件后缀名	主要作用
.shp	用于存储要素几何的主文件 ，必需文件
.shx	用于存储要素的几何索引文件，必需文件
.dbf	用于存储要素属性信息的的dBASE标，必需文件
.sbn和.sbx	用于存储要素空间索引的文件
.prj	用于存储坐标系信息的文件
.xml	用于存储shapefile的相关信息

        通常来说，我们简称的shp文件主要指前三种文件，最重要的就是shp和dbf文件，分别用于存储几何位置信息和属性信息。其他文件都是对其进行辅助完善的。（一般软件指定导出shapefile格式的时候，都会有这些文件，不用担心会少了哪个文件）

二、下载矢量数据文件

        首先，需要进入国家基础信息网站全国地理信息库，注册下载我们需要的数据。这里以安徽省为例，将I50和H50添加成果车，填写信息后可以免费下载。

三、Arcgis合并数据

        下载完成后，解压两个zip文件，得到两个gdb文件。打开Arcgis后，导入数据。数据类型如下图所示：

其中，BOUA为行政区面数据。

注意：部分沿海省份可能没有BOUA（例如山东省），此时可参照我这篇博客的方法进行获取Python一键获取市级行政单元shp文件

        导入方式：文件——添加数据——H50.gdb、I50.gdb——BOUA，导入结果如下：

此时，需要对数据进行合并处理：
中文版本：地理处理——ArcToolbox——数据管理工具——常规——合并
英文版本：ArcToolbox——Data Management Tools——General——Merge
添加数据后出现以下界面（注意一定要点右边的"文档"进行添加，否则两个BOUA名一样，不能合并！）：

点击确定后，需要等几分钟，合并完成后，如下图所示

两个区域的颜色保持一致了，但中间还有一条分割线需要进行处理，操作步骤如下：
中文版本：ArcToolbox——数据管理工具——制图综合——融合
英文版本：ArcToolbox——Data Management Tools——Generalization——Dissolve
勾选PAC和Name后，点击确定

四、筛选目标数据

合并数据完成后，点击右键打开属性表，观察可知，安徽省PAC的开头为34

以下，开始对数据进行筛选(这里就不写英文版本啦，因为翻译过来与软件差不多意思)：
分析工具——提取分析——筛选，在表达式输入SQL筛选语句：

"PAC">340000 AND "PAC"<350000

点击确定后就得到我们需要的目标省份的图层文件，最后一步导出shp文件即可：
右键数据——导出数据——保存类型(shapefile)
将生成的shp文件，导入Geoda中，显示如下：

如果要进行空间相关分析，还需在dbf文件中添加属性数据，可以不用Arcgis添加！只要用spss打开dbf即可！spss可保存为dbf文件

五、区县级shp转市级shp

可以先在excel中，打开生成的dbf文件，添加一列市编号town_id。即同一市级以下的区县记为同一编号，如下图所示：

由于excel没法保存为dbf文件，将该列数据复制到spss中，将其保存为dbf文件

导入Geoda，选择工具—dissolve（融合），将融合的变量设置为TOWN_ID，如下所示（低版本的Geoda可能没有这个功能，这里使用的是1.14版的）：

会生成新的shapefile文件，最终效果如下：

已将原来的104个区县单元转为16个市级单元。

以上就是本次分享的全部内容~