空间数据的处理讲稿.ppt

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

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

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

3、点、线、面和目标的编辑。辑。第五页，讲稿共八十六页哦6 62、空间数据编辑的内容、空间数据编辑的内容（1）修改错误数据）修改错误数据图形数据错误图形数据错误:点、线的丢失、点、线的丢失、重复、形状错误等；重复、形状错误等；属性数据错误属性数据错误;图形与属性数据配合错误等。图形与属性数据配合错误等。错误太多时，重新输入更省错误太多时，重新输入更省时间。时间。第六页，讲稿共八十六页哦7 7（2）空间数据的更新）空间数据的更新通过编辑更新数据，确保现势性。通过编辑更新数据，确保现势性。（3）空间数据的预处理和后处理）空间数据的预处理和后处理预处理：拓扑关系建立、曲线光滑、图预处理：拓扑关系建立、曲

4、线光滑、图幅拼接等；幅拼接等；后处理：图形整饰，如经纬线生成，图后处理：图形整饰，如经纬线生成，图例标注等。例标注等。第七页，讲稿共八十六页哦8 8（4）图幅的拼接图幅的拼接1 统一坐标，规范属性统一坐标，规范属性2 相邻图幅编号相邻图幅编号个位表示横向顺序，十位表示纵向顺序。个位表示横向顺序，十位表示纵向顺序。第八页，讲稿共八十六页哦9 9 3 图幅间坐标数据的匹配图幅间坐标数据的匹配同名边界点坐标应在允许范同名边界点坐标应在允许范围内，弧段左右多边形号相围内，弧段左右多边形号相同或相反。同或相反。不匹配时不匹配时:A.修改数据库中点的坐标；修改数据库中点的坐标；B.先对准两幅图的一条边缘先

5、对准两幅图的一条边缘线，再调整其它线段。线，再调整其它线段。4 删除公共边界删除公共边界返回第九页，讲稿共八十六页哦1010(1)平移变换平移变换 x*=x+dx y*=y+dy 其其中，中，dx、dy分别为分别为x,y方向平移量，相应的向方向平移量，相应的向量形式为：量形式为：x*,y*=x,y+dx,dy第三节第三节 空间数据的坐标变换空间数据的坐标变换1.1.图形的几何变换图形的几何变换图形的几何变换图形的几何变换第十页，讲稿共八十六页哦1111(2)比例变换比例变换 x*=x.Sx y*=y.Sy 其其中：中：Sx,Sy分别为分别为x,y方向的比例系数，矩方向的比例系数，矩阵为：阵为：

6、x*,y*=x,y.Sx 0 0 Sy第十一页，讲稿共八十六页哦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第十二页，讲稿共八十六页哦13132.投影变换投影变换投影变换会变形。投影投影变换会变形。投影变换是一复杂的数学运变换是一复杂的数学运算。可用软件转换。算。可用软件转换。返回第十三页，讲稿共

7、八十六页哦1414第四节第四节 空间数据的插值空间数据的插值空间插值（内插）是用已知点的数值估算空间插值（内插）是用已知点的数值估算未知点的数值的过程未知点的数值的过程。例如，在一个没有。例如，在一个没有记录数据的地点，其降水量可通过对附近记录数据的地点，其降水量可通过对附近气象站已知降水量记录的插值来估算出来。气象站已知降水量记录的插值来估算出来。在在GIS应用中，空间插值主要用于栅格数据，应用中，空间插值主要用于栅格数据，估算出格网中每个单元的值。因此，空间估算出格网中每个单元的值。因此，空间插值是将点数据转换成面数据的一种方法。插值是将点数据转换成面数据的一种方法。第十四页，讲稿共八十六

8、页哦1515连续空间连续空间:采用渐变模型采用渐变模型-点插值。点插值。离散空间离散空间:邻近元法邻近元法-区域区域插值。插值。第十五页，讲稿共八十六页哦1616一、一、点插值点插值整体拟合：模型由研究区域内所有采样点整体拟合：模型由研究区域内所有采样点上的全部特征观测值建立。（全局方法）上的全部特征观测值建立。（全局方法）局部拟合：用邻近的数据点来估计未知点局部拟合：用邻近的数据点来估计未知点的值。（局部方法）的值。（局部方法）返回第十六页，讲稿共八十六页哦17171、全局方法全局方法 利用每个可利用的控制点来构建一个方程利用每个可利用的控制点来构建一个方程或模型，该模型可用于估算未知数值。

9、或模型，该模型可用于估算未知数值。（1）趋势面分析趋势面分析 趋势面分析用多项式方程拟合已知数值的趋势面分析用多项式方程拟合已知数值的点点。该方程又称趋势面模型，线性趋势面。该方程又称趋势面模型，线性趋势面方程：方程：zx,yb0+b1x+b2y 特征值特征值z是是x和和y的函数。系数的函数。系数b由控制点估算。其拟合由控制点估算。其拟合程度可用相关系数确定。每个已知点观测值和估算值之间程度可用相关系数确定。每个已知点观测值和估算值之间的偏差或残差可以计算出来。的偏差或残差可以计算出来。第十七页，讲稿共八十六页哦1818拟合更复杂的面要求更高次的趋势面模型，拟合更复杂的面要求更高次的趋势面模型

10、，如三次趋势面方程：如三次趋势面方程：zx,yb0+b1x+b2y+b3x2+b4xy+b5y2+b6x3+b7x2y+b8xy2+b9y3三次面包含山和谷。统计检验能用于比较三次面包含山和谷。统计检验能用于比较几种可选模型的拟合程度。几种可选模型的拟合程度。GIS软件包如软件包如ARC/INFO提供了多达提供了多达12次趋势面模型。图次趋势面模型。图13.3是由是由105个数据点构建的爱达荷州年平均降水量的三阶趋势面等个数据点构建的爱达荷州年平均降水量的三阶趋势面等值线图。值线图。第十八页，讲稿共八十六页哦1919第十九页，讲稿共八十六页哦2020(2)回归模型回归模型 回归模型把方程中的一

11、个因变量与多个自变量联系回归模型把方程中的一个因变量与多个自变量联系起来，而后可用于预测或估算。起来，而后可用于预测或估算。例例:用雪水量用雪水量(SWE)作为作为因变量，位置和地形变量作为自变量，其流域模型的表达形式如因变量，位置和地形变量作为自变量，其流域模型的表达形式如下：下：SWEb0+blEASTING+b2SOUTHING+b3ELEV+b4PLANl000 这里这里EASTING和和SOUTHING分别对应一个高程格网的行数和列分别对应一个高程格网的行数和列数，数，ELEV是高程值，是高程值，PLANl000是表面曲率的测度。是表面曲率的测度。第二十页，讲稿共八十六页哦21212

12、 局部方法局部方法该方法用一组控制点样本来估算未知值，因此控制点对该该方法用一组控制点样本来估算未知值，因此控制点对该方法十分重要。第一要点是用于估算的控制点个数。方法十分重要。第一要点是用于估算的控制点个数。GIS软件包通常允许用户自己确定控制点的个数或用默认值软件包通常允许用户自己确定控制点的个数或用默认值(712个个)。你也许认为控制点越多，估算结果越精确。然。你也许认为控制点越多，估算结果越精确。然而这种设想的正确与否取决于控制点的分布与要估算单元而这种设想的正确与否取决于控制点的分布与要估算单元的关系以及空间自相关程度。控制点越多通常意味着估算的关系以及空间自相关程度。控制点越多通常

13、意味着估算越通用化。越通用化。第二十一页，讲稿共八十六页哦2222 控制点个数确定之后，下一步就是控制点选择。简单控制点个数确定之后，下一步就是控制点选择。简单办法是用离估算点最近的点作为控制点，另一种办法是办法是用离估算点最近的点作为控制点，另一种办法是用半径来选择控制点，半径的大小必须根据控制点的分用半径来选择控制点，半径的大小必须根据控制点的分布来调整。例如，布来调整。例如，ARC/INFO和和ArcView都提供了这两都提供了这两种方法。其他方法考虑了方向因素，例如，象限法或八方种方法。其他方法考虑了方向因素，例如，象限法或八方向法。象限法是从围绕每个要估算单元的四个方向选择控向法。象

14、限法是从围绕每个要估算单元的四个方向选择控制点，八方向法是从八个方向来选择控制点。制点，八方向法是从八个方向来选择控制点。第二十二页，讲稿共八十六页哦2323（1）泰森多边形泰森多边形 泰森泰森(Thiessen)多边形围绕已知点样本构建多边形围绕已知点样本构建而成，使得在泰森多边形内的任意点与多而成，使得在泰森多边形内的任意点与多边形内的已知点更接近，而不是与其他已边形内的已知点更接近，而不是与其他已知点接近。泰森多边形最初用于估算区域知点接近。泰森多边形最初用于估算区域降水量的平均值。降水量的平均值。第二十三页，讲稿共八十六页哦2424泰森多边形要求在已知点之间构建初始三角形，亦即泰森多边

15、形要求在已知点之间构建初始三角形，亦即连接已知点形成三角形。因为连接点的方法不同，形连接已知点形成三角形。因为连接点的方法不同，形成的三角形也不同。与构建不规则三角网相似，德劳成的三角形也不同。与构建不规则三角网相似，德劳内三角测量常用于构建泰森多边形。德劳内三角网测内三角测量常用于构建泰森多边形。德劳内三角网测量确保了每个已知点都与它最近的点相接，这样就使量确保了每个已知点都与它最近的点相接，这样就使得三角形尽量接近等边。经过三角形每条边的中点画得三角形尽量接近等边。经过三角形每条边的中点画垂线，连接起来就可以很容易地构建泰森多边形垂线，连接起来就可以很容易地构建泰森多边形(图图13.4)。

16、泰森多边形也称为冯罗诺。泰森多边形也称为冯罗诺(Voronoi)多边形。多边形。第二十四页，讲稿共八十六页哦2525第二十五页，讲稿共八十六页哦2626(2)密度估算密度估算密度估算基于点的分布及其已知值来量测密度估算基于点的分布及其已知值来量测格网的密度格网的密度。一种简单密度估算方法是将。一种简单密度估算方法是将格网置于点分布图上，将落在每个单元的格网置于点分布图上，将落在每个单元的点列表，将点值加和，将单元的点值总和点列表，将点值加和，将单元的点值总和除以单元大小，就估算得每个单元的密度。除以单元大小，就估算得每个单元的密度。第二十六页，讲稿共八十六页哦2727图图13.5显示一个简单密

17、度估显示一个简单密度估算的例子。输入是以算的例子。输入是以50m间间距标绘的看到鹿的地点的分布。距标绘的看到鹿的地点的分布。每个鹿的地点有一个计数值，每个鹿的地点有一个计数值，表示在该地一只鹿被看到多少表示在该地一只鹿被看到多少次。输出是密度格网：单元大次。输出是密度格网：单元大小是小是l hm2，密度用每公顷看，密度用每公顷看到的次数来表示。到的次数来表示。第二十七页，讲稿共八十六页哦2828核估算核估算(kernel estimation)是一种不同的密是一种不同的密度估算方法，它用核函数联系每一个点或度估算方法，它用核函数联系每一个点或观测点。表达为双变量概率密度函数，核观测点。表达为双

18、变量概率密度函数，核函数看起来像是一个隆起函数看起来像是一个隆起(bump)，以一个，以一个点为中心在一个定义的带宽或窗口范围内点为中心在一个定义的带宽或窗口范围内逐渐减小到逐渐减小到0(图图13.6)。核函数和带宽决定了。核函数和带宽决定了隆起的形状，其形状反过来决定了在估算隆起的形状，其形状反过来决定了在估算中的平滑量。在点中的平滑量。在点x上的核密度估算值是带上的核密度估算值是带宽范围内位于观测点宽范围内位于观测点xi上的隆起部分的总和。上的隆起部分的总和。第二十八页，讲稿共八十六页哦2929(x-xi)和和(y-yi)是带宽范围内的点是带宽范围内的点x和观测点和观测点xi之间的之间的x

19、、y坐标的偏差。坐标的偏差。第二十九页，讲稿共八十六页哦3030使用与简单估算方法相同使用与简单估算方法相同的输入数据，图的输入数据，图13.7显示显示核估算法构建的密度格网输核估算法构建的密度格网输出结果。格网上的密度值是出结果。格网上的密度值是期望值期望值(输出值的平均数输出值的平均数)而不而不是概率。虽然图是概率。虽然图13.7和图和图13.5间的差别不明显，但间的差别不明显，但核估算法通常能产生比简核估算法通常能产生比简单估算法更平滑的密度面。单估算法更平滑的密度面。第三十页，讲稿共八十六页哦3131(3)反距离权重插值反距离权重插值该方法该方法假设未知值的点受较近控制点的影假设未知值

20、的点受较近控制点的影响比较远控制点的影响更大响比较远控制点的影响更大。影响的程度影响的程度(或权重或权重)用点之间距离乘方的倒数表示用点之间距离乘方的倒数表示。乘。乘方为方为1.0意味着点之间数值变化率为恒定，意味着点之间数值变化率为恒定，该方法称为线性插值法。乘方为该方法称为线性插值法。乘方为2.0或更高或更高则意味着越靠近已知点，数值的变化率越则意味着越靠近已知点，数值的变化率越大，远离已知点趋于平稳。大，远离已知点趋于平稳。第三十一页，讲稿共八十六页哦3232式中：式中：z0是点是点0的估计值；的估计值；zi是控制点是控制点i的的z值；值；di是控制点是控制点i与点与点0间的距离；间的距

21、离；s是在估算中用到的控是在估算中用到的控制点的数目；制点的数目；k是指定的幂。是指定的幂。第三十二页，讲稿共八十六页哦3333图图13.8显示对爱达荷州显示对爱达荷州105个气象站点数据由反距离平个气象站点数据由反距离平方方法建立的年平均降水量曲面。图方方法建立的年平均降水量曲面。图13.9显示曲面的显示曲面的等雨量线图。等雨量线图。第三十三页，讲稿共八十六页哦3434(4)薄板样条函数法薄板样条函数法薄板样条函数薄板样条函数(Thin-plate splines)建立一个建立一个通过控制点的面，并使所有点的坡度变化通过控制点的面，并使所有点的坡度变化最小。换言之，最小。换言之，薄板样条函数

22、以最小曲率薄板样条函数以最小曲率面拟合控制点面拟合控制点。薄板样条函数的估计值由。薄板样条函数的估计值由下式计算：下式计算：第三十四页，讲稿共八十六页哦3535Q(x,y)Aidi2logdi+a+bx+cy 式中：式中：x和和y是要被插值的点的是要被插值的点的x、y坐标；坐标；di2(x-xi)2+(y-yi)2；xi和和yi是控制点是控制点i的的 x、y坐标。薄板样条函坐标。薄板样条函数包括二个部分：数包括二个部分：(a+bx+cy)表示局部趋势函数，表示局部趋势函数，di2logdi表示基本函数，可获得最小曲率的面。系数表示基本函数，可获得最小曲率的面。系数Ai、a、b和和c由以下线性方

23、程组决定：由以下线性方程组决定：式中：n为控制点的数目，fi为控制点i的已知值，系数的计算要求n+3联立方程.第三十五页，讲稿共八十六页哦3636 薄板样条函数的一个主要问题是在数据薄板样条函数的一个主要问题是在数据贫乏地区的坡度较大，常出现过伸贫乏地区的坡度较大，常出现过伸(overshoots)。用于校正过伸的方法包括薄。用于校正过伸的方法包括薄板张力样条板张力样条(thin-piate splines with tension)，规则样条，规则样条(regularized splines)和规则张力和规则张力样条样条(regularized splines with tension)等。

24、等。第三十六页，讲稿共八十六页哦3737规则样条函数的近似值规则样条函数的近似值与薄板样条函数有相同与薄板样条函数有相同的的局部趋势函数局部趋势函数，但是，但是基本函数基本函数取不同形式取不同形式：n n式中式中式中式中 是样条法中要用到的权重；是样条法中要用到的权重；是样条法中要用到的权重；是样条法中要用到的权重；d d是待定值的点和控制点是待定值的点和控制点是待定值的点和控制点是待定值的点和控制点i i间的距间的距间的距间的距离，离，离，离，c c是常数是常数是常数是常数0.5772150.577215；KK0 0(d/(d/)是修正的零次贝塞耳是修正的零次贝塞耳是修正的零次贝塞耳是修正的

25、零次贝塞耳(Besscl)(Besscl)函数函数函数函数它可由一个多项式方程估计。它可由一个多项式方程估计。它可由一个多项式方程估计。它可由一个多项式方程估计。值通常设为值通常设为值通常设为值通常设为0 00.50.5之间因为更大的之间因为更大的之间因为更大的之间因为更大的 值会导致数据贫乏地区趋于过伸。值会导致数据贫乏地区趋于过伸。值会导致数据贫乏地区趋于过伸。值会导致数据贫乏地区趋于过伸。ARC/INFOARC/INFO和和和和ArcViewArcView采用的采用的采用的采用的默认默认默认默认 值为值为值为值为0.10.1Q(x,y)Aidi2logdi+a+bx+cy第三十七页，讲稿

26、共八十六页哦3838薄板张力样条薄板张力样条(thin-plate splines with tension)法表达式：法表达式：式中：式中：式中：式中：是本张力法要用到的权重。如果是本张力法要用到的权重。如果是本张力法要用到的权重。如果是本张力法要用到的权重。如果 的权重被的权重被的权重被的权重被设为接近于设为接近于设为接近于设为接近于0 0，则用张力法与基本薄板样条法得到的估计，则用张力法与基本薄板样条法得到的估计，则用张力法与基本薄板样条法得到的估计，则用张力法与基本薄板样条法得到的估计值相似。值相似。值相似。值相似。ARC/INFOARC/INFO和和和和ArcViewArcView采

27、用的默认采用的默认采用的默认采用的默认 值为值为值为值为0.10.1。第三十八页，讲稿共八十六页哦3939 薄板样条函数及其变种适用于平滑和连续的面。样条薄板样条函数及其变种适用于平滑和连续的面。样条法也适用于气候数据的插值。图法也适用于气候数据的插值。图13.10，13.11分别显示分别显示用规则样条函数和张力样条函数法建立的年均降水量面。用规则样条函数和张力样条函数法建立的年均降水量面。两幅地图中的等雨量线都很平滑，比图两幅地图中的等雨量线都很平滑，比图13.9的等雨量线的等雨量线(用反距离平方插值法建立的用反距离平方插值法建立的)更为平滑。在数据贫乏地区更为平滑。在数据贫乏地区用样条法插

28、值的结果可能会比实际情况偏大，如图用样条法插值的结果可能会比实际情况偏大，如图13.10的北爱达荷的北爱达荷 45in等雨量线所例证的。等雨量线所例证的。第三十九页，讲稿共八十六页哦4040第四十页，讲稿共八十六页哦4141(5)克里金法克里金法克里金法克里金法(Kriging)是由南非采矿工程师是由南非采矿工程师D.G.Krigee创立的一种用于空间插值的地理统计创立的一种用于空间插值的地理统计方法。克里金法的原理是假设某种属性的空方法。克里金法的原理是假设某种属性的空间变化既不是完全随机也不是完全确定。空间变化既不是完全随机也不是完全确定。空间变化可能包括三种影响因素：空间相关因间变化可能

29、包括三种影响因素：空间相关因素，代表区域变量的变化；偏移或结构，代素，代表区域变量的变化；偏移或结构，代表趋势；还有随机误差。偏移出现与否和对表趋势；还有随机误差。偏移出现与否和对区域变量的解释导致了不同克里金法的出现。区域变量的解释导致了不同克里金法的出现。第四十一页，讲稿共八十六页哦4242 1.普通克里金法普通克里金法 假设不存在偏移，普通克里金法关注于空假设不存在偏移，普通克里金法关注于空间相关因素。衡量所选已知点之间空间相间相关因素。衡量所选已知点之间空间相关程度的测度是半方差关程度的测度是半方差(semivariance)，由，由下式计算：下式计算：h是已知点之间的距离，是已知点之

30、间的距离，n样本点的数量，样本点的数量，z是属性是属性值。值。第四十二页，讲稿共八十六页哦4343在不同距离的半方差值算出后，绘成半方差图，在不同距离的半方差值算出后，绘成半方差图，y轴代表偏差，轴代表偏差，x轴代表已知点之间的距离。半轴代表已知点之间的距离。半方差可分成三部分：熔核方差可分成三部分：熔核(nugget)、值域、值域(range)和基台和基台(sill)。熔核是在距离为。熔核是在距离为0处的半处的半方差，代表无关的空间噪音。值域是半方差的方差，代表无关的空间噪音。值域是半方差的空间相关部分，它显示半方差随着距离递增。空间相关部分，它显示半方差随着距离递增。超过值域范围，半方差趋

31、平于相对恒定值。达超过值域范围，半方差趋平于相对恒定值。达到恒定的半方差称为基台到恒定的半方差称为基台(sill)。第四十三页，讲稿共八十六页哦4444第四十四页，讲稿共八十六页哦4545半方差图半方差图(semivariogram)将半方差与距离关联起来。它将半方差与距离关联起来。它可单独用作空间相关的测度。但在克里金法中用作插值可单独用作空间相关的测度。但在克里金法中用作插值器器(interpolator)，半方差图必须与数学函数或模型拟合。例，半方差图必须与数学函数或模型拟合。例如，如，ARC/INFO提供了五种模型：高斯、线性、球面、圆提供了五种模型：高斯、线性、球面、圆形和指数模型。

32、拟合的半方差图便可用于估算任何给定距形和指数模型。拟合的半方差图便可用于估算任何给定距离的半方差。离的半方差。第四十五页，讲稿共八十六页哦4646普通克里金法在空间插值中直接使用拟合普通克里金法在空间插值中直接使用拟合半方差图。估算某点的半方差图。估算某点的z值的通用方程是：值的通用方程是：式中式中:z0是估计值，是估计值，zx是已知点的值；是已知点的值；Wx是是与每个已知点关联的权重；与每个已知点关联的权重；s是用于估算的是用于估算的已知点的数目。权重可由对一组联立方程已知点的数目。权重可由对一组联立方程的求解得到，例如，下列联立方程是由三的求解得到，例如，下列联立方程是由三个已知点估算一个

33、未知点的值所必须的：个已知点估算一个未知点的值所必须的：第四十六页，讲稿共八十六页哦4747 式中：式中：(hij)是已知点是已知点i和和j间的半方差，间的半方差，(hi0)是是已知点和未知点之间的半方差，已知点和未知点之间的半方差，是拉格朗日系是拉格朗日系数，它的加入是为了确保把估算误差降到最小。数，它的加入是为了确保把估算误差降到最小。上面的方程可改写成矩阵形式：上面的方程可改写成矩阵形式：第四十七页，讲稿共八十六页哦4848将左边的逆矩阵乘以右边的矩阵，可解将左边的逆矩阵乘以右边的矩阵，可解得权重向量。一旦知道了权重，方程得权重向量。一旦知道了权重，方程(13.17)便可用估算便可用估算

34、z0：z0z1W1+z2W2+z3W3第四十八页，讲稿共八十六页哦4949上例表明在上例表明在克里金法中用的权重不仅包括估算点和已知克里金法中用的权重不仅包括估算点和已知点之间的半方差，而且包括已知点之间的半方差点之间的半方差，而且包括已知点之间的半方差。这与。这与反距离权重插值法不同，后者只用适用于已知点和估反距离权重插值法不同，后者只用适用于已知点和估算点的权重。克里金法和其他局部方法之间的另一个算点的权重。克里金法和其他局部方法之间的另一个重要区别是：重要区别是：克里金法对每个估算点产生一个方差测克里金法对每个估算点产生一个方差测度来衡量估算值的可靠性度来衡量估算值的可靠性。上例中，方差

35、可由下式计。上例中，方差可由下式计算：算：第四十九页，讲稿共八十六页哦5050图图13.14显示由线性模型普通克里金法生成的年平均降水量曲面。显示由线性模型普通克里金法生成的年平均降水量曲面。图图 13.15显示克里金曲面的标准差分布。数据贫乏地区的标准差显示克里金曲面的标准差分布。数据贫乏地区的标准差最高，例如，爱达荷州的西南角。最高，例如，爱达荷州的西南角。第五十页，讲稿共八十六页哦51512通用克里金法通用克里金法通用克里金法通用克里金法(Universal Kriging)假设除了假设除了已知点之间的空间关系外，空间变量在已知点之间的空间关系外，空间变量在z值值上还有偏移或有结构因素。

36、一般地说，通上还有偏移或有结构因素。一般地说，通用克里金法具体表现为一个趋势，例如，用克里金法具体表现为一个趋势，例如，在克里金过程中的一个趋势面方程。在克里金过程中的一个趋势面方程。ARCINFO提供两种通用克里金法。提供两种通用克里金法。Universal 1使用一个平面，定义为如下一阶多项式：使用一个平面，定义为如下一阶多项式：第五十一页，讲稿共八十六页哦5252Mb1xi+b2yi式中：式中：M是偏移；是偏移；xi和和yi是已知点是已知点i的的x和和y坐坐标；标；b1和和b2是要估计的偏移系数。是要估计的偏移系数。Universal 2用二阶曲面，它被定义为一个二用二阶曲面，它被定义为

37、一个二阶多项式：阶多项式：M=b1xi+b2yi+b3x2+b4xiyi+b5y2系数系数bi必须按权重估算，这意味着通用克里必须按权重估算，这意味着通用克里金法比普通克里金法要求更多的联立方程金法比普通克里金法要求更多的联立方程用于估算未知值。用于估算未知值。第五十二页，讲稿共八十六页哦5353 图图13.17显示用通用克里金显示用通用克里金法与线性偏移建立的年均法与线性偏移建立的年均降水量曲面，图降水量曲面，图13.18显示显示在克里金曲面上标准差的在克里金曲面上标准差的分布。两幅地图与基于线分布。两幅地图与基于线性模型克里金法建立的地性模型克里金法建立的地图相似。主要的不同是在图相似。主

38、要的不同是在爱达荷州西南角出现了负爱达荷州西南角出现了负值，这是因为通用克里金值，这是因为通用克里金模型中线性偏移而导致插模型中线性偏移而导致插值不准的结果。值不准的结果。第五十三页，讲稿共八十六页哦5454 3其他克里金法其他克里金法除了普通克里金法和通用克里金法，还有其他的克里金法，除了普通克里金法和通用克里金法，还有其他的克里金法，例如例如:克里金块内插法克里金块内插法(block kriging)和联合克里金法和联合克里金法(co-kriging)。克里金块内插法估算某个小范围内或块内。克里金块内插法估算某个小范围内或块内的一个变量的平均值，联合克里金法在插值时增加的一个变量的平均值，

39、联合克里金法在插值时增加次要变量，这些次要变量与主要变量有相关关系。次要变量，这些次要变量与主要变量有相关关系。这种方法的原理是：变量之间的相关关系能用于提这种方法的原理是：变量之间的相关关系能用于提高主要变量的预测值精度。例如，用联合克里金法，高主要变量的预测值精度。例如，用联合克里金法，在降水量插值中，将高程作为一个附加变量会获得在降水量插值中，将高程作为一个附加变量会获得更好结果更好结果。第五十四页，讲稿共八十六页哦5555(6)各种局部方法的比较各种局部方法的比较用相同数据集，不同的局部方法会产生不同的结果。用相同数据集，不同的局部方法会产生不同的结果。图图13.19显示规则样条函数和

40、距离倒转平方法两种方显示规则样条函数和距离倒转平方法两种方法得出不同的年平均降水量估算结果。图法得出不同的年平均降水量估算结果。图13.20显示显示用线性模型普通克里金法和规则样条函数得出的不同用线性模型普通克里金法和规则样条函数得出的不同结果。虽然普通克里金法与规则样条函数两种方法间结果。虽然普通克里金法与规则样条函数两种方法间的差异比规则样条函数与倒转距离平方法间的差异更的差异比规则样条函数与倒转距离平方法间的差异更小，但这两幅地图有相似的总体模式。另外，正负误小，但这两幅地图有相似的总体模式。另外，正负误差多于差多于3in的，都位于数据贫乏地区。这说明不管用哪的，都位于数据贫乏地区。这说

41、明不管用哪种方法插值，都不能代替观测数据。种方法插值，都不能代替观测数据。第五十五页，讲稿共八十六页哦5656第五十六页，讲稿共八十六页哦5757 将估算值的差别绘成地图并不能确定哪种插值方法更将估算值的差别绘成地图并不能确定哪种插值方法更精确。交叉验证精确。交叉验证(cross-validation)分析通常用于确定不同分析通常用于确定不同插值方法的精度。对每一种插值方法，交叉验证分析重插值方法的精度。对每一种插值方法，交叉验证分析重复从已知数值集中删除一个已知点的过程，用剩下的已复从已知数值集中删除一个已知点的过程，用剩下的已知点估算被删除点的数值，并计算估算值与已知值的误知点估算被删除点

42、的数值，并计算估算值与已知值的误差。然后计算出误差统计值用于评价插值方法的精度。差。然后计算出误差统计值用于评价插值方法的精度。例如，交叉验证分析表明，在消除高程倾向克里金法例如，交叉验证分析表明，在消除高程倾向克里金法(elevation-detrended kriging)或联合克里金法或联合克里金法(co-kriging)中用高程作为中用高程作为x、y位置的附加变量，产生了比克里金位置的附加变量，产生了比克里金法好的结果。法好的结果。第五十七页，讲稿共八十六页哦5858 评定不同插值法精度的另一种方法是将已知点分成两部分样本：一部分样本用于给每种插值方法建模，另一部分样本用于检测各种方法

43、的精度。如果已知点的数目太少，则这种样本分离的方法不适用。第五十八页，讲稿共八十六页哦5959二、二、区域插值区域插值适用于离散空间数据。适用于离散空间数据。区域插值是从一种已知的分区数据中推断区域插值是从一种已知的分区数据中推断出同一地区的另一种分区数值的插值方法。出同一地区的另一种分区数值的插值方法。第五十九页，讲稿共八十六页哦6060 第六十页，讲稿共八十六页哦6161(1)叠置法叠置法认为源区和目标区内数据均匀分布。求解时，将认为源区和目标区内数据均匀分布。求解时，将目标区叠置到源区上，求出源区和目标区之间的目标区叠置到源区上，求出源区和目标区之间的交集。目标区计算公式：交集。目标区计

44、算公式：n UDi=USj Dj j=1 其中：其中：UDi为目标区第为目标区第i分区的值，分区的值，i=1,2.m；USj为源区第为源区第j个分区的值，个分区的值，j=1,2,.n；Dj 为第为第i个目标分区同有第个目标分区同有第j个源区的面积比。个源区的面积比。第六十一页，讲稿共八十六页哦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返回第六十二页，讲稿共八十六页哦6363(2)比重法比重法将源区的统计数据从同质性变成非同质性，将源区的统计数据从同质性变成非同质性，再

45、进行区域插值。再进行区域插值。1 叠置格网，确定精度。叠置格网，确定精度。2 将源区各分区的平均数赋将源区各分区的平均数赋予各格网点（图予各格网点（图a)。第六十三页，讲稿共八十六页哦64643 计算相邻的四计算相邻的四(八八)个格网点的平均值个格网点的平均值(图图b)。第六十四页，讲稿共八十六页哦65654 以分区为单位将以分区为单位将3的平均的平均值相加，得值为值相加，得值为Us，计算各分计算各分区系数区系数pi=Us/Us，Usi为源区为源区i分分区的人口数。并将各个格网点值区的人口数。并将各个格网点值乘以乘以pi，得调整后的格网值，得调整后的格网值(图图c)。第六十五页，讲稿共八十六页

46、哦66665不断循环不断循环2,3,4步，直到步，直到pi接近接近1(0.99)为止，为止，便便 可可计算目标区的内插值计算目标区的内插值(图图d)。例如，例如，Ud1区的内插值区的内插值=5.5+5.3+5.5+5.2+4.5=26返回第六十六页，讲稿共八十六页哦6767第五节第五节 矢量数据向栅格数据的转换矢量数据向栅格数据的转换先要确定分辨率。先要确定分辨率。矢量数据是直角坐标，原点在左下方；矢量数据是直角坐标，原点在左下方；栅格数据是行列坐标，原点在左上方。栅格数据是行列坐标，原点在左上方。第六十七页，讲稿共八十六页哦6868 其转换公式为：其转换公式为：X=(Xmax-Xmin)/J

47、 Y=(Ymax-Ymin)/I X，Y表示栅格边长。表示栅格边长。I、J为行、列数。为行、列数。第六十八页，讲稿共八十六页哦6969例如：已知某一地区例如：已知某一地区x方向为方向为15km，y方方向为向为30km，现要把该地区的地块图转换，现要把该地区的地块图转换成栅格数据。要求栅格的分辨率为成栅格数据。要求栅格的分辨率为30mx30m，则：，则：行数行数I=30km/30m=1000 列数列数J=15km/30m=500第六十九页，讲稿共八十六页哦7070一、边界代数法一、边界代数法转换的关键是对矢量表示的多边形边界转换的关键是对矢量表示的多边形边界内的所有栅格赋予多边形编码，形成栅内的

48、所有栅格赋予多边形编码，形成栅格数据阵列。格数据阵列。面的填充是根据边界的拓扑信息，通过面的填充是根据边界的拓扑信息，通过加减运算将边界位置信息动态地赋予各加减运算将边界位置信息动态地赋予各栅格。栅格。实现的前提是知道拓扑关系。实现的前提是知道拓扑关系。第七十页，讲稿共八十六页哦7171 填充过程：沿边界前进方向填充过程：沿边界前进方向y值上升时为值上升时为上行，上行，y值下降时为下行。上行时填充值值下降时为下行。上行时填充值为左多边形号减右多边形号，下行时填为左多边形号减右多边形号，下行时填充值为右多边形减左多边形号，将每次充值为右多边形减左多边形号，将每次填充值同该处的原始值作代数运算得到

49、填充值同该处的原始值作代数运算得到最终填充属性值。最终填充属性值。第七十一页，讲稿共八十六页哦7272N1N2弧上行，左弧上行，左-右为右为0-1=-1。弧左边栅。弧左边栅格值为格值为-1第七十二页，讲稿共八十六页哦7373N2N3弧下行，右弧下行，右-左为左为2-0=2，弧左边栅格，弧左边栅格值加值加2第七十三页，讲稿共八十六页哦7474N3N1弧下行，右弧下行，右-左左=3-0=3。弧左边栅格。弧左边栅格加加3第七十四页，讲稿共八十六页哦7575N1N4弧上行，左弧上行，左-右右=1-3=-2。弧左边加。弧左边加-2。第七十五页，讲稿共八十六页哦7676N4N2弧上行，左弧上行，左-右右=

50、1-2=-1，弧左边加，弧左边加-1第七十六页，讲稿共八十六页哦7777N4N3弧下行，右弧下行，右-左左=3-2=1，弧左边加，弧左边加1。返回第七十七页，讲稿共八十六页哦7878第六节第六节 栅格数据向矢量数据的转换栅格数据向矢量数据的转换转换的实质是将相同属性代码的栅格集合转转换的实质是将相同属性代码的栅格集合转变成由边界弧段以及拓扑关系。变成由边界弧段以及拓扑关系。1.栅格数据矢量化的过程栅格数据矢量化的过程(1)前处理。变成二值图形前处理。变成二值图形(0和和1)，线的细化。，线的细化。(2)矢量化。转变成直角坐标数据。矢量化。转变成直角坐标数据。(3)后处理。除去多余点；光滑处理（