DEM内插在软件开发中的应用.docVIP

数字高程模型内插在软件开发中的应用 摘 要 本文对数字高程模型（DEM）在地理信息、工程测绘软件开发中的数值内插、数据获取和数据分析应用等进行了介绍，最后对常用的计算机内插方法，如拉格朗日内插法、不规则三角网（TIN）内插、标准格网内插法进行了算法研究并编写计算机程序，最后根据工程实践，应用高程离散点建立数字高程模型对土方计算程序中的问题进行分析。 关键词 数字高程模型、内插方法、计算机编程 数字高程模型（DEM）的基本概念 数字高程模型（DEM）是指地面上按一定间距排列的规则格网上的高程数据集，它的数据标准包括头文件和数据集文件；它的主要应用是查询和分析地面高程、坡度、坡向和三维分析应用。数字高程模型一般是从野外实测地面特征点、特征线，按一定的内插方法计算出规则格网点的高程值。 从上世纪20年代起，随着计算机技术的发展，使得DEM内插的数值计算采用计算机编程求解，有些原本属于数学领域的内插方法借助计算机技术广泛应用工程实践中。 内插计算中常用的数值分析方法主要有按距离加权平均内插法、有限差分法、离散单元法和拉格朗日（Lagrange）插值等方法，这些方法都涉及较复杂的数学问题，应用于工程实践则需研究计算机算法，将这些数学模型用计算机能运行的程序进行解算，因此研究计算机算法是实现工程应用的必经之途径，因此本文除对内插的数值分析方法作了简单介绍以外，笔者还根据多年的工作实践，编写计算机程序对DEM方法进行土方计算和三维地面模型分析。 在求解变分问题中，物理量一般可表示为待定函数及其导数的积分式（如傅里叶函数和傅里叶积分），对其积分区域Ω，可仿照差分法的离散化方法（按一定方式划分网格并取结点），把Ω划分为有限个子区域Ωi，待定函数及其导数在子区域的结点上的值是待求的，这就是有限元法的基本理论，在各个Ωi内解析的插值函数（一般采用多项式插值，多项式系数一般由结点参数来确定，最常用的方法就是Lagrange插值法）来逼近。只要选取适当的插值幂次和结点参数的范围，一般就能满足对插值函数对精度和插值光滑性的要求，如后面还要讲到的标准格网内插采用的线性内插和五次幂内插法等均被广泛应用。总之，有限元法是以变分原理为基础吸取差分格式的思想而发展起来的一种比较好的数值解法，它把求解无限自由度的待定函数归结为求解有限单元的待定值问题，在数学上讲，具有按一定分布形式的结点及其一定类型的结点参数值的子区域Ωi称为单元。 常用的内插法及计算机编程在土工分析中的应用 2．1内插方法分析 不论采用何种方法在求解单元插值时，按距离加权平均内插法和拉格朗日内插法都是最常用和最简单的插值方法，也是目前计算机编程插值计算最常用的方法。 按距离平方加权平均内插法是根据已知的离散点在内插标准格网点的插值时，其计算公式为： ( 1 ) 式中：HP为内插点的内插值，Hi为各已知点的高程值，di为各已知点至内插点的距离，该距离是可以通过已知点至标准格网点的坐标计算出来。 拉格朗日内插法的基本原理是采用一个n次多项式（一元或二元）去拟合一个函数f(x)或f(x,y)，为计算方便和限制计算误差，一般n不大于5（阶数太大，除增加计算工作量外，有时还会出现不收敛情况和出现很大的计算误差），为保证内插的精度，必要时也采用分段插值，但由于离散点的特性，很容易在分段处其导数不连续，内插曲线不光滑，因此需要在插值结点的小区间内用低幂次插值，以保证内插曲线（或曲面）的光滑。 在数值内插中，最常用的方法是在现场采集数量有限的离散样本点，根据样本点求解拉格朗日内插多项式的系数，从而获得内插函数，求解任意点的解，在数学上则表示为：给定n个互异结点xi及对应的函数值yi = f(xi),(i = 1,2,···n)，写成多项式为： 对大多数问题是求解该多项式的系数ci，将n个离散点的（1）式写成矩阵方程： 求解的系数Ci的方程式为： 对于二元五次拉格朗日多项式来说（1）式可写为： 2．2 计算机编程内插算法分析 只要求得内插点至各已知的离散点距离，按（1）式就能很容易计算出内插点的Z值，但在实际工程中由于离散点太多，且离内插点太远的已知点，虽然其权倒数很小，但对于高程突变的点，其参与内插计算反而降低了内插精度，因此在实际编程时往往设定一个距离限值，在距离限值内的点才参与内插计算。 对于（4）式，要求解的系数Cki个数为： (个) 实际工程应用中，往往样本点（m个观测值）多于所要求解的变量Cn(mn），这样就需要对观测值z1、z2、···zm用最小二乘法进行高斯约化求解系数Cn： 首先求观测值与理论差的平方和： 再求σ对系数Ckj的偏导数，并令其