地理信息系统应用 空间数据插值处理 GIS数据空间插值处理--教案.doc

教 案 模块三 空间数据处理 知识点 空间数据的内插方法 空间数据的内插可以作如下简单的描述：设已知一组空间数据，它们可以是离散点的形式，也可以是分区数据的形式，现在要从这些离散数据中找到一个函数关系式。使该关系式最好地逼近这些己知的空间数据，并能根据该函数关系式推求出区域范围内其他任意点或任意分区的值。这种通过己知点或分区的数据，推求任意点或分区数据的方法就称为空间数据的内插。根据已知点和己知分区数据的不同，将空间数据内插分为点的内插和区域的内插。它们是地理信息系统数据处理常用的方法之一，广泛用于等值线自动制图、数字高程模型的建立、不同区域界线现象的相关分析和比较研究等。 一、点的内插 点的内插是研究具有连续变化特征现象(例如地面高程等)的数值内插方法。根据内插精度的不同，可以将点内插法分为精确和概略的两类(图3.93)。但是，不论采用何种内插方法，建立按一定网络(例如三角网、格网等)形式排列的地面点高程数据，一般需要经过数据取样、数据处理和数据记录三个过程。数据取样是指数据点的选取和坐标的确定；数据处理是以数据点作为控制基础，用某一数学模型来模拟地表面，进行内插加密计算，确定三角网或格网节点处的特征值；数据记录是将建立的节点特征值记录于存储器内，以供分析应用. （一）数据取样 建立按一定格网形式排列的地面点高程，一般先从现有地图的等高线上进行数据取样。取样点可以按地性线(山脊线、山谷线、坡度变换线)、沿等高线或沿断面线布设。将数据点选择在地性线坡度改变处，或沿等高线在方向改变的地点，即根据地形变化的点。这样，数据点落在地形特征点上，能很好地控制地表面形态(图3.94)，这种方法称为随机取样方案。该方案任意两点间直线部分相对于地面高程上的最大差距，应按系统的功能、地面形态待征和对地形分析的精度要求来确定。以保证具有足够的地面信息，作为数据处理的控制依据。 点内插法 点内插法 精确的 粗略的 按距离加权内插法 多项式内插法 样条函数内插法 有限差分法 克里格内插法 指数级数趋势面内插法 傅里叶级数内插法 最小二乘样条函数内插法 按距离加权最小二乘内插法 图3.93 点内插法的分类 当按一定间隔的格网取样时，取一个二次曲线来代表地表面相应部分的曲线，而用间隔中点的线性内插值与二次曲线内插值所得到的断面高程的差异来判断格网密度是否恰当。如图3.95所示，设拟定的取样点1、2和3，相应的高程为h1，h2和h3，取等距高间隔为单位数，则间隔中点的高程为： 线性内插值： （4-14） 二次线性内插值 （4-15） 其高程差异 （4-16） 图3.94 沿等高线和地性线布设取样点 此值应在数字地面模型精度要求的限差之内，否则要缩短相同间距：对地形变化显著的地方，要辅以地形特征点为数据点，以提高数据点控制地面形态特征的能力。 12 1 2 3 二次曲线 图3.95 高程内插间距的确定方法 图3.96移动拟合方法 对采集的离散高程点，在插值运算之前，须进行认真的检查。检查其点位与记录的高程值是否匹配一致，点位是否有偏离现象，特征点高程是否有遗漏，有否重复记录的情况等等。检查工作可以在图形终端上进行，与原因进行核对，也可以通过打表输出检查，直至获得净化的取样数据文件，便可以转入插值运算。 （二）数据内插 由于取样的数据点呈离散分布形式，或者数据点虽按格网排列，但格网的密度不能满足使用的要求，这样就需要以数据点为基础进行插值运算。插值运算的手段是选择一个合理的数学模型，利用已知点上的信息求出函数的待定系数。由于地表面形态千变万化，既无规律性又无重复性，取用一个低次多项式来拟合整个地表面形态是不切合实际的。若采用高次多项式模拟地表面，又会出现函数的不稳定。因此，通常采用局部函数内插和逐点内插两种方式。 (1)逐点内插法以插值点位中心，定义一个局部函数去拟合周围的数据点，数据点的范围随插值点位置的变化而变化。因此又成移动曲面法。 = 1 \* GB3 ①移动拟合法。移动拟合法是典型的逐点内插法。对每一个待定点取用一个多项式曲面拟合该点附近的地表面，此时取待定点作平面坐标的原点，并用待定点为圆心，以R为半径的圆内各数据点来定义函数的待定系数，如图3.96所示。 设取二次多项式来拟合，则待求点的高程可写成一般式为 （4-17） 将坐标原点平移到待定点处，即 （4-18） 代入得到移动拟合法二次多项式插值公式