一种数字河网提取时集水面积阈值的确定方法技术

本发明专利技术公开了一种数字河网提取时集水面积阈值的确定方法，利用流域遥感影像数据进行流域水体识别，栅格化水体识别结果；利用流域DEM数据计算出流域中每一个栅格单元的汇流累积值，生成汇流累积栅格；按汇流累积值的最小阈值提取流域中的最密可能河道；以最密可能河道上的栅格单元为中心，根据栅格化的水体识别结果，逐河道栅格计算一定范围内的河道可靠性；根据栅格化的河道可靠性，建立阈值与河道可靠性均值的关系曲线，确定临界阈值点，该临界阈值点即为集水面积阈值。本发明专利技术具有数据来源稳定可靠、计算效率高、结果客观合理等优点，有利于流域集水面积阈值的确定及数字河网的提取。

A Method for Determining the Threshold of Catchment Area in Digital River Network Extraction

【技术实现步骤摘要】
一种数字河网提取时集水面积阈值的确定方法
本专利技术涉及水文
，特别是一种数字河网提取时集水面积阈值的确定方法。
技术介绍
河网特征反映了流域的综合水文特征。随着遥感、地理信息以及数字流域等技术的发展，基于栅格数字高程模型(DEM，DigitalElevationModel)的分布式水文模型以其充分考虑降雨和下垫面条件空间变化的特点，现已成为流域水文模型的发展趋势。河网的准确提取为模型计算提供保障，为了进一步促进流域水文模型的发展，需要更深入研究集水面积阈值的确定方法。集水面积阈值是提取河网中的关键参数，集水面积小于阈值，则该栅格判定为坡地栅格，集水面积大于等于阈值则判定为河道栅格。阈值越小，提取的水系越密，河长越长；阈值越大，提取的水系越疏，河长越短。从技术角度讲，阈值可以设定为1至流域栅格个数之间的任意一个值，从理论角度讲，阈值只能是一个值，或者一个无限小的范围。阈值的选择主观性强，阈值的选取不同，往往会得到不同的水系结构，随意性很大。对集水面积阈值估算的主观性限制了分布式建模的应用，不利于国内分布式水文模型的发展。遥感是通过传感器接收来自陆地表面的电磁辐射，进而反演地表参数的过程。太阳辐射被反射、吸收和透射的比例随着地物类型和条件的不同而变化，不同地物形成的不同反射光谱反射率。通过遥感影像数据进行水体识别的方法较为成熟，有单波段法、比值法、谱间关系法、监督分类法等。理论上讲，遥感影像数据水体识别结果是现实水体的数据化反映，但是因为基于遥感影像研究实现水域面积提取的模型算法无法识别任一地物，也无法完美区分水体与其他地物，也由于遥感影像资料自身因云量影像产生误差，或者是流域本身存在小水体，水体识别结果在干流不连续，在支流表现为散乱的点，无法直接由遥感影像水体识别结果提取流域水系。针对以上不足，如何结合DEM与遥感影像数据，确定客观合理的集水面积阈值，正是专利技术人需要解决的问题。
技术实现思路
本专利技术所要解决的技术问题是克服现有技术的不足而提供一种数字河网提取时集水面积阈值的确定方法，本专利技术提高了确定流域集水面积阈值的精度与可靠性。本专利技术为解决上述技术问题采用以下技术方案：根据本专利技术提出的一种数字河网提取时集水面积阈值的确定方法，包括以下步骤：步骤1、利用流域遥感影像数据进行流域水体识别，栅格化水体识别结果；步骤2、利用流域DEM数据计算出流域中每一个栅格单元的汇流累积值，生成汇流累积栅格；步骤3、按汇流累积值的最小阈值提取流域中的最密可能河道；步骤4、以最密可能河道上的栅格单元为中心，根据栅格化的水体识别结果，逐河道栅格计算一定范围内的河道可靠性；步骤5、根据栅格化的河道可靠性，建立阈值与河道可靠性均值的关系曲线，确定临界阈值点，该临界阈值点即为集水面积阈值。作为本专利技术所述的一种数字河网提取时集水面积阈值的确定方法进一步优化方案，步骤1具体如下：利用不同地物光谱特征，采用基于遥感影像的水体指数的计算方法，得到栅格化的水体识别结果。作为本专利技术所述的一种数字河网提取时集水面积阈值的确定方法进一步优化方案，所述步骤2中利用流域DEM数据计算出流域中每一个栅格单元的汇流累积值，具体如下：步骤2.1、初始化流域中每一个栅格单元的汇流累积值，取值为1；步骤2.2、以栅格单元Cell为中心，通过周围栅格的高程值与该栅格单元的高程值的对比，找出与其相比最低的栅格单元CellD；步骤2.3、将Cell作为出流栅格，CellD作为入流栅格，入流栅格的汇流累积值在自身原有汇流累积值上加出流栅格的汇流累积值；步骤2.4、逐栅格循环，计算得到每一个栅格单元中的汇流累积值A。作为本专利技术所述的一种数字河网提取时集水面积阈值的确定方法进一步优化方案，所述步骤3具体如下：步骤3.1、结合流域实际自然地理情况，设置汇流累积值的最小阈值Min；步骤3.2、利用Min对步骤2中计算得到的汇流累积栅格进行重分类，生成最密可能河道，汇流累积栅格中汇流累积值高于Min的判定为河道栅格，汇流累积值低于Min的判定为坡地栅格。作为本专利技术所述的一种数字河网提取时集水面积阈值的确定方法进一步优化方案，所述步骤4具体如下：步骤4.1、结合遥感影像与DEM信息，设置河道栅格与水体栅格容许偏差X；水体栅格是：步骤1中得到栅格化的水体识别结果，标记出识别为水体的栅格，生成水体栅格；步骤4.2、逐次访问河道栅格中的每一个栅格单元，寻找以该单元为中心，X*X范围内的水体单元；步骤4.3、根据4.2中已选择的水体单元，计算河道栅格单元的河道可靠性；其中，C为河道栅格单元的河道可靠性，n为X*X范围内水体单元的总数，i为X*X范围内水体单元的编号，xi为第i个水体单元的行号，yi为第i个水体单元的列号，x为河道栅格的行号，y为河道栅格的列号。作为本专利技术所述的一种数字河网提取时集水面积阈值的确定方法进一步优化方案，所述步骤5具体如下：步骤5.1、结合流域实际自然地理情况，设置汇流累积值的最大阈值Max；步骤5.2、在阈值范围[Min,Max]内，逐阈值进行河道判断，统计河道栅格单元的个数num及该num个河道栅格的河道可靠性，计算对应阈值下的河道可靠性均值Cmean；其中，j为某一阈值对应的num个河道栅格中的河道栅格编号,Cj为该阈值下的第j个河道栅格对应的河道可靠性；步骤5.3、建立阈值与河道可靠性均值的关系曲线，曲线斜率转折点即为临界阈值点，该临界阈值点即为集水面积阈值。本专利技术采用以上技术方案与现有技术相比，具有以下技术效果：(1)本专利技术以DEM为基础，结合遥感影像数据，建立阈值与遥感影像数据的函数关系，确定集水面积阈值，这样既保证了计算结果的精度与可靠性，同时解决了人为主观因素对阈值确定的影响问题；(2)本专利技术设计的集水面积阈值，有利于分布式水文模型的直接调用，用于分布式水文模型计算单元的类型划分，可以进一步促进数字水文学以及分布式模型的深入发展，为流域洪水预报提供科学依据，为河道防洪和水库调度等提供了决策支持；(3)本方法主要应用流域数字高程模型与遥感影像数据，数据来源稳定可靠，方法中变量之间的函数关系明确，有利于集水面积阈值的自动化生成，同时保证了结果的客观合理性。附图说明图1是本专利技术的计算流程示意图。图2为本专利技术中计算出的流域水体栅格示意图。图3为本专利技术中计算出的流域汇流累积栅格示意图。图4为本专利技术中计算出的流域最密可能河道栅格示意图。图5为本专利技术中计算出的流域最疏可能河道栅格示意图。图6为本专利技术中计算出的河道单元的河道可靠性栅格分布示意图。图7为本专利技术中计算出的阈值与河道可靠性均值的关系曲线。图8为本专利技术屯溪流域河道水系结果示意图。图9为本专利技术屯溪流域河道计算结果和实测结果对比示意图。具体实施方式为了使本专利技术的目的、技术方案和优点更加清楚，下面将结合附图及具体实施例对本专利技术进行详细描述。如图1所示，本专利技术一种数字河网提取时集水面积阈值的确定方法，包括以下步骤：S1、利用流域遥感影像数据进行流域水体识别，栅格化水体识别结果(如图2)，包括以下步骤：利用不同地物光谱特征，采用基于遥感影像的水体指数的计算方法，得到栅格化的水体识别结果，标记出识别为水体的栅格，生成水体栅格。S2、利用流域DEM数据计算出流域中每一个栅格单元的汇流累积值，生成汇流累积栅本文档来自技高网...

【技术保护点】
1.一种数字河网提取时集水面积阈值的确定方法，其特征在于，包括以下步骤：步骤1、利用流域遥感影像数据进行流域水体识别，栅格化水体识别结果；步骤2、利用流域DEM数据计算出流域中每一个栅格单元的汇流累积值，生成汇流累积栅格；步骤3、按汇流累积值的最小阈值提取流域中的最密可能河道；步骤4、以最密可能河道上的栅格单元为中心，根据栅格化的水体识别结果，逐河道栅格计算一定范围内的河道可靠性；步骤5、根据栅格化的河道可靠性，建立阈值与河道可靠性均值的关系曲线，确定临界阈值点，该临界阈值点即为集水面积阈值。

【技术特征摘要】
1.一种数字河网提取时集水面积阈值的确定方法，其特征在于，包括以下步骤：步骤1、利用流域遥感影像数据进行流域水体识别，栅格化水体识别结果；步骤2、利用流域DEM数据计算出流域中每一个栅格单元的汇流累积值，生成汇流累积栅格；步骤3、按汇流累积值的最小阈值提取流域中的最密可能河道；步骤4、以最密可能河道上的栅格单元为中心，根据栅格化的水体识别结果，逐河道栅格计算一定范围内的河道可靠性；步骤5、根据栅格化的河道可靠性，建立阈值与河道可靠性均值的关系曲线，确定临界阈值点，该临界阈值点即为集水面积阈值。2.根据权利要求1所述的一种数字河网提取时集水面积阈值的确定方法，其特征在于，步骤1具体如下：利用不同地物光谱特征，采用基于遥感影像的水体指数的计算方法，得到栅格化的水体识别结果。3.根据权利要求1所述的一种数字河网提取时集水面积阈值的确定方法，其特征在于，所述步骤2中利用流域DEM数据计算出流域中每一个栅格单元的汇流累积值，具体如下：步骤2.1、初始化流域中每一个栅格单元的汇流累积值，取值为1；步骤2.2、以栅格单元Cell为中心，通过周围栅格的高程值与该栅格单元的高程值的对比，找出与其相比最低的栅格单元CellD；步骤2.3、将Cell作为出流栅格，CellD作为入流栅格，入流栅格的汇流累积值在自身原有汇流累积值上加出流栅格的汇流累积值；步骤2.4、逐栅格循环，计算得到每一个栅格单元中的汇流累积值A。4.根据权利要求1所述的一种数字河网提取时集水面积阈值的确定方法，其特征在于，所述步骤3具体如下：步骤3.1、结合流域实际自然地理情况，设置汇流累...

【专利技术属性】
技术研发人员：徐杰，马亚楠，李致家，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32