arcgis10多边形叠加分析与重分类.docx
《arcgis10多边形叠加分析与重分类.docx》由会员分享，可在线阅读，更多相关《arcgis10多边形叠加分析与重分类.docx（12页珍藏版）》请在冰点文库上搜索。
arcgis10多边形叠加分析与重分类
3.1、线和面的叠合
启用\ex16\ex16.mxd，激活dataframe1，有3个图层，“区界”仅用于显示，不参与计算。
本练习要求得到“计算范围”内的道路网密度，该范围和区界多边形不重合。
右击dataframe1，选择“属性”，进入常规选项，确认“地图”和“显示”单位均为米，主菜单中选用“地理处理-环境…”，进一步设置：
工作空间-当前工作空间：
\ex16\temp
工作空间-临时工作空间：
\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:
\GIS\ex16\temp\Union1.shp
JoinAttributes（Optional）：
ALL所有字段都合并
按OK键继续，软件出现处理过程提示窗，右下角出现Completed提示，处理完毕，按Close键关闭。
叠合后的图层Union1出现在目录表、地图窗口内。
3计算叠合后的多边形面积
打开图层属性表“AttributesofUnionl”，在表的由下侧选用菜单Options/AddField…，在随后的对话框中定义所增加的字段：
Name：
New_Area
Type：
double
Precision：
10
Scale：
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，调出Layer
Properties对话框，选择DefinitionQuery标签，点击按钮QueryBuilder…，输入组合查询条件（单双引号必须都是英文字符）：
"Union1.HIGHT"<=505AND"Union1.LANDUSE"LIKE'R%'
按OK键返回，图层Unionl中只有高程小于或等于505，土地使用为住宅，其土地使用（LANDUSE）属性为R开头的字符串，才进入选择集，按“确定”键关闭LayerProperties对话框，可以看到Unionl中的要素被过滤。
过滤成功后，图层的要素明显减少。
再打开Unionl的LayerProperties窗口，选择Symbology标签，定义损失密度图层的图例，数量-分级色彩，中部字段框内，选择“值”；Lossden；
按Classify…按钮
Classification/Method:
NaturalBreaks（Jenks）
Classification/Classes:
3
Classification/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，勾选“总和”
制定输出表：
\ex16\temp\Sum_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）。
两个图层叠合的结果是生成一个新的图层，除图形数据来自被叠合的两个图层外，还保持了原来两个图层各自的属性数据。
INPUTPolygonUNIONPolygonOUTPUTPolygon
INPUTPolygonINTERSECTPolygonOUTPUTPolygon
图11-8二种典型叠合的示意
3.3、多边形归并
ArcMap可以按要素的属性进行合并，如将相同类型的点、线、面合并为点簇、线簇、多边形簇。
本练习对乡进行合并，使乡合并为县。
图11-9dataframe1的显示
启动地图文档\ex16\ex16.mxd，激活dataframe3（见图11-9），仅有一个面状图层“乡”。
地图窗口在按钮条点击图标
，出现ArcToolbox窗口，为了有效利用显示屏，可将ArcToolbox窗口拖动到目录表下侧，展开ArcToolbox/DataManagementTools/Generalize制图综合，鼠标双击Dissolve融合，出现Dissolve对话框：
（合并操作）
InputFeatures：
乡
OutputFeatureClass：
d:
\GIS\ex16\temp\dissolve1.shp
DissolveField（s）（optional）：
合并的字段名
√CO_NAME勾选，按每个乡所在县的名称相同合并
StatisticsField（s）（optional）：
需计算的字段
FieldStatisticType
AREA
SUM
下拉添加多边形面积字段名，下拉选择计算方法为累加
POP
SUM
下拉添加人口字段名，下拉选择计算方法为累加
图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.722261461844
Polygon广宁县665.428857468291
Polygon志远县245.959860334346
Polygon罗丰县510.328586618648
图11-11按乡合并成县
升级会员 