利用ARCGIS进行地类计算与统计

本文格式为Word版，下载可任意编辑——利用ARCGIS进行地类计算与统计地类计算与统计

一、数据准备。

应用到的数据包括社界（DWG文件）和所求年份的现状图（Shape或GeoDatabase）DWG文件的注记的插入点（Text为左下角点）要落在对应的社界面上，这样才能保证数据转换后注记和面层的一一对应。

二、数据处理。

数据的处理主要包括数据的转换、点面属性连接、数据相交三个部分。

2.1数据转换

2.1.1建立数据集

（1）点击Geodatabase：

开启ArcCatalog程序，找适合的路径然后→新建→Personal

-1-

▽

（2）新建数据库后，进入数据库，→新建数据集→输入名称外，其余使用默认设置，不用修改，直接，直到：

-2-

▽

输入数据集名称

▽

坐标系统选Unknown(未知),或者用选DWG文件的投影

▽

-3-

容限使用默认即可，点击

▽

2.1.2将DWG数据导入数据集

！！首先必需确定DWG文件的路径没有中文名（D:\\pssj\\sj.dwg），否则导入会一片空白（1）进入数据集aaa，→→

（2）在中添加DWG文件的注记层（Annotation）和面层（Polygon）

▽

-4-

添加注记和面层，然后点导入

▽

2.1.3检查修改面层的拓扑

！！面层可能存在裂缝和重叠错误，这样会带来计算面积的错误，因此要进行拓扑的检查（1）数据集aaa中，→→

（2）前面两步使用默认直接，选择要素的时候勾选面层然后

▽

-5-

使用默认等级，然后

▽

添加拓扑规则，规则选择和，分两次添加，然后

▽

点击，然后选验证拓扑

▽

-6-

（3）修改拓扑错误

开启ArcMap，添加aaa_Topology，即可看到拓扑检查结果（红色部分）

可以看出，面层存在重叠，不存在裂缝，修要修改重叠部分，采用挖空的方法

▽→→用（

）选中重叠部分：

▽

-7-

→（clip，可以将与选择部分有重叠的所有面擦除）

▽

直接，对所有重叠部分重复以上步骤

▽

在ArcMap工具栏位置→调出拓扑工具栏→

重新验证拓扑看还有没有拓扑错误

▽

修改完没有拓扑问题后，→→

▽

-8-

拓扑错误已经消除，可以进行下一步操作

！！假使导入DWG文件的线层，然后用线层构面的话，可以省略掉拓扑检查和修改这一步，操作会相对简单些。但拓扑检查是保证数据质量最常用的检查方式，因此这里将这部分内容写上，大家可以对比一下两种方法，以达到一方面简化操作步骤增加效率，一方面可以增加对拓扑的理解，保证数据质量的效果。！！重叠部分的挖空和缝隙的填补操作是比较常用的编辑方式。挖空用命令，填补缝隙则通过修改节点或者完成，在此不写具体的操作步骤。

2.2点面属性连接

!!点面属性连接是根据点落在面内部的特性，通过空间连接运算，将点的属性添加到面的属性里面。具体来说就是标注点（Annotation）的Refname字段存储了社名，要将社名这个属性添加到对应的社界面上。

2.2.1添加数据，检查位置，保存字段

（1）在ArcMap中添加标注层（sj_dwg_Annotation）和面层（sj_dwg_Polygon）

-9-

原CAD图对应ARCMAP图

检查一下点的位置，确定每个面上对应唯一的注记点，排除异常状况。（2）保存字段

保存字段是保存有用的字段参与空间连接，去除不需要的字段，使结果简单明白。

开启相应图层的

▽

-10-

切换到这一栏

▽

对于注记层，点，然后勾选，保存Refname字段（存储社名的字段）对于面层，点即可，不需要保存属性字段

2.2.2进行空间相交运算

点击工具栏的按钮，开启ArcToolbox工具箱

▽

▽

-11-

依照→→→（空间连接）

开启空间连接对话框

！！早期计算地类采用（相交）工具，一样可以得到结果且有助于理解连接的过程，但操作相对繁杂，应用工具操作比较简单，对结果检查比较简单

▽

对话框中，下拉选择面层，下拉选择注记层，则自己选择路径并为输出的要素类命名（在此命名为SJ）。其余使用默认。点击

▽

ArcMap自动加载空间连接生成的图层（SJ），开启其

-12-

▽

找到相应图层→→

▽

存储社名，存储面积，字段是空间连接过程中自动生成的，表示连接到这个面的注记点的个数，假使为0则表示这个面上没有点，假使为2或3则表示这个面上有2个或者3个点，这样就不符合一一对应原则，要检查修改。由上图可知全部为1，点与面一一对应，可以进行下一步运算。！！到了这一步，社界面积就计算出来了，可以通过Access统计每个社的面积。

2.3图层相交求地类面积

（1）在上节开启的ArcMap中添加2023年的现状地类Shape数据（Geodatabase数据也可，假使求2023年地类则添加2023年现状图，以此类推）

（2）由图上可以看出社界和现状图重叠在一起，相交运算是通过现状图的界限将村界切开，并且把现状图的属性放到切开的小块社界中（可以参看Intersect工具的帮助）

-13-

（3）开启相交工具对话框→→→（相交）

在中下拉选择SJ（社界）和07MZDL（07现状图），选择路径并输入名称（在此命名为SJDL）。其余使用默认，点击开始相交运算。运算完成后会自动将SDDL加载到ArcMap中。

▽

上图中红色是原先SJ层边界，淡绿填充层是相交后生成的SJDL层，可以看出SJDL实际上

-14-

就是07地类界限将SJ层切开，或者SJ层界限将07地类切开，亦即两者的公共部分保存。（4）开启SJDL的图层属性，在字段对话框中保存（社名）、（地类代码）、(地类名称)和（面积）四个字段（参见2.2.1）。（5）开启SJDL的属性表（参见2.2.2）

属性表里面每一条记录对应SJDL中的一个面，说明白这块地所属的社、地类名称和代码以及面积这些信息，譬如说腊圃村村委这块地被07现状边界切成了四块，对应前面四条记录。这个表是最基本的表，后面的数据都将由这个表整理得出。（6）导出属性表

在属性表界面中，右下角的菜单，→，开启导出表对话框

▽

-15-

导出数据对话框中，栏选allrecords（所有数据），输出表里选路径并命名dbf文件（这里命名为DLSJ.dbf）,然后

▽

由上图可以看出，导出的dbf文件包含两个：dbf后缀和xml后缀文件。dbf后缀文件存储所有数据；xml是元数据文档，起说明作用，可以忽略。

▽

用Exceldbf后缀文件，xls文件（DLSJ.xls）。至此，地类面积已经求出并导出为Excel文件，后续统计以导出的Excel文件为基础。

三、地类数据统计

（1）开启Access（→→下图）

(2)用Access开启DLSJ.xls

-16-

(3)上图中，点击后，会出现导入对话框

这个不修改，直接

▽

这里要勾选，然后。这里直接点击则使用默认名为表名，假使想修改表名则点击，然后再。

▽

-17-

如上图所示，在数据库的表视图里可以看到导入的表（DLSJ）

▽

双击表DLSJ可以查看导入内容是否正确。

（4）地类统计，合并一致地类的面积。

由上图可以看到，每个社都存在有同一地类有多块地的状况，如村委有两块1110的地，布尾社有多块1310的地。一般状况我们制表只需知道红线范围内每个社每个地类有多少面积，所以要进行统计将同一个社的同一地类的所有地块的面积加起来。

-18-

点击→双击

▽

添加要进行查询的表，选中DLSJ然后→

▽

在查询对话框中，分别RefName、DLDM、DLMC和Shape_Area，添加四个字段到查询字段中。

▽

-19-

在字段栏内部→，调出选项

▽

▽

在选项中，除了Shape_Area对应栏选外，其余栏选。其他选项不用修改！！这个查询实际上是将同一个分组的一致项进行总计，首先是合并同一个社的一致地类，然后合并同一地类的不同地块面积。

▽

选项设立好后，点击工具栏的（运行），则得到查询结果。

▽

-20-

ArcMap自动把转化的注记层（SJAnno）加载数据库中生成了SJAnno要素类

到这一步，注记要素在ArcGIS中已经生成，剩下的就是将其和面层导出为CAD格式

▽

开启FMEWorkbench（→→）！！FME是加拿大Safe软件公司开发的一款商用数据转换软件，在数据格式转换方面较为常用，功能十分强大。

▽

根据向导进行数据转换

直接

▽

-31-

这一步是选择源数据的格式▽

选中ESRIGeodatabase（MDB），然后

▽

-32-

格式选好后，进行下一步

▽

这一步是对MDB进行定位，点击“?〞浏览数据库并将其添加，结果如图，点击

到下一步▽

-33-

这一步是选择目标格式▽

选择ACAD格式，然后，然后

▽

-34-

这一步是对目标格式进行设置，譬如版本、属性等。点击

▽

在Release中选择2000，其他不变，然后，然后

！！这个选项规定了输出的CAD版本，默认是2023，由于CAD2023打不开2023版本的dwg文件，因此选择比较通用的2000版本。

▽

-35-

直接点击（完成）

▽

此对话框是选择参与转换的要素类，只勾选面层（SJ）和注记层（SJAnno）即可，然后

▽

-36