mapgis实验三 空间分析

实验报告实验项目3：空间分析综合实验课程：_GIS软件及其应用

实验指导

空间分析综合实验

所属课程名称：GIS软件及其应用

实验属性：综合

实验学时：6

（1）指导思想

在具体工作中会遇到空间分析的问题，本实验拟解决两个实际的问题，通过问题的解决，让同学熟练掌握软件的相关功能，开阔解决问题的思路。

（2）实验目的和要求

目的和要求：

1、掌握空间分析的基本方法

2、掌握图件中相关属性的转移办法

3、掌握图件空间分析求的相关面积和影响测算方法

4、解决实际工作中容易碰到的问题

一、实验步骤

(一)矢量化房屋与湖面。

将数据文件导入新建的文件工程中，在“工程文件”列表中单击右键，选择“新建区”，使用【区编辑】|【输入弧段】工具，照着已有的房屋矢量化。完成后选择【区编辑】|【区编辑】|【输入区】，再单击各个房屋即可得到区文件，命名为“房屋”。

使用相同的办法新建一个区文件，矢量化湖泊，并且命名为“湖泊”。结果如图1.

图1

(二)将居民数、楼房号变成相应属性

在区文件“房屋”的属性中添加字段“楼房号”与“居民数”，并输入如图2的数据。其中居民数应该使用整型字段而不是双精度，因为人数不可能出现小数位，只是当时没有在意这点，好在也不影响结果。

图2

将点文件的楼房号与居民数分别输入到区文件“房屋”的相应字段中。打开【库管理】|【属性库管理】，在菜单栏中点击【文件】|【装区文件】，将“房屋”打开，选择【属性】|【编辑单个属性】|【编辑单个区属性】，在弹出的对话框中一个个输入对应值。如图3. 输入完成后最终结果如图4.

图3 图4

输入属性这个步骤，本来应该使用属性表连接的方式比较快，不过在这个实验中，两个点文件本身都没有自带的相应属性，所以用不成。如果两个点文件自带属性，那么属性表连接具体步骤如下：打开点文件，点击【属性】|【属性输出】，选择

id、居民点、楼房号属性。输出类型为默认的mapgis表格（*.wb），输出文件则

自己设置，再打开需要连接这两种属性的区文件“房屋”，点击【属性】|【属性

连接】，在弹出的对话框，“连接文件”选择区文件“房屋”，“被连文件”则选择刚刚的表格文件，根据字段连接，在“连入字段”框框中选择需要的字段，确定即可。如下图。

(三)将事故点和湖泊进行BUFFer

打开【空间分析】|【空间分析】，弹出“空间分析子系统”界面，将“事故点”

文件装入。在菜单栏上选择【空间分析】|【缓冲区分析】|【输入缓冲区半径】，在弹出的对话框中输入数值100如图5.下一步点击【空间分析】|【缓冲区分析】|【求一个点缓冲区】，点击相应的点即可。结果如图6.保存文件为“事故点缓冲”。

缓冲区的半径根据题目给出的1:10000比例尺求出。

图5

图6

以同样的方法步骤将湖泊进行缓冲分析，只是在输入缓冲半径时输入15.结果如图7.保存文件为“湖泊缓冲”。

图7

(四)根据问题进行空间分析

1)事故点发生有毒物质泄漏，要求1公里内的人员都需撤离，哪些楼房受影响？共需

撤离多少人？若一辆车44人，需组织多少车辆帮助撤离？

将区文件“事故点缓冲”与“房屋”打开。选择【空间分析】|【区空间分析】|【区对区相交分析】，在弹出的对话框中分别选择叠加文件1：“事故点缓冲”、叠加文件2：“房屋”，点击确定，在弹出的模糊半径对话框中不改变默认值，直接确定，保存文件为“影响范围”，相交的结果如图8. 选择【属性分析】|【浏览属性数据】，可以得到受波及的楼房信息，如图9.

图8 图9

保存的结果会直接在界面中显示，无需再打开。选择【属性分析】|【单属性统计】，在弹出的“选择文件类型和统计属性”对话框中，选择文件“影响范围”，“选择属性”为“区属性”，“选择分析属性”为“居民数”，点击确定即可，结果如图

10.

图10

得出的结果为728人，728/44≈17辆车，需要17辆。

2)图左下角有一湖，现沿湖边要建设一绿化带，宽度为150米，绿化带面积有多大？

若每2平方米种植4棵树，略需多少树苗？

将区文件“湖泊”与“湖泊缓冲”打开，选择【空间分析】|【区空间分析】|【区对区相减分析】，在弹出的对话框中，叠加文件1为“湖泊缓冲”，叠加文件2为“湖泊”，点击确认即可，模糊半径不予理会。得到的结果如图11.

图11

绿化带的面积为8024.523是图上的单位，即8024.523mm²，将2m²化为图上单位，则是：2000000mm²/10000=200mm²。则每200平方毫米种植四棵树，最终结果

为：8024.523/200*4≈160棵。

二、总结

本实验我遇到的两个问题，一是属性的连接，二是面积的换算。

首先在属性连接过程中，我先尝试了属性表连接的方法，不过因为点文件的属性为空值，没有成功。于是用了一个一个值输入的方法，很繁琐，而且容易出错，需要有正确的图对照。其次面积换算时，我想将面积8024.523mm²化成实地面积，不过涉及的到0很多，一直算错，到现在为止我才搞明白，8024.523mm²=0.008024523m²*10000.还是将实际的2m²化为图上单位比较容易算。

总的来说两个问题都不是大问题，还是能完美解决的。

实验中使用的区文件，是我自己手动矢量化，因为当时没想到还能使用拓扑重建。

现在记住了可在检查完线之后将其转为弧段，再进行区拓扑错误检查，拓扑重建。这样得出的结果才是没有误差的正确的，最终面积为8045平米。