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

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

《arcgis10多边形叠加分析与重分类.docx》由会员分享，可在线阅读，更多相关《arcgis10多边形叠加分析与重分类.docx（11页珍藏版）》请在冰豆网上搜索。

arcgis10多边形叠加分析与重分类多边形叠加分析与重分类3.1、线和面的叠合启用ex16ex16.mxd，激活dataframe1，有3个图层，“区界”仅用于显示，不参与计算。

本练习要求得到“计算范围”内的道路网密度，该范围和区界多边形不重合。

右击dataframe1，选择“属性”，进入常规选项，确认“地图”和“显示”单位均为米，主菜单中选用“地理处理-环境”，进一步设置：

工作空间-当前工作空间：

ex16temp工作空间-临时工作空间：

ex16/temp确定返回。

ArcToolbox-AnalysisTools-Overlay-Intersect，出现以下对话框：

输出类型（可选）：

LINE，下拉选项，输出线要素类。

图11-1叠加分析对话框确定，经计算，产生要素类intersect1，自动加载，可看出原来的“道路”被切割后限定在“计算范围”内。

打开intersect1属性表，添加字段，进一步输入：

名称：

R_length类型：

浮点型精度：

8（字段宽度）小数位数：

2图11-2叠加产生需要计算的道路线段按确定，右击R_length，选择计算几何，进一步设置：

属性：

长度坐标系：

使用数据源的坐标系单位：

米按确定结束，R_length的取值为相交叠加处理后线要素的几何长度，还可以看到该属性表有一个字段Area，所有记录中取值均为1579952.9，这是“计算范围”多边形的面积，叠加分析过程中被带入了intersect1，右击R_length-Statistics统计，可以看到叠加分析处理后有55个路段，长度总和为11169.67，该范围的道路网密度如下：

11169.67/1579952.9*1000=7.07（km/km2）3.2、多边形叠合练习简介激活dataframe2，可以看到“高程”、“地块”二个多边形图层（见图11-3）。

“高层”多边形是由地形等高线组合而成。

打开图层属性表，可以看到该属性表有字段Hight，表示该多边形的最大高程。

打开图层属性表“Attributeof地块”，可以看到该属性表有Landuse，Value,Class等字段，分别表示土地使用、估计财产、地基类型。

在目录表的下部点击标签Source，可以看到，图层数据的存放路径。

点击按钮AddData，在d:

/GIS/ex16/路径下，选择独立属性表found.dbf，用Add按钮，该表加载。

鼠标右键点击该表，选Open打开，可以看到该表有Class，Para二个属性，表示地基类型和损失系数（图11-4），关闭属性表窗口。

图11-3dataframe2的显示图11-4地块属性表、地基损失参数表的逻辑关系本例为一个假设的洪水淹没损失估计，损失与如下因素有关：

（1）洪水水位的高程为505m，淹没范围按等高线考虑；

（2）住宅用地为财产损失的估计对象；（3）淹没引起的损失除了与居民的财产有关，也与场地的稳定性有关；（4）计算被淹没的面积，估计损失的财产。

1计算地块财产密度用鼠标右键打开图层“地块”的要素属性表，添加字段，进一步输入：

名称：

V_A类型：

浮点型精度：

7（字段宽度）小数位数：

6鼠标右键点击字段名V_A，选用字段计算器FieldCalculator，出现提示：

V_A=，用鼠标点击输入：

Value/P_area按确定键结束，可以看到字段V_A被赋值。

属性表关闭。

2空间叠合如果ArcToolbox窗口未显示，在地图窗口点击图标，调出ArcToolbox，展开ArcToolbox/AnalysisTools/Overlay，鼠标双击Union，出现Union对话框，在第一行InputFeatures提示下，展开选择“高程”，再展开选择“地块”，可以看到这二个图层名出现在Features列表中。

OutputFeatureClass：

d:

GISex16tempUnion1.shpJoinAttributes（Optional）：

ALL所有字段都合并按OK键继续，软件出现处理过程提示窗，右下角出现Completed提示，处理完毕，按Close键关闭。

叠合后的图层Union1出现在目录表、地图窗口内。

3计算叠合后的多边形面积打开图层属性表“AttributesofUnionl”，在表的由下侧选用菜单Options/AddField，在随后的对话框中定义所增加的字段：

Name：

New_AreaType：

doublePrecision：

10Scale：

2图11-5叠合后的多边形专题图层按OK键返回，再为字段New_Area赋值。

右击-计算几何，出现计算几何对话框：

属性是面积，单位是平方米，按“确定”，可以观察到字段New_Area中的计算结果为叠合后多边形的面积。

4计算地块估计损失、地块损失密度对union1，选用菜单Options/AddFied。

，在对话框内输入：

Name：

Estloss（地块的估计损失）Type:

double（数据类型为数值型）Precision：

10（宽度为10）Scale：

2（小数点保留2位）按OK键后，再选用菜单Options/AddField：

Name:

Lossden（单位面积的损失密度）Type:

double（数据类型为数值型）Precision：

8（宽度为8）Scale：

3（小数点保留3位）按OK键确认。

下一步将地基类型损失系数表连接到叠合多边形属性表。

缩小当前的属性表窗口，到目录表中点击Source标签（如果没有表名found，点击AddData，选择/ex16/found.dbf，加载）。

鼠标右键点击Union1.shp，选用快捷菜单的JoinsandRelates/Join，执行表和表的合并连接操作，在JoinData对话框中输入表连接的条件：

Whatdoyouwanttojointo：

Joinattributesfromatable下拉选择，执行表和表连接1.Choosethefieldinthislayerthatthejoinwilltakeplace：

CLASS下拉选择Union1属性表中的关键字段2.Choosethetabletojointothislayerorloadthetable：

found下拉选择被连表名3.Choosethefieldinthetabletobasethejoinonly：

CLASS下拉选择found表中的关键字段按OK键确认，提示是否增加属性索引，回答“No”，完成合并连接，可以看到属性表“AttributeofUnionl.shp”多了字段found.PARA，即损失系数。

鼠标右键点击表中的字段名Union1.Estloss，选用FieldCalculator，不使用VBA程序，在Union1.Estloss=的提示下，借助鼠标在文本框输入：

Union1.New_Area*Union1.V_A*found.PARA先按回车键表示输入完毕，再按OK键确认，字段Union1.Estloss被赋值，即：

地块估计损失=地块财产密度叠合后的多边形面积损失系数再用鼠标右键点击字段名Union1.Lossden，选用FieldCalculator，不使用VBA程序，在Union1.Lossden=的提示下，用鼠标输入：

Union1.V_A*found.PARA先按回车键表示输入完毕，再按OK键确认，字段Union1.lossden被赋值，即：

地块的损失密度=地块财产密度损失系数5地块过滤，分析结果表达关闭属性表窗口，返回dataframe2，过滤图层的要素，双击图层名Unionl，调出LayerProperties对话框，选择DefinitionQuery标签，点击按钮QueryBuilder，输入组合查询条件（单双引号必须都是英文字符）：

Union1.HIGHT=505ANDUnion1.LANDUSELIKER%按OK键返回，图层Unionl中只有高程小于或等于505，土地使用为住宅，其土地使用（LANDUSE）属性为R开头的字符串，才进入选择集，按“确定”键关闭LayerProperties对话框，可以看到Unionl中的要素被过滤。

过滤成功后，图层的要素明显减少。

再打开Unionl的LayerProperties窗口，选择Symbology标签，定义损失密度图层的图例，数量-分级色彩，中部字段框内，选择“值”；Lossden；按Classify按钮Classification/Method:

NaturalBreaks（Jenks）Classification/Classes:

3Classification/Method:

Manual下拉选择，手动方式分类在右侧BreakValues对话框中，输入：

.030.060.090按OK键返回，可以看到右侧的图例表：

Symbol（符号）Range（分类范围）Label（图例标识，中文输入）对应的颜色符号.005-.030低对应的颜色符号.031-.060中对应的颜色符号.061-.090高图11-6分级设置根据需要调整多边形的填充符号、颜色，按“应用”按钮，观察地图显示效果，满意后按“确定”键，关闭LayerProperties窗口。

可以关闭地块图层的显示状态，再打开，观察显示效果（参见图11-7）。

最后一步为汇总损失值，用鼠标右键打开图层Unionl的属性表“AttributeofUnionl”，鼠标右键点击该表的字段名LAND_ID，选汇总Summarize：

选择汇总字段：

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原来的值是地块的原始面积，经计算后，变成该地块的淹没比例。

关闭属性表窗口。

图11-7分析结果专题图显示练习结束，选File/Exit，关闭MapDocument，退出ArcMap。

6叠合功能简介空间数据的叠合（Overlay）是GIS的典型功能。

ArcMap可作点和面、线和面、面和面的叠合处理。

图11-8为二种典型的面和面相叠合的示意图（Union和Intersect）。

两个图层叠合的结果是生成一个新的图层，除图形数据来自被叠合的两个图层外，还保持了原来两个图层各自的属性数据。

INPUTPolygonUNIONPolygonOUTPUTPolygonINPUTPolygonINTERSECTPolygonOUTPUTPolygon图11-8二种典型叠合的示意3.3、多边形归并ArcMap可以按要素的属性进行合并，如将相同类型的点、线、面合并为点簇、线簇、多边形簇。

本练习对乡进行合并，使乡合并为县。

图11-9dataframe1的显示启动地图文档ex16ex16.mxd，激活dataframe3（见图11-9），仅有一个面状图层“乡”。

地图窗口在按钮条点击图标，出现ArcToolbox窗口，为了有效利用显示屏，可将ArcToolbox窗口拖动到目录表下侧，展开ArcToolbox/DataManagementTools/Generalize制图综合，鼠标双击Dissolve融合，出现Dissolve对话框：

（合并操作）InputFeatures：

乡OutputFeatureClass：

d:

GISex16tempdissolve1.shpDissolveField（s）（optional）：

合并的字段名CO_NAME勾选，按每个乡所在县的名称相同合并StatisticsField（s）（optional）：

需计算的字段FieldStatisticTypeAREASUM下拉添加多边形面积字段名，下拉选择计算方法为累加POPSUM下拉添加人口字段名，下拉选择计算方法为累加图11-10融合设置按OK键继续，软件出现处理过程提示窗，右下角出现Completed提示，表示处理完毕，按Close键关闭。

叠合后的图层dissolve1出现在目录表、地图窗口内。

每个乡按所在县合并，空间处理的结果是取消了多边形乡的边界，保留了县的边界，对AREA和POP二个字段作了累加计算。

用鼠标右键打开图层dissolve1的Properties，选择Symbology标签，为每个县定一种颜色（类别-唯一值UniqueValues为图例类型，CO_NAME为ValueField），点击Symbology标签左下角的Addallvalues添加分类的值，按“确定”键，可以看到志远县在几何上由三个互不相邻的多边形组成（地图右上角），但是每个县在空间上是一个整体，在属性表中也只有一条记录，打开要素属性表显示结果如下：

Dissolve_Shape*CO_NAMESUM_AREASUM_POP94（合并后的要素类型）（县名）（合并后的面积）（合并后的人口）Polygon兴益县437.722261461844Polygon广宁县665.428857468291Polygon志远县245.959860334346Polygon罗丰县510.328586618648图11-11按乡合并成县