[发明专利]一种估算无资料地区库塘坝蓄量的方法

权利要求书

1.一种估算无资料地区库塘坝蓄量的方法，其特征在于，包括以下步骤：

步骤1，利用流域DEM数据计算流域出口点以上的每一个栅格单元的汇流累积值，生成汇流累积栅格；

步骤2，根据汇流累积栅格提取流域中的河道栅格；

步骤3，利用流域遥感影像数据进行流域水体识别，栅格化水体识别结果，生成水体栅格；

步骤4，以每一水体栅格为中心，统计到该水体栅格的距离在聚集容许值范围内的水体栅格，计算该中心水体栅格的水体聚集度；

步骤5，根据水体栅格的水体聚集度、水体栅格到河道栅格的最短距离计算水体栅格的水体可靠度，剔除水体可靠度低于可靠度阈值的水体栅格；

步骤6，选取若干水体栅格作为样本，在数字地图中量取或估算样本水体栅格的测量水深；

步骤7，基于选取的样本水体栅格的测量水深与水体聚集度，建立水深与水体聚集度的关系曲线，根据水体聚集度，逐水体栅格估算对应的计算水深；

步骤8，判定相互距离在一定范围内的水体栅格属于同一水利工程，进行水体划分，根据计算水深求得各水利工程的蓄量；

所述的步骤4具体如下：

步骤4.1，假定以某一水体栅格为中心，一定范围内的水体栅格距该水体栅格越近，数量越多，则该水体栅格越可能属于较大的水利工程，该水体栅格的水深越可能较大，称该水体栅格为中心水体栅格，周围一定范围内的水体栅格的聚集程度为水体聚集度，设定该范围的最大值为聚集容许值R；

步骤4.2，逐次访问每一个水体栅格，寻找以该水体栅格为中心，聚集容许值R为半径范围内的水体栅格集合；

步骤4.3，根据步骤4.2中已选择的水体栅格集合，计算中心水体栅格的水体聚集度；

其中，C为中心水体栅格的水体聚集度，n为水体栅格集合内水体栅格的总数，i为水体栅格集合内水体栅格的编号，x_i为第i个水体栅格的行号，y_i为第i个水体栅格的列号，x为中心水体栅格的行号，y为中心水体栅格的列号，α为水体栅格距离权重；

所述的步骤5具体如下：

步骤5.1，找寻距离水体栅格最近的河道栅格，最短距离为Dmin；

步骤5.2，计算水体栅格的水体可靠度P；

P＝C*Dmin^β；其中，β为河道栅格距离权重；

步骤5.3，设定水体栅格的水体可靠度的阈值为可靠度容许值Pmax，若水体栅格的水体可靠度小于可靠度容许值Pmax，则剔除该水体栅格。

2.根据权利要求1所述的一种估算无资料地区库塘坝蓄量的方法，其特征在于，所述的步骤1中，利用流域DEM数据计算流域出口点以上的每一个栅格单元的汇流累积值，生成汇流累积栅格，具体如下：

步骤1.1，初始化流域中每一个栅格单元的汇流累积值，取值为1；

步骤1.2，以栅格单元Cell为中心，计算周围栅格与该栅格单元的高程差与距离的比值，找出比值最大值对应的栅格单元Cell_D；

步骤1.3，将Cell作出流栅格，Cell_D作为入流栅格，方向由Cell指向Cell_D，入流栅格的汇流累积值在自身原有汇流累积值上加出流栅格的汇流累积值；

步骤1.4，逐栅格循环，计算得到每一个栅格单元的汇流累积值。

3.根据权利要求1所述的一种估算无资料地区库塘坝蓄量的方法，其特征在于，所述的步骤2具体如下：

步骤2.1，结合流域实际自然地理情况，设置汇流累积值的阈值T；

步骤2.2，利用阈值T对步骤1中计算得到的汇流累积栅格进行重分类，汇流累积栅格中汇流累积值大于等于T的判定为河道栅格，汇流累积值小于T的判定为坡地栅格。

4.根据权利要求1所述的一种估算无资料地区库塘坝蓄量的方法，其特征在于，所述的步骤3具体如下：利用不同地物光谱特征，采用基于遥感影像的水体指数的计算方法，得到栅格化的水体识别结果，被识别为水体的栅格记为水体栅格。