[发明专利]一种基于反地形汇流特征的喀斯特地貌单元提取方法

权利要求书

1.一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，包括如下步骤：

S1. 基于研究区DEM投影获取反地形DEM，计算地表汇点，基于距离阈值聚合汇点，小于距离阈值的地表汇点被分配给同一个喀斯特地貌单元；

S2.以属于同一个喀斯特地貌单元中的地表汇点作为倾泻点提取反地形流域边界，作为实际地形中扩展的喀斯特地貌单元边界；

S3.使用正地形DEM，基于汇流阈值提取汇流网络，基于汇流网络和坡度成本栅格，计算从汇流网络到区域内其他位置的坡度成本距离；

S4.基于坡度成本距离结果，设置坡度成本距离阈值，小于坡度成本距离阈值的区域为非山体区域，使用扩展的喀斯特地貌单元擦除非山体区域，得到喀斯特地貌单元。

2.根据权利要求1所述的一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，所述S1中研究区DEM投影时，采用双线性内插法进行投影栅格重采样。

3.根据权利要求1所述的一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，所述S1中反地形DEM的获取方法如下：使用一个大于研究区最大高程值的值减去研究区DEM。

4.根据权利要求1所述的一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，所述S1中地表汇点的计算方法如下：使用D8算法计算反地形流向，根据反地形流向计算研究区反地形地表汇点。

5.根据权利要求1所述的一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，所述S1中聚合汇点时，选取的距离阈值为60m。

6.根据权利要求1所述的一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，所述S3中汇流网络的提取方法如下：使用优先漫水算法计算正地形流向，根据正地形流向计算正地形汇流累计量，并将大于汇流阈值的区域提取为汇流网络。

7.根据权利要求1所述的一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，所述S3中汇流阈值为100。

8.根据权利要求1所述的一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，所述S3中坡度成本栅格的计算方法如下：使用二阶差分法计算正地形地表坡度，再加上一个极小值即可获得坡度成本栅格。

9.根据权利要求1所述的一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，所述S4中非山体区域的提取方法具体如下：基于坡度成本距离阈值将坡度成本距离分类为二值化栅格，小于坡度成本距离阈值的设置为1，大于等于坡度成本距离阈值的设置为NODATA；将上述二值化栅格转化为矢量面，删除二值化栅格中GRIDCODE值不为1的矢量面要素，得到非山体区域。

10.根据权利要求1所述的一种基于反地形汇流特征的喀斯特地貌单元提取方法，其特征在于，所述S4中坡度成本距离阈值为500。