[发明专利]一种考虑山区平原地貌差异的子流域单元划分方法

权利要求书

1.一种考虑山区平原地貌差异的子流域单元划分方法，其特征在于：包括如下步骤，

S1、以DEM数据为基础，计算流域内各栅格流向和汇流累积量，并基于设定的河网阈值提取模拟河网；

S2、设定平原单元长度，对平原区及其周边栅格进行N*N大小聚合，并对每片N*N区域设置平原单元标记序号；

S3、从流域出口开始，以模拟河网栅格为主干，向上溯源遍历，根据遇到的栅格类型进行计算单元划分，并对划分后的各单元赋值一个唯一的子流域编号，该编号为从1开始的自然数；划分出的单元包括平原单元和山区单元；

S4、利用栅格流向信息，确定相邻平原单元和山区单元之间的流向；

S5、根据平原单元和山区单元之间的流向关系，从河流出口单元开始溯源遍历，生成子流域上下游拓扑关系表；

步骤S2具体包括如下内容，

S21、将平原单元长度设定为栅格长度的整数倍，以确保划分的平原单元和山区单元边界吻合；

S22、从栅格图层起始行列开始，逐行逐列按照步长进行累加遍历，每次对N行N列共N*N个相邻栅格进行检查，如果该区域内的平原栅格的个数大于或等于山区栅格的个数，则将该N*N区域标记为平原单元，并对区域内所有栅格设置相同的平原单元标记序号，所述序号为从1开始的自然数；若平原栅格的个数小于山区栅格的个数，则不做处理。

2.根据权利要求1所述的考虑山区平原地貌差异的子流域单元划分方法，其特征在于：栅格流向采用D8算法计算得出，即将该栅格流向周边8个栅格中坡度最陡的栅格的流向作为该栅格的流向；栅格的汇流累积量即为上游所有流入当前栅格的栅格数量之和。

3.根据权利要求1所述的考虑山区平原地貌差异的子流域单元划分方法，其特征在于：河网阈值用于定义各个栅格是属于河网栅格或坡面栅格，当栅格的汇流累积量大于河网阈值，则将该栅格定义为河网栅格，否则将该栅格定义为坡面栅格；所有的所述河网栅格构成模拟河网，所述河网阈值须保证提取的模拟河网在源头区与实际河网源头一致。

4.根据权利要求1所述的考虑山区平原地貌差异的子流域单元划分方法，其特征在于：所述平原单元为N*N栅格型单元，山区单元为根据汇流路径确定的子流域型单元，所述平原单元包括单元内有河道的平原河道单元和单元内无河道的平原坡面单元；所述山区单元包括单元内有河道的山区河道单元和单元内无河道的山区坡面单元。