基于ArcGIS的泰森多边形法计算区域平均雨量.docxVIP

ArcGIS计算区域平均雨量王玉德（辽宁省水文水资源勘测局朝阳分局， 辽宁 朝阳 122000）［摘要］ 区域平均降雨量是防汛部门进行洪涝 ArcGIS 计算区域平均雨量 王玉德 （辽宁省水文水资源勘测局朝阳分局， 辽宁 朝阳 122000） ［摘要］ 区域平均降雨量是防汛部门进行洪涝预报与评估的重要依据， 实时、 准确的区域平均降雨量计算结果， 对防汛预报的准确性起着至关重要的作用。 在众多区域平均降雨量计算方法中， 泰森多边形法以适合区域内雨 量站及降雨量分布不均匀的情况而被广泛应用。 利用基于 ArcGIS 平台的泰森多边形法计算区域平均雨量， 解 决了传统方法计算泰森多边形面积以及面积权重系数的难题， 提高了区域平均降雨量的时效性和准确性。 ［关键词］ 区域平均雨量； 泰森多边形； ArcGIS ［中图分类号］ TV 125 ［文献标识码］ B 时使得计算工作量和工作难度加大， 很难满足时 引言 1 效性准确性的要求。 为解决这一难题，笔者在多年 防汛、 水文预报工作中总结出一种较传统的方法 更高效、更准确、更易操作的基于 ArcGIS 平台的新 方法。 区域平均降雨量是防汛部门进行洪涝预报与 评估的重要依据，实时、准确的区域平均降雨量计 算结果， 对防汛预报的准确性起着至关重要的作 用。 通常雨量站所观测到的降雨量，只能代表该雨 量站周围较小范围的降雨情况， 不能代表全流域 或某一区域的平均降雨量，因此，单独某个雨量站 泰森多边形法的原理 2 泰森多边形是荷兰气象学家 A·H·Thiessen 提 的降雨数据不能作为洪涝预报与评估的依据。 正 出的， 它是一种通过离散雨量点的雨量计算区域 因如此， 我们需要采用全流域或某一区域的所有 雨量站的降雨数据来计算区域平均降雨量。 平均雨量的方法。 即将相邻的离散雨量点连成三 角形，再做出三角形各边的垂线，于是每个离散点 目前计算区域平均降雨量的方法很多， 常用 周围被垂线围成一个多边形 （这个多边形就叫泰 的有算术平均法、数值法、等值线法、泰森多边形 法等。 在这些方法中，泰森多边形法最适合流域或 区域内雨量站或降雨量分布不均匀的情况， 是我 国计算区域平均雨量常用的方法 ，被水利、气象、 环境等部门广泛应用。 然而，在传统的泰森多边形法中，确定多边形 面积和面积权重系数非常困难， 以往的方法是利 用求积仪或者高级语言编程来获得多边形面积 。 森多边形），用这个泰森多边形中的唯一雨量观测 的降雨数据来代表这个多边形范围内的降雨量 。 按照为此做法， 在全流域或某一区域内就构成了 一个泰森多边形网。 在这个泰森多边形网内，每一 个泰森多边形都对应一个雨量观测点， 这个雨量 观测点的降雨数据就代表对应泰森多边形内的降 雨数据。 于是，通过泰森多边形网内所有雨量观测 点的降雨数据就可以算出区域平均降雨量。 设每个雨量观测点的降雨量为 xi，每个对应的 泰森多边形的面积为 fi，则区域平均降雨量可按下 式求得[1]： 当雨量站增加或减少时， 各站的权重系数还得重 新计算，可能还会出现雨量站迟报 、漏报、错报等 情况，使得雨量的权重发生变化，当以上情况出现 ［收稿日期］ 2014-04-24 ［作者简介］ 王玉德（1978-），男，辽宁朝阳人，工程师，现从事水文软件开发、防汛抗旱、水文水资源工作。 吉林水利基于 ArcGIS 的泰森多边形法计算区域平均雨量王玉德2014 年 06 月泰森多边形的原理和特点， 通过 ArcGIS 需要完成如下工作。3.1 ArcGIS 对三角网格的优化构建三角形网格是构建泰森多边形的先决条 件，如图 1（2）所示。 ArcGIS 在绘制三角形时，要通 过程序判断得出最优三角形。所谓最优三角形就是满足三角形各边最短的 三角形，即：从所有雨量站点任取一点作为第一个 nn 吉林水利 基于 ArcGIS 的泰森多边形法计算区域平均雨量 王玉德 2014 年 06 月 泰森多边形的原理和特点， 通过 ArcGIS 需要完成 如下工作。 3.1 ArcGIS 对三角网格的优化 构建三角形网格是构建泰森多边形的先决条 件，如图 1（2）所示。 ArcGIS 在绘制三角形时，要通 过程序判断得出最优三角形。 所谓最优三角形就是满足三角形各边最短的 三角形，即：从所有雨量站点任取一点作为第一个 n n x軃= f1 x1 +f2 x2 +…+fn xn = 1 Σf x =ΣA x i i i i f1 +f2 +…+fn F i = 1 i = 1 式中：xi 为雨量观测点的降雨量，单位为 mm；fi 为泰森多边形的面积，单位为 km2；n 为区域内雨量 观测点或泰森多边形的个数；F 为区域的总面积单 位为 km2；A