《实习3-MAPGIS空间分析》由会员分享，可在线阅读，更多相关《实习3-MAPGIS空间分析（11页珍藏版）》请在金锄头文库上搜索。
1、实习3 MAPGIS空间分析与专题图制作实习目的： 掌握MAPGIS空间分析子系统的基本功能，以云南维西地质空间数据为例，掌握区空间分析方法,完成矿产分布与地层关系的专题图制作。一、空间分析方法这是一种将两个地图要素迭加产生一个新的要素层的操作，其结果是原来的要素被分割、剪断、套合，然后生成新的要素，新要素综合了原来两层要素所具有的属性。也就是说，空间迭加，不仅产生新的空间特征，还将输入特征的属性联系起来，产生新的属性。空间迭加分为矢量数据和栅格数据两种类型。对于矢量数据，采用矢量迭加方法，该方法对矢量的空间数据进行分割、剪断、套合等操作，对和矢量相关的属性进行连接，迭加结果是新的矢量数据和属
2、性数据。对于栅格数据则采用栅格加权迭加方法，该方法将两个栅格文件的对应元素加权相加，作为迭加结果的对应元素。比如，为了分析矿点和落雪组（Pt1l）地层的相关关系，在检索出 Pt1l地层数据后，就需要将Pt1l 地层数据和矿点数据进行迭加，得到地层和矿点之间的套合关系和两层要素的属性连接表，以便进一步作定量统计分析。矢量空间分析包括以下三种类型：1）迭加分析。包括区对区迭加分析、线对区迭加分析、点对区迭加分析、区对点迭加分析和点对线迭加分析。2）缓冲区分析。包括点BUFFER分析、线BUFFER分析、区BUFFER分析。3）多层立体迭置。设有原文件A和B，迭加结果为文件C，其中:A文件属性为：
3、缺省字段， f1等B文件属性为： 缺省字段， f2等缺省字段随着图形类型的不同而不同，分别是：点文件：标志码。线文件：标志码、线长度。区文件：标志码、面积、周长。迭加过程如下：A + B = C其中文件的图形类型和文件相同，而属性则是文件与文件属性连接的结果。下面主要介绍空间分析中的叠加分析（下一次实习有缓冲区分析和属性检索）1.区对区迭加分析包括合并、相交、相减、判别四种方式。迭加结果用限影表示，迭加结果的属性为：标志码、面积、周长，f1、区号、f2其中区号为第二个文件的区号。合并：属于A或属于B的区域标志码 面积 周长 F1 标志码 面积 周长 f2 标志码 面积 周长 f1 区号 F21
4、 320.5 61.2 a 1 280.7 50.1 b 1 198.2 51.3 a2 122.3 42.1 a 1 b+ =3 158.4 53.4 1 b相交：属于A且属于B的区域标志码 面积 周长 f1 标志码 面积 周长 f2 标志码 面积 周长 F1 区号 f21 320.5 61.2 a + 1 280.7 50.1 b = 1 122.3 42.1 a 1 b相减：属于A不属于B的区域标志码 面积 周长 f1 标志码 面积 周长 f2 标志码 面积 周长 F1 区号 F21 320.5 61.2 a + 1 280.7 50.1 b = 1 198.2 51.3 a 判别：属于
5、A的区域标志码 面积 周长 F1 标志码 面积 周长 f2 标志码 面积 周长 f1 区号 F21 320.5 61.2 a 1 280.7 50.1 b 1 198.2 51.3 a+ =2 122.3 42.1 a 1 b2.线对区迭加分析包括相交、判别、相减三两种方式，迭加结果文件仍然是线文件，迭加结果的属性为：标志码、线长度、f1、区号、f2线图元用标号表示。相交：穿过区域的线段部分标志码 线长度 f1 标志码 面积 周长 f2 标志码 线长度 f1 区号 f21 167.0 a + 1 320.5 61.2 b = 1 80.8 a 1 b判别：所有线图元标志码 线长度 f1 标志码
6、 面积 周长 f2 标志码 线长度 F1 区号 f21 167.0 a 1 320.5 61.2 b 1 32.2 a2 80.8 a 1 b+ =3 44.0 a相减：区域以外的线段标志码 线长度 f1 标志码 面积 周长 f2 标志码 线长度 f11 167.0 a 1 320.5 61.2 b 1 32.2 a+ =3 44.0 a3.点对区迭加分析包括相交、判别、相减三种方式，迭加结果文件仍然是点文件，结果属性为：标志码、f1、区号、f2相交：落在区域上的点标志码 f1 标志码 面积 周长 f2 标志码 f1 区号 f21 p1 1 320.5 61.2 b 1 p1 1 b2 p2
7、2 p2 1 b3 p3= =判别：所有点图元标志码 f1 标志码 面积 周长 f2 标志码 f1 区号 f21 p1 1 320.5 61.2 b 1 p1 1 b2 p2 2 p2 1 b3 p3+ =3 p3相减：区域以外的点图元标志码 f1 标志码 面积 周长 f2 标志码 f11 p1 1 320.5 61.2 b 1 P32 p23 p3+ =4.区对点迭加分析分为相减和相交两种方式。迭加结果为区文件，结果属性和原始区文件相同。相交：保留那些有点落在上面的区域标志码 面积 周长 F1 标志码 f2 标志码 面积 周长 f11 320.5 61.2 b 1 p1 1 320.5 61
8、.2 b2 150.7 45.1 c 2 p2+3 p3=相减：保留那些没有点落在上面的区域标志码 面积 周长 f1 标志码 f2 标志码 面积 周长 f11 320.5 61.2 b 1 p1 2 150.7 45.1 c2 150.7 45.1 c 2 p2+3 p3=5.点对线迭加分析迭加结果为点文件，该方法保留所有点，找到距离某点最近的线并计算出点线之间的距离，然后将线号和点线距离记录到属性中。点线距离定义如下： 对任意点D和曲线L，假设L由n个离散点d0、d1、d2dn构成，则D到d0、d1dn的距离分别为S0、S1Sn；D到直线段（d0，d1）、（d1，d2）（dn-1，dn）的法
9、线距离分别为：lidii到 的 法 线 距 离(,)1若 法 线 距 离 存 在若 法 线 距 离 不 存 在那么点D到曲线L的距离S=min（S0，S1，S2Sn，l1，l2ln）。迭加结果的属性为： 标志码、f1、线号、点线距离、f2标志码 f1 标志码 线长度 f2 标志码 f1 线号 点线距离 f21 p1 1 23.6 l1 1 p1 1 5.2 l12 p2 2 67.0 l2 2 p2 1 4.4 l13 p3+ =3 p3 2 3.8 l2二、实习内容与步骤1.实习资料云南维西地质空间数据工程文件：维西地质图.mpj，包括矿产.wt、断层.wl、地层界线.wl、岩浆界线.wl、
10、地层.wp、岩浆.wp、脉岩.wp、以及图例.wt、图例.wl、图例.wp。维西幅地质图比例尺为1:20万，图幅号为 G-47-(10)经纬度范围： 99:00:00-100:00:00E，26:40:00-27:20:00N2.MAPGIS系统设置在MAPGIS的主界面中选定系统设置或在编辑等子系统中选定设置修改目录环境，对工作目录、系统库目录等参数进行修改,如下图：3、维西地层区文件与矿产点文件的空间相交分析选定空间分析，进入空间分析子系统；选定文件，分别装入矿产点文件和地层区文件；选择空间分析区空间分析区对点相交分析，进行地层区文件与矿产点文件的空间相交分析，如右图； 保存分析结果的区文
11、件。4、专题图底制作根据维西幅地质图图幅号G-47-(10)和起点经纬度（纬度2640，经度99），选择实用服务投影变换系列标准图框根据图幅号生成图框，生成1：20万图框，如下图；保存1：20万图框的点、线、区文件，作为专题图底图。5、整图变换由于底图坐标与空间分析结果图坐标存在平移现象，需要进行整图变换。以底图内框的左下角坐标为原点，对空间分析结果文件和矿产点文件、图例点线面文件进行整图平移。在编辑子系统中，打开或装入空间分析结果文件和矿产点文件、图例点文件、图例线文件和图例面文件；选择其它整图变换键盘输入参数，输入X、Y位移量（X=-2500,Y=-14753.75），如上图。换名保存平移后的点、线、区文件。6、制作专题图创建新的专题工程文件，逐次添加底图和平移后的空间分析结果、图例区、线、点文件，形成矿产分布与地层关系的专题图，如下图。7、作业提交云南维西矿产分布与地层关系的专题图工程文件和点、线、面文件。将工程文件和相关点、线、面文件放在一个文件夹，双击工程文件，要能打开所有文件；压缩文件夹名为：实习3-班号-学号-姓名，下星期二（4月10日）由每个班的学习委员收齐后，拷贝给我。
