[分享]多层面叠合分析

楼主^#

更多发布于：2004-08-05 17:28

一、 目的
多层面叠合分析是空间数据分析的基本方法，包括栅格数据的信息复合分析与矢量数据的叠置分析，本实验的目的在于：
1、加深对多层面叠合分析基本原理、方法的认识；
2、熟练掌握ARCVIEW多层面叠合分析的技术方法。
3、结合实际、掌握利用多层面叠合分析方法解决地学空间分析问题的能力。
二、 实验准备
1、软件准备：Arcview
2、数据准备：
矢量文件：文件point.dbf，point.shp，point.shx，文件river.dbf，river.shp，river.shx，polygon.dbf，文件polygon..shp，polygon.shx，文件polygon1.dbf，polygon1.shp，polygon1.shx，文件polygon2.dbf，polygon2.shp，polygon2.shx，文件overlay.dbf，overlay.shp，overlay.shx，文件roads.dbf，roads.shp，roads.shx，文件roads1.dbf，roads1.shp，roads1.shx。
栅格文件：文件夹gridline1和info，文件夹gridline2和info，文件夹query6和info，文件夹query5和info。
三、 实验内容
1、基于栅格数据的复合分析
(1) 原理验证
1) 添加gridline1层面和gridline2层面，并利用【Analysis】菜单中的【Map Calculator】命令进行栅格层面的叠加，观察其结果；
2) 应用【Analysis】菜单中的【Reclassify】命令，分别将gridline1和gridline2层面进行重分类，使得原有的No data数据为0，再进行叠加运算，观察其结果与上一步的结果有何不同，考虑其原因。
(2) 应用分析
query5层面和query6层面分别表示的是一些小工业镇区的城市化范围及其附近一个自然保护区的保护范围，要求利用栅格数据的复合分析，计算出这些小工业镇区的城市化范围。
方法：添加query5层面和query6层面，从【Analysis】菜单中选择【Map Query】命令，执行not(query6)andquery5运算，得出所求的范围。
2、基于矢量数据的叠置分析
Arcview中矢量数据的叠置分析是在Geoprocessing模块的支持下进行的，首先在【File】菜单中选择【Extension】命令，添加Geoprocessing模块（如图1）。
<v:group><v:shapetype><v:stroke joinstyle="miter"></v:stroke><v:path connecttype="rect" gradientshapeok="t"></v:path></v:shapetype><v:shape><v:textbox style="MARGIN-TOP: 4.187pt; LEFT: auto; MARGIN-LEFT: 7.218pt; WIDTH: 147.625pt; TOP: auto; HEIGHT: 16.218pt; mso-next-textbox: #_x0000_s1071">
<TABLE cellSpacing=0 cellPadding=0 width="100%">

<TR>
<TD>
<DIV>
图1. 添加Geoprocessing模块</DIV></TD></TR></TABLE></v:textbox></v:shape><v:shapetype><v:stroke joinstyle="miter"></v:stroke><v:formulas><v:f eqn="if lineDrawn pixelLineWidth 0 "></v:f><v:f eqn="sum @0 1 0 "></v:f><v:f eqn="sum 0 0 @1 "></v:f><v:f eqn="prod @2 1 2 "></v:f><v:f eqn="prod @3 21600 pixelWidth "></v:f><v:f eqn="prod @3 21600 pixelHeight "></v:f><v:f eqn="sum @0 0 1 "></v:f><v:f eqn="prod @6 1 2 "></v:f><v:f eqn="prod @7 21600 pixelWidth "></v:f><v:f eqn="sum @8 21600 0 "></v:f><v:f eqn="prod @7 21600 pixelHeight "></v:f><v:f eqn="sum @10 21600 0 "></v:f></v:formulas><v:path connecttype="rect" gradientshapeok="t" extrusionok="f"></v:path><lock v:ext="edit" aspectratio="t"></lock></v:shapetype><v:shape><v:imagedata src="5.files/image001.jpg" title="1"></v:imagedata></v:shape><w:wrap type="none"></w:wrap><w:anchorlock></w:anchorlock></v:group>
(1) Dissolve融合分析
1) 添加polygon层面并激活，从【Analysis】菜单中选择【Geoprocessing Wizard】命令，选择Dissolve features based on anattribute命令，执行下一步（如图2）。
<v:group><v:shape><v:textbox style="MARGIN-TOP: 3.718pt; LEFT: auto; MARGIN-LEFT: 7.187pt; WIDTH: 147.625pt; TOP: auto; HEIGHT: 16.187pt; mso-next-textbox: #_x0000_s1080">
<TABLE cellSpacing=0 cellPadding=0 width="100%">

<TR>
<TD>
<DIV>
图2. 选择根据属性值融合命令</DIV></TD></TR></TABLE></v:textbox></v:shape><v:shape><v:imagedata src="5.files/image003.jpg" title="2"></v:imagedata></v:shape><w:wrap type="none"></w:wrap><w:anchorlock></w:anchorlock></v:group>
2) 在随后出现的对话框中选择要进行融合的层面，要进行融合的属性字段以及确定融合后的结果文件的名称与路径，完成后按Finish键（如图3）。
<v:group><v:shape><v:textbox style="MARGIN-TOP: 3.812pt; LEFT: auto; MARGIN-LEFT: 7.187pt; WIDTH: 111.625pt; TOP: auto; HEIGHT: 16.218pt; mso-next-textbox: #_x0000_s1083">
<TABLE cellSpacing=0 cellPadding=0 width="100%">

<TR>
<TD>
<DIV>
图3. 设定融合的参数</DIV></TD></TR></TABLE></v:textbox></v:shape><v:shape><v:imagedata src="5.files/image004.jpg" title="3"></v:imagedata></v:shape><w:wrap type="none"></w:wrap><w:anchorlock></w:anchorlock></v:group>
3) 点击融合后的结果Disslv1-polygon层面，观察其属性，与polygon层面进行比较。可以看到，Disslv1-polygon层面是将polygon层面中属性字段type值一样的要素进行合并。(如图4)
<v:group><v:shape><v:textbox style="MARGIN-TOP: 4.281pt; LEFT: auto; MARGIN-LEFT: 7.187pt; WIDTH: 120.625pt; TOP: auto; HEIGHT: 16.187pt; mso-next-textbox: #_x0000_s1086">
<TABLE cellSpacing=0 cellPadding=0 width="100%">

<TR>
<TD>
<DIV>
图4. 融合后的结果层面</DIV></TD></TR></TABLE></v:textbox></v:shape><v:shape><v:imagedata src="5.files/image006.jpg" title="4"></v:imagedata></v:shape><w:wrap type="none"></w:wrap><w:anchorlock></w:anchorlock></v:group>
(2) Merge合并分析
1) 添加polygon1、polygon2层面并激活，执行【Geoprocessing Wizard】命令，选择Merge themes together命令，执行下一步。
2) 从列表中选择要进行合并的层面，确定输出的文件名称、路径，执行操作。
3) 由结果层面及其属性可以看出，进行合并后，它们的属性及图形都进行了合并。
(3) Clip裁剪分析
1) 添加并激活overlay、polygon层面，执行【Geoprocessing Wizard】命令，选择Clip one theme based on another命令，执行下一步。
2) 从列表中选择要进行裁剪的层面（input theme）以及用来裁剪的层面，确定输出的文件名称、路径，执行操作。
<wrapblock><v:group><v:shape><v:textbox style="MARGIN-TOP: 4.031pt; LEFT: auto; MARGIN-LEFT: 7.281pt; WIDTH: 129.625pt; TOP: auto; HEIGHT: 16.218pt; mso-next-textbox: #_x0000_s1040">
<TABLE cellSpacing=0 cellPadding=0 width="100%">

<TR>
<TD>
<DIV>
图5. Clip操作的参数设定</DIV></TD></TR></TABLE></v:textbox></v:shape><v:shape><v:imagedata src="5.files/image008.jpg" title="5"></v:imagedata></v:shape><w:wrap type="topandbottom"></w:wrap></v:group></wrapblock> 注意进行裁剪的层面与用来裁剪的层面的不同，若只想对进行裁剪的层面或用来裁剪的层面中的某一个或几个要素进行操作，可将其选中后再执行Geoprocessing Wizard操作。在这一步时，就可选择Use Selected features only选项（如图5）。
3) 由结果层面及其属性可以看出，裁剪分析是以用来裁剪的整个层面与被裁剪的层面进行运算的，输出结果将继承进行裁剪的层面的所有属性。
(4) Intersect相交分析
1) 添加并激活overlay、polygon层面，执行【Geoprocessing Wizard】命令，选择Intersect two themes命令，执行下一步。
2) 从列表中选择要进行相交分析的层面（input theme）以及用来相交的层面，确定输出的文件名称、路径，执行操作。
3) 由结果层面及其属性可以看出，相交分析是以用来相交的层面的中的要素（全部或一部分）与被相交的层面进行相交运算，输出结果将继承两层面的所有属性。
由相同的两层面进行不同的分析（Clip和Intersect），比较它们的结果有什么不同，体会它们的原理及应用。
(5) Union合并分析
1) 添加并激活overlay、polygon层面，执行【Geoprocessing Wizard】命令，选择Union two themes命令，执行下一步。
2) 从列表中选择要进行合并分析的层面（input theme）以及用来合并的上覆层面，确定输出的文件名称、路径，执行操作。
3) 由结果层面及其属性可以看出，合并分析是以用来合并的层面的中的要素（全部或一部分）与被合并的层面中的要素（全部或一部分）进行合并运算，输出结果将合并两层面的所有要素，并继承两层面的所有属性。
将Union分析的结果与前两种分析的结果进行比较，注意观察Union分析的特点。
(6) Assign data by location按空间位置的属性联接分析
联接分析是将两个层面按相同的位置进行属性联接，根据数据类型的空间关系不同，可分为三种：最近距离、是否在内部和是否是其一部分。其中，进行联接的层面的属性表称为目标表，用来联接的属性表称为源表。
1) Nearest最近距离
当一个点层面（源层面）与一个点层面或线层面（目标层面）进行联接分析时，进行最近距离的联接分析，同时增加一个距离字段到目标表中，用来记录目标层面中的要素与最近点的距离。
I、添加并激活point、river层面，执行【Geoprocessing Wizard】命令，选择Assign data by location命令，执行下一步。
<v:group><v:shape><v:textbox style="MARGIN-TOP: 4.312pt; LEFT: auto; MARGIN-LEFT: 7.781pt; WIDTH: 147.625pt; TOP: auto; HEIGHT: 16.187pt; mso-next-textbox: #_x0000_s1074">
<TABLE cellSpacing=0 cellPadding=0 width="100%">

<TR>
<TD>
<DIV>
图6. Nearest类型的联接分析</DIV></TD></TR></TABLE></v:textbox></v:shape><v:shape><v:imagedata src="5.files/image010.jpg" title="6"></v:imagedata></v:shape><w:wrap type="none"></w:wrap><w:anchorlock></w:anchorlock></v:group>
II、从列表中选择目标层面（the theme to assign data to）以及源层面（the theme to assign data from），确定输出的文件名称、路径，执行操作。（结果见图6）
2) Inside是否在内部
当一个面层面（源层面）与一个点（或线、面）层面（目标层面）进行联接分析时，进行是否在内部的联接分析。
I、添加并激活point、river层面，执行【Geoprocessing Wizard】命令，选择Assign data by location命令，执行下一步。
<v:group><v:shape><v:textbox style="MARGIN-TOP: 4.218pt; LEFT: auto; MARGIN-LEFT: 7.781pt; WIDTH: 138.625pt; TOP: auto; HEIGHT: 16.218pt; mso-next-textbox: #_x0000_s1077">
<TABLE cellSpacing=0 cellPadding=0 width="100%">

<TR>
<TD>
<DIV>
图7. Inside类型的联接分析</DIV></TD></TR></TABLE></v:textbox></v:shape><v:shape><v:imagedata src="5.files/image012.jpg" title="7"></v:imagedata></v:shape><w:wrap type="none"></w:wrap><w:anchorlock></w:anchorlock></v:group>
II、从列表中选择目标层面（the theme to assign data to）以及源层面（the theme to assign data from），确定输出的文件名称、路径，执行操作。（结果见图7）
3) Part of是否是其一部分
当两个线层面进行联接分析时，属于是否是其一部分的联接分析。
以下两个线层面roads，roads1分别表示的是某一公路网及这一公路网中需要维修的路段，现需要将需维修路段的属性添加到公路网中，方法如下：
I、添加并激活roads、roads1层面，执行【Geoprocessing Wizard】命令，选择Assign data by location命令，执行下一步。
II、 <wrapblock><v:group><v:shape><v:textbox style="MARGIN-TOP: 4.031pt; LEFT: auto; MARGIN-LEFT: 7.187pt; WIDTH: 165.625pt; TOP: auto; HEIGHT: 16.218pt; mso-next-textbox: #_x0000_s1049">
<TABLE cellSpacing=0 cellPadding=0 width="100%">

<TR>
<TD>
<DIV>
图8. 公路网的需维修信息的联接</DIV></TD></TR></TABLE></v:textbox></v:shape><v:shape><v:imagedata src="5.files/image014.jpg" title="8"></v:imagedata></v:shape><w:wrap type="topandbottom"></w:wrap></v:group></wrapblock> 从列表中选择目标层面（the theme to assign data to）以及源层面（the theme to assign data from），确定输出的文件名称、路径，执行操作。（结果见图8）
四、 实习报告要求
将所做工作以幻灯片形式做以汇报，内容包括原理、过程、结果。

喜欢0 评分0

GIS麦田守望者，期待与您交流。

举报回复

cheng_312 路人甲注册日期2006-06-06 发帖数48 QQ 铜币273枚威望0点贡献值0点银元0个加关注写私信	1楼^# 发布于：2007-05-14 10:53 <P>(6) Assign data by location按空间位置的属性联接分析</P> <P>这个功能在arcgis9.0里到哪去了哦？</P>
	举报回复(0) 喜欢(0) 评分

ynkm 路人甲注册日期2004-05-26 发帖数264 QQ 铜币27枚威望0点贡献值0点银元0个加关注写私信	2楼^# 发布于：2005-01-06 11:48 <img src="images/post/smile/dvbbs/em01.gif" /><img src="images/post/smile/dvbbs/em02.gif" />
	举报回复(0) 喜欢(0) 评分

superfond 路人甲注册日期2004-11-11 发帖数101 QQ 铜币280枚威望0点贡献值0点银元0个加关注写私信	3楼^# 发布于：2005-01-06 01:44 <img src="images/post/smile/dvbbs/em01.gif" /><img src="images/post/smile/dvbbs/em02.gif" />
	举报回复(0) 喜欢(0) 评分