arcgis10多边形叠加分析与重分类.docx

1、arcgis10 多边形叠加分析与重分类多边形叠加分析与重分类 3.1、线和面的叠合 启用ex16ex16.mxd，激活 data frame1，有 3 个图层，“区界”仅用于显示，不参与计算。本练习要求得到“计算范围”内的道路网密度，该范围和区界多边形不重合。右击 data frame1，选择“属性”，进入常规选项，确认“地图”和“显示”单位均为米，主菜单中选用“地理处理-环境”，进一步设置：工作空间-当前工作空间：ex16temp 工作空间-临时工作空间：ex16/temp 确定返回。ArcToolbox-Analysis Tools-Overlay-Intersect，出现以下对话框：输

2、出类型（可选）：LINE，下拉选项，输出线要素类。图 11-1 叠加分析对话框 确定，经计算，产生要素类 intersect1，自动加载，可看出原来的“道路”被切割后限定在“计算范围”内。打开 intersect1 属性表，添加字段，进一步输入：名称：R_length 类型：浮点型 精度：8（字段宽度）小数位数：2 图 11-2 叠加产生需要计算的道路线段 按确定，右击 R_length，选择计算几何，进一步设置：属性：长度 坐标系：使用数据源的坐标系 单位：米 按确定结束，R_length 的取值为相交叠加处理后线要素的几何长度，还可以看到该属性表有一个字段 Area，所有记录中取值均为 1

3、579952.9，这是“计算范围”多边形的面积，叠加分析过程中被带入了 intersect1，右击 R_length-Statistics 统计，可以看到叠加分析处理后有 55个路段，长度总和为 11169.67，该范围的道路网密度如下：11169.67/1579952.9*1000=7.07(km/km2)3.2、多边形叠合练习简介 激活 data frame2，可以看到“高程”、“地块”二个多边形图层（见图 11-3）。“高层”多边形是由地形等高线组合而成。打开图层属性表，可以看到该属性表有字段 Hight，表示该多边形的最大高程。打开图层属性表“Attribute of 地块”，可以看到

4、该属性表有Landuse，Value,Class 等字段，分别表示土地使用、估计财产、地基类型。在目录表的下部点击标签 Source，可以看到，图层数据的存放路径。点击按钮 Add Data，在 d:/GIS/ex16/路径下，选择独立属性表 found.dbf，用 Add 按钮，该表加载。鼠标右键点击该表，选 Open 打 开，可以看到该表有 Class，Para 二个属性，表示地基类型和损失系数（图 11-4），关闭属性表窗口。图 11-3 data frame2 的显示 图 11-4 地块属性表、地基损失参数表的逻辑关系 本例为一个假设的洪水淹没损失估计，损失与如下因素有关：（1）洪水水

5、位的高程为 505m，淹没范围按等高线考虑；（2）住宅用地为财产损失的估计对象；（3）淹没引起的损失除了与居民的财产有关，也与场地的稳定性有关；（4）计算被淹没的面积，估计损失的财产。1 计算地块财产密度 用鼠标右键打开图层“地块”的要素属性表，添加字段，进一步输入：名称：V_A 类型：浮点型 精度：7（字段宽度）小数位数：6 鼠标右键点击字段名 V_A，选用 字段计算器 Field Calculator，出现提示：V_A=，用鼠标点击输入：Value/P_area 按确定键结束，可以看到字段 V_A 被赋值。属性表关闭。2 空间叠合 如果 ArcToolbox 窗口未显示，在地图窗口点击 图

6、标，调出 ArcToolbox，展开 ArcToolbox/Analysis Tools/Overlay，鼠标双击Union，出现 Union 对话框，在第一行 Input Features 提示下，展开选择“高程”，再展开选择“地块”，可以看到这二个图层名出现在 Features 列表中。Output Feature Class：d:GISex16tempUnion1.shp JoinAttributes(Optional)：ALL 所有字段都合并 按 OK 键继续，软件出现处理过程提示窗，右下角出现 Completed 提示，处理完毕，按 Close 键关闭。叠合后的图层 Union1 出

7、现在目录表、地图窗口内。3 计算叠合后的多边形面积 打开图层属性表“Attributes of Unionl”，在表的由下侧选用菜单 Options/Add Field，在随后的对话框中定义所增加的字段：Name：New_Area Type：double Precision：10 Scale：2 图 11-5 叠合后的多边形专题图层 按 OK 键返回，再为字段 New_Area 赋值。右击-计算几何，出现计算几何对话框：属性是面积，单位是平方米，按“确定”，可以观察到字段 New_Area 中的计算结果为叠合后多边形的面积。4 计算地块估计损失、地块损失密度 对 union1，选用菜单 Opt

8、ions/Add Fied。，在对话框内输入：Name：Estloss（地块的估计损失）Type:double（数据类型为数值型）Precision：10（宽度为 10）Scale：2 （小数点保留 2 位）按 OK 键后，再选用菜单 Options/Add Field：Name:Lossden（单位面积的损失密度）Type:double（数据类型为数值型）Precision：8（宽度为 8）Scale：3（小数点保留 3 位）按 OK 键确认。下一步将地基类型损失系数表连接到叠合多边形属性表。缩小当前的属性表窗口，到目录表中点击 Source 标签（如果没有表名 found，点击 Add D

9、ata，选择/ex16/found.dbf，加载）。鼠标右键点击 Union1.shp，选用快捷菜单的 Joins and Relates/Join，执行表和表的合并连接操作，在 Join Data 对话框中输入表连接的条件：What do you want to join to：Join attributes from a table 下拉选择，执行表和表连接 1.Choose the field in this layer that the join will take place：CLASS 下拉选择 Union1 属性表中的关键字段 2.Choose the table to join

10、 to this layer or load the table：found 下拉选择被连表名 3.Choose the field in the table to base the join only：CLASS 下拉选择 found 表中的关键字段 按 OK 键确认，提示是否增加属性索引，回答“No”，完成合并连接，可以看到属性表“Attribute of Unionl.shp”多了字段 found.PARA，即损失系数。鼠标右键点击表中的字段名 Union1.Estloss，选用 Field Calculator，不使用 VBA 程序，在 Union1.Estloss=的提示下，借助鼠标

11、在文本框输入：Union1.New_Area*Union1.V_A*found.PARA 先按回车键表示输入完毕，再按 OK 键确认，字段 Union1.Estloss 被赋值，即：地块估计损失=地块财产密度 叠合后的多边形面积 损失系数 再用鼠标右键点击字段名 Union1.Lossden，选用 Field Calculator，不使用 VBA 程序，在 Union1.Lossden=的提示下，用鼠标输入：Union1.V_A*found.PARA 先按回车键表示输入完毕，再按 OK 键确认，字段 Union1.lossden 被赋值，即：地块的损失密度=地块财产密度 损失系数 5 地块过滤

12、，分析结果表达 关闭属性表窗口，返回 data frame2，过滤图层的要素，双击图层名 Unionl，调出 Layer Properties 对话框，选择 Definition Query 标签，点击按钮 Query Builder，输入组合查询条件（单双引号必须都是英文字符）：Union1.HIGHT=505 AND Union1.LANDUSE LIKE R%按 OK 键返回，图层 Unionl 中只有高程小于或等于 505，土地使用为住宅，其土地使用（LANDUSE）属性为 R 开头的字符串，才进入选择集，按“确定”键关闭 Layer Properties 对话框，可以看到 Union

13、l 中的要素被过滤。过滤成功后，图层的要素明显减少。再打开 Unionl 的 Layer Properties 窗口，选择 Symbology 标签，定义损失密度图层的图例，数量-分级色彩，中部字段框内，选择“值”；Lossden；按 Classify 按钮 Classification/Method:Natural Breaks（Jenks）Classification/Classes:3 Classification/Method:Manual 下拉选择，手动方式分类 在右侧 Break Values 对话框中，输入：.030.060.090 按 OK 键返回，可以看到右侧的图例表：Sym

14、bol（符号）Range（分类范围）Label（图例标识，中文输入）对应的颜色符号.005-.030 低 对应的颜色符号.031-.060 中 对应的颜色符号.061-.090 高 图 11-6 分级设置 根据需要调整多边形的填充符号、颜色，按“应用”按钮，观察地图显示效果，满意后按“确定”键，关闭 Layer Properties 窗口。可以关闭地块图层的显示状态，再打开，观察显示效果（参见图 11-7）。最后一步为汇总损失值，用鼠标右键打开图层 Unionl 的属性表“Attribute of Unionl”，鼠标右键点击该表的字段名 LAND_ID，选 汇总 Summarize：选择汇总

15、字段：union1.LAND_ID 选择一个或多个要包括在输出表中的汇总统计信息：展开 union1.P_area，勾选“平均”展开 union1.NEW_Area，勾选“总和”展开 union1.Estloss，勾选“总和”制定输出表：ex16tempSum_Output.dbf 保存类型为 dBASE 表 确定，添加到地图文档，打开 Sum_Output，在右键点击 Average_P_area，选用“字段计算器”，在 Average_P_area=的提示下，利用鼠标双击字段名，在文本框实现输入：Sum_NEW_Ar/Ave_P_area*100 确定。Average_P_area 原来的

16、值是地块的原始面积，经计算后，变成该地块的淹没比例。关闭属性表窗口。图 11-7 分析结果专题图显示 练习结束，选 File/Exit，关闭 Map Document，退出 ArcMap。6 叠合功能简介 空间数据的叠合（Overlay）是 GIS 的典型功能。ArcMap 可作点和面、线和面、面和面的 叠合处理。图 11-8 为二种典型的面和面相叠合的示意图（Union 和 Intersect）。两个图层叠合的 结果是生成一个新的图层，除图形数据来自被叠合的两个图层外，还保持了原来两个图层各自的属性数据。INPUT Polygon UNION Polygon OUTPUT Polygon I

17、NPUT Polygon INTERSECT Polygon OUTPUT Polygon 图 11-8 二种典型叠合的示意 3.3、多边形归并 ArcMap可以按要素的属性进行合并，如将相同类型的点、线、面合并为点簇、线簇、多边形簇。本练习对乡进行合并，使乡合并为县。图 11-9 data frame1 的显示 启动地图文档ex16ex16.mxd，激活 data frame3（见图 11-9），仅有一个面状图层“乡”。地图窗口在按钮条点击图标，出现 ArcToolbox 窗口，为了有效利用显示屏，可将 ArcToolbox 窗口拖动到目录表下侧，展开 ArcToolbox/Data Man

18、agement Tools/Generalize制图综合，鼠标双击 Dissolve融合，出现 Dissolve 对话框：（合并操作）Input Features：乡 Output Feature Class：d:GISex16tempdissolve1.shp Dissolve Field(s)(optional)：合并的字段名 CO_NAME 勾选，按每个乡所在县的名称相同合并 Statistics Field(s)(optional)：需计算的字段 Field Statistic Type AREA SUM 下拉添加多边形面积字段名，下拉选择计算方法为累加 POP SUM 下拉添加人口字

19、段名，下拉选择计算方法为累加 图 11-10 融合设置 按 OK 键继续，软件出现处理过程提示窗，右下角出现 Completed 提示，表示处理完毕，按 Close 键关闭。叠合后的图层 dissolve1 出现在目录表、地图窗口内。每个乡按所在县合并，空间处理的结果是取消了多边形乡的边界，保留了县的边界，对 AREA 和 POP 二个字段作了累加计算。用鼠标右键打开图层 dissolve1 的 Properties，选择 Symbology 标签，为每个县定一种颜色（类别-唯一值 Unique Values为图例类型，CO_NAME 为 Value Field），点击 Symbology 标

20、签左下角的 Add all values 添加分类的值，按“确定”键，可以看到志远县在几何上由三个互不相邻的多边形组成（地图右上角），但是每个县在空间上是一个整体，在属性表中也只有一条记录，打开要素属性表显示结果如下：Dissolve_Shape*CO_NAME SUM_AREA SUM_POP94（合并后的要素类型）（县名）（合并后的面积）（合并后的人口）Polygon 兴益县 437.722261 461844 Polygon 广宁县 665.428857 468291 Polygon 志远县 245.959860 334346 Polygon 罗丰县 510.328586 618648 图 11-11 按乡合并成县