[发明专利]一种自动识别流域数字化过程中出口点的方法

说明书

技术领域

本发明涉及流域数字化，尤其是识别流域数字化过程中出口点的方法。

背景技术

为了更好地描述水文现象、探索和揭示水文规律，水文学研究正逐步迈向数字化，近年来，随着大范围、高精度数字高程模型（DEM）的出现，数字水文学的研究得到了进一步发展。利用DEM对流域进行数字化也就成为了构建流域水文模型、进行水文模拟与研究的前提。

DEM原始数据是进行流域数字化的基础数据，在利用该基础数据进行流域数字化时，其关键是如何精确、便捷地识别出流域出口点。流域出口点的准确定位对数字化结果有着重要影响，在流域数字化过程中，即使是位置相近的几个点，也会生成完全不同的研究流域边界和集水区域。因此，如果流域出口点定位不准确，不仅限制了水文模拟的精度，也不利于真实掌握研究流域的降雨-径流响应。

目前，流域数字化过程中对于流域出口点的识别方法主要是依据水文站点的经纬度坐标，但是由于水文站点坐标定位时本身存在误差，导致了站点所对应位置往往并不是流域的真实出口点，现有方法则需要多次人为干预对站点位置进行修正，才能确定出流域出口，生成流域集水区域。

但是，人为干预修正流域出口点非常繁琐复杂，尤其当研究流域较多时，数字化过程也将更加的费时费力，不利于数字水文学的发展。

发明内容

发明目的：提供一种自动识别流域数字化过程中出口点的方法，以解决现有技术需要人为修正出口点，效率低下的问题。

技术方案：一种自动识别流域数字化过程中出口点的方法，包括如下步骤：

S1、读取研究区域的数字高程模型的原始数据并进行预处理，依次获得流向矩阵、累积汇水面积矩阵、栅格控制面积矩阵与河道栅格矩阵；

S2、基于水文站点的经纬度坐标，对水文站点在河道栅格矩阵内的位置进行定位，并判断水文站点是否位于河道栅格上，若不位于河道栅格上，则执行S3；若位于河道栅格上，该河道栅格即为流域潜在出口点，则执行S4；

S3、根据流向矩阵和水文站点，识别流域潜在出口点，并进行标记；

S4、根据栅格控制面积矩阵，提取流域潜在出口点的控制面积，并计算该面积与水文站集水面积之间的相对误差δ；

S5、根据设定的面积许可误差，判断δ的绝对值是否小于该面积许可误差，若判断结果为否，则执行S6；若判断结果为是，则执行S7；

S6、根据流向矩阵，对流域的潜在出口点进行修正，然后重复S4和S5；

S7、自动识别过程结束，当前的流域潜在出口点即为流域出口点。

所述数字高程模型的原始数据采用ArcInfo ASCII格式，地理坐标系统采用的是WGS84经纬度坐标。

栅格控制面积矩阵的提取步骤包括：

1）由数字高程模型分辨率计算得到单个栅格单元的面积；

2）以累积汇水面积矩阵为基础，将每个栅格单元的上游累积栅格数乘以单个栅格单元的面积，得到栅格控制面积矩阵。

河道栅格矩阵的提取步骤包括：

1）设定一个生成河道栅格的阈值；

2）利用栅格控制面积矩阵，按照给定的阈值对每个栅格单元进行标记，若栅格单元的控制面积小于给定的阈值，将该单元标记为无数据区；若栅格单元的控制面积大于给定的阈值，将该单元标记为河道栅格，得到河道栅格矩阵。

对水文站点在河道栅格矩阵内的位置进行定位包括如下步骤：

1）由水文站点的纬度坐标计算该站点在河道栅格矩阵内所处的行数：

SNx＝nrows-INT[(SLat-yllcorner)×60×(60/(3600×cellsize))]；

式中，SNx为水文站点在河道栅格矩阵内所处的行数；SLat为水文站点的纬度坐标；nrows为河道栅格矩阵的总行数；yllcorner为河道栅格矩阵西南角栅格单元对应的纬度坐标；cellsize为DEM分辨率；INT[]为取整函数；

2）由水文站点的经度坐标计算该站点在河道栅格矩阵内所处的列数：

SNy＝INT[(SLon-xllcorner)×60×(60/(3600×cellsize))]+1；

式中，SNy为水文站点在河道栅格矩阵内所处的列数；SLon为水文站点的经度坐标；xllcorner为河道栅格矩阵西南角栅格单元对应的经度坐标；