GIS+中格网DEM+生成等高线算法的一种新改进.pdfVIP

GIS 中格网DEM 生成等高线算法的一种新改进 王欢 魏海平 解放军信息工程大学测绘学院三系 摘要：等高线和数字高程模型(DEM)都是用来表示地面的高程信息的。目前，国家高程库的 建立是以等高线作为存储格式的，因此，将从格网DEM 中快速提取精度高的等高线数据是十分有 必要的。本文在格网 DEM 生成等高线的传统算法的基础上，根据其所产生的搜索速度与精度等问 题，引入状态矩阵的概念。通过判断单位格网的每条棱与要追踪等高线的相交情况获得该格网 DEM 的状态矩阵，为了提高精度，对状态矩阵进行修正，然后根据修正后的状态矩阵的值判断等 高线追踪的起点并依据一定的法则进行追踪。对于生成的等高线存在的问题，提出了一系列的解 决方法。 关键字：数字高程模型 等高线 状态矩阵 格网 DEM 1 引言 等高线和数字高程模型(DEM)都是用来表示地面的高程信息的。前者是地面高程相等点的连 线在水平面上的投影，用以来表达地面的起伏形态。后者是以数字的形式来表示某区域地面的高 (x , y ) z 程信息，它由该区域的一系列地面点坐标 及其相应的高程值 所组成的。目前，国家高程 库的建立是以等高线作为存储格式的，因此，将从 DEM 中快速提取精度高的等高线数据是十分 有必要的。 从 DEM 的结构上可以分为规则格网模型和不规则三角网模型（TIN ）两种类型。它们各有 其优缺点，规则格网数据结构简单，易于构网，缺点是当地面较平坦时，数据冗余大。TIN 能够 充分地表达地形的结构特征，具有可变的分辨率,数据冗余量小，但它导致了数据存储和操作的 复杂性。在由DEM 提取等高线过程中，由于TIN 的定位点具有不确定性，从而确定了从中插值 的等高点的不确定性，因此无法很好的把等高线追踪出来。所以提取等高线一般是通过是规则格 网DEM 来进行的。 规则格网 DEM 提取等高线的一般流程如下：对于要提取的高程为 H 的等高线，首先根据各 个格网点上的高程值确定等高线与之相交的棱，然后插值计算出等高线上各坐标点的坐标，依据 一定的法则将这些坐标点连起来组成等高线，最后使用光滑算法使等高线平滑。这种算法搜索速 度慢，并且精度不高。因此本文在传统算法的基础上，提出了一种改进方法，以增加搜索速度， 提高精度。 2 规则格网 DEM 提取等高线算法的改进 对于格网 DEM 提取等高线的一般算法所存在的问题，本文引入状态矩阵来加以解决。状态 矩阵分为横向和纵向两种类型，它们用来表示格网 DEM 中各条棱与等高线的相互关系。通过对 状态矩阵的修正，我们可以提高等高线的搜索速度和精度。 1 2.1 创建状态矩阵 提取等高线时，要获得等高线所经过的坐标点，该坐标点是格网 DEM 中格网的棱与等高线 的交点来确定的。状态矩阵就是用于判断格网 DEM 中各条格网的棱与某高程的等高线是否有交 点的矩阵。它分为横向状态矩阵和纵向状态矩阵，分别用来表示格网的水平方向与竖直方向棱与 等高线的相交情况。若相交，则矩阵相应的值为 1，否则，值为 0 。因此，状态矩阵形式如下： ⎡ h(0,0) h(0,1) Κ h(0, nCol −2) ⎤ ⎢ ⎥ StatusH ⎢ h(1,0) h(1,1) Κ h(1, nCol −2) ⎥ ⎢ Μ Μ Μ Μ ⎥ ⎢ ⎥ ⎣h(nRow −1,0) h(nRow −1,1) Κ h(nRow −1,