空间数据的处理精选PPT.ppt

《空间数据的处理精选PPT.ppt》由会员分享，可在线阅读，更多相关《空间数据的处理精选PPT.ppt（86页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、关于空间数据的处理1 1第1页，讲稿共86张，创作于星期二2 2内容提要内容提要第一节第一节 空间数据处理的目的和内容空间数据处理的目的和内容第二节第二节 空间数据的编辑空间数据的编辑第三节第三节 空间数据的坐标变换空间数据的坐标变换第四节第四节 空间数据的插值空间数据的插值第五节第五节 矢量数据向栅格数据的转换矢量数据向栅格数据的转换第六节第六节 栅格数据向矢量数据的转换栅格数据向矢量数据的转换第2页，讲稿共86张，创作于星期二3 3第一节第一节 空间数据处理的目的和内容空间数据处理的目的和内容1、目的、目的(1)数据的净化数据的净化;(2)数据的规范化。数据的规范化。第3页，讲稿共86张，

2、创作于星期二4 42、内容、内容(1)编辑处理；编辑处理；(2)变换处理；变换处理；(3)编码和压缩处理编码和压缩处理;(4)数据的插值；数据的插值；(5)数据类型的转换。数据类型的转换。返回第4页，讲稿共86张，创作于星期二5 5第二节第二节 空间数据的编辑空间数据的编辑1、概述、概述编辑是对输入的数据进行检查、改错、更编辑是对输入的数据进行检查、改错、更新及加工，以得到净化的数据。新及加工，以得到净化的数据。图形数据的编辑分两种：图形数据的编辑分两种：1图形参数编辑：线型、宽、色图形参数编辑：线型、宽、色,符号尺寸、符号尺寸、颜色颜色,面域图案及颜色等设定。面域图案及颜色等设定。2几何数据

3、编辑：点、线、面和目标的编几何数据编辑：点、线、面和目标的编辑。辑。第5页，讲稿共86张，创作于星期二6 62、空间数据编辑的内容、空间数据编辑的内容（1）修改错误数据）修改错误数据图形数据错误图形数据错误:点、线的丢失、点、线的丢失、重复、形状错误等；重复、形状错误等；属性数据错误属性数据错误;图形与属性数据配合错误等。图形与属性数据配合错误等。错误太多时，重新输入更省错误太多时，重新输入更省时间。时间。第6页，讲稿共86张，创作于星期二7 7（2）空间数据的更新）空间数据的更新通过编辑更新数据，确保现势性。通过编辑更新数据，确保现势性。（3）空间数据的预处理和后处理）空间数据的预处理和后处

4、理预处理：拓扑关系建立、曲线光滑、图预处理：拓扑关系建立、曲线光滑、图幅拼接等；幅拼接等；后处理：图形整饰，如经纬线生成，图后处理：图形整饰，如经纬线生成，图例标注等。例标注等。第7页，讲稿共86张，创作于星期二8 8（4）图幅的拼接图幅的拼接1 统一坐标，规范属性统一坐标，规范属性2 相邻图幅编号相邻图幅编号个位表示横向顺序，十位表示纵向顺序。个位表示横向顺序，十位表示纵向顺序。第8页，讲稿共86张，创作于星期二9 9 3 图幅间坐标数据的匹配图幅间坐标数据的匹配同名边界点坐标应在允许范同名边界点坐标应在允许范围内，弧段左右多边形号相围内，弧段左右多边形号相同或相反。同或相反。不匹配时不匹配

5、时:A.修改数据库中点的坐标；修改数据库中点的坐标；B.先对准两幅图的一条边缘先对准两幅图的一条边缘线，再调整其它线段。线，再调整其它线段。4 删除公共边界删除公共边界返回第9页，讲稿共86张，创作于星期二1010(1)平移变换平移变换 x*=x+dx y*=y+dy 其其中，中，dx、dy分别为分别为x,y方向平移量，相应的向方向平移量，相应的向量形式为：量形式为：x*,y*=x,y+dx,dy第三节第三节 空间数据的坐标变换空间数据的坐标变换1.1.图形的几何变换图形的几何变换图形的几何变换图形的几何变换第10页，讲稿共86张，创作于星期二1111(2)比例变换比例变换 x*=x.Sx y

6、*=y.Sy 其其中：中：Sx,Sy分别为分别为x,y方向的比例系数，矩方向的比例系数，矩阵为：阵为：x*,y*=x,y.Sx 0 0 Sy第11页，讲稿共86张，创作于星期二1212(3)旋转变换旋转变换 将点将点(x,y)旋转旋转 角角 x=A.cos,y=A.sin x*=A.cos(+)=A.(cos.cos -sin.sin)=x.con -y.sin y*=A.sin(+)=A(sin.con +con.sin)=x.sin +y.con 矩阵为：矩阵为：x*,y*=x,y.con sin-sin con第12页，讲稿共86张，创作于星期二13132.投影变换投影变换投影变换会变形

7、。投影投影变换会变形。投影变换是一复杂的数学运变换是一复杂的数学运算。可用软件转换。算。可用软件转换。返回第13页，讲稿共86张，创作于星期二1414第四节第四节 空间数据的插值空间数据的插值空间插值（内插）是用已知点的数值估算空间插值（内插）是用已知点的数值估算未知点的数值的过程未知点的数值的过程。例如，在一个没有。例如，在一个没有记录数据的地点，其降水量可通过对附近记录数据的地点，其降水量可通过对附近气象站已知降水量记录的插值来估算出来。气象站已知降水量记录的插值来估算出来。在在GIS应用中，空间插值主要用于栅格数据，应用中，空间插值主要用于栅格数据，估算出格网中每个单元的值。因此，空间估

8、算出格网中每个单元的值。因此，空间插值是将点数据转换成面数据的一种方法。插值是将点数据转换成面数据的一种方法。第14页，讲稿共86张，创作于星期二1515连续空间连续空间:采用渐变模型采用渐变模型-点插值。点插值。离散空间离散空间:邻近元法邻近元法-区域区域插值。插值。第15页，讲稿共86张，创作于星期二1616一、一、点插值点插值整体拟合：模型由研究区域内所有采样点整体拟合：模型由研究区域内所有采样点上的全部特征观测值建立。（全局方法）上的全部特征观测值建立。（全局方法）局部拟合：用邻近的数据点来估计未知点局部拟合：用邻近的数据点来估计未知点的值。（局部方法）的值。（局部方法）返回第16页，

9、讲稿共86张，创作于星期二17171、全局方法全局方法 利用每个可利用的控制点来构建一个方程利用每个可利用的控制点来构建一个方程或模型，该模型可用于估算未知数值。或模型，该模型可用于估算未知数值。（1）趋势面分析趋势面分析 趋势面分析用多项式方程拟合已知数值的趋势面分析用多项式方程拟合已知数值的点点。该方程又称趋势面模型，线性趋势面。该方程又称趋势面模型，线性趋势面方程：方程：zx,yb0+b1x+b2y 特征值特征值z是是x和和y的函数。系数的函数。系数b由控制点估算。其拟合由控制点估算。其拟合程度可用相关系数确定。每个已知点观测值和估算值程度可用相关系数确定。每个已知点观测值和估算值之间的

10、偏差或残差可以计算出来。之间的偏差或残差可以计算出来。第17页，讲稿共86张，创作于星期二1818拟合更复杂的面要求更高次的趋势面模型，拟合更复杂的面要求更高次的趋势面模型，如三次趋势面方程：如三次趋势面方程：zx,yb0+b1x+b2y+b3x2+b4xy+b5y2+b6x3+b7x2y+b8xy2+b9y3三次面包含山和谷。统计检验能用于比较三次面包含山和谷。统计检验能用于比较几种可选模型的拟合程度。几种可选模型的拟合程度。GIS软件包如软件包如ARC/INFO提供了多达提供了多达12次趋势面模型。图次趋势面模型。图13.3是由是由105个数据点构建的爱达荷州年平均降水量的三阶趋势个数据点

11、构建的爱达荷州年平均降水量的三阶趋势面等值线图。面等值线图。第18页，讲稿共86张，创作于星期二1919第19页，讲稿共86张，创作于星期二2020(2)回归模型回归模型 回归模型把方程中的一个因变量与多个自变量联系起回归模型把方程中的一个因变量与多个自变量联系起来，而后可用于预测或估算。来，而后可用于预测或估算。例例:用雪水量用雪水量(SWE)作为因作为因变量，位置和地形变量作为自变量，其流域模型的表达形式如下：变量，位置和地形变量作为自变量，其流域模型的表达形式如下：SWEb0+blEASTING+b2SOUTHING+b3ELEV+b4PLANl000 这里这里EASTING和和SOUT

12、HING分别对应一个高程格网的行数和列数，分别对应一个高程格网的行数和列数，ELEV是高程值，是高程值，PLANl000是表面曲率的测度。是表面曲率的测度。第20页，讲稿共86张，创作于星期二21212 局部方法局部方法该方法用一组控制点样本来估算未知值，因此控制点对该方法用一组控制点样本来估算未知值，因此控制点对该方法十分重要。第一要点是用于估算的控制点个数。该方法十分重要。第一要点是用于估算的控制点个数。GIS软件包通常允许用户自己确定控制点的个数或软件包通常允许用户自己确定控制点的个数或用默认值用默认值(712个个)。你也许认为控制点越多，估算结。你也许认为控制点越多，估算结果越精确。然

13、而这种设想的正确与否取决于控制点的分果越精确。然而这种设想的正确与否取决于控制点的分布与要估算单元的关系以及空间自相关程度。控制点越布与要估算单元的关系以及空间自相关程度。控制点越多通常意味着估算越通用化。多通常意味着估算越通用化。第21页，讲稿共86张，创作于星期二2222 控制点个数确定之后，下一步就是控制点选择。简单控制点个数确定之后，下一步就是控制点选择。简单办法是用离估算点最近的点作为控制点，另一种办法是办法是用离估算点最近的点作为控制点，另一种办法是用半径来选择控制点，半径的大小必须根据控制点的分用半径来选择控制点，半径的大小必须根据控制点的分布来调整。例如，布来调整。例如，ARC

14、/INFO和和ArcView都提供了这两都提供了这两种方法。其他方法考虑了方向因素，例如，象限法或八方种方法。其他方法考虑了方向因素，例如，象限法或八方向法。象限法是从围绕每个要估算单元的四个方向选择控向法。象限法是从围绕每个要估算单元的四个方向选择控制点，八方向法是从八个方向来选择控制点。制点，八方向法是从八个方向来选择控制点。第22页，讲稿共86张，创作于星期二2323（1）泰森多边形泰森多边形 泰森泰森(Thiessen)多边形围绕已知点样本构建多边形围绕已知点样本构建而成，使得在泰森多边形内的任意点与多而成，使得在泰森多边形内的任意点与多边形内的已知点更接近，而不是与其他已边形内的已知

15、点更接近，而不是与其他已知点接近。泰森多边形最初用于估算区域知点接近。泰森多边形最初用于估算区域降水量的平均值。降水量的平均值。第23页，讲稿共86张，创作于星期二2424泰森多边形要求在已知点之间构建初始三角形，亦即连接泰森多边形要求在已知点之间构建初始三角形，亦即连接已知点形成三角形。因为连接点的方法不同，形成的三角已知点形成三角形。因为连接点的方法不同，形成的三角形也不同。与构建不规则三角网相似，德劳内三角测量常形也不同。与构建不规则三角网相似，德劳内三角测量常用于构建泰森多边形。德劳内三角网测量确保了每个已知用于构建泰森多边形。德劳内三角网测量确保了每个已知点都与它最近的点相接，这样就

16、使得三角形尽量接近等边。点都与它最近的点相接，这样就使得三角形尽量接近等边。经过三角形每条边的中点画垂线，连接起来就可以很容易经过三角形每条边的中点画垂线，连接起来就可以很容易地构建泰森多边形地构建泰森多边形(图图13.4)。泰森多边形也称为冯罗诺。泰森多边形也称为冯罗诺(Voronoi)多边形。多边形。第24页，讲稿共86张，创作于星期二2525第25页，讲稿共86张，创作于星期二2626(2)密度估算密度估算密度估算基于点的分布及其已知值来量测密度估算基于点的分布及其已知值来量测格网的密度格网的密度。一种简单密度估算方法是将。一种简单密度估算方法是将格网置于点分布图上，将落在每个单元的格网

17、置于点分布图上，将落在每个单元的点列表，将点值加和，将单元的点值总和点列表，将点值加和，将单元的点值总和除以单元大小，就估算得每个单元的密度。除以单元大小，就估算得每个单元的密度。第26页，讲稿共86张，创作于星期二2727图图13.5显示一个简单密度估显示一个简单密度估算的例子。输入是以算的例子。输入是以50m间间距标绘的看到鹿的地点的分距标绘的看到鹿的地点的分布。每个鹿的地点有一个计布。每个鹿的地点有一个计数值，表示在该地一只鹿被数值，表示在该地一只鹿被看到多少次。输出是密度格看到多少次。输出是密度格网：单元大小是网：单元大小是l hm2，密度，密度用每公顷看到的次数来表示。用每公顷看到的

18、次数来表示。第27页，讲稿共86张，创作于星期二2828核估算核估算(kernel estimation)是一种不同的密是一种不同的密度估算方法，它用核函数联系每一个点或度估算方法，它用核函数联系每一个点或观测点。表达为双变量概率密度函数，核观测点。表达为双变量概率密度函数，核函数看起来像是一个隆起函数看起来像是一个隆起(bump)，以一个，以一个点为中心在一个定义的带宽或窗口范围内点为中心在一个定义的带宽或窗口范围内逐渐减小到逐渐减小到0(图图13.6)。核函数和带宽决定了。核函数和带宽决定了隆起的形状，其形状反过来决定了在估算隆起的形状，其形状反过来决定了在估算中的平滑量。在点中的平滑量。

19、在点x上的核密度估算值是带上的核密度估算值是带宽范围内位于观测点宽范围内位于观测点xi上的隆起部分的总和。上的隆起部分的总和。第28页，讲稿共86张，创作于星期二2929(x-xi)和和(y-yi)是带宽范围内的点是带宽范围内的点x和观测点和观测点xi之间的之间的x、y坐标的偏差。坐标的偏差。第29页，讲稿共86张，创作于星期二3030使用与简单估算方法相同使用与简单估算方法相同的输入数据，图的输入数据，图13.7显示显示核估算法构建的密度格网核估算法构建的密度格网输出结果。格网上的密度输出结果。格网上的密度值是期望值值是期望值(输出值的平均数输出值的平均数)而不是概率。虽然图而不是概率。虽然

20、图13.7和和图图13.5间的差别不明显，间的差别不明显，但核估算法通常能产生比但核估算法通常能产生比简单估算法更平滑的密度简单估算法更平滑的密度面。面。第30页，讲稿共86张，创作于星期二3131(3)反距离权重插值反距离权重插值该方法该方法假设未知值的点受较近控制点的影假设未知值的点受较近控制点的影响比较远控制点的影响更大响比较远控制点的影响更大。影响的程度影响的程度(或权重或权重)用点之间距离乘方的倒数表示用点之间距离乘方的倒数表示。乘。乘方为方为1.0意味着点之间数值变化率为恒定，意味着点之间数值变化率为恒定，该方法称为线性插值法。乘方为该方法称为线性插值法。乘方为2.0或更高或更高则

21、意味着越靠近已知点，数值的变化率越则意味着越靠近已知点，数值的变化率越大，远离已知点趋于平稳。大，远离已知点趋于平稳。第31页，讲稿共86张，创作于星期二3232式中：式中：z0是点是点0的估计值；的估计值；zi是控制点是控制点i的的z值；值；di是是控制点控制点i与点与点0间的距离；间的距离；s是在估算中用到的是在估算中用到的控制点的数目；控制点的数目；k是指定的幂。是指定的幂。第32页，讲稿共86张，创作于星期二3333图图13.8显示对爱达荷州显示对爱达荷州105个气象站点数据由反距离平个气象站点数据由反距离平方方法建立的年平均降水量曲面。图方方法建立的年平均降水量曲面。图13.9显示曲

22、面的等显示曲面的等雨量线图。雨量线图。第33页，讲稿共86张，创作于星期二3434(4)薄板样条函数法薄板样条函数法薄板样条函数薄板样条函数(Thin-plate splines)建立一个建立一个通过控制点的面，并使所有点的坡度变化通过控制点的面，并使所有点的坡度变化最小。换言之，最小。换言之，薄板样条函数以最小曲率薄板样条函数以最小曲率面拟合控制点面拟合控制点。薄板样条函数的估计值由。薄板样条函数的估计值由下式计算：下式计算：第34页，讲稿共86张，创作于星期二3535Q(x,y)Aidi2logdi+a+bx+cy 式中：式中：x和和y是要被插值的点的是要被插值的点的x、y坐标；坐标；di

23、2(x-xi)2+(y-yi)2；xi和和yi是控制点是控制点i的的 x、y坐标。薄板样条函数包括坐标。薄板样条函数包括二个部分：二个部分：(a+bx+cy)表示局部趋势函数，表示局部趋势函数，di2logdi表表示基本函数，可获得最小曲率的面。系数示基本函数，可获得最小曲率的面。系数Ai、a、b和和c由以下线性方程组决定：由以下线性方程组决定：式中：n为控制点的数目，fi为控制点i的已知值，系数的计算要求n+3联立方程.第35页，讲稿共86张，创作于星期二3636 薄板样条函数的一个主要问题是在数据薄板样条函数的一个主要问题是在数据贫乏地区的坡度较大，常出现过伸贫乏地区的坡度较大，常出现过伸

24、(overshoots)。用于校正过伸的方法包括薄。用于校正过伸的方法包括薄板张力样条板张力样条(thin-piate splines with tension)，规则样条，规则样条(regularized splines)和规则张力和规则张力样条样条(regularized splines with tension)等。等。第36页，讲稿共86张，创作于星期二3737规则样条函数的近似值规则样条函数的近似值与薄板样条函数有相同与薄板样条函数有相同的的局部趋势函数局部趋势函数，但是，但是基本函数基本函数取不同形式取不同形式：n n式中式中式中式中 是样条法中要用到的权重；是样条法中要用到的权重

25、；是样条法中要用到的权重；是样条法中要用到的权重；d d是待定值的点和控制点是待定值的点和控制点是待定值的点和控制点是待定值的点和控制点i i间的距间的距间的距间的距离，离，离，离，c c是常数是常数是常数是常数0.5772150.577215；KK0 0(d/(d/)是修正的零次贝塞耳是修正的零次贝塞耳是修正的零次贝塞耳是修正的零次贝塞耳(Besscl)(Besscl)函数函数函数函数它可由一个多项式方程估计。它可由一个多项式方程估计。它可由一个多项式方程估计。它可由一个多项式方程估计。值通常设为值通常设为值通常设为值通常设为0 00.50.5之间因为更之间因为更之间因为更之间因为更大的大的

26、大的大的 值会导致数据贫乏地区趋于过伸。值会导致数据贫乏地区趋于过伸。值会导致数据贫乏地区趋于过伸。值会导致数据贫乏地区趋于过伸。ARC/INFOARC/INFO和和和和ArcViewArcView采用的默认采用的默认采用的默认采用的默认 值为值为值为值为0.10.1Q(x,y)Aidi2logdi+a+bx+cy第37页，讲稿共86张，创作于星期二3838薄板张力样条薄板张力样条(thin-plate splines with tension)法表达式：法表达式：式中：式中：是本张力法要用到的权重。如果是本张力法要用到的权重。如果是本张力法要用到的权重。如果是本张力法要用到的权重。如果 的权

27、重被设的权重被设为接近于为接近于0 0，则用张力法与基本薄板样条法得到的估，则用张力法与基本薄板样条法得到的估，则用张力法与基本薄板样条法得到的估，则用张力法与基本薄板样条法得到的估计值相似。计值相似。计值相似。计值相似。ARC/INFO和和ArcViewArcView采用的默认采用的默认采用的默认采用的默认 值为值为值为值为0.1。第38页，讲稿共86张，创作于星期二3939 薄板样条函数及其变种适用于平滑和连续的面。样条薄板样条函数及其变种适用于平滑和连续的面。样条法也适用于气候数据的插值。图法也适用于气候数据的插值。图13.10，13.11分别显示分别显示用规则样条函数和张力样条函数法建

28、立的年均降水量用规则样条函数和张力样条函数法建立的年均降水量面。两幅地图中的等雨量线都很平滑，比图面。两幅地图中的等雨量线都很平滑，比图13.9的等的等雨量线雨量线(用反距离平方插值法建立的用反距离平方插值法建立的)更为平滑。在数据贫更为平滑。在数据贫乏地区用样条法插值的结果可能会比实际情况偏大，如图乏地区用样条法插值的结果可能会比实际情况偏大，如图13.10的北爱达荷的北爱达荷 45in等雨量线所例证的。等雨量线所例证的。第39页，讲稿共86张，创作于星期二4040第40页，讲稿共86张，创作于星期二4141(5)克里金法克里金法克里金法克里金法(Kriging)是由南非采矿工程师是由南非采

29、矿工程师D.G.Krigee创立的一种用于空间插值的地理统计创立的一种用于空间插值的地理统计方法。克里金法的原理是假设某种属性的空方法。克里金法的原理是假设某种属性的空间变化既不是完全随机也不是完全确定。空间变化既不是完全随机也不是完全确定。空间变化可能包括三种影响因素：空间相关因间变化可能包括三种影响因素：空间相关因素，代表区域变量的变化；偏移或结构，代素，代表区域变量的变化；偏移或结构，代表趋势；还有随机误差。偏移出现与否和对表趋势；还有随机误差。偏移出现与否和对区域变量的解释导致了不同克里金法的出现。区域变量的解释导致了不同克里金法的出现。第41页，讲稿共86张，创作于星期二4242 1

30、.普通克里金法普通克里金法 假设不存在偏移，普通克里金法关注于空假设不存在偏移，普通克里金法关注于空间相关因素。衡量所选已知点之间空间相间相关因素。衡量所选已知点之间空间相关程度的测度是半方差关程度的测度是半方差(semivariance)，由，由下式计算：下式计算：h是已知点之间的距离，是已知点之间的距离，n样本点的数量，样本点的数量，z是属性是属性值。值。第42页，讲稿共86张，创作于星期二4343在不同距离的半方差值算出后，绘成半方差图，在不同距离的半方差值算出后，绘成半方差图，y轴代表偏差，轴代表偏差，x轴代表已知点之间的距离。半轴代表已知点之间的距离。半方差可分成三部分：熔核方差可分

31、成三部分：熔核(nugget)、值域、值域(range)和基台和基台(sill)。熔核是在距离为。熔核是在距离为0处的半处的半方差，代表无关的空间噪音。值域是半方差的方差，代表无关的空间噪音。值域是半方差的空间相关部分，它显示半方差随着距离递增。空间相关部分，它显示半方差随着距离递增。超过值域范围，半方差趋平于相对恒定值。达超过值域范围，半方差趋平于相对恒定值。达到恒定的半方差称为基台到恒定的半方差称为基台(sill)。第43页，讲稿共86张，创作于星期二4444第44页，讲稿共86张，创作于星期二4545半方差图半方差图(semivariogram)将半方差与距离关联起来。它将半方差与距离关

32、联起来。它可单独用作空间相关的测度。但在克里金法中用作插可单独用作空间相关的测度。但在克里金法中用作插值器值器(interpolator)，半方差图必须与数学函数或模型拟合。，半方差图必须与数学函数或模型拟合。例如，例如，ARC/INFO提供了五种模型：高斯、线性、球面、提供了五种模型：高斯、线性、球面、圆形和指数模型。拟合的半方差图便可用于估算任何给圆形和指数模型。拟合的半方差图便可用于估算任何给定距离的半方差。定距离的半方差。第45页，讲稿共86张，创作于星期二4646普通克里金法在空间插值中直接使用拟合普通克里金法在空间插值中直接使用拟合半方差图。估算某点的半方差图。估算某点的z值的通用

33、方程是：值的通用方程是：式中式中:z0是估计值，是估计值，zx是已知点的值；是已知点的值；Wx是是与每个已知点关联的权重；与每个已知点关联的权重；s是用于估算的是用于估算的已知点的数目。权重可由对一组联立方程已知点的数目。权重可由对一组联立方程的求解得到，例如，下列联立方程是由三的求解得到，例如，下列联立方程是由三个已知点估算一个未知点的值所必须的：个已知点估算一个未知点的值所必须的：第46页，讲稿共86张，创作于星期二4747 式中：式中：(hij)是已知点是已知点i和和j间的半方差，间的半方差，(hi0)是已知点和未知点之间的半方差，是已知点和未知点之间的半方差，是拉格朗日是拉格朗日系数，

34、它的加入是为了确保把估算误差降到最系数，它的加入是为了确保把估算误差降到最小。上面的方程可改写成矩阵形式：小。上面的方程可改写成矩阵形式：第47页，讲稿共86张，创作于星期二4848将左边的逆矩阵乘以右边的矩阵，可解得将左边的逆矩阵乘以右边的矩阵，可解得权重向量。一旦知道了权重，方程权重向量。一旦知道了权重，方程(13.17)便可用估算便可用估算z0：z0z1W1+z2W2+z3W3第48页，讲稿共86张，创作于星期二4949上例表明在上例表明在克里金法中用的权重不仅包括估算点和已知克里金法中用的权重不仅包括估算点和已知点之间的半方差，而且包括已知点之间的半方差点之间的半方差，而且包括已知点之

35、间的半方差。这与。这与反距离权重插值法不同，后者只用适用于已知点和估算反距离权重插值法不同，后者只用适用于已知点和估算点的权重。克里金法和其他局部方法之间的另一个重要点的权重。克里金法和其他局部方法之间的另一个重要区别是：区别是：克里金法对每个估算点产生一个方差测度克里金法对每个估算点产生一个方差测度来衡量估算值的可靠性来衡量估算值的可靠性。上例中，方差可由下式计。上例中，方差可由下式计算：算：第49页，讲稿共86张，创作于星期二5050图图13.14显示由线性模型普通克里金法生成的年平均降水量曲面。图显示由线性模型普通克里金法生成的年平均降水量曲面。图 13.15显示克里金曲面的标准差分布。

36、数据贫乏地区的标准差显示克里金曲面的标准差分布。数据贫乏地区的标准差最高，例如，爱达荷州的西南角。最高，例如，爱达荷州的西南角。第50页，讲稿共86张，创作于星期二51512通用克里金法通用克里金法通用克里金法通用克里金法(Universal Kriging)假设除了假设除了已知点之间的空间关系外，空间变量在已知点之间的空间关系外，空间变量在z值值上还有偏移或有结构因素。一般地说，通上还有偏移或有结构因素。一般地说，通用克里金法具体表现为一个趋势，例如，用克里金法具体表现为一个趋势，例如，在克里金过程中的一个趋势面方程。在克里金过程中的一个趋势面方程。ARCINFO提供两种通用克里金法。提供两

37、种通用克里金法。Universal 1使用一个平面，定义为如下一阶多项式：使用一个平面，定义为如下一阶多项式：第51页，讲稿共86张，创作于星期二5252Mb1xi+b2yi式中：式中：M是偏移；是偏移；xi和和yi是已知点是已知点i的的x和和y坐坐标；标；b1和和b2是要估计的偏移系数。是要估计的偏移系数。Universal 2用二阶曲面，它被定义为一个二用二阶曲面，它被定义为一个二阶多项式：阶多项式：M=b1xi+b2yi+b3x2+b4xiyi+b5y2系数系数bi必须按权重估算，这意味着通用克里必须按权重估算，这意味着通用克里金法比普通克里金法要求更多的联立方程金法比普通克里金法要求更

38、多的联立方程用于估算未知值。用于估算未知值。第52页，讲稿共86张，创作于星期二5353 图图13.17显示用通用克里金显示用通用克里金法与线性偏移建立的年均降法与线性偏移建立的年均降水量曲面，图水量曲面，图13.18显示在克显示在克里金曲面上标准差的分布。里金曲面上标准差的分布。两幅地图与基于线性模型克两幅地图与基于线性模型克里金法建立的地图相似。主里金法建立的地图相似。主要的不同是在爱达荷州西南要的不同是在爱达荷州西南角出现了负值，这是因为通角出现了负值，这是因为通用克里金模型中线性偏移而用克里金模型中线性偏移而导致插值不准的结果。导致插值不准的结果。第53页，讲稿共86张，创作于星期二5

39、454 3其他克里金法其他克里金法除了普通克里金法和通用克里金法，还有其他的克除了普通克里金法和通用克里金法，还有其他的克里金法，例如里金法，例如:克里金块内插法克里金块内插法(block kriging)和联合克和联合克里金法里金法(co-kriging)。克里金块内插法估算某个小范围内或。克里金块内插法估算某个小范围内或块内的一个变量的平均值，联合克里金法在插值时增加次块内的一个变量的平均值，联合克里金法在插值时增加次要变量，这些次要变量与主要变量有相关关系。这种方法要变量，这些次要变量与主要变量有相关关系。这种方法的原理是：变量之间的相关关系能用于提高主要变量的预的原理是：变量之间的相关

40、关系能用于提高主要变量的预测值精度。例如，用联合克里金法，在降水量插值中，将测值精度。例如，用联合克里金法，在降水量插值中，将高程作为一个附加变量会获得更好结果高程作为一个附加变量会获得更好结果。第54页，讲稿共86张，创作于星期二5555(6)各种局部方法的比较各种局部方法的比较用相同数据集，不同的局部方法会产生不同的结果。图用相同数据集，不同的局部方法会产生不同的结果。图13.19显示规则样条函数和距离倒转平方法两种方法得出显示规则样条函数和距离倒转平方法两种方法得出不同的年平均降水量估算结果。图不同的年平均降水量估算结果。图13.20显示用线性模型显示用线性模型普通克里金法和规则样条函数

41、得出的不同结果。虽然普普通克里金法和规则样条函数得出的不同结果。虽然普通克里金法与规则样条函数两种方法间的差异比规则样通克里金法与规则样条函数两种方法间的差异比规则样条函数与倒转距离平方法间的差异更小，但这两幅地图条函数与倒转距离平方法间的差异更小，但这两幅地图有相似的总体模式。另外，正负误差多于有相似的总体模式。另外，正负误差多于3in的，都位的，都位于数据贫乏地区。这说明不管用哪种方法插值，都不于数据贫乏地区。这说明不管用哪种方法插值，都不能代替观测数据。能代替观测数据。第55页，讲稿共86张，创作于星期二5656第56页，讲稿共86张，创作于星期二5757 将估算值的差别绘成地图并不能确

42、定哪种插值方法更精确。将估算值的差别绘成地图并不能确定哪种插值方法更精确。交叉验证交叉验证(cross-validation)分析通常用于确定不同插值方分析通常用于确定不同插值方法的精度。对每一种插值方法，交叉验证分析重复从已知法的精度。对每一种插值方法，交叉验证分析重复从已知数值集中删除一个已知点的过程，用剩下的已知点估算被数值集中删除一个已知点的过程，用剩下的已知点估算被删除点的数值，并计算估算值与已知值的误差。然后计算删除点的数值，并计算估算值与已知值的误差。然后计算出误差统计值用于评价插值方法的精度。例如，交叉验证出误差统计值用于评价插值方法的精度。例如，交叉验证分析表明，在消除高程倾

43、向克里金法分析表明，在消除高程倾向克里金法(elevation-detrended kriging)或联合克里金法或联合克里金法(co-kriging)中用高程作为中用高程作为x、y位置的附加变量，产生了比克里金法好的结果。位置的附加变量，产生了比克里金法好的结果。第57页，讲稿共86张，创作于星期二5858 评定不同插值法精度的另一种方法是将已知点分成两部分样本：一部分样本用于给每种插值方法建模，另一部分样本用于检测各种方法的精度。如果已知点的数目太少，则这种样本分离的方法不适用。第58页，讲稿共86张，创作于星期二5959二、二、区域插值区域插值适用于离散空间数据。适用于离散空间数据。区域

44、插值是从一种已知的分区数据中推断区域插值是从一种已知的分区数据中推断出同一地区的另一种分区数值的插值方法。出同一地区的另一种分区数值的插值方法。第59页，讲稿共86张，创作于星期二6060 第60页，讲稿共86张，创作于星期二6161(1)叠置法叠置法认为源区和目标区内数据均匀分布。求解时，将认为源区和目标区内数据均匀分布。求解时，将目标区叠置到源区上，求出源区和目标区之间的目标区叠置到源区上，求出源区和目标区之间的交集。目标区计算公式：交集。目标区计算公式：n UDi=USj Dj j=1 其中：其中：UDi为目标区第为目标区第i分区的值，分区的值，i=1,2.m；USj为源区第为源区第j个

45、分区的值，个分区的值，j=1,2,.n；Dj 为第为第i个目标分区同有第个目标分区同有第j个源区的面积比。个源区的面积比。第61页，讲稿共86张，创作于星期二6262 UD1=35x(3/7)+30 x(2/6)=25UD2=35x(4/7)+10 x(1/3)=23.3UD3=30 x(4/6)+10 x(2/3)=26.6返回第62页，讲稿共86张，创作于星期二6363(2)比重法比重法将源区的统计数据从同质性变成非同质性，将源区的统计数据从同质性变成非同质性，再进行区域插值。再进行区域插值。1 叠置格网，确定精度。叠置格网，确定精度。2 将源区各分区的平均数赋将源区各分区的平均数赋予各格

46、网点（图予各格网点（图a)。第63页，讲稿共86张，创作于星期二64643 计算相邻的四计算相邻的四(八八)个格网点的平均值个格网点的平均值(图图b)。第64页，讲稿共86张，创作于星期二65654 以分区为单位将以分区为单位将3的平均的平均值相加，得值为值相加，得值为Us，计算各分计算各分区系数区系数pi=Us/Us，Usi为源区为源区i分分区的人口数。并将各个格网点值区的人口数。并将各个格网点值乘以乘以pi，得调整后的格网值，得调整后的格网值(图图c)。第65页，讲稿共86张，创作于星期二66665不断循环不断循环2,3,4步，直到步，直到pi接近接近1(0.99)为止，为止，便便 可计可

47、计算目标区的内插值算目标区的内插值(图图d)。例如，例如，Ud1区的内插值区的内插值=5.5+5.3+5.5+5.2+4.5=26返回第66页，讲稿共86张，创作于星期二6767第五节第五节 矢量数据向栅格数据的转换矢量数据向栅格数据的转换先要确定分辨率。先要确定分辨率。矢量数据是直角坐标，原点在左下方；矢量数据是直角坐标，原点在左下方；栅格数据是行列坐标，原点在左上方。栅格数据是行列坐标，原点在左上方。第67页，讲稿共86张，创作于星期二6868 其转换公式为：其转换公式为：X=(Xmax-Xmin)/J Y=(Ymax-Ymin)/I X，Y表示栅格边表示栅格边长。长。I、J为行、列数。为

48、行、列数。第68页，讲稿共86张，创作于星期二6969例如：已知某一地区例如：已知某一地区x方向为方向为15km，y方方向为向为30km，现要把该地区的地块图转换，现要把该地区的地块图转换成栅格数据。要求栅格的分辨率为成栅格数据。要求栅格的分辨率为30mx30m，则：，则：行数行数I=30km/30m=1000 列数列数J=15km/30m=500第69页，讲稿共86张，创作于星期二7070一、边界代数法一、边界代数法转换的关键是对矢量表示的多边形边界转换的关键是对矢量表示的多边形边界内的所有栅格赋予多边形编码，形成栅内的所有栅格赋予多边形编码，形成栅格数据阵列。格数据阵列。面的填充是根据边界

49、的拓扑信息，通过面的填充是根据边界的拓扑信息，通过加减运算将边界位置信息动态地赋予各加减运算将边界位置信息动态地赋予各栅格。栅格。实现的前提是知道拓扑关系。实现的前提是知道拓扑关系。第70页，讲稿共86张，创作于星期二7171 填充过程：沿边界前进方向填充过程：沿边界前进方向y值上升时为值上升时为上行，上行，y值下降时为下行。上行时填充值值下降时为下行。上行时填充值为左多边形号减右多边形号，下行时填为左多边形号减右多边形号，下行时填充值为右多边形减左多边形号，将每次充值为右多边形减左多边形号，将每次填充值同该处的原始值作代数运算得到填充值同该处的原始值作代数运算得到最终填充属性值。最终填充属性

50、值。第71页，讲稿共86张，创作于星期二7272N1N2弧上行，左弧上行，左-右为右为0-1=-1。弧左边栅。弧左边栅格值为格值为-1第72页，讲稿共86张，创作于星期二7373N2N3弧下行，右弧下行，右-左为左为2-0=2，弧左边栅格，弧左边栅格值加值加2第73页，讲稿共86张，创作于星期二7474N3N1弧下行，右弧下行，右-左左=3-0=3。弧左边栅格。弧左边栅格加加3第74页，讲稿共86张，创作于星期二7575N1N4弧上行，左弧上行，左-右右=1-3=-2。弧左边加。弧左边加-2。第75页，讲稿共86张，创作于星期二7676N4N2弧上行，左弧上行，左-右右=1-2=-1，弧左边加