一种利用数字高程模型计算地貌单位线初始概率的方法技术

本发明专利技术提供了一种利用数字高程模型计算地貌单位线初始概率的方法，包括以下主要步骤：提取流域数字高程模型(DEM)数据；利用流域DEM数据提取流域分级河道栅格；基于分级河道栅格提取河流交汇点；依据河流交汇点提取各子流域出口点；运用子流域出口点提取各子流域面积；以各子流域面积为基础分析计算得到地貌单位线初始概率。本发明专利技术公开的一种利用数字高程模型计算地貌单位线初始概率的方法具有数据来源稳定可靠、计算效率高、结果客观合理等优点，有利于地貌单位线初始概率的直接求取。

A method for computing initial probability of geomorphologic unit hydrograph using digital elevation model

The present invention provides a method for calculating the initial probability of geomorphologic unit hydrograph using digital elevation model, includes the following main steps: extraction of digital elevation model (DEM) data using DEM data extraction; watershed basin classification of river river river grid; Grid Classification extraction flow intersection based on intersection point extraction based on the river; each sub basin exit point the outlet of the watershed; extraction area of each sub basin; in the area of each sub basin is calculated by analysis of geomorphological unit hydrograph based initial probability. A method of calculating the initial probability of geomorphologic unit hydrograph using digital elevation model method disclosed by the invention has the advantages of stable and reliable data sources, high computation efficiency and the result is objective and reasonable, is conducive to the geomorphological unit hydrograph initial probability is obtained.

【技术实现步骤摘要】
一种利用数字高程模型计算地貌单位线初始概率的方法
本专利技术属于水文
，具体涉及一种基于数字高程模型计算地貌单位线初始概率的方法。
技术介绍
现代研究认为，流域水文响应是地貌扩散和水动力扩散对降落在流域上具有一定时空分布的净雨共同作用的结果。这种对于流域汇流机理的揭示，使得水文学家和自然地理学家试图从理论上建立流域地形地貌特征与流域水文响应之间的联系成为可能。其中的一种方法就是以Horton-Strahler河流分级理论为基础，借用统计物理学中大量“粒子”运动宏观表现的方法来建立地貌单位线。在建立基于Horton-Strahler河流分级的地貌瞬时单位线中重要的一个步骤就是计算水滴随机降落到流域上各级坡面的概率(即初始概率)，以及水滴从低级别河流到高级别河流的转移的概率(即转移概率)。在1979年Rodiguze-Iturbe等人指出初始概率均与水系结构有关，并在Horton地貌参数的基础上运用Smart定律对其进行估算，从而得到水滴从降落在流域面后选择某一条汇流路径运动至流域出口的概率(即路径概率)，进而计算得到流域汇流时间的概率分布密度，即地貌瞬时单位线。但是Smart定律本身就是一种基于数据统计而的得到的规律，运用其估算初始概率时往往并不能得到很好地精度，从而限制了地貌瞬时单位线的应用，不利于地貌瞬时单位线的发展。
技术实现思路
为了解决现有技术中存在的不足，本专利技术提供了一种基于数字高程模型计算地貌单位线初始概率的方法，具有数据来源稳定可靠、计算效率高、结果客观合理等优点，有利于地貌单位线初始概率的直接求取。为解决上述问题，本专利技术具体采用以下技术方案：一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，包括以下步骤：步骤1，提取流域数字高程模型(DEM)数据；步骤2，利用流域DEM数据提取流域分级河道栅格；步骤3，基于分级河道栅格提取河流交汇点；步骤4，依据河流交汇点提取各子流域出口点；步骤5，运用子流域出口点提取各子流域面积；步骤6，以各子流域面积为基础分析计算得到地貌单位线初始概率。前述的一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，所述步骤1中利用地形信息系统软件ArcGIS构建地理处理工作流，获取流域数字高程模型(DEM)数据，包括以下步骤：1.1，填洼，得到填洼之后的栅格Fill；1.2，计算流向，得到流向栅格Dir；1.3，计算汇流流量，得到累计流量栅格Acc；1.4，确定流域出口站点StationPiont；1.5，提取目标流域AimWatershed。前述的一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，所述步骤2中利用流域DEM数据提取流域分级河道栅格的步骤包括：2.1，由流域DEM数据经过填洼，计算流向，计算汇流流量得到累计流量栅格Acc；2.2，设置汇流阈值threshold，依据汇流阈值threshold对累计流量栅格Acc进行重分类，低于汇流阈值threshold的栅格赋空值NODATA，高于汇流阈值threshold的栅格赋值为1，提取出流域河道栅格RiverRaster；2.3，在2.2中提取的流域河道栅格RiverRaster的基础上进行河流分级，将流域河道栅格RiverRaster依据Strahler分级法进行河流分级，得到分级河道栅格。前述的一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，所述步骤3中基于分级河道栅格提取河流交汇点的步骤包括：3.1，将分级河道栅格转换成矢量格式的线型河道LineRiver；3.2，将分级河道栅格转换成矢量格式的点型河道PointRiver；3.3，对于点型河道PointRiver，按线型河道LineRiver中的链与其空间位置相交的原则，逐一提取出点型河道PointRiver中与线型河道LineRiver中的水系的各条链相交的点集RallyPoint：RallyPoint＝LineRiver*PointRiver(式1)其中：*为相交运算；3.4，将点集RallyPoint之间作相交运算，提取出点集RallyPoint之间的交点，即为河网中各条河流链之间的交汇点IntersectPointij：IntersectPointij＝RallyPointi*RallyPointj(i，j＝1,2……n，且i≠j)(式2)其中：*为相交运算，n为点集RallyPoint的总个数，i、j为RallyPoint的序列号，且i≠j。前述的一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，所述步骤4中依据河流交汇点提取各子流域出口点的步骤包括：4.1，提取所用栅格的分辨率a，依据分辨率a确定搜索距离d：4.2，对步骤3中提取出来的各条河流链之间的交汇点进行邻域分析，提取出距每一个交汇点在搜索距离d之内的栅格单元RasterUnit；4.3，对4.2中提取出的每一个的栅格单元RasterUnit进行分析，判断其是否为子流域的出口点，找出栅格单元中的分支河流链的注入点即各子流域出口点InjectPoint：(InjectPoint＝Con(RasterUnit＞threshold)and(RasterUnit＜max(RasterUnit)))(式4)其中：Con为ArcGIS中条件判断工具函数，threshold为汇流阈值。前述的一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，所述步骤5中运用子流域出口点提取各子流域面积的步骤包括：5.1，在步骤2的2.2中提取出的流域河道栅格RiverRaster基础上，将流域河道栅格RiverRaster转换为矢量格式Rivershp：Rivershp＝Rastertoshp(RiverRaster)(式5)其中：Rastertoshp为ArcGIS中栅格转换为矢量的工具函数；5.2，将4.3中提取出来的各子流域出口点InjectPoint与5.1中提取出来的矢量格式Rivershp做相交计算，提取出各子流域出口点InjectPoint所在的矢量格式Rivershp的小分块即子流域出口分块InjectChunk，并读取子流域出口分块InjectChunk的值，然后乘以每一个栅格单元RasterUnit的面积a2，得到相对应的各子流域的面积AChunk：InjectChunk＝InjectPoint*Rivershp(式6)其中：*为相交运算；AChunk＝Value(InjectChunk)×a2(式7)其中：Value(InjectChunk)为子流域出口分块InjectChunk小块中的值。前述的一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，所述步骤6中以各子流域面积为基础分析计算得到地貌单位线初始概率的步骤包括：6.1，将子流域出口分块InjectChunk与经过其中的河流做连接汇总分析，将河流的级别属性赋予InjectChunk；6.2，将同一级的子流域出口分块InjectChunk的值相加得到areaAw，并用高一级的河流的汇流面积areaAw+1减去低一级的河流的汇流面积areaAw计算得到每一级河流的汇水面积Aw：Aw＝areaAw+1-areaAw(式8)6.3，基于每一级河流的汇水面本文档来自技高网...

【技术保护点】
一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，包括以下步骤：步骤1，提取流域数字高程模型数据；步骤2，利用流域数字高程模型数据提取流域分级河道栅格；步骤3，基于分级河道栅格提取河流交汇点；步骤4，依据河流交汇点提取各子流域出口点；步骤5，运用子流域出口点提取各子流域面积；步骤6，以各子流域面积为基础分析计算得到地貌单位线初始概率。

【技术特征摘要】
1.一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，包括以下步骤：步骤1，提取流域数字高程模型数据；步骤2，利用流域数字高程模型数据提取流域分级河道栅格；步骤3，基于分级河道栅格提取河流交汇点；步骤4，依据河流交汇点提取各子流域出口点；步骤5，运用子流域出口点提取各子流域面积；步骤6，以各子流域面积为基础分析计算得到地貌单位线初始概率。2.根据权利要求1所述的一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，步骤1中利用地形信息系统软件ArcGIS构建地理处理工作流，获取流域数字高程模型数据，包括以下步骤：1.1，填洼，得到填洼之后的栅格Fill；1.2，计算流向，得到流向栅格Dir；1.3，计算汇流流量，得到累计流量栅格Acc；1.4，确定流域出口站点StationPiont；1.5，提取目标流域AimWatershed。3.根据权利要求2所述的一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，步骤2中利用流域数字高程模型数据提取流域分级河道栅格的步骤包括：2.1，由流域数字高程模型数据经过填洼、计算流向、计算汇流流量得到累计流量栅格Acc；2.2，设置汇流阈值threshold，依据汇流阈值threshold对累计流量栅格Acc进行重分类，低于汇流阈值threshold的累计流量栅格赋空值NODATA，高于汇流阈值threshold的累计流量栅格赋值为1，提取出流域河道栅格RiverRaster；2.3，在2.2中提取的流域河道栅格RiverRaster的基础上进行河流分级，将流域河道栅格RiverRaster依据Strahler分级法进行河流分级，得到流域分级河道栅格。4.根据权利要求3所述的一种利用数字高程模型计算地貌单位线初始概率的方法，其特征在于，步骤3中基于分级河道栅格提取河流交汇点的步骤包括：3.1，将分级河道栅格转换成矢量格式的线型河道LineRiver；3.2，将分级河道栅格转换成矢量格式的点型河道PointRiver；3.3，对于点型河道PointRiver，按线型河道LineRiver中的链与其空间位置相交的原则，逐一提取出点型河道PointRiver中与线型河道LineRiver中的水系的各条链相交的点集RallyPoint：RallyPoint＝LineRiver*PointRiver(式1)其中：*为相交运算；3.4，将点集RallyPoint之间作相交运算，提取出点集RallyPoint之间的交点，即为河网中各条河流链之间的交汇点IntersectPointij：IntersectPointij＝RallyPointi*RallyPointj(i，j＝1,2……n，且i≠j)(式2)其中：*为相交运算，n为点集RallyPoint的总个数，i、j为RallyPoint的序列号，且i≠j。5....

【专利技术属性】
技术研发人员：童冰星，李致家，邓元倩，
申请(专利权)人：河海大学，
类型：发明
国别省市：江苏,32