第七章 空间数据分析PPT讲稿.ppt

《第七章 空间数据分析PPT讲稿.ppt》由会员分享，可在线阅读，更多相关《第七章 空间数据分析PPT讲稿.ppt（97页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、第七章 空间数据分析第1页，共97页，编辑于2022年，星期一分析技术分析技术:空间图形数据的拓扑运算；非空间属性数据运算；空间和非空间数据的联合运算。空空间间分分析析是综合分析空间数据的技术的通称。空间分析有着十分丰富的内涵，它是构成地理信息系统的核心部分之一，在整个地理数据的应用中发挥着举足轻重的作用，也是GIS区别与其它信息系统的一个显著标志。第2页，共97页，编辑于2022年，星期一空间分析主要获得空间位置及其属性描述两方面的信息。主要目的：以空间数据库为基础解决用户的实际问题运用各种几何逻辑运算手段第3页，共97页，编辑于2022年，星期一7.1 空间索引在介绍空间索引之前，先谈谈什

2、么叫“索引“。对一个数据集做”索引“，是为了提高对这个数据集检索的效率。书的”目录“就是这本书内容的”索引“，当我们拿到一本新书，想查看感兴趣内容的时候，我们会先查看目录，确定感兴趣的内容会在哪些页里，直接翻到那些页，就OK了，而不是从第一章节开始翻，一个字一个字地找我们感兴趣的内容，直到找到为止，这种检索内容的效率也太低了，如果一本书没有目录，可以想象有多么不方便可见书的目录有多重要，索引有多重要。现在大家对索引有了感性认识，那什么是“空间索引“呢？”空间索引“也是”索引“，是对空间图形集合做的一个”目录“，提高在这个图形集合中查找某个图形对象的效率。目前，常见空间索引类型有BSP树树，KD

3、B树，树，R树树，R+树和树和CELL树树，空间索引的性能的优越直接影响空间数据库和地理信息系统的整体性能 第4页，共97页，编辑于2022年，星期一空间索引类型R 树是基于磁盘的索引结构，是 B 树(一维)在高维空间的自然扩展，易于与现有数据库系统集成，能够支持各种类型的空间查询处理操作，在实践中 得到了广泛的应用，是目前最流行的空间索引方法之一。R 树索引比较适合静态数据，例如用作底图的数据和不经常编辑的数据（属性数据的编辑除外）推荐建立 R 树索引。R 树空间索引方法是通过设计虚拟的矩形目标，将一些空间位置相近的目标对象，包含在这个矩形内，如下面左图所示，把这些矩形作为空间索引的节点，它

4、 含有所包含的空间对象的指针。例如，A 中包含空间对象 D、E、F，B 中包含空间对象 H、I、J、K，C 中包含空间对象 L、M。如下面右图所示为 R 树索引的例子。第5页，共97页，编辑于2022年，星期一在进行空间检索的时候，首先判断哪些矩形落在检索窗口内，再进一步判断哪些目标是被检索的内容。这样可以提高检索速度。特点特点（1）R树索引的空间检索效率非常高，检索精度也很高。（2）R树索引的索引数据和地理数据不存储在一起，且R树索引的更新涉及到R树的平衡和分裂，比较复杂，更新速度较慢，大规模编辑后的查询效率降低，需要重建空间索引，适用于只读数据或者不经常进行空间数据编辑的数据。适用情况适用

5、情况R树索引比较适合静态数据，例如用作底图的数据和不常被编辑的数据（属性数据的编辑除外）推荐建立R树索引。第6页，共97页，编辑于2022年，星期一四叉树是一种重要的层次化数据集结构，主要用来表达二维坐标下空间层次关系，实际上它是一维二叉树在二维空间的扩展。其基本思路是：把整个工作空间编码为0，用水平线和 垂直线将工作空间划分为相等的四部分，各每一部分按顺序分别编码为00、01、02和03，左上角为00，右下角为03；对划分后的格子按照同样的方式继续划分，比如对第00格继续划分，划分后的四个小格子分别编码为000、001、002和003。根据需要，还可以继续划分，比如第031格子继续划分为03

6、10、0311、0312和0313，如下图所示。以此类推，可以一直划分下去，划分的层数取决于研究区域对象的数量和平均大小。第7页，共97页，编辑于2022年，星期一假设上图所示的对象为整个工作空间的对象，则根据这个工作空间建立的四叉树结构如下图所示：现在在 SuperMap 中，四叉树最多允许分成13层。基于希尔伯（Hilbert）编码的排序规则，从四叉树中可确定索引类中每个对象实例的被索引属性值是 属于哪个最小范围。从而提高了检索效率。四叉树索引适用于小数据量的高并发编辑。第8页，共97页，编辑于2022年，星期一根据数据集的某一属性字段或根据给定的一个范围，将空间对象进行分类，通过索引进行

7、管理已分类的空间对象，以此提高查询检索速度。根据上述原理，创建图库索引有两种方式，根据字段索引或者根据范围索引。字段索引：即根据数据集的某一属性字段将空间对象进行分类，通过索引进行管理已分类的空间对象，以此提高查询检索速度。建议使用与位置信息相关的字段进行图库索引的 创建，例如对于全国县级行政区域图，可以使用表示行政区域代码的字段进行图库索引的创建。第9页，共97页，编辑于2022年，星期一范围索引：即根据给定的一个范围（图幅的长和宽）将空间对象进行分类，通过索引进行管理已分类的空间对象，以此提高查询检索速度。对于按标准比例尺分幅存储(如1:25万 数据、1:10万数据、1:5万数据等)的数据

8、合并到数据库中后生成的数据集，范围索引有着优异的效果，可以提供非常好的查询性能。在海量数据（即对象数超过百万个的数据）的显示和查询时优势尤其明显。建立图库索引后，在显示的时候，只显示你需要查看的那一幅或几幅地图，大大缩短显示时间。如下图所示为根据给定范围进行图库索引的示例：第10页，共97页，编辑于2022年，星期一图库索引示意图 特别地，图库索引还可以被缓存到本地（以二进制文件形式体现），这样可更进一步提高查询检索的速度。第11页，共97页，编辑于2022年，星期一动态索引结合了 R 树索引与四叉树索引的优点，提供非常好的并发编辑支持，具有很好的普适性。若不能确定数据适用于哪种空间索引，推荐

9、建立动态索引。动态索引采用划分多层网格的方式来组织管理数据，它的基本方法是将数据集按照一定的规则划分成相等或不相等的网格，记录每一个地理对象所占的网格 位置。在GIS中常用的是规则网格。当用户进行空间查询时，首先计算出用户查询对象所在的网格，通过该网格快速查询所选地理对象。可以优化查询操作。第12页，共97页，编辑于2022年，星期一动态索引示意图 第13页，共97页，编辑于2022年，星期一特点特点当数据集在进行浏览操作的时候，动态索引方式的速度会比较快；索引更新和并发的能力好；索引的空间检索精度高，准确度高。注意：此索引类型支持数据集的动态并发编辑。适用情况适用情况动态索引结合了R 树索引

10、与四叉树索引的优点，提供非常好的并发编辑支持，具有很好的普适性。若不能确定数据适用于哪种空间索引，推荐建立动态索引。补充说明补充说明建立动态索引后，用户可以按图幅进行本地缓存，可以大大提高查询和浏览的速度。第14页，共97页，编辑于2022年，星期一查询和定位空间对象，并对空间对象进行量算是查询和定位空间对象，并对空间对象进行量算是GIS的基本功能之一的基本功能之一，它是它是GIS进行高层次分析的基础。进行高层次分析的基础。在在GIS中，为进行高层次分析，往往需要查询定位中，为进行高层次分析，往往需要查询定位空间对象，并用一些简单的量测值对地理分布或现空间对象，并用一些简单的量测值对地理分布或

11、现象进行描述，如长度、面积、距离等。实际上象进行描述，如长度、面积、距离等。实际上，空空间分析首先始于空间查询和量算，它是空间分析的间分析首先始于空间查询和量算，它是空间分析的定量基础。定量基础。7.2 空间信息查询第15页，共97页，编辑于2022年，星期一1、概念空间查询的的一般问题是“有没有”，“是什么”，“在什么位置”，“怎么到达”2、查询对象图形中的信息属性表中的信息其他信息（例如关系等）第16页，共97页，编辑于2022年，星期一3、意义信息管理：通过查询可以获取数据信息，进行信息管理和数据更新特定信息提取：提取需要的信息，去除无用信息，便于使用空间分析基础：查询结果一般是对查找信

12、息或数据的报告，研究对于这些信息单独提出进行相关分析。第17页，共97页，编辑于2022年，星期一图形和属性的互查是最常用的查询，主要有类：图形和属性的互查是最常用的查询，主要有类：1、按属性信息的要求来查询定位空间位置、按属性信息的要求来查询定位空间位置，称为称为“属性查图形属性查图形”。如在中。如在中国行政区划图上查询人口大于国行政区划图上查询人口大于4000万且城市人口大于万且城市人口大于1000万的省有哪些万的省有哪些?称为称为SQL查询查询.2、根据对象的空间位置查询有关的属性信息、根据对象的空间位置查询有关的属性信息，称为称为“图形查属性图形查属性”。如一。如一般的般的GIS软件都

13、提供一个软件都提供一个“INFO”工具，让用户利用鼠标，用点选、画工具，让用户利用鼠标，用点选、画线、矩形、圆、不规则多边形等工具选中地物，并显示所查询对象的线、矩形、圆、不规则多边形等工具选中地物，并显示所查询对象的属性列表，可进行有关统计分析。属性列表，可进行有关统计分析。第18页，共97页，编辑于2022年，星期一1、基于空间关系查询、基于空间关系查询 空间实体间存在多种空间关系，包括拓扑、距离、方位等。空间实体间存在多种空间关系，包括拓扑、距离、方位等。如查找满足下列条件的城市：如查找满足下列条件的城市：在京沪线的东部；距离京沪线不超过在京沪线的东部；距离京沪线不超过50公里公里;城市

14、人口大于城市人口大于100万万;城市区域面积城市区域面积5000平方公里平方公里.第19页，共97页，编辑于2022年，星期一简单的点线面相互关系拓扑查询包括：简单的点线面相互关系拓扑查询包括：面面查询面面查询:如与某个多边形相邻的多边形有哪些；如与某个多边形相邻的多边形有哪些；面线查询面线查询:如某个多边形内包含哪些线；如某个多边形内包含哪些线；面点查询面点查询:如某个多边形内有哪些点状地物；如某个多边形内有哪些点状地物；线面查询线面查询:如某条线经过的多边形有哪些；如某条线经过的多边形有哪些；线线查询线线查询:如与某条河流相连的支流有哪些；如与某条河流相连的支流有哪些；线点查询线点查询:如

15、某条道路上有哪些桥梁，某条输电线上有哪些变电站；点面查询如某条道路上有哪些桥梁，某条输电线上有哪些变电站；点面查询:如某个点落如某个点落在那个多边形内；在那个多边形内；点线查询点线查询:如某个结点由哪些线相交而成；如某个结点由哪些线相交而成；第20页，共97页，编辑于2022年，星期一2、基于空间关系和属性特征查询、基于空间关系和属性特征查询传统的传统的SQL并不能处理空间查询，对并不能处理空间查询，对GIS而言，需要而言，需要对对SQL进行扩展进行扩展,主要包括空间数据与属性数据的匹配等主要包括空间数据与属性数据的匹配等如地址匹配查询如地址匹配查询 根据街道的地址来查询事物的空间位置和属性信

16、息是根据街道的地址来查询事物的空间位置和属性信息是GIS特有的一种查询特有的一种查询功能，功能，这种查询利用地理编码，输入街道的门牌号，就可以知道大致的位置和所在的街区。这种查询利用地理编码，输入街道的门牌号，就可以知道大致的位置和所在的街区。它它对空间分布的社会、经济调查和统计很有帮助，只要在调查表中添加了地址，对空间分布的社会、经济调查和统计很有帮助，只要在调查表中添加了地址，GIS就可就可以自动地从空间位置的角度来统计分析各种经济社会调查资料。另外，这种查询也经常以自动地从空间位置的角度来统计分析各种经济社会调查资料。另外，这种查询也经常用于公用事业管理，事故分析等方面，如邮政、通讯、供

17、水、供电、治安、消防、医疗用于公用事业管理，事故分析等方面，如邮政、通讯、供水、供电、治安、消防、医疗等领域。等领域。第21页，共97页，编辑于2022年，星期一7.3 空间量算空间数量量算主要是指对空间对象的几何参数进行量算，以获得关于这些几何体的属性信息。在地理数据库中，主要的空间数量量算有：Distance(geometry,geometry):double：两个几何体的距离测量 Length(geometry):double：几何测量 Area(geometry):double：面积测量 Buffer(double):Geometry：缓冲区运算，返回与当前几何体距离小于或等于某个距离

18、的点的集合。ConvexHull():Geometry：凸包运算，返回一个包含当前几何体所有点的凸包。由于地球是一个椭球体，因此在空间数据库之中，所有的几何体都是带有地图投影与空间参照系的。空间数量量算都是在一定的投影与参照系下进行。对于若干几何体，由于不同参照系的椭球参数不同，如果把他们放在不同的空间参照系下进行空间数量量算，会得出不同的结果。第22页，共97页，编辑于2022年，星期一7.3.1 几何量算1长度长度 线状物体的长度是最基本的形态参数之一，在矢量线状物体的长度是最基本的形态参数之一，在矢量数据格式下，线由点组成，线状物体表示为一个坐标数据格式下，线由点组成，线状物体表示为一个

19、坐标串串(Xi,Yi)，而线长度可由两点间直线距离相加得到。，而线长度可由两点间直线距离相加得到。则线状物体长度的计算公式为：则线状物体长度的计算公式为：第23页，共97页，编辑于2022年，星期一空间信息量算空间信息量算2面积面积 多多边边形形的的面面积积是是一一个个重重要要指指标标。多多边边形形边边界界可可以以分分解解为为上上下下两两半半，其其面面积积就就是是上上半半边边界界下下的的积积分分值值与与下下半半边边界界下下的的积积分分值值之之差差。设设面面状状物物体体的的轮轮廓廓边边界界由由一一个个点点的的序序列列P1(x1,y1),P2(x2,y2),，Pn(xn,yn)表示，其面积为表示，

20、其面积为:第24页，共97页，编辑于2022年，星期一2面积面积YoXSS1S2S=S2-S1第25页，共97页，编辑于2022年，星期一YoXS1(X1,y1)(X2,y2)(X3,y3)(X4,y4)(X5,y5)S1=(x2-x1)(y1+y2)/2+(x3-x2)(y2+y3)/2+(x4-x3)(y3+y4)/2+(x5-x4)(y4+y5)/2第26页，共97页，编辑于2022年，星期一空间信息量算空间信息量算3弯曲度弯曲度弯弯曲曲度度是是描描述述线线状状物物体体弯弯曲曲程程度度的的一一个个重重要要参参数数，它它定定义义为为曲曲线线长长度度与与曲曲线线的的两两个个端端点点之之间间长

21、长度的比值，即度的比值，即:第27页，共97页，编辑于2022年，星期一 4.4.质心量算质心量算 质心是描述地理现象空间分布的一个重要指标。质质心心可可简简单单地地描描述述为为地地理理目目标标保保持持均均匀匀分分布布的的平平衡衡点点。质心通常定义为一个多边形或面的几何中心，当多边形比较简单时，计算很容易。当多边形形状复杂时，计算也更加复杂。质心量算可用于对地理分布变化的跟踪；计算目标物对周围地区的经济辐射范围。如应用质心量测分析人口变迁、土地类型变化等。第28页，共97页，编辑于2022年，星期一分两种情况：分两种情况：1）面面状状目目标标的的重重心心。可可以以理理解解为为多多边边形形内内的

22、的平平衡衡点点，正正如如一一块块均均质质木木块块被被悬挂起来的平衡点。悬挂起来的平衡点。空间信息量算空间信息量算第29页，共97页，编辑于2022年，星期一面面状状目目标标重重心心可可以以通通过过计计算算梯梯形形重重心心的的平平均均值值而而得得到到。将将多多边边形形的的各各个个顶顶点点投投影影到到x轴轴上上，就就得得到到一一系系列列梯梯形形（如如图图），所所有有梯梯形形重重心心的的联联合合就就确确定定了了整整个个多多边边形形的的重心。重心。空间信息量算空间信息量算按梯形计算重心位置按梯形计算重心位置第30页，共97页，编辑于2022年，星期一空间信息量算空间信息量算按梯形计算重心位置按梯形计算

23、重心位置设多边形的顶点序列设多边形的顶点序列设多边形的顶点序列设多边形的顶点序列(x(xi i,y,yi i)按顺时针编按顺时针编按顺时针编按顺时针编 码，则其重心的计算公式为：码，则其重心的计算公式为：码，则其重心的计算公式为：码，则其重心的计算公式为：其中，其中，其中，其中，和和和和 是第是第是第是第i i个梯形的重心的个梯形的重心的个梯形的重心的个梯形的重心的x x坐标和坐标和坐标和坐标和y y坐标，是梯形的面积。坐标，是梯形的面积。坐标，是梯形的面积。坐标，是梯形的面积。它们由下式得到：它们由下式得到：它们由下式得到：它们由下式得到：第31页，共97页，编辑于2022年，星期一按梯形计

24、算重心位置按梯形计算重心位置第32页，共97页，编辑于2022年，星期一空间信息量算空间信息量算可可理理解解为为其其分分布布中中心心。其其重重心心计计算算方方法法是是取取离离散散目目标标的的加加权权平平均均中中心心，它它是是离离散散目目标标保保持持均均匀匀分分布布的的平平衡衡点点。计计算算公式为：公式为：2）面状分布离散目标的重心）面状分布离散目标的重心其中，其中，i为离散目标物，为离散目标物，Wi为该目标物权重。为该目标物权重。Xi与与Yi为其坐标。为其坐标。第33页，共97页，编辑于2022年，星期一空间信息量算空间信息量算5.形状量算形状量算 当当把把城城市市作作为为单单个个面面状状目目

25、标标看看待待时时，可可以以直直接接使使用用面面状状目目标标的的形形状状系系数数，如如形形状状率率、圆圆形形率率、紧紧凑凑度度等等,这这些些指指标标计计算算较较简简单单，但但只只反反映映一一个个抽抽象的形状；象的形状；当当把把城城市市作作为为面面状状目目标标的的集集合合看看待待时时，可可以以使使用用放放射射状状指指数数、标标准准面面积积指指数数等等形形状状系系数数，这这些些指指标标计计算算较较复复杂杂，但但反反映映了了城城市市内内部部的的具具体体联联系系。在在多多数数指指标标中中，都都以以圆圆形形作作为为城城市市的的标标准准形形状。状。第34页，共97页，编辑于2022年，星期一一、空间信息分类

26、 空间信息分类方法是地理信息系统功能组成的重要组成部分。与地图上所负载的已经专门分类和处理的数据相比较，地理信息系统存储的数据具有原始数据的性质，这样用户就可以根据不同的使用目的对数据进行任意提取和分析。7.3 空间信息分类与统计分析第35页，共97页，编辑于2022年，星期一空间信息分类与统计分析1.主成分分析法 地理问题往往涉及大量的相互关联的自然和社会要素，太多的变量，无疑增加了分析问题的难度与复杂性。主成分分析就是把原来多个变量化为少数几个综合指标的一种统计分析方法。它克服了变量选择时的冗余和相关。要使这些较少的综合指标能够尽量多地反映原来较多指标的信息，那么它们的形式应当是原来变量指

27、标的线性组合，且它们之间必须是独立的。第36页，共97页，编辑于2022年，星期一空间信息分类与统计分析设有n个样本，p个变量指标 x1，x2，xp，它们的综合指标为 Z1，Z2，Zp，（mp）,则应有Lij是据相关系数矩阵求出的相应的特征向量。第37页，共97页，编辑于2022年，星期一空间信息分类与统计分析2.层次分析法 层次分析（AHP）是系统分析的数学工具之一，它把人的思维过程层次化、数量化，并用数学方法为分析、决策、预报、控制提供定量的依据。AHP方法把相互关联的要素按隶属关系分为若干层次，请有经验的专家对各层次各要素的相对重要性给出定量指标，利用数学方法综合专家意见给出各层次各要素

28、的相对重要性权值，作为综合分析的基础。第38页，共97页，编辑于2022年，星期一空间信息分类与统计分析选择旅游地的层次结构选择旅游地的层次结构 P3选择旅游地选择旅游地景景色色费费用用居居住住饮饮食食旅旅途途P1P2目标层 准则层 方案层 第39页，共97页，编辑于2022年，星期一空间信息分类与统计分析3.系统聚类分析法 聚类分析是根据地理变量的属性或特征的相似性、亲疏程度，用数学方法逐步地化型分类的一种方法。聚类分析应有如下几步操作：数据处理聚类分析统计量的选择和计算选择合适的聚类方法进行聚类分析第40页，共97页，编辑于2022年，星期一 一、空间聚类一、空间聚类空空间间聚聚类类是是根

29、根据据预预先先设设定定的的聚聚类类条条件件，使使符符合合条条件件的的区区域域输输出出在在图图上上，不不符符合合条条件件的的区区域域为为空空白白。在空间聚类中常用的是逻辑运算，用逻辑表达式来分析处理非几何特性之间的逻辑关系。常用逻辑运算包括：（1）逻辑交运算。若子集为A,B，对其进行逻辑交运算，得交集C=AB,如图6-1（a）所示。（2）逻辑并运算。若子集为A,B，对其进行逻辑并运算，得并集C=AB,如图6-1（b）所示。（3）逻辑非运算。若子集为A,B，对其进行逻辑非运算，得并集C=A-B,如图6-1（c）所示。（a）(b)(c)第41页，共97页，编辑于2022年，星期一 逻辑设定条件可以按

30、照需要的属性设定，也可按照空间要素的集合条件设定。如左下图所示，是一幅某镇土地利用现状栅格图（1.耕地，2.园地，3.林地，4.居民点，5.独立工矿，6.水域，7未利用地），设定条件可以是：E=（属性=“水域”）（面积1公顷）（水域邻接居民地）其输出的结果图如右下图所示。这类聚类条件的设定常用于位址规划。第42页，共97页，编辑于2022年，星期一 二、空间聚合二、空间聚合空空间间聚聚合合是是根根据据预预先先设设定定的的聚聚合合条条件件，在在同同一一图图层层上上进进行行数数据据类类别别的的合合并并或或转转换换，以以实实现现空空间间地地域域的的兼兼并并，从从而而将将复复杂杂的的空空间间数数据据合

31、合并并成成预预定定的的类类别别。空间聚合的结果往往是将复杂的类别转换为较简单的类别，大多数以小比例尺图形输出。当从大比例尺图形向小比例尺图形转换时，常使用这种方法。空间聚合分析第43页，共97页，编辑于2022年，星期一 三、区域提取三、区域提取区域提取是以某一区域为提取范围，提取另一图层中相应地区的属性数据，从而得到一新的目标图。第44页，共97页，编辑于2022年，星期一空间信息分类与统计分析4.判别分析法 判别分析是根据某一地理对象的各种特征指标来分辨或判断其类型归属问题的一种统计方法。它对于地理类型的划分和边界判断具有重要的理论和实践意义。包括：（1）两类判别；（2）多类判别；（3）逐

32、步判别。第45页，共97页，编辑于2022年，星期一空间信息分类与统计分析两类判别:设A、B两个地区，有m个地理要素，在这些地理要素中，A、B分别有nA、nB组地理特征值，依据费歇准则，建立一个判别函数式，即：判别函数应满足两类地区之间的均值差应尽可能地大，各类内部离差平方和应尽可能地小，即二者I达到最大，才能将A和B两个地区清楚地分开。第46页，共97页，编辑于2022年，星期一空间信息分类与统计分析二、空间统计分析 1.1.常规统计分析常规统计分析 常规统计分析主要完成对数据集合的均值、总和、方差、频数、峰度系数等参数的统计分析。第47页，共97页，编辑于2022年，星期一空间信息分类与统

33、计分析2.空间自相关分析 空间自相关分析是认识空间分布特征、选择适宜的空间尺度来完成空间分析的最常用方法。目前普遍使用空间自相关系数MoranI指数，其计算公式如下：N表示空间实体数目；xi表示空间实体的属性值；是xi的平均值；Wij1表示空间实体i与 j相邻，Wij 0表示空间实体i与j不相邻。第48页，共97页，编辑于2022年，星期一空间信息分类与统计分析3.回归分析回归分析（回归分析（regression analysis)是确定两种或两种以上变是确定两种或两种以上变数间相互依赖的定量关系的一种统计分析方法。运用十分广泛，数间相互依赖的定量关系的一种统计分析方法。运用十分广泛，回归分析

34、按照涉及的自变量的多少，可分为一元回归分析和多回归分析按照涉及的自变量的多少，可分为一元回归分析和多元回归分析；按照自变量和因变量之间的关系类型，可分为线元回归分析；按照自变量和因变量之间的关系类型，可分为线性回归分析和性回归分析和非线性回归非线性回归分析。如果在回归分析中，只包括一个分析。如果在回归分析中，只包括一个自变量和一个因变量，且二者的关系可用一条直线近似表示，这自变量和一个因变量，且二者的关系可用一条直线近似表示，这种回归分析称为一元线性回归分析。如果回归分析中包括两个或种回归分析称为一元线性回归分析。如果回归分析中包括两个或两个以上的自变量，且因变量和自变量之间是线性关系，则称为

35、两个以上的自变量，且因变量和自变量之间是线性关系，则称为多元线性回归多元线性回归分析。分析。第49页，共97页，编辑于2022年，星期一空间信息分类与统计分析一元线性回归分析步骤：（1）建立模型；（2）求解模型。（3）对回归模型的可信度进行检验。第50页，共97页，编辑于2022年，星期一空间信息分类与统计分析4.趋势分析 通过数学模型模拟地理特征的空间分布与时间过程，把地理要素时空分布的实测时间点之间的不足部分内插或预测出来。第51页，共97页，编辑于2022年，星期一一、缓冲区的概念和作用 邻近度描述了地理空间中两个地物距离相近的程度，其确定是空间分析的一个重要手段。缓冲区分析是解决邻近度

36、问题的空间分析工具之一。7.5 缓冲区分析第52页，共97页，编辑于2022年，星期一缓冲区分析缓冲区就是根据点、线、面地理实体，建立其周围一定宽度范围内的扩展距离图。缓冲区实际上是一个独立的多边形区域，它的形态和位置与原来因素有关。第53页，共97页，编辑于2022年，星期一缓冲区分析二、缓冲区的建立 点的缓冲区的生成比较简单，是以点实体为圆心，以测定的距离为半径绘圆，这个圆形区域即为缓冲区。线和面的缓冲区生成，实质上是求折线段的平行线。第54页，共97页，编辑于2022年，星期一缓冲区分析（1）折线曲率过大时的相应处理（a）(b)缓冲区转角的处理第55页，共97页，编辑于2022年，星期一

37、缓冲区分析（2）重叠缓冲区的处理第56页，共97页，编辑于2022年，星期一缓冲区分析（3）不同缓冲区宽度时的处理第57页，共97页，编辑于2022年，星期一缓冲区分析三、缓冲区分析应用实例本例以道路拓宽中的拆迁指标计算，说明如何利用缓冲区分析，计算一条道路拓宽过程中的拆迁指标。操作步骤如下：（1）明确目的和标准：目的是计算由于道路拓宽而需拆迁的建筑物的建筑面积和房产价值。道路拓宽标准：道路从原有的20m拓宽到60m；道路尽量保持直线；部分位于拆迁区内的10层以上的建筑不拆除。（2）数据准备：现状道路数字地图；区域内建筑物分布图及其相关信息。第58页，共97页，编辑于2022年，星期一缓冲区分

38、析（3）空间操作：以拓宽半径建立道路缓冲区；将缓冲区与建筑物层数据进行拓扑叠加，形成包括部分或全部位于道路拓宽区内建筑物信息的新图。（4）统计分析：根据标准对部分落入拆迁区的10层以上高楼，从选择组中去除，并调整局部拓宽道路边界；对需拆迁建筑物进行拆迁指标计算。（5）结果输出：以地图或表格形式将分析结果输出。第59页，共97页，编辑于2022年，星期一 叠置分析是GIS最常用的提取空间信息的手段之一。GIS的叠置分析是将有关主题层组成的数据层面，进行叠加产生一个新数据层面的操作，其结果综合了原来两层或多层要素所具有的属性。叠置分析不仅包含空间关系的比较，还包括属性关系的比较。7.6 叠置分析第

39、60页，共97页，编辑于2022年，星期一叠置分析 叠置分析的分类：（1）从叠置条件上分条件叠置无条件叠置（2）从数据结构上分矢量叠置分析栅格叠置分析以特定的逻辑、算术表达式为条件，对两组或两组以上的图件中相关要素进行叠置。也称全叠置，将同一地区、同一比例尺的两图层或多图层进行叠合，得到该地区多因素组成的新分区图。用来求解两层或两层以上矢量数据的某种集合，实现拓扑叠置，得到新的空间特性和属性关系。用来求解两层或两层以上栅格数据的某种集合，得到的是新的栅格属性。第61页，共97页，编辑于2022年，星期一叠置分析 一、矢量数据的叠置分析 叠置分析是指在统一空间参照系统条件下，把两层或两层以上的专

40、题要素图层进行叠置，以产生空间区域的多重属性特征，或建立地理要素之间的空间对应关系。根据叠置目的的不同可分为合成叠置与统计叠置两类。第62页，共97页，编辑于2022年，星期一叠置分析 合成叠置与统计叠置（1）合成叠置 （2）统计叠置第63页，共97页，编辑于2022年，星期一叠置分析 空间叠置分析根据叠置对象图形特征的不同，分为点与多边形的叠置、线与多边形的叠置、多边形与多边形的叠置三种类型。第64页，共97页，编辑于2022年，星期一叠置分析 1.点与多边形叠置点与多边形叠置点与多边形点与多边形叠置分析（据黄杏元等，2001年）第65页，共97页，编辑于2022年，星期一叠置分析 2.线与

41、多边形叠置线与多边形叠置 线与多边形的叠置是确定一图层上的弧段落在另一图层的哪个多边形内，以便为图层的每条弧段建立新的属性。计算过程通常是计算线与多边形的交点，只要相交，就产生一个结点，将原线打断成多条弧段，并将原线和多边形的属性信息一起赋给新弧段。第66页，共97页，编辑于2022年，星期一叠置分析 2.线与多边形叠置线与多边形叠置 线与多边形叠置分析（据黄杏元等，2001年）第67页，共97页，编辑于2022年，星期一叠置分析 3.多边形与多边形叠置分析多边形与多边形叠置分析 多边形与多边形叠置分析是将两个或两个以上多边形图层进行叠加产生一个新的多边形图层的操作，其结果将原来多边形要素分割

42、成新多边形，新多边形要素综合了原来所有叠置图层的属性，用于解决地理变量的多准则分析、区域多重属性的模拟分析、地理特征的动态变化分析，以及区域信息提取等。第68页，共97页，编辑于2022年，星期一叠置分析 多边形的叠置分析基本的步骤多边形的叠置分析基本的步骤（1）对原始多边形数据形成拓扑关系；（2）多层多边形数据的空间叠置，形成新的 层；（3）对新层中的多边形重新进行拓扑组建；（4）剔除多余的多边形，提取感兴趣的部 分。第69页，共97页，编辑于2022年，星期一叠置分析 二、栅格数据的叠置分析栅格数据的叠置是一个比较简单的过程，层间叠置可通过像元之间的各种运算来实现。第70页，共97页，编辑

43、于2022年，星期一叠置分析 栅格叠置分析示例土地利用土地利用地形坡度地形坡度人口密度人口密度交通便利度交通便利度得得 分分得得 分分得得 分分得得 分分方面方面1方面方面2方面方面3方面方面4自然保护区自然保护区组组 合合适宜性适宜性第71页，共97页，编辑于2022年，星期一一、泰森多边形的概念和生成 设平面上有 n 个互不重叠的离散数据点，则其中的任意一个离散数据点Pi都有一个临近范围Bi，在Bi中的任一个点同Pi点之间距离小于它同其他离散数据点之间距离。这里的Bi域是一个不规则多边形，该多边形称为泰森多边形。泰森多边形分析第72页，共97页，编辑于2022年，星期一泰森多边形分析泰森多

44、边形的生成 泰森多边形的生成是将所有相邻离散数据点连成三角形，作这些三角形各边的垂直平分线，于是每一个离散数据点周围的若干条垂置平分线相交组成一个多边形。第73页，共97页，编辑于2022年，星期一泰森多边形分析泰森多边形和泰森三角形泰森多边形泰森多边形泰森三角形泰森三角形第74页，共97页，编辑于2022年，星期一泰森多边形分析二、泰森多边形的特点（1）每个多边形内含有且仅含有一个离散数据点。（2）若区域内的任意一点 K(xk,yk)位于包含离散点 Pi(xi,yi)的多边形内，则 K(xk,yk)到 Pi(xi,yi)之间的距离总小于它到其他离散点 Pj(xj,yj)之间距离。即：(xk-

45、xi)2+(yk-yi)21/2(xk-xi)2+(yk-yi)21/2第75页，共97页，编辑于2022年，星期一泰森多边形分析二、泰森多边形的特点（3）若点M(xm,ym)位于包含离散点Pi(xi,yi)与含Pj(xj,yj)的两个多边形的公共边上，则M(xm,ym)到Pi(xi,yi)之间的距离等于它到离散点Pj(xj,yj)之间距离。即：(xm-xi)2+(ym-yi)21/2=(xm-xi)2+(ym-yi)21/2（4）泰森多边形的任意一个顶点周围存在三个离散点，将其连成三角形后其外接圆的圆心即为该顶点，该三角形称泰森三角形。第76页，共97页，编辑于2022年，星期一泰森多边形分

46、析三、泰森多边形分析的应用 泰森多边形的分析方法及构成的多边形和三角网在区域地学分析中，具有广泛的意义。利用泰森多边形可以确定一些商业中心、工厂或其他的经济活动点的影响范围。泰森多边形分析并不只限于经济地理方面的应用，还可用于近邻关系分析，区域专题现象的分级统计分析，绘制区域等值线图等。第77页，共97页，编辑于2022年，星期一7.7数字高程模型分析第78页，共97页，编辑于2022年，星期一 空间数据往往是根据用户要求获取的采样观测值，其分布往往是不规则的。一般来讲，在已存在观测点的区域内估计未观测点的特征值的过程称数据内插，在已存在观测点的区域范围之外估计未观测点的特征值的过程称推估。它

47、们是GIS数据处理常用的方法之一，广泛应用于等值线自动制图、数字高程模型的建立、不同区域界线现象的相关分析和比较研究等。数字高程模型 空间插值第79页，共97页，编辑于2022年，星期一空间插值 对于二维的更复杂的数据，需用二次或高次多项式。二次线性变化曲面数学模型：z=b0+b1x+b2x2 二次趋势面数学模型：z=b0+b1x+b2y+b3x2+b4xy+b5y2 趋势面是平滑函数，很难正好通过原始数据点，除非数据点少而且曲面的次数高才能使曲面正好通过原始数据点。第80页，共97页，编辑于2022年，星期一空间插值线性回归分析线性回归分析 高次多项式高次多项式 第81页，共97页，编辑于2

48、022年，星期一空间插值二、局部内插法 1.双线性多项式内插如果数据是按正方形格网节如果数据是按正方形格网节点布置如左图所示，这时可点布置如左图所示，这时可先用点先用点A A和和B B及及C C和和D D两对点的两对点的高程，线性内插出点高程，线性内插出点MM和和NN的的高程（高程（ZMZM和和ZNZN），然后再由），然后再由ZMZM和和ZNZN直线内插待求点直线内插待求点P P的的高程高程ZPPZPP。内插公式为：。内插公式为：第82页，共97页，编辑于2022年，星期一空间插值2.双三次多项式（样条函数）在分块插值区内，每一分块可以定义出一个不同的多项式曲面，当 n 次多项式与其相邻分块的

49、边界上所有 n-1次导数都连续时，称为样条函数。在数据点为方格网的情况下，采用三次曲面来描述格网内的内插值时，则待定点内插值ZP为：ZP=a1x3y3+a2x2y3+a3xy3+a4y3+a5x3y2+a6x2y2+a7xy2+a8y2+a9x3y+a10 x2y+a11xy+a12y+a13x3+a14x2+a15x+a16。第83页，共97页，编辑于2022年，星期一空间插值3.移动拟合法 它是待定点为中心的典型的逐点内插法。其方法是：定义一个合适的局部函数去拟合周围的数据点，通过求解拟合函数，求出待定点的内插值。通常做法是将坐标点放置在待定点上，而采用的数据点应落在半径为R的圆内。第84

50、页，共97页，编辑于2022年，星期一空间插值移动拟合方法 第85页，共97页，编辑于2022年，星期一一、网络的概念和构成网络是一个由点和线的二元关系构成的系统，通常用来描述某种资源或物质在空间上的运动。网络分析的理论基础是图论，所用数据结构为非线性图数据结构。网络分析是运筹学模型中的一个基本模型，它的根本目的是研究、筹划一项网络工程如何安排，并使其运行效果最好，如一定资源的最佳分配，从一地到另一地的运输费用最低等。7.8 网络分析第86页，共97页，编辑于2022年，星期一第八节 网络分析 网络数据模型是真实世界中网络系统的抽象表示。网络是有若干线性实体互连而成的一个系统，资源经由网络来传