[发明专利]基于六边形格网的多流向流域特征提取方法

说明书

本发明公开了一种基于六边形格网的多流向流域特征提取方法，该方法首先按双线性插法将规则格网DEM转换成球面六边形格网表示，分别计算得到六边形格网的顶点高程和中心点高程；然后依次连接顶点和中心点，将每个六边形格网划分成三角形，计算每个三角形中心点高程，根据三角形与邻近六边形格网高程差确定水流流向；根据与邻近六边形格网高程差比例，计算格网流量；最后遍历格网，跟踪流域路径。本发明利用六边形格网作为基础格网，解决了基于规则格网DEM流域特征提取中，因边邻近、角邻近格网距离不一致引起的地形表达失衡性问题；利用三角形与邻近六边形格网高程差比例，确定流向和流量，解决了原六边形单流向算法无法真实描述水流特征的问题。

技术领域

本发明涉及流域水文特征提取技术领域，特别涉及一种基于六边形格网的多流向流域特征提取方法。

背景技术

流域水文特征提取与分析是水文水资源的重要组成部分，为水文建模、地形地貌建模、水土侵蚀评估、水资源管理等提供重要的基础信息。平向汇流的转移路径通常以水流方向、径流节点、径流网络和流域边界等水文要素表达。其中，水流方向是流域水文特征提取的基础。

随着DEM数据获取能力的提升，基于DEM的流域水文特征提取得到广泛研究。其中基于规则格网的DEM因其数据结构简单、计算处理方便等特点成为流域水文特征提取的主要研究对象。基于规则格网DEM的水流方向算法主要的方法包括单流向算法、多流向算法和混合流向算法。

单流向算法假设每个格网有且仅有一个水流方向。典型的算法有D8算法：根据与8个邻近格网间的高程差确定水流方向，水流向高程差最大的格网。Rho8算法：引入一个均匀分布的随机变量计算坡度权重，在D8算法的基础上，增加了水流方向选择的随机性。

多流向算法假设每个格网有多个水流方向，即水流可以流向比它低的所有相邻格网。典型的算法有FD8算法：在D8算法基础上，根据高程差和比例分配流向和水流量。FRho8算法：在Rho8算法基础上添加多方向权重因子。多流向算法避免了水流偏向某格网，能较好的模拟处于坡面的水域特征。

混合流向算法将单流向和多流向算法结合，根据流域地形特征，分别采用不同的流向算法。在水流辐合区域，使用单流向算法，在坡面-水流发散区域，使用多流向算法。

基于六边形格网的算法：因为六边形格网只有边邻近格网，与邻近格网中心点距离相等，在处理最近邻域、移动路径等邻域处理问题中具有一定的优势。将基于规则格网的DEM数据转换至六边形格网，利用与邻近格网间的高程差，计算水流方向。

现有方法存在如下不足：

(1)正方形格网与其边邻近和角邻近格网间距离之比为流域水文特征提取中水流量的计算依赖于格网间距离，不一致的格网间距导致边邻近方向和角邻近方向上地形表达的失衡性，影响流域水文特征提取精度。

(2)基于六边形格网的地形特征提取研究中，水流方向基于单流向计算，而真实流域水流特征较复杂，在水流辐合区，单个格网水流全部流向高程更低格网，而在坡面，水流呈发散状。因此，单流向算法不能反映坡面的发散水流状况，无法真实描述水流特征。

发明内容

本发明针对现有技术的缺陷，提供了一种基于六边形格网的多流向流域特征提取方法。

为了实现以上发明目的，本发明采取的技术方案如下：

一种基于六边形格网的多流向流域特征提取方法，以球面六边形格网作为基础格网表达地形，分解每个六边形格网为三角形，根据与邻近六边形格网的高程差确定水流流向，根据高程差比例计算流量。

进一步地，基于六边形格网的多流向流域特征提取方法，包括以下步骤：

S1：基于六边形格网的DEM构建，包括以下子步骤：