GIS7第八章空间分析

第八章空间分析课题

第八章空间分析目的要求

通过本章的学习，了解GIS中基本空间分析的功能及其实现算法。理解叠置分析的概念和类型，掌握多边形叠置分析的步骤和方法。理解缓冲区的概念和作用，并能说明其应用方法，理解生成线缓冲区的算法。教学重点

空间查询与量算、缓冲区分析、叠置分析、网络分析、空间插值、分类统计教学难点

网络分析、空间插值、分类统计教学课时

8课时教学方法

课堂实例+讲授空间分析的目的空间分析:指用于分析空间目标的一系列技术处理,其目的是：(1)描述与认知空间数据分布特征，如点线面的空间分异状况；(2)解释空间现象与空间模式的形成机理，如城市土地利用变化研究；(3)调控在地理空间上发生的事件，如水资源的合理配置；(4)预测预报，如洪水的预测预报。空间分析的主要方法(1)基于地图的空间图形分析，如GIS中的缓冲区、叠加分析、数字高程模型,数字地面模型等；(2)

空间动力学分析，有城市扩张模型(驱动力等)、空间价格竞争模型(区位优势)、空间择位模型（中心地等）；(3)空间信息分析，是指根据数据或统计方法建立的模型，如空间聚类、空间自相关、回归模型等。实际上，同一种空间分析方法和模型可以归属于不同的类型。例如，中心地属于空间信息分析模型，但中心地的形成又追求服务距离最短的动力学机制，可归属于空间动力学模型。1．空间信息查询与量算查询和定位空间对象，并对空间对象进行量算是GIS的基本功能之一，它是GIS进行高层次分析的基础。在GIS中，为进行高层次分析，往往需要查询定位空间对象，并用一些简单的量测值对地理分布或现象进行描述，如长度、面积、距离等。实际上，空间分析首先始于空间查询和量算，它是空间分析的定量基础。（1）空间信息查询图形和属性的互查是最常用的查询，主要有两类：1、按属性信息的要求来查询定位空间位置，称为“属性查图形”。如在中国行政区划图上查询人口大于4000万且城市人口大于1000万的省有哪些?称为SQL查询.2、根据对象的空间位置查询有关的属性信息，称为“图形查属性”。如一般的GIS软件都提供一个“INFO”工具，让用户利用鼠标，用点选、画线、矩形、圆、不规则多边形等工具选中地物，并显示所查询对象的属性列表，可进行有关统计分析。（1）空间信息查询1、基于空间关系查询

空间实体间存在多种空间关系，包括拓扑、距离、方位等。如查找满足下列条件的城市： 在京沪线的东部；距离京沪线不超过50公里;

城市人口大于100万;城市区域面积5000平方公里.

简单的点线面相互关系拓扑查询包括：面面查询:如与某个多边形相邻的多边形有哪些；面线查询:如某个多边形内包含哪些线；面点查询:如某个多边形内有哪些点状地物；线面查询:如某条线经过的多边形有哪些；线线查询:如与某条河流相连的支流有哪些；线点查询:如某条道路上有哪些桥梁，某条输电线上有哪些变电站点面查询:如某个点落在那个多边形内；点线查询:如某个结点由哪些线相交而成；（1）空间信息查询

2、基于空间关系和属性特征查询

传统的SQL并不能处理空间查询，对GIS而言，需要对SQL进行扩展,主要包括空间数据与属性数据的匹配等3、地址匹配查询根据街道的地址来查询事物的空间位置和属性信息是GIS特有的一种查询功能，这种查询利用地理编码，输入街道的门牌号，就可以知道大致的位置和所在的街区。（1）空间信息查询几何量算

1．长度

线状物体的长度是最基本的形态参数之一，在矢量数据格式下，线由点组成，线状物体表示为一个坐标串(Xi,Yi)，而线长度可由两点间直线距离相加得到。则线状物体长度的计算公式为：

（2）、空间信息量算-12．面积

多边形的面积是一个重要指标。多边形边界可以分解为上下两半，其面积就是上半边界下的积分值与下半边界下的积分值之差。设面状物体的轮廓边界由一个点的序列P1(x1,y1),P2(x2,y2),…，Pn

(xn,yn)表示，其面积为:几何量算

（2）、空间信息量算-12．面积YoXSS1S2S=S2-S1å=++=niiiiiyxyxS11121几何量算（2）、空间信息量算-1YoXS1(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（2）、空间信息量算-1形状量算

面状地物形状量测包括空间一致性问题，即有孔多边形和破碎多边形的处理；多边形边界特征描述问题当把城市作为单个面状目标看待时，可以直接使用面状目标的形状系数，如形状率、圆形率、紧凑度等,这些指标计算较简单，但只反映一个抽象的形状；当把城市作为面状目标的集合看待时，可以使用放射状指数、标准面积指数等形状系数，这些指标计算较复杂，但反映了城市内部的具体联系。在多数指标中，都以圆形作为城市的标准形状。

（2）、空间信息量算-1欧拉数=（孔数）-（碎片数-1）欧拉数空间一致性最常用的指标是欧拉函数，用来计算多边形的破碎程度和孔的数目。（2）、空间信息量算-1

1）形状比(FORMRATIO)

该指标能反映城市的带状特征，城市的带状特征越明显则形状比越小。显然，如果城市为狭长带状分布，其长轴两端的联系是不便捷的。形状比=A/L2

其中，A为区域面积，L为区域最长轴的长度。（2）、空间信息量算-1

伸延率=L/L’

式中，L为区域最长轴长度，L’为区域最短轴长度。

2）伸延率(ELONGATIONRATIO)

该指标反映城市的带状延伸程度，带状延伸越明显则延伸率越大，反映城市的离散程度越大。（2）、空间信息量算-1

公式：紧凑度=P

其中，A为面积，P为周长。该指标反映城市的紧凑程度，其中圆形区域被认为最紧凑，紧凑度为1。其它形状的区域，其离散程度越大则紧凑度越低。3）紧凑度(COMPACTNESSRATIO)A/

2p（2）、空间信息量算-1质心量算质心是描述地理对象空间分布的一个重要指标。例如要得到一个全国的人口分布等值线图，而人口数据只能到县级，所以必须在每个县域里定义一个点作为质心，代表该县的数值，然后进行插值计算全国人口等值线质心通常定义为一个多边形或面的几何中心，当多边形比较简单，比如矩形，计算很容易。但当多边形形状复杂时，计算很复杂（2）、空间信息量算-2（2）、空间信息量算-2分两种情况：

1）面状目标的重心。可以理解为多边形内的平衡点，正如一块均质木块被悬挂起来的平衡点。（2）、空间信息量算-2面状目标重心可以通过计算梯形重心的平均值而得到。将多边形的各个顶点投影到x轴上，就得到一系列梯形（如图），所有梯形重心的联合就确定了整个多边形的重心。

按梯形计算重心位置

按梯形计算重心位置设多边形的顶点序列(xi,yi)按顺时针编码，则其重心的计算公式为：其中，和是第i个梯形的重心的x坐标和y坐标，是梯形的面积。它们由下式得到：åå=iiiGAAXX/∑∑/=iiiGAAYYiY

按梯形计算重心位置ïîïíì-++=-++=-+=++++++++6/))((6/))((2/))((121211212111iiiiiiiiiiiiiiiiiiiiixxyyyyAYyyxxxxAXxxyyA

可理解为其分布中心。其重心计算方法是取离散目标的加权平均中心，它是离散目标保持均匀分布的平衡点。计算公式为：2）面状分布离散目标的重心

其中，i为离散目标物，Wi为该目标物权重。Xi与Yi为其坐标。∑∑∑∑=,=iiiiiGiiiiiGWYWYWXWX（2）、空间信息量算-3（2）、空间信息量算-3距离量算“距离”是人们日常生活中经常涉及到的概念，它描述了两个事物或实体之间的远近程度。最常用的距离概念是欧氏距离，无论是矢量结构，还是栅格结构都很容易实现在GIS中，距离通常是两个地点之间的计算，但有时人们想知道一个地点到所有其它地点的距离，这时得到的距离是一个距离表面

（2）、空间信息量算-3各向同性和各向异性的距离表面欧氏距离、曼哈顿距离和一种非欧氏距离2．空间变换

为了满足特定空间分析的需要，需对原始图层及其属性进行一系列的逻辑或代数运算，以产生新的具有特殊意义的地理图层及其属性，这个过程称为空间变换。矢量数据结构空间变换复杂繁琐，栅格数据结构空间变换容易。基于栅格结构的空间变换可分为三种方式：单点变换：假定运算不受其他邻点属性影响。邻域变换：考虑邻域单元的影响，平滑，坡度，坡向分析等。区域变换：考虑整个变换区域的属性影响，如整体插值，求和，归组3．再分类

再分类（Reclassification）：根据不同的需要对原始数据再次进行分类和提取的过程。3．再分类基于非空间属性分类可用经典的数理统计算法，如主成分分析、层次分析、聚类分析、判别分析，不改变已有属性，仅根据属性，划分到相应的类别中。矢量数据结构：点、线地物直接修改属性表中的数值实现。面状地物属性修改需同时修改几何形状和属性（去公共边界，属性统一）。栅格数据结构：点线面均通过修改属性值并改变图例表示。归组：最常用和最简单的再分类。ArcGIS中3D分析和空间分析工具条中都有再分类工具4．缓冲区分析

1、概念邻近度：描述了地理空间中两个地物距离相近的程度，其确定是空间分析的一个重要手段。交通沿线或河流沿线的地物有其独特的重要性，公共设施的服务半径，大型水库建设引起的搬迁，铁路、公路以及航运河道对其所穿过区域经济的发展的重要性等，均是一个邻近度问题。缓冲区分析是解决邻近度问题的空间分析工具之一。4．缓冲区分析缓冲区：地理空间目标的一种影响范围或服务范围。缓冲区分析：是指根据分析对象的点、线、面实体，自动建立其周围一定距离的带状区，用以识别这些实体或者主体对邻近对象的辐射范围或者影响程度，是解决临近度问题的空间分析工具之一。它在交通、林业、资源管理、城市规划中有着广泛的应用。例如：湖泊和河流周围的保护区的定界；汽车服务区的选择；民宅区远离街道网络的缓冲区的建立等。4．缓冲区分析2缓冲区主要的类型（1）基于点要素的缓冲区：通常以点为圆心、以一定距离为半径的圆（2）基于线要素的缓冲区：通常是以线为中心轴线，距中心轴线一定距离的平行条带多边形。（3）基于面要素的缓冲区：向外或向内扩展一定距离以生成新的多边形。点、线、多边形的缓冲区另外，还有特殊形态的缓冲区，如点对象有三角形、矩形、圆形；线对象有双侧对称，双侧不对称或是单侧缓冲区，对于面对象有内侧和外侧缓冲区。第七节缓冲区分析

不同缓冲区宽度时的处理图示4．缓冲区分析3空间缓冲区分析过程(1)建立缓冲区以图形元素为基础，拓宽或紧缩一定宽度而形成的区域。这个宽度通常是等距的，也可以是不等距的缓冲区。

(2)缓冲区分析根据建立的缓冲区，对缓冲区内的空间信息形态、特征、分布作进一步分析。栅格缓冲区的建立将栅格数据表示为一个二值(0,1)矩阵(M×N)，其中“0”像元为空白位置，“1”元素为空间物体所占据的位置。经过距离变换，计算出每个“0”元素与最近的“1”元素的距离，即背景像元与空间物体的最小距离。假设缓冲区的宽度为d，则缓冲区边界就是距离为d的各个背景像元的集合。4．缓冲区分析矢量缓冲区的建立矢量缓冲区常见的有角平分法。角平分法由三步组成，即逐个线段计算简单平行线，尖角光滑矫正和自相交处理。尖角光滑矫正除角平分线法之外，还可采取圆弧法，但矫正过程都很复杂，难以完备地实现。4．缓冲区分析难以保证双线的等宽性角分线法矢量缓冲区计算的基本问题是双线问题。凸角圆弧法BdRd=r/sin(B/2)保证了平行曲线的等宽性折点凸凹性的判断把相邻两个线段看成两个向量，其方向取坐标点序方向。若前一个向量以最小角度转向第二个向量时呈逆时钟方向，则为凸顶点。反之为凹顶点折点凸凹性的判断baABCAB=(XB-XA,YB-YA)BC=(XC-XB,YC-YB)S=AB×

BC=(XB-XA)(YC-YB)-(XC-XB)(YB-YA)S>0,ABC呈逆时钟，顶点为凸。S<0,ABC呈顺时钟，顶点为凹。S=0,ABC共线。缓冲取边线的自相交复杂对象或对象集合的缓冲区边线可能自相交。右图就是一个例子。当轴线的弯曲空间不允许双线的边线无压的通过时，就会产生若干个自相交多边形。重叠多边形与岛屿多边形岛屿多边形是缓冲区边线的有效组成部分；重叠多边形不是缓冲区边线的有效组成部分，不参与缓冲区的最终构建。岛屿重叠多边形岛屿重叠多边形与岛屿多边形的自动判别首先定义轴线坐标点序为其方向，缓冲区双线分成左右边线，左右边线自相交多边形的判别情形恰好相反。对于左边线，岛屿自相交多边形呈逆时钟方向。重叠自相交多边形呈顺时钟方向。重叠缓冲区的处理图示4．缓冲区分析4、缓冲区分析的三要素在进行空间缓冲区分析时，通常要将研究的问题抽象为以下三类要素:①主体

表示分析的主要目标，一般分为点源、线源和面源三种类。②邻近对象表示受主体影响的客体，例如行政界线变更时所涉及的居民区、森林遭砍伐时所影响的水土流失范围等。③对象的作用条件表示主体对邻近对象施加作用的影响条件或强度。4．缓冲区分析解题流程解题过程

首先要以区域的道路分布图、河流分布图、森林分布图为数据源。解题流程见图所示。道路分布图森林分布图河流分布图结束生成道路周围5km缓冲区叠置生成河流周围1km缓冲区叠置5空间缓冲区分析实例缓冲区分析实例1

已知一伐木公司，获准在某林区采伐，为防止水土流失，规定不得在河流周围1km内采伐林木。另外，为便于运输，决定将采伐区定在道路周围5km之内。请找出符合上述条件的采伐区，输出森林采伐图。4．缓冲区分析将该地区具有相同比例尺且进行配准的道路分布图、河流分布图、森林分布图，进行预处理和数字化；利用河流分布图生成1km的等距离缓冲区；利用道路分布图生成5km的等距离缓冲区；森林分布图中可采伐林地、道路缓冲区及河流缓冲区图进行叠置，叠置条件表达式为:

采伐区＝森林分布图中可伐林地∩道路周围5km缓冲区

∩非河流周围1km缓冲区将上述3张图进行两两叠置，所得结果即为森林采伐图。4．缓冲区分析缓冲区分析实例2如已知一湖泊，要求在它周围5000m内必需禁止任何污染性工业企业存在，在它周围500m内必需禁止建筑任何永久性建筑物。解题步骤：（1）先建立缓冲区；（2）同现有污染性工业企业图叠置，显示在范围内应禁止的污染性工业企业；（3）同现有永久性建筑物图叠置，显示在范围内应禁止的永久性建筑物。4．缓冲区分析缓冲区分析实例4如已知流域水系图，要求计算在它周围5公里内所有建筑物,以便进行分析。解题步骤：（1）先建立缓冲区；（2）同现有污染性工业企业图叠置，显示在范围内应禁止的污染性工业企业；（3）同现有永久性建筑物图叠置，显示在范围内应禁止的永久性建筑物。4．缓冲区分析ARC/INFO命令:Buffer<in_cover><out_cover>{buffer_item}{buffer_table}{buffer_distance}{fuzzy_tolerance}{LINE|POLY|POINT|NODE}{ROUND|FLAT}{FULL|LEFT|RIGHT}对于本例:Bufferbufferbufferline##50001lineflat5．叠加分析大部分GIS软件以分层的方式组织地理景观，将地理景观按主题分层提取，同一地区的整个数据层集表达了该地区地理景观的内容。每个主题层可以用矢量结构的点线面图层文件表达，也可以用栅格结构的图层文件进行表达。叠加分析是GIS最常用的提取空间隐含信息的手段之一。该方法来源于传统的透明材料叠加。既将来自不同数据源的图纸绘于透明纸上，在透光桌上将其叠放在一起，然后用笔勾出感兴趣的部分。5．叠加分析叠加分析：将有关主题层组成的数据层面，进行叠加产生一个新数据层面的操作，其结果综合了原来两层或多层要素所具有的属性。也就是说，叠加分析不仅生成了新的空间关系，还将输入数据层的属性联系起来产生了新的属性关系。叠加分析是对新要素的属性按一定的数学模型进行计算分析，进而产生用户需要的结果或回答用户提出的问题。5．叠加分析分为以下三类视觉信息叠加矢量数据叠加点与多边形叠加线与多边形叠加多边形叠加栅格图层叠加5．叠加分析一、视觉信息叠加视觉信息叠加是将不同侧面的信息内容叠加显示在结果图件或屏幕上，以便研究者判断其相互空间关系，获得更为丰富的空间信息。地理信息系统中视觉信息叠加包括以下几类：点状图，线状图和面状图之间的叠加显示面状图区域边界之间或一个面状图与其他专题区域边界之间的叠加遥感影象与专题地图的叠加专题地图与数字高程模型（DEM）叠加显示立体专题图视觉信息叠加不产生新的数据层面，只是将多层信息复合显示，便于分析5．叠加分析二、基于矢量数据的叠置分析基于矢量数据的叠置分析是参与分析的两个图层的要素均为矢量数据。点面叠置分析

线面叠置分析面面叠置分析5．叠加分析矢量数据——点与多边形叠加

点与多边形的叠加分析，实质是计算包含关系（包含分析），判断各个点的归属（落在哪个多边形内）。叠置的结果是为每点产生一个新的属性。例如，井位与规划区叠加，可找到包含每个井的区域。5．叠加分析

点与多边形的叠加通常不产生新的数据层，只是把属性信息叠加到原图层中，然后通过属性查询间接获得点与多边形叠加的需要信息。

在完成点与多边形几何关系计算后，还要进行属性信息处理。最简单的方式是将多边形属性信息叠加到其中的点上（或反之）。如有多点分布在同一多边形，则要采取一些特殊的规则。矢量数据——点与多边形叠加矢量数据——点与多边形叠加例如：将水井与规划区图层相叠置，可确定每口井所属的规划区范围。*1*3*4*2ADBC点属性1234多边形属性1属性2ABCD点多边形点属性面属性1面属性21A2D3C4B5．叠加分析线与多边形的叠置是将线状要素层或网状要素层和多边形叠置，对线和多边形求交运算，根据每个线要素同多边形的关系，以形成新的空间目标集、新的属性表，得到线与多边性联合的属性表。线与多边形的叠合的目的是确定某一线状图层上的弧段落在另一多边形图层上的哪个多边形内，以便为图层的每条弧段建立新的属性。矢量数据——线与多边形叠加5．叠加分析例如：当确定某一行政区内各种等级道路的里程数时，就需要将道路图与境界图相叠置，计算弧段与多边形边界的交点，在交点处截断弧段，并对弧段重新编号，建立弧段与多边形的归属关系。线与多边形叠置分析多边形叠加分析也称为Polygon-on-polygon叠置，它是指同一地区、同一比例尺的两组或两组以上的多边形要素的数据文件进行叠置。参加叠置分析的两个图层应都是矢量数据结构。若需进行多层叠置，也是两两叠置后再与第三层叠置，依次类推。其中被叠置的多边形为本底多边形，用来叠置的多边形为上覆多边形，叠置后产生具有多重属性的新多边形。矢量数据——多边形与多边形的叠加分析5．叠加分析矢量数据的多边形叠置分析5．叠加分析

多边形与多边形的叠置是指将两个不同图层的多边形要素相叠加，根据两组多边形的交点来建立多重属性的多边形或进行多边形范围内的属性特征的统计分析。

多边形叠置过程分几何求交过程和属性确定过程，算法的核心是多边形求交。①对两个多边形进行边界求交和弧段分割运算，并以新弧段为单位重建拓扑关系；②判断重建多边形落在原始多边形层的哪个多边形内，从而建立新叠置多边形与原始多边形的关系，并抽取属性。多边形与多边形的叠置原理多边形叠加分析矢量数据的多边形叠置分析5．叠加分析取本底多边形取上叠多边形两个多边形相交？两个多边形是包含关系？多边形求交，重建拓扑关系属性抽取，建立关联连接上叠多边形结束？本底多边形结束？ynnn建立包含关系属性抽取yyy结束多边形的叠置流程图矢量数据的多边形叠置分析多边形叠置就是地理数据层的叠加，分为合成叠置与统计叠置。合成叠置是将同一地区、同一比例尺的两组或更多的多边形要素的数据文件进行叠置，根据两组多边形边界的交点来建立具有多重属性的多边形，如图(a)。统计叠置是将多边形数据层叠加，进行多边形范围的属性特征的统计分析，如图(b)。5．叠加分析合成叠置统计叠置(a)合成叠置

(b)统计叠置

多边形叠加产生碎屑多边形多边形叠加产生碎屑多边形并(Union)：保留两个输入层中所有多边形；交(Intersect)：保留两个输入层中的公共区域判断(Identity)：输出以其中一输入图层为界的边界内的所有多边形多边形与多边形的叠置分析功能删除(Erase)：输出层为保留以其中一输入图层为控制边界之外的所有多边形。即：在将更新的特征加入之前，须将控制边界之内的内容删除。更新(Update)：输出层为一个经删除处理后的图层与一个新特征图层进行合并后的结果

切割(Clip)：输出层为按一个图层的边界，对另一个图层的内容要素进行截取后的结果。

参考：黄杏元.地理信息系统概论（修订版），第五章第2节，空间叠合分析相关内容。逻辑判断叠置法ABCABCABCABCABCABCA.AND.B.AND.CA.NOT.(B.AND.C)A.AND.B.OR.CA.OR.B.OR.CA.XOR.B.XOR.CA.AND.(B.OR.C)5.空间叠加日喀则地区沙漠化类型分布图日喀则地区行政区界线图叠置分析实例：有日喀则地区的行政界线图层和通过遥感技术提取的该区沙漠化类型分布图层，求日喀则地区各县沙漠化类型统计数据。图层1图层25.空间叠加行政区界图属性表图层1：图形，属性表日喀则地区行政区界线图5.空间叠加日喀则地区沙漠化类型分布图属性表图层2：图形，属性表日喀则地区沙漠化类型分布图5.空间叠加空间叠置图叠置图属性表使用IDENTITY命令，进行两图层空间叠置，得到叠置图5.空间叠加叠置图层属性表图层1属性表图层2属性表5.空间叠加数学运算叠置法：算术运算、函数运算算术运算指两层以上的对应网格值经加、减运算，而得到新的栅格数据系统的方法。三、基于栅格数据的叠置分析11111111111111111111111111111111111131111111ABCD=A+B+CE=|A-B|F=D-E111111111112222315.空间叠加三、基于栅格数据的叠置分析

特点：栅格数据的叠置算法，虽然数据存贮量比较大，但运算过程比较简单。变换方法：（1）点变换（2）区域变换方法（3）邻域变换方法基于栅格数据的叠置分析是参与分析的两个图层的要素均为栅格数据。5.空间叠加（1）点变换点变换只依据参与叠置图层相应点的属性值进行新的运算，既与各图层的邻域点的属性无关，也不受区域内一般特征的影响。运算方法包括：算术运算，指数运算，三角函数运算等5.空间叠加244633463377225423643366-=0-21-2-102-200-1-1点变换示意结果特征：运算后得到的新属性值可能与原图层的属性意义完全不同。5.空间叠加(2)区域变换方法指在计算新图层相应的属性值时，不仅与原图层对应的栅格的属性值有关，而且要顾及原图层所在区域的集合特征（区域长度、面积、周长等）。(3)邻域变换方法指在计算新图层相应的属性值时，不仅考虑原图层对应的栅格及其属性，而且还应顾及与该栅格相关联的邻域或者影响半径内的栅格属性值的影响。5.空间叠加101060601010206030303060101060601010306030303060-=0000001000000点变换实例1：土地利用变化区域探测80年遥感影像90年遥感影像点变换后影像点变换后影像分析通过80和90年两期影像的相减运算后得到变换影像，如果：变换影像值＝0；说明该区未发生变化变换影像值≠0；说明该区已发生变化10耕地20居民点30水域40草地50未利用地60林地Legend注意：此处的遥感影像可以是分类结果，也可以是原始的遥感影像。在一般应用中，多使用原始的遥感影像，可提高变化探测速度。5.空间叠加信息获取研究方法变化过程分析机制分析点变换实例2：土地利用变化分析更新土地利用图社经统计数据动态变化研究变化信息的解译与制图变化数据管理原始遥感信息专题图件数据预处理空间操作与分析图像预处理影像分类变化信息探测结果比较变化信息分类高分辨率图像融合土地利用/覆盖动态变化研究框架5.空间叠加1985年(TM)1995年(TM)1998年(SPOT)2000年(TM)

不同时期土地利用图信息获取5.空间叠加土地类型转移矩阵:根据两个不同时间（t和t+△t）的土地利用图计算从一种类型到另一种类型的转换概率，来分析土地利用变化过程。它依靠GIS技术，将两个不同时间的土地利用图进行删格化处理，计算t时刻上A类有多少格网点转换成在t+△t时刻的B，C，D等类型，转换点数占该类型总数的百分比可称为转移概率。研究方法使用景观生态学方法，计算土地类型转移矩阵。5.空间叠加对任意两期土地利用类型图和，按照下式的地图代数方法，可以求得：由k时期到k+1

时期的土地利用变化图Ci×j，它表现了土地利用变化的类型及其空间分布。

5.空间叠加

12┄nai+1a11a12┄a1na1+2a21a22┄a2na2+┆┆┆┆

Nan1a32┄annan+a+ja+1a+2┄a+na注：表中ai+代表土地类型I在转化前的总量；a+j代表土地类型I在转化后的总量；a是土地总面积

类型转移距阵示例

5.空间叠加变化过程分析5.空间叠加5.空间叠加土地利用地形坡度人口密度交通便利度得分得分得分得分方面1方面2方面3方面4自然保护区组合适宜性栅格叠置分析示例 6.网络分析网络数据结构的基本组成部分和属性1、链（Link）网络中流动的管线如街道、河流、水管，其状态属性包括阻力和需求。2、结点（Node）网络中链的结点，如港口、车站等，其状态属性包括阻力和需求等。结点中的特殊类型障碍（Barrier），禁止网络上流动的点。拐点（Turn），出现在网络中的分割点上，其状态有属性和阻力，如拐弯的时间和限制（如在8点到18点不允许左拐）。中心（Center），是接受或分配资源的位置，如水库、商业中心，电站等，其状态包括资源容量（如总量），阻力限额（中心到链的最大距离或时间）。站点（Stop），在路径选择中资源增减的结点，如库房、车站等，其状态属性有资源需求，如产品数量。

路径分析静态求最佳路径：在给定每条链上的属性后，求最佳路径。N条最佳路径分析：确定起或终点，求代价最小的N条路径，因为在实际中最佳路径的选择只是理想情况，由于种种要素而要选择近似最佳路径。最短路径或最佳耗费路径：确定起点终点和要经过的中间点、中间连线，求最短路径或最佳耗费路径。动态最佳路径分析：实际网络中权值是随权值关系式变化的，可能还会临时出现一些障碍点，需要动态的计算最佳路径。

最小值问题求解如：有数组intArray[n]（Array[i]<1000,i=0,n）,如何找出它们之中的最小值？int

FindMin(int*array,intbound){

intmin=1000;

for(inti=0;i<bound;i++) if(array[i]<=min) min=array[i]; returnmin;}……

intarray[7]={789,33,7898,7565,76,22,88};

intresult=FindMin(array,7);ASSERT(result==22)……最佳路径求解最佳路径求解有多种不同的方法，其中Dijkstra算法适合于求解某个起点（源点）到网络中的其它各个结点的最佳路径。Dijkstra算法（1）1、引进一个辅助向量D，它的每个分量D[i]表示当前所找到的从起点vm

到每个终点vi的最短路径的长度。假设用带权的邻接矩阵arcs来表示带权有向图，arcs[i][j]表示弧<vi,vj

>上的权值。若<vi,vj

>不连通，则arcs[i][j]=∞。那么D[i]的初值为：D[i]=arcs[m][i]vi∈V此外，将已找到的从vm

出发的最短路径的终点的集合记为S，它的初始状态为空集。Dijkstra算法（2）2、选择vj

使得D[j]=Min{D[i]|vi∈V-S}Vj就是当前求得的一条从vm出发的最短路径的终点。其中j为这条最短路径的终点，将其加入到终点集合S，令

S=S∪{j}Dijkstra算法（3）3、修改辅助向量D，即修改从vm出发到集合V-S上任一顶点vk可达的最短路径长度。显然，若D[j]+arcs[j][k]<D[k]，则表明从vm出发，经过vj到达vk的路径更短。因此，如果D[j]+arcs[j][k]<D[k]，则修改D[k]为：

D[k]=D[j]+arcs[j][k]VmVjVkarcs[j][k]=8D[k]=15D[j]=5V-SSDijkstra算法（4）4、重复操作第二步、第三步共n-1次。由此求得从v到图上其余各顶点的最短路径是依路径长度递增的序列。

例子20V5V0V4V1V3V210060301010505带权有向图邻接矩阵例子(思路)ACiBi 如图所示，A为所求最短距离的起点，其他Bi,Ci

为终点。 我们要求的是一系列最短距离。我们先假定这些最短距离互不相等。那么我们可以把这些最短距离按升序（从小到大）排列。 我们把所有顶点分为两类C和B.令A到Bi这些顶点的最短距离不为无穷大。A到Ci这些顶点的最短距离为无穷大。

这就说明A到Ci中的点要么不通，要么通过Bi中的点与之连接。

例子(思路)ACiBi

这样，对于A到Ci中任何一个点的最小距离，我们总可以在Bi中找到一点，使得A到这一点的最小距离小于前一个距离。（因为：A到Ci中的点要么不通，要么通过Bi中的点与之连通。）

因此，我们可以先不考虑Ci中的点。例子(思路) 于是，对于右图，我们第一步只考虑下图：V5V0V4V21003010Bi={v2,v4,v5}V5V0V4V1V3V21006030101050520例子(思路)我们用mindist[]这个数组来保存由v0到其它顶点的最小距离，这些距离按升序排列。考虑右图：第一步，通过比较，我们知道mindistance[v0][v2]=mindist[0]=10,（v0-v2)，这是我们求出的第一个最小距离。一旦我们得到v2，我们就可以知道：V5V0V4V21003010V5V0V4V1V3V21006030101050520例子(思路)V0跟v2直接连通到的点v3

之间的最小距离不再是无穷大，它应当是mindistance[v0][v2]+dis[v2][v3],这样我们应当把v3放入前面的集合Bi中。（注意：有多少v2直接连通到的点都应当考虑进来。）V5V0V4V3V2100301050Bi={v2,v4,v5,v3}V5V0V4V1V3V21006030101050520例子(思路)第二步，我们把与v2直接连通的点v3考虑进来。dis[0][5]=100;dis[0][4]=30;dis[0][2]=10;dis[0][3]=60;除10以外，30是最小的。我们可以证明30是v0到其它顶点除10以外最小的。V5V0V4V3V2100301050V5V0V4V1V3V21006030101050520例子(思路)不可能存在这样一个点Vn，使得10<mindistance[0][n]<30.原因如前所述。VnV5V0V4V3V2100301050BiV5V0V4V1V3V21006030101050520例子(思路)不可能存在这样一个点Vn，使得10<mindistance[0][n]<30.原因如前所述。VnV5V0V4V3V2100301050BiV5V0V4V1V3V21006030101050520例子(思路)这样我们得到我们的第二个最小距离：Mindistance[v0][v4]=mindist[1]=30,(v0-v4)接下来，我们把v4与之直接连通的点考虑进来。。。V5V0V4V3V2100301050BiV5V0V4V3V2100301050V5V0V4V1V3V21006030101050520例子 以v0为起点，计算它到其它各顶点的最短路径，计算过程中最短路径长度向量D的变化见D0-D4，计算出的各条最短路径见表4-4。úúúúúúûùêêêêêêë饥=10030100D例子终点从v0到其它各结点的最短路径v1∞∞∞∞∞v210(v0,v2)v3∞60(v0,v2,v3)50(v0,v4,v3)v430(v0,v4)30(v0,v4)v5100(v0,v5)100(v0,v5)90(v0,v4,v5)60(v0,v4,v3,v5)vjv2v4v3v5V5V0V4V1V3V21006030101050520例子起点终点最短路径路径长度v0v1无

v2(v0,v2)10

v3(v0,v4,v3)50

v4(v0,v4)30

v5(v0,v4,v3,v5)60V5V0V4V1V3V21006030101050520求最短路径的方法第六节空间网络分析（2）地址匹配实质是对地理位置的查询，它涉及到地址的编码(Geocode)。地址匹配与其他网络分析功能结合起来，可以满足实际工作中非常复杂的分析要求。所需输入的数据，包括地址表和含地址范围的街道网络及待查询地址的属性值。2、地址匹配第六节空间网络分析（3）3、资源分配资源分配网络模型由中心点(分配中心)及其状态属性和网络组成。资源分配的两种方式：（1）由分配中心向四周输出这种分配功能可以解决资源的有效流动和合理分配，其在地理网络中的应用与区位论中的中心地理论类似。（2）由四周向中心集中涉及资源分配的中心选址问题中心点选址问题中，最佳选址位置的判定标准，是使其所在的顶点与图中其它顶点之间的最大距离达到最小。这个选址问题实际上就是求网络图的中心点问题。这类选址问题适宜于医院、消防站、超市、学校等服务设施的布局问题。

中心选址问题的图论描述

设G=(V,E)是一个无向赋权连通图，其中V={v1,v2,…,vn}，E={e1,e2,…,en}。连接两个顶点的边的权值代表该两顶点之间的距离。对于每个顶点vi，它与各顶点之间的最短路径长度为di1,di2,…,din。顶点vi的最大服务距离是这几个最短路径长度中的最大值，记为e(vi0)。e(vi0)=max(di1,di2,…,din)那么，中心点选址问题，就是求图G的中点vi0，使得该顶点的最大服务距离达到最小，即

e(vi0)=min{e(vi)}中心选址问题的实例例如某县要在其所辖的8个乡镇之一修建一个消防站，为8个乡镇服务，要求消防站至最远乡镇的距离达到最小。假设该8个乡镇之间的交通网络被抽象为图所示的无向赋权连通图，权值为乡镇之间的距离。下面求解消防站应设在哪个乡镇，即哪个顶点？中心选址问题的实例8v6v8v1v5v4v2v39363253757v7中心选址问题的实例首先，用Dijkstra算法计算出每一个顶点vi至其它各顶点vj的最短路径长度dij(i,j=1,2,…,6)，写出距离矩阵：

中心选址问题的实例其次，求距离矩阵中每行的最大值，即各个顶点的最大服务距离，得

e(v1)=14,e(v2)=15,e(v3)=20,e(v4)=12,e(v5)=15,e(v6)=17,e(v7)=12,e(v8)=20

最后计算最大服务距离的最小值。显然，e(v4)=e(v7)=min{e(vi)}。所以，消防站应建在v4或v7点所在的乡镇即可。7．空间插值概念空间插值空间数据的插值即对一组已知空间数据,可以是离散点的形式,也可以是分区数据的形式,要从这些数据中找到一个函数关系式,使该关系式最好地逼近已知的空间数据,并能根据该函数关系式推求出区域范围内其它任意点或任意分区的值。用于将离散点的测量数据转换为连续的数据曲面。广泛应用于等值线自动制图、数字高程模型的建立、不同区域界线现象的相关分析和比较研究等。空间插值理论假设空间位置上越靠近的点，越可能具有相似的特征值；而距离越远的点，其特征值相似的可能性越小。7．空间插值需要进行空间插值的情况：现有的离散曲面的分辨率、像元大小或方向与所要求的不符现有的连续曲面的数据模型与所需的数据模型不符现有的数据不能完全覆盖所要求的区域范围7．空间插值(1)对不足或缺失数据的估计.由于观测台站分布的密度及分布位置的原因,不可能任何空间地点的数据都能实测得到,需要用到插值,以了解区域内观测变量的完整空间分布;(2)数据的网格化.规则格网能够更好地反映连续分布的空间现象,并对他们的变化作出模拟.对已知观测台站的观测数据进行空间内插,便可得到格网化数据;(3)内插等值线.以等值线的形式直观地显示数据的空间分布;(4)对不同分区未知数据的推求空间插值方法的主要目标摘自：朱求安,张万昌,余钧辉.基于GIS的空间插值方法研究[J].江西师范大学学报(自然科学版),2004,28(3),183-188.7．空间插值进行空间插值时,一般包括以下过程:(1)空间插值数据源的获取;(2)对数据进行分析,找出源数据的分布特性、统计特性,便于选择最恰当的插值方法;(3)插值方法的选择并进行插值计算;(4)对插值结果的评价;(5)运用多种插值方法进行计算,对各种方法的插值结果进行比较、分析并选择最佳的插值方法.摘自：朱求安,张万昌,余钧辉.基于GIS的空间插值方法研究[J].江西师范大学学报(自然科学版),2004,28(3),183-188.7．空间插值空间插值的应用：在对自然科学及社会学科的数据处理中,空间插值往往是一个必不可少的工具,尤其在水文、气象、生态、环境、地质等地学学科以及社会、经济等领域得到了广泛的应用.数字地面模型是地学研究中重要的分析工具之一,它的生成便依赖于空间插值方法;在气象气候学中,气象要素(如降水、温度、太阳辐射等)在空间尺度上连续分布的数据,对各类模型的研究有着重要意义，因为各种气象要素的观测站分布往往是稀疏而不均匀的,在各个台站观测的点数据基础上,要推算出空间面上气象要素的分布,空间插值方法是有力的工具;摘自：朱求安,张万昌,余钧辉.基于GIS的空间插值方法研究[J].江西师范大学学报(自然科学版),2004,28(3),183-188.7．空间插值空间插值的应用：在生态学上,空间插值可以用来估算生物种群的分布情况,为病虫害的防治与综合治理等方面提供了理论和实践基础;在地质学上,空间插值方法是探矿及矿藏分布情况的有力分析工具,可以对土壤及所含各种元素的分布情况作出估算，为生产实践提供支持;可以为环境污染情况及扩散的分析提供辅助手段；除了地学上,空间插值方法在社会、经济等方面也得到了很好的应用,也渗透到了其它的学科,如医学上通过空间插值构建复杂性动态脑地形图,对一些脑疾病患者脑电数据进行分析等.摘自：朱求安,张万昌,余钧辉.基于GIS的空间插值方法研究[J].江西师范大学学报(自然科学版),2004,28(3),183-188.7．空间插值空间插值的数据源

摄影测量得到的正射航片或卫星影像卫星或航天飞机的扫描影像野外测量采样数据，采样点随机分布或有规律的线性分布（沿剖面线或沿等高线）数字化的多边形图、等值线图7．空间插值空间插值的数据通常是复杂空间变化有限的采样点的测量数据，这些已知的测量数据称为“硬数据”。如果采样点数据比较少的情况下，可以根据已知的导致某种空间变化的自然过程或现象的信息机理，辅助进行空间插值，这种已知的信息机理，称为“软信息”。但通常情况下，由于不清楚这种自然过程机理，往往不得不对该问题的属性在空间的变化作一些假设，例如假设采样点之间的数据变化是平滑变化，并假设服从某种分布概率和统计稳定性关系。7．空间插值完全规则采样当区域景观大量存在有规律的空间分布模式时,会得到片面的结果，如有规律间隔的数或沟渠空间点采样的方法：7．空间插值完全随机采样随机采样点之间的位置是不相关的完全随机采样会导致采样点的分布不均，一些点的数据密集，另一点的数据缺少。空间点采样的方法：7．空间插值成层随机采样：随机采样和规则采样的结合，即单个的点随机的分布于规则的格网内空间点采样的方法：聚集采样：用于分析不同尺度的空间变化7．空间插值规则断面采样：常用于河流、山坡剖面的测量等值线采样：是数字化等高线图插值数字高程模型最常用的方法。空间点采样的方法：7．空间插值各种不同的采样方式7．空间插值空间插值方法整体插值方法边界内插方法

趋势面分析

变换函数插值多元回归分析其他局部插值方法最近邻点法：泰森多边形方法

移动平均插值方法：距离倒数插值

样条函数插值方法

空间自协方差最佳插值方法：克里金插值

其中最重要的有线性内插法、双线性插值、趋势面插值、样条插值及克里金(Kriging)插值，这些方法有的是整体拟合，有的是局部拟合。用研究区域所有的采样点的数据进行全区特征拟合用邻近的数据点来估计未知点的值7．空间插值空间插值方法的分类空间插值方法依据不同的标准,有多种不同的分类方法.[1]黄杏元等依据已知点和已知分区数据的不同,可将空间数据插值分为点的内插和区域的内插。[2]邬伦等则对应分为空间内插和外推两种:空间内插算法是一种通过已知点的数据推求同一区域其它未知点数据的计算方法;空间外推算法则是通过已知区域的数据,推求其它区域数据的方法.摘自：朱求安,张万昌,余钧辉.基于GIS的空间插值方法研究[J].江西师范大学学报(自然科学版),2004,28(3),183-188.7．空间插值空间插值方法的分类[3]李新等依据空间插值的基本假设和数学本质将空间内插分类为:几何方法、统计方法、空间统计方法、函数方法、随机模拟方法、物理模型模拟方法和综合方法。[4]Nigel.M.Waters在他的文章中对各种分类方法进行了概括,包括

:(1)点的插值和面的插值;(2)全局性插值与局部性插值;(3)精确插值和近似插值;(4)统计性插值和确定性插值;(5)渐变插值和突变插值.7．空间插值空间插值方法的分类[4]Nigel.M.Waters的文章还指出：对于点的插值包括:最近距离法插值,如泰森多边形法,B样条插值,Kriging插值,在经验知识基础上的手工目视插值,趋势面分析法,傅立叶变换及移动距离权重平均法等.对于面的插值(区域插值)计算,根据在插值过程中是否采用辅助数据将其分为两类:无辅助数据的面插值法和有辅助数据的面插值法.7．空间插值——整体插值方法边界内插法假设任何重要的变换发生在边界上，边界内的变化是均匀的，同质的，即在各方向都是相同的。这种概念模型经常用于土壤和景观制图，可以通过定义“均质的”土壤单元、景观图斑来表达其他的土壤、景观特征属性7．空间插值——整体插值方法在边界内插的实际应用中，应考虑数据源是否符合这些理论假设：属性值在图斑或景观单元内是随机变化的，不是有规律的同一类别的所有图斑存在同样的类方差（噪声）所有的属性值都呈正态分布所有的空间变化发生在边界上，是突变而不是渐变7．空间插值——整体插值方法趋势面分析：某种地理属性在空间的连续变化，可以用一个平滑的数学平面加以描述。这种根据采样点的属性数据与地理坐标的关系，进行多元回归分析得到平滑数学平面方程的方法即趋势面分析。思路：先用已知采样点数据拟合出一个平滑的数学平面方程，再根据该方程计算无测量值的点上的数据。该理论假设：地理坐标是独立变量、属性值是独立变量且正态分布，同样回归误差也是与位置无关的独立变量整体内插主要采用趋势面分析技术。7．空间插值——整体插值方法趋势面分析方法适用于:能以空间的视点诠释趋势和残差;观测有限,内插也基于有限的数据.趋势面是个平滑函数,很难正好通过原始数据点,除非是数据点少且趋势面次数高才能使曲面正好通过原始数据点,所以趋势面分析是一个近似插值方法7．空间插值——整体插值方法变换函数插值法根据一个或多个空间参量的经验方程进行整体空间插值，这种经验方程称为变换函数.变换函数插值法属于近似的空间插值。7．空间插值——整体插值方法多元回归分析多元回归在各种统计方法中,使用较多的是回归分析,其特点是不需要分布的先验知识.多元回归在数学形式上与趋势面很相似,但是它们又有着显著的不同.在多元回归中,存在多重共线性,但它并非内在的,可以通过逐步回归解决7．空间插值——整体插值方法其他方法傅里叶级数法,对于周期性的数据序列,如海浪,可以利用傅里叶级数将它们分解为正弦波和余弦波.小波变换方法,小波分析是傅立叶分析发展的重大突破,并已经从傅立叶分析独立出来.傅里叶级数和小波变换方法它们特别是在遥感影象分析方面应用广泛,但它们需要的数据量大.7．空间插值——局部插值方法局部插值步骤：定义一个邻域或搜索范围搜索落在此邻域范围的数据点选择表达这有限个点的空间变化的数学函数为落在规则格网单元上的数据点赋值。重复这个步骤直到格网上的所有点赋值完毕7．空间插值——局部插值方法最近邻点法Thiessen多边形或Voronoi多边形采用一种极端的边界内插方法，只用最近的单个点进行区域插值。泰森多边形按数据点位置将区域分割成子区域，每个子区域包含一个数据点，各子区域到其内数据点的距离小于任何到其他数据点的距离，并用其内数据点进行赋值。连接所有数据点的连线形成Delaunay三角形，与不规则三角网TIN具有相同的拓扑结构。7．空间插值——局部插值方法移动平均插值距离倒数插值综合了泰森多边形的邻近点方法和趋势面分析的渐变方法的长处，假设未知点X0属性值是在局部领域内中所有数据点的距离加权平均值。是加权移动平均方法的一种距离倒数插值方法是GIS软件根据点数据生成栅格图层的最常见方法。7．空间插值——局部插值方法样条函数插值样条函数是一个分段函数，进行一次拟合只有与少数点拟合，同时保证曲线段连接处连续。因此样条函数可以修改少数数据点配准而不必重新计算整条曲线，插值速度快，能够保留局部的变化特征。缺点是其误差不能直接估算。并且在实际中还要解决样条块的定义及如何在三维空间中将这些块拼成复杂曲面，又不引入原始曲面中所没有的异常现象等问题。7．空间插值——局部插值方法样条函数插值常用样条函数是三次样条函数，用于三维情况，此时进行的是曲面内插。采用三次曲面来描述格网内的内插值时，则待定点内插值ZP为：ZP=a1x3y3+a2x2y3+a3xy3+a4y3+a5x3y2 +a6x2y2+a7xy2+a8y2+a9x3y+a10x2y+ a11xy+a12y+a13x3+a14x2+a15x+a16。7．空间插值——局部插值方法样条函数插值B样条函数是感兴趣区间以外均为零的其他样条的和，可按简单的方法用地磁多项式进行局部拟合。常用于数字化的线划在显示之前进行平滑处理，例如土壤、地质图上的各种边界，传统的制图总希望绘出较平滑的曲线。存在的问题是平滑之后周长和面积的计算与平滑前不同。7．空间插值——局部插值方法克里金插值又叫空间自协方差最佳插值方法该方法充分吸收了地理统计的思想，认为任何在空间连续性变化的属性是非常不规则的，不能用简单的平滑数学函数进行模拟，可以用随机表面给予较恰当的描述。这些连续性变化的空间属性称为“区域性变量”，如像气压、高程、及其他连续性变化的描述指标变量。这种应用地理统计方法进行空间插值的方法被称为克里金（Kriging）插值。7．空间插值——局部插值方法克里金插值Kriging插值方法在空间相关范围分析的基础上，用相关范围内的采样点来估计待插点属性值。（1）数据检验与分析，删去明显偏离实际的采样数据点。（2）数据预处理。（3）绘制方差图，了解空间变量的集聚范围与方向。（4）克里金插值估计。相对来说，克里金插值则能较好地反映各种地形变化，但克里金方法的计算量很大，因此在对大面积区域大数据量内插时，这是一个不能不考虑的因素。

线性内插法用线性曲面方程Zp=b0+b1x+b2y拟合待定点附近的地形面（例如TIN模型中的三角平面可视为线性曲面）。可用3个数据点求解待定系数，然后根据线性曲面方程计算待定点的高程值。

7．空间插值——局部插值方法双线性插值法用曲面方程Zp=a0+a1x+a2y+a3xy拟合待定点附近的地形面，也可以表示为 Zp=(ax+b)(cy+d)。即认为待定点附近的地形面的高程Zp在x轴方向呈线性变化，同时在y轴方向上也呈线性变化。可用待定点附近的4个数据点来计算双线性曲面函数的待定系数，再根据曲面函数计算待定点的高程值。

7．空间插值——局部插值方法移动拟合法（1）移动拟合法是用二次多项式来拟合地面高程，即Zp=Ax2+Bxy+Cy2+Dx+Ey+F，通常把坐标原点平移到待定点之后利用上式，采用待定点为圆心的圆周内6个数据点求解函数的待定系数，然后根据曲面方程计算Zp

7．空间插值——局部插值方法移动拟合法（2）

当数据点多于6个时，可按最小二乘法求解待定系数。根据数据点距待定点远近不同，利用“距离越远对插值点的影响越小”的思想，可赋以适当的权值，权重可采用

W=1/di2

W=[(R-di)/di]2

W=exp(-di2/R2)等形式（其中di为距离，R为采样圆半径），该方法称为按距离加权最小二乘内插法。7．空间插值——局部插值方法移动拟合法（3）

首先建立误差方程式：再根据最小二乘法原理，按vwvT=min的方法，建立方程式求解待定系数。 移动拟合法的结果显示误差分布和地形坡度大小没有特别显著的关系，在地性线地区无法显示出这些转折地地形。7．空间插值——局部插值方法8．空间统计分析常规统计分析主要完成对数据集合的均值、总和、方差、频数、峰度系数等参数的统计分析。空间自相关分析是认识空间分布持征、选择适宜的空间尺度来完成空间分析的最常用的方法。2、空间自相关分析1、常规统计分析8.空间统计分析对于空间自相关分析，目前普遍使用空间自相关系数表示，其计算公式如下I

=N—Wij∑∑Wij

(xi–x)(xj–x)xi–x*其中:N表示空间实体数目；xi表示空间实体的属性值；x是xi的平均值；Wij

表示实体i与j的空间关系，它通过拓扑关系获得。

Wij＝1表示空间实体i与j相邻；

Wij

＝0表示空间实体I与j不相邻。I的值介于–1与1之间：

I＝1表示空间自正相关，空间实体呈聚合分布；

I＝-1表示空间自负相关，空间实体呈离散分布；

I＝0表