[ppt模板]空间数据的分析.ppt

1,如何为一公园选择合理位置？,空间数据分析,2,What is“Spatial Analysis”？,3,An Early Example:,1854年8月到9月英国伦敦霍乱流行时，当局始终找不到发病的原因，后来医生约翰斯诺(John Snow)参与调查。,4,他在绘有霍乱流行地区所有道路、房屋、饮用水机井等内容的1：6500比例尺地图上，标出了每个霍乱病死者的居住位置，得到了霍乱病死者居住分布图。,5,霍乱病死者居住分布图(John Snow,1854),6,霍乱死亡病例与水源的关系,7,斯诺博士分析了这张分布图，马上想到霍乱病源之所在死者住家都集中于饮用“布洛多斯托”井水的地方及周围。,8,在这个例子中，患者的居住地与饮用水井之间的空间位置关系提示了霍乱病的发病根源。,9,定 义,空间分析是基于空间数据的分析技术，它以地学原理为依托，通过分析算法，从空间数据中获取有关地理对象的空间位置、空间分布、空间形态、空间形成、空间演变等信息。根本目的在于：通过对空间数据的深加工和分析获取新的信息。,10,空间分析的主要内容,11,1、空间位置：借助于空间坐标系传递空间对象的定位信息，是空间对象表述的研究基础，即投影与转换理论。2、空间分布：同类空间对象的群体定位信息，包括分布、趋势、对比等内容。3、空间形态：空间对象的几何形态。4、空间距离：空间物体的接近程度。5、空间关系：空间对象的相关关系，包括拓扑、方位、相似、相关等。,各种空间分析方法,12,基本的空间分析包括以下方面：空间索引 7.1 空间查询 7.1 空间量算 7.2空间统计分析 7.3缓冲区分析7.4叠加分析 7.5数字高程模型7.6网络分析7.7 空间建模与空间决策支持系统 7.8空间分析与空间动态建模 7.10,面向应用的分析,简单的空间分析,复杂的空间分析,7.1 空间索引,空间索引就是根据空间对象的位置和形状或空间对象之间的某种空间关系按一定的顺序排列的一种数据结构，其中包含空间对象的概要信息，如对象的标识、外接矩形及指向空间对象实体的指针。,13,索引类型主要有：一、格网型空间索引二、BSP树空间索引三、KDB树空间索引四、R树和R+树五、CELL树,14,一、格网型空间索引,将研究区域用横竖线条划分大小相等和不等的格网，记录每个格网所包含的空间实体。当用户进行查询时，首先计算出用户查询对象所在格网，然后再在该格网中快速查询所选空间实体。,15,7.1 空间信息查询,16,空间数据的查询就是依靠数据库所储存的空间与属性信息来回答现实世界中一些应用问题。该查询工作可分为两步：首先，借助于空间索引，在空间数据库中快 速检索出被选空间实体；然后，根据空间数据和属性数据的连接即可 得到该空间实体的属性列表。,基于属性特征查询 一般来说，基于属性信息的查询操作主要是在属性数据库中完成的。基于空间关系和属性特征的查询（SQL）一种空间扩展SQL查询语言GeoSQL,17,18,GeoSQL的实现过程,19,查询过程分为三种类型：根据数据库中的数据及信息，直接回答人们的问题；通过逻辑表达式完成查询；根据现有数据模型，构造复杂模型，回答更为“复杂”的问题。,20,应用举例:为一家新开设的银行选址。要求：1、远离目前存在的银行;2、附近有大量的人口数量;3、结果以矢量数据的文件形式输出。,21,人口密度专题图(popden),银行位置专题图(bank.shp),街区专题图(street.shp),22,1、查询提取私人储蓄额高于10000000的银行。,2、计算生成距离远近分布图。,23,3、根据距离和人中密度专题图层，查询提取远离已有银行500且人口密度大于3000的区域。,4、将提取的栅格图层转为矢量图层，并与街区专题图进行叠加显示。,24,定位查询分层查询区域查询条件查询空间关系查询,最基本的查询功能，可查询图形与属性数据,一、查询内容,25,查询分层存放的图形与属性数据,定位查询分层查询区域查询条件查询空间关系查询,26,定位查询分层查询区域查询条件查询空间关系查询,查询区域内的图形与属性数据,27,定位查询分层查询区域查询条件查询空间关系查询,根据条件表达或查询图形与属性数据,28,定位查询分层查询区域查询条件查询空间关系查询,又称拓扑查询，面与面，线与线，点与点，点与线，点与面，线与面,二、GIS查询系统具备的功能,（1）空间数据处理的功能。查询处理器针对一定的数据模型和实际应用问题，必须设计一些程序模块和过程，对空间数据进行基本运算。这些基本运算和操作，除了常规的算数、统计及逻辑运算外，还应具有空间搜索、再分类、叠加、邻域、连通等空间数据的分析操作。,29,（2）空间数据处理的控制功能。这种控制是对计算机指令进行具体的空间数据运算。空间查询处理器的这种控制规范是以一种高级语言的形式来表达，称为查询语言，过去表现为按一定的词法和句法进行人机交互，现在表现为过程、控件和菜单等形式。,30,（3）构造应用模型的能力。为了面向专业领域的应用，GIS 的查询处理器必须提供一种开发语言和接口，使其具有使用和组合空间数据的基本操作，构造应用模型的能力。,31,7.2 空间量算,32,空间信息的自动化量算是地理信息系统所具有的重要功能，也是进行空间分析的定量化基础。,33,其中的主要量算有：一、几何量算 二、形状量算三、质心量算四、距离量算,一、几何量算,几何量算对点、线、面、体4类目标物而言，其含义不同的：点状目标：坐标；线状目标：长度、曲率、方向；面状目标：面积、周长等；体状目标：表面积、体积等。,34,1.线的长度计算,线长度可通过求算两点间直线距离得到。矢量数据结构中的求算：栅格数据结构中的求算：,35,2.面状地物的面积,36,矢量数据结构：面积和周长的计算。在平面直角坐标系中，计算面积时，对于每条边，计算y值以下面积，求其代数和，便是多边形面积值，周长则是线段之和。,栅格数据结构：统计具有相同属性值的格网数目。,37,二、形状量算,地物外形是影像处理中模式识别的一个重要部分。目标物的外观是多变的，很难找到一个准确的量对其进行描述。,38,面状地物形状量测两个基本考虑：空间一致性问题：即有孔多边形和破碎多 边形的处理；多边形边界特征的描述问题。,39,度量空间一致性最常用的指标是欧拉函数，用来计算多边形的破碎程度和孔的数目。欧拉函数的结果是一个数，称为欧拉数。欧拉函数的计算公式为：欧拉数=（孔数）（碎片数1）,40,41,关于多边形边界描述问题，最常用的指标包括多边形长短轴之比，周长面积比，面积长度比。,42,43,其中，P为目标物周长，A为目标物面积。,r1，目标物为紧凑型；r1，目标物为标准型；r1，目标物为膨胀型。,44,45,三、质心量算,46,质心是描述地理对象分布的一个重要指标。质心通常定义为多边形的几何中心或重心。,47,48,在有些情况下，质心描述的是分布中心，而不是绝对几何中心。,如果考虑其他一些因素，对其赋予一定的权重系数，称为加权平均中心。式中，Wi第i个离散目标物权重；Xi，Yi第i个离散目标物的坐标。,49,四、距离量算,50,距离描述了两个事物或实体之间的远近程度。各向同性表面简单距离；欧式距离各向异性表面耗费距离。非标准欧式距离,7.3 空间数据统计分析,51,52,常规统计分析主要完成对数据集合的均值、总和、方差、频度、峰度系数等的统计分析。其他的统计分析往往对空间数据及其相关信息进行抽象，采用更加复杂的数据模型进行分析，主要包括：回归分析 趋势分析 层次分析,主成分分析 聚类分析 判别分析,数据统计模型,从人类认识的角度来看有精确的和模糊的两种类型，因为绝大多数地理现象难以用精确的定量关系划分和表示，因此模糊的模型更为实用，结果也往往更接近实际，模糊评价一般经过四个过程：评价因子的选择与简化；多因子重要性指标(权重)的确定；因子内各类别对评价目标的隶属度确定；选用某种方法进行多因子综合。,53,一、主成分分析二、层次分析法三、系统聚类分析四、判别分析五、趋势面分析六、最优分割分级法,54,一、主成分分析,55,56,地理系统是多要素的复杂系统。在地理学研究中，多变量问题是经常会遇到的。变量太多，无疑会增加分析问题的难度与复杂性，而且在许多实际问题中，多个变量之间是具有一定的相关关系的。,因此，人们会很自然地想到，能否在相关分析的基础上，用较少的新变量代替原来较多的旧变量，而且使这些较少的新变量尽可能多地保留原来变量所反映的信息？,57,58,事实上，这种想法是可以实现的，主成分分析方法就是综合处理这种问题的一种强有力的工具。主成分分析是把原来多个变量划为少数几个综合指标的一种统计分析方法。从数学角度来看，这是一种降维处理技术。,例如：100个学生的数学、物理、化学、语文、历史、英语的成绩如下表（部分）。,59,二、层次分析法(AHP),60,日常生活中有许多决策问题。决策是指在面临多种方案时需要依据一定的标准选择某一种方案。例1 购物 买饭，则要依据色、香、味、价格等方面的因素选择某种饭菜。例2 择业面临毕业，可能有高校、科研单位、企业等单位可以去选择，一般依据工作环境、工资待遇、发展前途、住房条件等因素择业。,61,层次分析法是美国运筹学家Saaty教授于二十世纪80年代提出的一种实用的多方案或多目标的决策方法。其主要特征是，它合理地将定性与定量的决策结合起来，按照思维、心理的规律把决策过程层次化、数量化。该方法自1982年被介绍到我国以来，以其定性与定量相结合地处理各种决策因素的特点，以及其系统灵活简洁的优点，迅速地在我国社会经济各个领域内，如能源系统分析、城市规划、经济管理、科研评价等，得到了广泛的重视和应用。,62,层次分析法简介,一、层次分析法基本原理,63,分解,建立,确定,计算,判断,实际问题,层次结构,多个因素,诸因素的相 对重要性,权向量,综合决策,二 层次分析法的基本步骤,64,1 建立层次结构模型,2 构造成对比较矩阵,3 层次单排序及一致性检验,4 层次总排序及其一致性检验,65,尺度,第 个因素与第 个因素的影响相同,第 个因素比第 个因素的影响稍强,第 个因素比第 个因素的影响强,第 个因素比第 个因素的影响明强,第 个因素比第 个因素的影响绝对地强,含义,比较尺度：（19尺度的含义）,2,4,6,8表示第个因素相对于第个因素的影响介于上述两个相邻等级之间。不难定义以上各尺度倒数的含义，根据。,66,三 层次分析法的优点和局限性,1 系统性 层次分析法把研究对象作为一个系统，按照分解、比较判断、综合的思维方式进行决策，成为系统分析的重要工具。,2 实用性 层次分析法把定性和定量方法结合起来，能处理许多用传统的最优化技术无法着手的实际问题，应用范围很广，同时，这种方法使得决策者与决策分析者能够相互沟通，决策者甚至可以直接应用它，这就增加了决策的有效性。,67,3 简洁性 具有中等文化程度的人即可以了解层次分析法的基本原理并掌握该法的基本步骤，计算也非常简便，并且所得结果简单明确，容易被决策者了解和掌握。,该法的局限性主要表现在以下几个方面：,第一 只能从原有的方案中优选一个出来，没有办法得出更好的新方案。,68,第二 该法中的比较、判断以及结果的计算过程都是粗糙 的，不适用于精度较高的问题。第三 从建立层次结构模型到给出成对比较矩阵，人主观 因素对整个过程的影响很大，这就使得结果难以让 所有的决策者接受。当然采取专家群体判断的办法 是克服这个缺点的一种途径。,三、系统聚类分析,69,首先将n个样本各自成一类，然后规定类与类之间的距离，选择距离最近的两类合并成一个新类，计算新类与其他类的距离，再将距离最小的两类进行合并。这样每次减少一类，直到达到所需的分类数或所有样本都归为一类为止。,70,一、聚类要素的数据处理,在聚类分析中，聚类要素的选择是十分重要的，它直接影响分类结果的准确性和可靠性。在地理分类和分区研究中，被聚类的对象常常是多个要素构成的。不同要素的数据往往具有不同的单位和量纲，其数值的变异可能是很大的，这就会对分类结果产生影响。因此当分类要素的对象确定之后，在进行聚类分析之前，首先要对聚类要素进行数据处理。,71,假设有m 个聚类的对象，每一个聚类对象都有n个要素构成。它们所对应的要素数据可用表3.4.1给出。,72,表3.4.1 聚类对象与要素数据,二、直接聚类法,原理 先把各个分类对象单独视为一类，然后根据距离最小的原则，依次选出一对分类对象，并成新类。每经过m-1次就可以把全部分类对象归为一类，这样就可以根据归并的先后顺序作出聚类谱系图。,73,图3.4.1 直接聚类谱系图,74,实例分析,下表给出了某农业生态经济系统各个区域单元的有关数据，下面我们运用系统聚类法，对该农业生态经济系统进行聚类分析，步骤如下：(1)采用欧氏距离测度21个区域单元之间的距离；(2)选用组平均法，计算类间的距离，依据不同的聚类标准（距离），对各样本（各区域单元）进行聚类，并作出聚类谱系图。,75,76,某农业生态经济系统各区域单元的有关数据,77,78,某农业生态经济系统区域单元的系统聚类（组平均法）谱系图,从聚类分析谱系图（图3.4.5）可以看出，在不同的聚类标准（距离）下，聚类结果不同，当距离标准逐渐放大时，21个区域单元被依次聚类。当距离为0时，每个样本为单独的一类；当距离为5，则21个区域单元被聚为16类；当距离为10，则21个区域单元被聚为9类；当距离为15，则21个区域单元被聚为5类；当距离为20，则21个区域单元被聚为3类；最终，当聚类标准（距离）扩大到25时，21个区域单元被聚为1类。,79,四、判别分析,80,有一些昆虫的性别很难看出，只有通过解剖才能够判别；但是雄性和雌性昆虫在若干体表度量上有些综合的差异。于是统计学家就根据已知雌雄的昆虫体表度量（这些用作度量的变量亦称为预测变量）得到一个标准，并利用这个标准来判别其他未知性别的昆虫。这样的判别虽然不能保证百分之百准确，但至少大部分判别都是对的，而且用不着杀死昆虫来进行判别了。,81,判别分析和前面的聚类分析主要不同点就是，在聚类分析中一般人们事先并不知道或一定要明确应该分成几类，完全根据数据来确定。而在判别分析中，至少有一个已经明确知道类别的“训练样本”，利用这个数据，就可以建立判别准则，并通过预测变量来为未知类别的观测值进行判别了。,82,五、趋势面分析,83,一、空间趋势面分析概述趋势面分析 用数学的方法，以数学模型来模拟（或拟合）地理数据的空间分布及其区域性变化趋势的方法。,84,趋势面的性质与特点是一种光滑的数学曲面，它能集中地代表地理数据在大范围内的空间分布变化趋势与实际上的地理曲面不同，它只是实际曲面的一种近似值。实际曲面包括趋势面和剩余（或离差）曲面两部分，即实际曲面趋势面+剩余曲面,85,二、趋势面分析的数学模型,（一）趋势面分析的数学原理设以Zi(xi，yi)表示某一地理特征值在空间上的分布。其中(xi，yi)为平面上点的坐标。任一观测点Zi可分解为两个部分，即,86,趋势面参数的确定（最小二乘法）使每一个观测值与趋势值的残差平方和为最小，即,87,（二）多项式趋势面的数学模型,88,89,多项式方程作为趋势面方程因为任何函数在一定范围内总可以用多项式来逼近，并可调整多项式的次数来满足趋势面分析的需要，一般来说，多项式的次数越高则趋势值越接近于观测值，而剩余值越小,90,7.4 缓冲区分析,91,邻近度（Proximity）描述了地理空间中两个地物距离相近的程度，其确定是空间分析的一个重要手段。缓冲区分析是解决临近度问题的空间分析工具之一。,92,所谓缓冲区就是地理空间目标的一种影响范围或服务范围。从数学的角度看，缓冲区分析的基本思想是给定一个空间对象或集合，确定它们的邻域，邻域的大小由邻域半径R决定。,93,94,对于一个空间实体Oi，其缓冲区定义为：,对于空间实体集合：,其缓冲区定义为：,95,96,缓冲区类型：一、基于矢量结构的缓冲区分析二、基于栅格结构的缓冲区分析,1.矢量点的缓冲区：如建立污染源缓冲区，该区不能有饮用水源通过,97,98,2.矢量线的缓冲区：,线的窄缓冲区线的宽缓冲区,99,公路噪声污染，在公路两侧建立缓冲区，该区内不建立居民区；为防止水土流失，河流两侧建立缓冲区，该区内森林不许砍伐。,100,3.矢量面的缓冲区：,101,湖泊周围一定范围内为水源涵养林，限砍。,缓冲区计算的基本问题是双线问题。双线问题有很多另外的名称，如图形加粗，加宽线，中心线扩张等，它们指的都是相同的操作,102,103,角分线法（简单平行线法）,角分线法的缺点是难以最大限度保证双线的等宽性，尤其是在凸侧角点在进一步变锐时，将远离轴线顶点。,104,凸角圆弧法,105,对于简单情形，缓冲区是一个简单多边形，但当计算形状比较复杂的对象或多个对象集合的缓冲区时，就会产生若干个自相交多边形。,106,107,自相交多边形分为两种情况：岛屿多边形和重叠多边形。,108,算法比较简单，核心问题是距离变换。,栅格数据,缓冲区,二、基于栅格结构的缓冲区分析,109,应用举例：土地适宜性评价,1）目的：找出适宜耕种的地区2）标准：在对现有林区中开发 道路沿线300m范围内不能种值；河流沿线500m范围内不能种植。3）准备数据：道路分布图 河流分布图 森林分布图,110,4）空间操作,111,112,113,114,115,7.5 叠加分析,116,地理信息系统的叠加分析，是将有关主题层组成的数据层面进行叠加，产生一个新数据层面的操作，其结果综合了原来两层或多层要素所具有的属性。叠加分析不仅包含空间关系的分析，还包含属性关系的比较。,117,GIS叠加分析可以分为以下几类：一、视觉信息叠加二、点与多边形叠加三、线与多边形叠加四、多边形叠加五、栅格图层叠加,118,一、视觉信息叠加,视觉信息叠加是将不同侧面的信息内容叠加显示在结果图件或屏幕上，以便研究者判断其相互空间关系，获得更为丰富的空间信息。,119,120,视觉信息叠加包括以下几类：1.点状图，线状图和面状图之间的叠加显示；2.面状图区域边界之间或一个面状图与其他专 题区域边界之间的叠加；3.遥感影像与专题地图的叠加4.专题地图与数字高程模型（DEM）叠加,视觉信息叠加不产生新的数据层面，只是将多层信息复合显示，便于分析。,121,二、点与多边形的叠加,点与多边形叠加，实际上是计算多边形对点的包含关系。除了进行点与多边形的几何关系计算，还要进行属性信息处理。,122,123,1.点对多边形的叠加：,相交：落在区域上的点；相减：区域以外的点；判别：所有的点图元。,124,2.多边形对点的叠加,叠加结果属性与原始区文件相同相交：保留那些有点落在上面的区域；相减：保留那些没有点落在上面的区域。,125,三、线对多边形的叠加,线对多边形的叠加是确定一图层上的弧段落在另一图层的哪个多边形内，以便为图层的每条新弧段建立新的属性。,126,127,相交：穿过区域的线段部分；相减：区域以外的线段部分；判别：所有线图元。,128,四、多边形与多边形的叠加,多边形与多边形的叠加是指同一地区、同一比例尺的两组或两组以上的多边形要素的数据文件进行叠加。叠加后产生是具有多重属性的新多边形。,129,130,根据叠加结果最后欲保留的空间特征的不同要求，一般的GIS软件都提供了以下四种类型的多边形叠加操作。,131,132,合并,判别,相交,相减,？,合并:属于A或属于B的区域;相交:属于A且属于B的区域;相减:属于A不属于B的区域;判别:属于A的区域.,133,由于矢量结构的精度有限，几何对象不可能完全匹配，叠加结果可能会出现一些碎屑多边形。,134,135,五、栅格数据叠加分析,基于栅格数据的叠加分析是参与分析的两个或多个图层的要素均为栅格数据。栅格数据的叠加算法虽然数据存储量比较大，但是运算过程比较简单。这个特点和前面讲的矢量数据刚好相反。,136,通用土壤流失方程 E=RKLSCP E:平均土壤流失量；R:降雨强度；K:土壤可蚀性；L:坡长；S:坡度；C:耕作因子；P:水土保持措施因子,137,栅格图层叠加的形式：A.基于数学运算的叠加运算B.二值逻辑叠加,138,A.基于数学运算的叠加运算,森林地区融雪经验模型 M=0.19T+0.17D M是融雪速度（cm/d）；T是空气温度；D是露点温度,139,B.二值逻辑叠加,常作为数据库查询的工具。这种数据库查询通常分为两步：首先，进行再分类操作，为每个条件创建一个 新图层，通常是二值图层，1代表符合 条件，0表示不符合条件。第二，进行二值逻辑叠加操作，得到想查询的 结果。,140,7.6 数字高程模型分析,141,数字地面模型与数字高程模型,数字地面（地形）模型（DTM）：是通过地表点集的空间坐标及其属性数据表示表面特征的地学模型。是带有空间位置特征和地面属性特征的数字描述。（Digital Terrain Model）数字高程模型（DEM）：DTM中属性为高程的要素。（Digital Elevation Model）,142,143,一、DEM的数据采集,以航空或航天遥感图像为数据源；以地形图为数据源；以地面实测记录为数据源。,144,这种方法是由航空或航天遥感立体像对，用摄影测量的方法建立空间地形立体模型，量取密集数字高程数据，建立DEM。,1.以航空或航天遥感图像为数据源；,145,2.以地形图为数据源:,主要以比例尺不大于1：1万的国家近期地形图为数据源，从中量取中等密度地面点集的高程数据，建立DEM。,146,3.以地面实测记录为数据源,用电子速测仪和电子手簿或测距经纬仪配合PC1500等袖珍计算机，在已知点位的测站上，观测到目标点的方向、距离和高差三个要素。计算出目标点的x、y、z三维坐标，存储于电子手簿或袖珍计算机中，成为建立DEM的原始数据。这种方法一般用于建立小范围大比例尺(比例尺大于1：5000)区域的DEM，对高程的精度要求较高。,二、DEM的数据表示,1.数学分块曲面表示法 2.规则格网表示法 3.不规则三角网(TIN)表示法,147,148,1.数学分块曲面表示法 这种方法把地面分成若干个块，每块用一种数学函数，以连续的三维函数高平滑度地表示复杂曲面，并使函数曲面通过离散采样点。,2.规则网格模型,将区域空间分为规则的网格单元（可以是正方形或三角形），每个单元对应一个数值。,149,优点：计算机处理以栅格为基础的矩阵很方便，使高程矩阵成为最常见的DEM；缺点：在平坦地区出现大量数据冗余；若不改变格网大小，就不能适应不同的地形条件。,150,151,3.不规则三角网(TIN)将区域划分为相邻的三角面网络，区域中任意点落在三角面顶点、线或三角形内，落在顶点其高程与顶点相同，落在线上则由两个顶点线性插值得到，落在三角形内则由三个顶点插值得到。,152,TIN表示法利用所有采样点取得的离散数据，按照优化组合的原则，把这些离散点(各三角形的顶点)连接成相互连续的三角面(在连接时，尽可能地确保每个三角形都是锐角三角形或是三边的长度近似相等)。,优点：采用不规则三角网减少网格方法的数据冗余。采用不规则三角网可根据情况减少野外作业量。相对平坦的地方采集点少,地形变化剧烈的地方采集点多.,153,154,三、DEM建立的主要方法,1.手工数字化等高线,2.从遥感影像上自动提取,155,由*.shp线文件生成的江西潋水河流域小尺度区域TIN模型,手工数字化等高线,156,汉源流沙河流域三维景观,从遥感影像上自动提取,四、DEM模型在GIS中的应用,1.坡度的计算2.坡向的计算3.剖面分析4.可视性分析,157,1.坡度的计算,在流域提取、泥石流分析和植物生长环境研究中都需要坡度分析。,158,规则网格模型中地表基本单元的坡度等于其法向量N与Z轴之夹角。,159,=v=,160,2.坡向分析,坡向是地表单元的法向量在OXY平面上的投影与X轴之间的夹角，坡向通常要换算成正北方向起算的角度。,161,=v=,162,3.剖面分析,163,剖面线的确定：既可以人工输入，也可以利用鼠标实时确定。,剖面图的绘制应在格网DEM或三角网DEM上进行。已知两点的坐标A(X1,X2)，B(x2,y2)，则可用求出两点连线与格网或三角网的夹角，以及各交点之间的距离，用内插法求出交点的高程.,164,线性内插法求某点高程,Z代表高程，则，Zb=(Zc-Za)/S2*S1+Za,165,A,B,C,S1,S2,然后按选定的垂直比例尺和水平比例尺，按距离和高程绘制出剖面图。,166,4.可视性分析,167,可视性分析的基本因子有两个:一个是两点之间的通视性(Intervisibility);另一个是可视域(ViewShed)，即对于给定的观察点所覆盖的区域。,168,169,O,A,B,C,P,五、其他方面的应用,1.由TIN获取任意点P的高程已知Q1、Q2、Q3三个顶点，高程线性内插 求P的Z,170,平面方程,Z=Z1-（x-x1)(y21z31-y31z21)+(y-y1)(z21x31-z31x21)/(x21y31-x31y21),2.由TIN进行曲面拟合由于在TIN中可以获取任意点高程,进行密集插值可获得拟合的地形曲面.,171,3.土石方计算在DEM基础上进行土石方计算非常方便。,172,4.利用DEM绘制等高线图：如图所示，利用DEM绘制等高线图，是以格网点高程数据或者将离散的高程数据由栅格追踪法原理转换为矢量等值线所产生的。,173,5.利用DEM绘制地面晕渲图：晕渲图是以通过模拟实际地面本影与落影的方法有效反映地形起伏的重要的地图制图学方法。在各种小比例尺地形图、地理图，以及各类有关专题地图上得到非常广泛的应用。,174,6.透视立体图的绘制：立体图是表现物体三维模型最直观形象的图形，它可以生动逼真地描述制图对象在平 面和空间上分布的形态特征和构造关系。,175,7.7 网络分析,176,177,网络由一组线状要素（一系列联结的弧段）相互联结组成的，是物质、信息流通的通道，非计算机网络。,网络模型是对现实世界网络的抽象。在模型中，网络由链(Link)、结点(Node)、站点(Stop)、中心(Center)和转向点(Turn)组成。建立一个好的网络模型的关键是清楚地认识现实网络的各种特性与以网络模型的要素(Link,Node,Stop,Center,Turn)表示的特性之间的关系。,178,网络组成要素,结点(Node)：网络中任意两条线段的交点链(Link)：连接两个结点的弧段。中心(Center)：网络中位于结点处，具有沿着链收集和发放资源能力的设施，如邮局、电站、水库等 站点(Stop)：资源沿着网络路径流动时被分配或收集的位置转向点(拐点，Turn)：链路相交处，资源流向发生改变的点,179,180,一、网络图论基础,1.网络图：由一些点及点之间的线段所组成的图形。这些图形与地图，与各种函数图形有较大的区别。,181,182,2.路与回路 图中的一条路，就是由图中的一个定点、一条边，再一个定点、一条边排列而成，而且要求排在它前面的定点和排在它后面的定点都是它的端点。对于有向路来说，要求排在每一条边之前和之后的顶点分别是这条边的起点和终点。,183,184,3.连通性：,185,186,(1),(2),二、网络分析的应用,1.服务点的最优区位问题2.最小生成树3.最短路径分析,187,1.服务点的最优区位问题：在城市管理中，利用GIS技术确定服务点的最优区位问题十分重要，如确定幼儿园、商场、消防队、医院等的最优位置，以达到服务、资源的最优配置。,188,算法A,设G是一个有n个顶点，即V=v1,v2,vn,m条边,即E=e1,e2,en,的无向连通图，那么对于每一个顶点vi，它与各顶点间的最短路径的长度为di1,di2,din，上式中，最大数称为顶点vi的最大服务距离，用e(vi)表示。为了得到服务点的最优区位，需要解决如下问题，即求出一个点vio，使得e(vio)具有最小的值。,189,190,v1,v2,v3,v4,v5,v6,3,3,4,4,3,3,2,2,1）计算G的距离表,2）计算每一行的最大值，得e(v1)=6,e(v2)=7,e(v3)=6,e(v4)=7,e(v5)=6,e(v6)=7;3）求得mine(vi)=6,从而，定出v1,v3,v5均是G的中心。,191,算法B,设G是一个有n个顶点，即V=v1,v2,vn,m条边,即E=e1,e2,en,的无向连通图，那么对于每一个顶点vi，它与各顶点间的最短路径的长度为di1,di2,din，设每个顶点有一个正负荷a(vi),i=1,2,n,现求出一个顶点vi，使得S(vi)=最小，此点被认为是G的中央点。,192,193,1）计算G的距离方阵,194,2）求算：,S(v1)=122.2，S(v2)=71.3，S(v3)=69.5，S(v4)=69.5，S(v5)=108.5，S(v6)=72.8，S(v7)=95.3，3）得到G的中央点是v3和v4。,195,2.最小生成树,生成树是图的极小连通子图。一个连通的赋权图G可能有很多的生成树。设T为图G的一个生成树，若把T中各边的权数相加，则这个和数称为生成树T的权树。在G的所有生成树中，权数最小的生成树称为G的最小生成树。,196,构造最小生成树的依据有两个：1）在网中选择n-1条边，连接网中的n个顶点。2）尽可能选取权值为最小的边。,197,最小生成树的算法克罗斯克法,1）先把图G中各边按权数从小到大重新排列，并取权数最小的一条边为T中的边；2）在剩下的边中，按顺序取下权数最小一条边，若该边与T中已有的边构成回路，则舍去，否则选进T中；3）重复2操作，直到有m-1条边被选进T中，这m-1条边就是G的最小生成树。,198,199,1,2,6,5,4,3,19,33,18,21,16,11,14,6,6,5,3.路径分析,在空间网络分析中，路径问题占有重要位置。人们常想知道在地理空间网络中，两指定结点间是否存在路径，如果有则特别希望找出其中的最短路径。,200,迪克斯查算法,为了求出最短路径，需先计算网络任意两点间的距离，并形成nn阶距离矩阵。W=Wi,j式中，Wi,j为网络中边ei,j的距离。在矩阵W中：当i,j间有边相连接时，wi,j0；当i,j间无边相连接时，wi,j=当i=j时，wi,j=0,201,202,V1,V2,V3,V4,V5,V6,V7,V8,1,2,3,3,3,2,3,3,3,1,1,4,2,最短路径搜索步骤如下：1）令起点标号为0，其他结点暂时标号为Wij。2）对未被定标的结点全部给出暂时标号，其值为min j的旧标号，（i的旧标号+wj）,这里i是前一步刚被标定的结点，wij是边eij的权，如果结点i和j不相邻接，wij=。,203,3）找出所有暂时标号的最小值，用它作为响应结点的固定标号。如果存在几个有统一最小标号值的结点，则可任取一个加以定标。4）重复进行2和3，直至指定的终点L被定标时为止。用此法可直接得到由起点K到其他结点的最短路径的长度，那就是该结点的定标数值。,204,上图中，从结点v1到v7的最短路径的标号过程如下：,205,其中括号内是暂时标号，没有括号的为定标。距离起点越近的定顶点，越先得到固定标号。采用回溯的方法，可以得到从起点K到其他结点的最短路径经由的结点。如：从结点v1到v7，最短路径长度为7.,206,7.8 空间分析模型与空间决策支持,207,208,空间分析与应用模型的关系,空间分析是基本的、解决一般问题的理论和方法，而一般应用模型是不可枚举的，是复合、复杂的，解决专门问题的理论和方法，它的解决应以空间分析的基本方法和算法模型为基础。（局部与整体的关系）,209,空间分析模型是指用于GIS空间分析的数学模型，是在GIS空间数据基础上建立起来的模型，是通过作用于原始数据和派生数据的一组顺序的、交互的空间分析操作命令，对一个空间决策过程进行的模拟。,一、空间分析模型的有关知识,210,1.GIS常用的分析模型,211,1)相关分析模型-用来分析研究各种地理要素数据之间相互关系的一种有效手段。2)趋势面分析模型-用来将现象的空间分布特征及其区域变化趋势模拟出来。3)预测模型-反映地理要素的动态发展规律，并用于预测分析。-常用回归模型。4）聚类模型描述各种地理要素数据之间的近似程度，相似的可合并。,212,地理信息系统不仅要完成管理大量复杂的地理数据的任务，更为重要的是要完成地理分析、评价、预测和辅助决策的任务，必须发展广泛的适用于地理信息系统的地理分析模型，这是地理信息系统走向实用的关键。模型：所谓模型，就是将系统的各个要素，通过适当的筛选，用一定的表现规则描写出来的简明的映象。通常表达了某个系统的发展过程或发展结果。,2.GIS空间分析模型分类,地学模型：是用来描述地理系统各个要素之问相互关系和客观规律的，它用信息的、语言的、数学的或其它表达形式，通常反映地学过程及其发展趋势或结果。是在对系统所描述的具体对象与过程，进行大量专业研究的基础上，总结出来的客观规律的抽象或模拟。地学模型也称为专题分析模型。,213,214,数学模型：是应用数学的语言和工具，对部分现实世界的信息(现象、数据)加以翻译、归纳的产物。数学模型经过演绎、推导，给出数学上的分析、预报、决策或控制，再经过解释回到现实世界，完成实践理论实践这一循环。,3.模型的作用和特点,215,应用模型是联系GIS应用系统与常规专业研究的纽带 应用模型是综合利用GIS应用系统中大量数据的工具 应用模型是GIS应用系统解决各种实际问题的武器 应用模型是GIS应用系统向更高技术水平发展的基础,216,（1）空间分析模型是联系GIS应用系统与专业领域的纽带，必须以广泛、深入的专业研究为基础；（2）空间分析模型是综合利用GIS中大量数据的工具，数据的综合分析和应用主要通过模型来实现；（3）空间分析模型是分析型和辅助决策型GIS区别于管理型GIS的一个重要特征，是解决空间分析和辅助决策问题的核心。,4.空间分析模型的意义,217,5.地理信息建模系统（GIMS）,是研究如何根据给定条件（如已知数据和约束条件）自动生成解决问题（如确定候选地址）的整个操作过程。它能支持面向用户的空间分析模型的定义、生成和检验的环境，支持与用户交互式的基于GIS的分析、建模和决策。GIMS是目前GIS研究的热点问题之一。,二、空间分析过程及其模型1.空间分析过程：空间分析的目的是解决某类与地理空间有关的问题，通常涉及多种空间分析操作的组合。空间分析的一般步骤是：评价分析的目的和评价准则；准备分析数据；进行空间分析操作；进行结果分析；解释、评价结果；结果数尺。,218,219,例1：道路拓宽改建过程中的拆迁指标计算,1）明确分析的目的和标准目的：计算由于道路拓宽而需拆迁的建筑物的面积和房产价值；道路拓宽改建的标准是：a)道路从原有的20m拓宽至60m；b)拓宽道路应尽量保持直线；c)部分位于拆迁区内的10层以上的建筑不拆除。2）准备进行分析的数据涉及两类信息：一类是现状道路图；另一类是分析区域内建筑物分布图及相关的信息；,3）进行空间操作a)选择拟拓宽的道路，根据拓宽半径，建立道路的缓冲区。b)将此缓冲区与建筑物层数据进行拓扑叠加，产生一幅新图，此图包括所有部分或全部位于缓冲区内的建筑物信息。4）进行统计分析a)对全部或部分位于拆迁区内的建筑物进行选择，凡部分落入拆迁区且楼层高于10层以上的建筑物，将其从选择组中去除，并对道路的拓宽边界进行局部调整。b)对所有需拆迁的建筑物进行拆迁指标计算。5）将分析结果以地图或表格的形式打印输出。,220,221,2.空间分析建模 地图分析方法是一种分析和综合空间数据（地图数据）的方法，其目的是有助于对空间数据的解释。要求地图数据结构、数据处理方法、过程控制要具有一定的规范和统一标准。地图数据规范要求空间数据具有一定的数据模型和数据结构，有一定的区域框架结构。,222,数据处理控制规范 系统对数据运算进行指令控制，包括识别图层，陈述执行函数