CN109815611A - 一种基于数字流域的流域边界生成方法 - Google Patents
一种基于数字流域的流域边界生成方法 Download PDFInfo
- Publication number
- CN109815611A CN109815611A CN201910102879.3A CN201910102879A CN109815611A CN 109815611 A CN109815611 A CN 109815611A CN 201910102879 A CN201910102879 A CN 201910102879A CN 109815611 A CN109815611 A CN 109815611A
- Authority
- CN
- China
- Prior art keywords
- pixel
- water
- basin
- flow
- neighborhood
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Alarm Systems (AREA)
- Image Processing (AREA)
Abstract
本发明公开了一种基于数字流域的流域边界生成方法,首先基于数字高程模型数据DEM计算水流方向矩阵;根据水流方向矩阵,计算集水面积;生成水系和单元;生成流域边界;流域边界与降雨预报图叠加得到流域降雨量预报数值。本发明通过DEM模型方式提取流域边界,精确绘制了水电站所在流域,使天气预报图中流域范围可视化,实现流域区域内来水的准确预测。
Description
技术领域
本发明属于水文技术领域,具体涉及一种基于数字流域的流域边界生成方法。
背景技术
降雨是水电企业产生效益的重要影响因素之一,而未来降雨预报是梯级水库优化调度决策的主要依据。中央气象台每日滚动发布的未来24~168小时降水预报是面向公众的社会服务产品,其作为中国最高水平的气象产品,对水库调度具有很强指导意义,但并没有针对水电企业的专业服务,没有针对流域进行边界划分,流域区域(边界)不明确导致专业人员、非专业人员只能凭经验判断降雨区域,即无法准确判断降雨是否在流域区域内。流域区域降雨预报不准确会影响到来水预测,引起发电方式安排不科学合理等一系列问题,无法满足水电精益化调度分析的需求。
发明内容
为解决上述问题,本发明提出一种基于数字流域的流域边界生成方法,实现流域区域内来水的准确预测。
本发明采用如下技术方案,一种基于数字流域的流域边界生成方法,具体步骤如下:
1)基于数字高程模型数据(Digital Elevation Model,DEM)计算水流方向矩阵;
2)根据水流方向矩阵,计算集水面积;
3)生成水系和单元;
4)生成流域边界;
5)流域边界与降雨预报图叠加得到流域降雨量预报数值。
优选地,在所述步骤1)之前对数字高程模型数据DEM进行预处理,具体为预先对数字高程模型数据DEM进行填洼处理,得到无洼地的数字高程模型数据DEM。
优选地,所述步骤1)计算水流方向矩阵的方法为,利用D8算法计算水流的路径及方向,得到与DEM维数相同的水流方向矩阵,具体包括以下步骤:
11)将中心像元的八邻域像元编码,所述八邻域像元编码用于表示水流方向;
12)根据中心像元及八邻域像元的高程关系,计算下降方向,公式为为:
其中,z为邻域像元的高程;value为中心像元的高程;distance表示像元中心之间的距离;
下降方向最大值即为最陡下降方向;
若八邻域像元中存在多个最陡下降方向,则扩大相邻像元范围,直至得到唯一最陡下降方向为止;
13)利用最陡下降方向对应的像元中的编码对输出像元进行赋值,最终得到与DEM维数相同的水流方向矩阵。
优选地,所述步骤2)中计算集水面积的方法为:根据水流方向矩阵计算得出汇流累积量,即上游累积集水面积,每个像元的汇流累积量表示最终汇流经过像元的上游像元的水量,具体方法如下:
依次遍历水流方向矩阵的像元,若当前像元的八邻域像元中没有流入的水流,则当前像元的汇流累积量为零,即当前像元的输出值为零;如果当前像元的八邻域像元中存在汇入方向的水流,则根据分配的水量权重叠加计算汇流累积量,作为当前像元的输出值;
最终获得当前区域中每个像元的汇流累积量,形成汇流累积量矩阵,即为集水面积矩阵。
优选地,所述步骤3)中生成水系和单元的方法为:
在集水面积矩阵中,大于预设的面积阈值的像元为河道,在集水面积矩阵上标注出河道的位置,生成河流网络栅格;
细化河流网络栅格,直到留下彼此连通的由单个栅格点组成的图形;
细化后的河流网络栅格即为骨架栅格,从骨架栅格的西北角开始,根据八邻域进行搜索,依次跟踪相邻栅格,将栅格点转换为矢量点坐标序列并连线,得到矢量化的水系,水系间的分水岭构成单元。
优选地,所述步骤4)中生成流域边界的方法为:确定一个流域出水口,即集水区的最低点,结合水流方向矩阵,从出水口开始遍历每个像元的八邻域,寻找水流方向流向当前像元的邻域单元,记为当前像元的上游像元;
以上游像元继续向上游逐层搜索,直到某一层上游像元的八邻域均不流向当前像元时停止;
最终生成一个栅格流域的包围圈,即流域边界。
优选地,所述步骤5)的方法具体为,设定流域边界投影方式与降雨预报图的投影方式相同,将流域边界与降雨预报图叠加,确定落在流域边界内的格点,获取格点降雨预报数值,从而计算出流域降雨量预报数值。
发明所达到的有益效果:本发明是一种基于数字流域的流域边界生成方法,实现流域区域内来水的准确预测。本发明采用DEM模型方式提取流域边界,精确绘制了水电站所在流域,使天气预报图中流域范围可视化,提高精益调度气象支持信息的准确性,来水预测精度提高,能够用于指导水库调度,为水电节水增发创造经济效益。
附图说明
图1是本发明实施例中一种中心像元的八邻域像元编码方式示意图;
图2是本发明实施例中一种集水面积计算示意图;
图3是本发明实施例中一种流域边界图;
图4是本发明实施例中一种流域边界与天气预报叠加图。
具体实施方式
下面根据附图并结合实施例对本发明的技术方案作进一步阐述。
实施例1:
一种基于数字流域的流域边界生成方法,具体步骤如下:
1)基于数字高程模型数据DEM计算水流方向矩阵;
2)根据水流方向矩阵,计算集水面积;
3)生成水系和单元;
4)生成流域边界;
5)流域边界与降雨预报图叠加得到流域降雨量预报数值。
实施例2:
在实施例1的基础上,在所述步骤1)之前对数字高程模型数据DEM进行预处理,具体为预先对数字高程模型数据DEM进行填洼处理,得到无洼地的数字高程模型数据DEM。
DEM在离散化过程中的插值误差和采样误差,造成许多洼地(凹陷型洼地和阻挡型洼地),这些洼地将在水流方向计算时,造成有些水不能流出流域边界,从而产生很大的误差或不能计算出合理的结果。流域信息的提取需要高质量的DEM,要求没有坑和坝,需要对DEM进行预处理,即对原始的DEM进行填洼处理,得到无洼地的DEM。
填洼过程是对DEM进行地表面水流动的模拟,假设有无限量的降水,雨水不断在地表形成径流,由于洼地是在周边地形中无水流方向的像元,则可以根据水流方向判断出洼地的位置,继而将凹陷处填平。通过不断地迭代,即填平之后再次判断是否存在洼地,如果存在则继续填洼处理,直到设定高程阈值(判断洼地的高程差)内所有的洼地填充完成。
作为一种较佳的实施例,所述步骤1)计算水流方向矩阵的方法为,利用D8算法计算水流的路径及方向,得到与DEM维数相同的水流方向矩阵,具体包括以下步骤:
11)将中心像元的八邻域像元编码,所述八邻域像元编码用于表示水流方向,如图1所示;每个水流方向(水流离开中心像元的指向)都有固定的编码,以编码矩阵表示该区域水流的方向;
12)根据中心像元及八邻域像元的高程关系,计算下降方向,公式为为:
其中,z为邻域像元的高程;value为中心像元的高程;distance表示像元中心之间的距离;如果像元大小为1,则两个正交像元之间的距离为1,两个对角线像元之间的距离为
下降方向最大值即为最陡下降方向;
若八邻域像元中存在多个最陡下降方向,则扩大相邻像元范围,直至得到唯一最陡下降方向为止;
水流方向的连线即为水流路径;
13)利用最陡下降方向对应的像元中的编码对输出像元进行赋值,最终得到与DEM维数相同的水流方向矩阵。
优选地,所述步骤2)中计算集水面积的方法为:根据水流方向矩阵计算得出汇流累积量,即上游累积集水面积,每个像元的汇流累积量表示最终汇流经过像元的上游像元的水量,具体方法如下:
依次遍历水流方向矩阵的像元,若当前像元的八邻域像元中没有流入的水流,则当前像元的汇流累积量为零,即当前像元的输出值为零;如果当前像元的八邻域像元中存在汇入方向的水流,则根据分配的水量权重叠加计算汇流累积量,作为当前像元的输出值;所述水量权重可以缺省,缺省值为1,即认为以规则格网表示的数字高程模型的每个单元都有一个单位的水量。
最终获得当前区域中每个像元的汇流累积量,形成汇流累积量矩阵,即为集水面积矩阵,如图2所示。
高流量的单元是集中流动区域,汇流累积的数值越大,该区域越容易形成地表径流。因此集水面积的计算可用于识别河道。
优选地,所述步骤3)中生成水系和单元的方法为:
在集水面积矩阵中,大于预设的面积阈值的像元为河道,在集水面积矩阵上标注出河道的位置,生成河流网络栅格;本实施例中面积阈值设置为38841。
细化河流网络栅格,直到留下彼此连通的由单个栅格点组成的图形;
所述细化河流网络栅格具体为从曲线的边缘开始,每次剥掉一个栅格宽的一层,直到留下彼此连通的由单个栅格点组成的图形。
细化后的河流网络栅格即为骨架栅格,从骨架栅格的西北角开始,根据八邻域进行搜索,依次跟踪相邻栅格,将栅格点转换为矢量点坐标序列并连线,得到矢量化的水系,水系间的分水岭构成单元。
优选地,所述步骤4)中生成流域边界的方法为:确定一个流域出水口,即集水区的最低点,结合水流方向矩阵,从出水口开始遍历每个像元的八邻域,寻找水流方向流向当前像元的邻域单元,记为当前像元的上游像元;
以上游像元继续向上游逐层搜索,直到某一层上游像元的八邻域均不流向当前像元时停止;
最终生成一个栅格流域的包围圈,即流域边界,如图3所示。
优选地,所述步骤5)的方法具体为,设定流域边界投影方式与降雨预报图的投影方式相同,将流域边界与降雨预报图叠加,确定落在流域边界内的格点,获取格点降雨预报数值,从而计算出流域降雨量预报数值。
中央气象台降雨量预报图每日8点、20点各更新一次。采用自动下载程序每隔半小时刷新、下载新的预报图至本地文件夹。
将生成好的流域边界与降雨量预报图叠加,设定流域边界投影方式(与中央气象台降雨量预报图采用同样的投影方式),将流域边界与降雨预报图叠加,如图4所示,确定落在流域边界内的格点,获取格点降雨预报数值,从而计算出流域降雨量预报数值。
以上所述仅为本发明的优选实施方式,用以说明本发明的技术方案,而非对其限制;应当指出:对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (7)
1.一种基于数字流域的流域边界生成方法,其特征在于,包括以下步骤:
1)基于数字高程模型数据DEM计算水流方向矩阵;
2)根据水流方向矩阵,计算集水面积;
3)生成水系和单元;
4)生成流域边界;
5)流域边界与降雨预报图叠加得到流域降雨量预报数值。
2.根据权利要求1所述的一种基于数字流域的流域边界生成方法,其特征在于,在所述步骤1)之前对数字高程模型数据DEM进行预处理,具体为预先对数字高程模型数据DEM进行填洼处理,得到无洼地的数字高程模型数据DEM。
3.根据权利要求1所述的一种基于数字流域的流域边界生成方法,其特征在于,所述步骤1)计算水流方向矩阵的方法为,利用D8算法计算水流的路径及方向,得到与DEM维数相同的水流方向矩阵,具体包括以下步骤:
11)将中心像元的八邻域像元编码,所述八邻域像元编码用于表示水流方向;
12)根据中心像元及八邻域像元的高程关系,计算下降方向,公式为为:
其中,z为邻域像元的高程;value为中心像元的高程;distance表示像元中心之间的距离;
下降方向最大值即为最陡下降方向;
若八邻域像元中存在多个最陡下降方向,则扩大相邻像元范围,直至得到唯一最陡下降方向为止;
13)利用最陡下降方向对应的像元中的编码对输出像元进行赋值,最终得到与DEM维数相同的水流方向矩阵。
4.根据权利要求1所述的一种基于数字流域的流域边界生成方法,其特征在于,所述步骤2)中计算集水面积的方法为:根据水流方向矩阵计算得出汇流累积量,即上游累积集水面积,每个像元的汇流累积量表示最终汇流经过像元的上游像元的水量,具体方法如下:
依次遍历水流方向矩阵的像元,若当前像元的八邻域像元中没有流入的水流,则当前像元的汇流累积量为零,即当前像元的输出值为零;如果当前像元的八邻域像元中存在汇入方向的水流,则根据分配的水量权重叠加计算汇流累积量,作为当前像元的输出值;
最终获得当前区域中每个像元的汇流累积量,形成汇流累积量矩阵,即为集水面积矩阵。
5.根据权利要求1所述的一种基于数字流域的流域边界生成方法,其特征在于,所述步骤3)中生成水系和单元的方法为:
在集水面积矩阵中,大于预设的面积阈值的像元为河道,在集水面积矩阵上标注出河道的位置,生成河流网络栅格;
细化河流网络栅格,直到留下彼此连通的由单个栅格点组成的图形;
细化后的河流网络栅格即为骨架栅格,从骨架栅格的西北角开始,根据八邻域进行搜索,依次跟踪相邻栅格,将栅格点转换为矢量点坐标序列并连线,得到矢量化的水系,水系间的分水岭构成单元。
6.根据权利要求1所述的一种基于数字流域的流域边界生成方法,其特征在于,所述步骤4)中生成流域边界的方法为:确定一个流域出水口,即集水区的最低点,结合水流方向矩阵,从出水口开始遍历每个像元的八邻域,寻找水流方向流向当前像元的邻域单元,记为当前像元的上游像元;
以上游像元继续向上游逐层搜索,直到某一层上游像元的八邻域均不流向当前像元时停止;
最终生成一个栅格流域的包围圈,即流域边界。
7.根据权利要求1所述的一种基于数字流域的流域边界生成方法,其特征在于,所述步骤5)的方法具体为,设定流域边界投影方式与降雨预报图的投影方式相同,将流域边界与降雨预报图叠加,确定落在流域边界内的格点,获取格点降雨预报数值,从而计算出流域降雨量预报数值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910102879.3A CN109815611B (zh) | 2019-02-01 | 2019-02-01 | 一种基于数字流域的流域边界生成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910102879.3A CN109815611B (zh) | 2019-02-01 | 2019-02-01 | 一种基于数字流域的流域边界生成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109815611A true CN109815611A (zh) | 2019-05-28 |
CN109815611B CN109815611B (zh) | 2020-01-14 |
Family
ID=66605131
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910102879.3A Active CN109815611B (zh) | 2019-02-01 | 2019-02-01 | 一种基于数字流域的流域边界生成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109815611B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110737931A (zh) * | 2019-08-28 | 2020-01-31 | 中国铁路设计集团有限公司 | 一种基于ArcGIS的铁路桥渡水文关键参数提取方法 |
CN111008259A (zh) * | 2019-12-06 | 2020-04-14 | 中国水利水电科学研究院 | 一种流域降雨相似性搜索方法 |
CN113094916A (zh) * | 2021-04-21 | 2021-07-09 | 东莞理工学院 | 一种通过dem数据精准模拟生成河网的方法 |
CN117171287A (zh) * | 2023-11-02 | 2023-12-05 | 航天宏图信息技术股份有限公司 | 地形流域的提取划分方法、装置、电子设备和存储介质 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902844A (zh) * | 2012-09-03 | 2013-01-30 | 南京师范大学 | 基于大数据量dem数据的子流域划分方法 |
CN107220754A (zh) * | 2017-05-18 | 2017-09-29 | 长江水利委员会长江科学院 | 一种县域尺度山洪灾害风险评估方法 |
CN108399309A (zh) * | 2018-03-16 | 2018-08-14 | 中国水利水电科学研究院 | 一种大尺度复杂地形区分布式水文模型的子流域划分方法 |
-
2019
- 2019-02-01 CN CN201910102879.3A patent/CN109815611B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102902844A (zh) * | 2012-09-03 | 2013-01-30 | 南京师范大学 | 基于大数据量dem数据的子流域划分方法 |
CN107220754A (zh) * | 2017-05-18 | 2017-09-29 | 长江水利委员会长江科学院 | 一种县域尺度山洪灾害风险评估方法 |
CN108399309A (zh) * | 2018-03-16 | 2018-08-14 | 中国水利水电科学研究院 | 一种大尺度复杂地形区分布式水文模型的子流域划分方法 |
Non-Patent Citations (2)
Title |
---|
江岭: ""基于DEM的流域地形分析并行算法关键技术研究"", 《中国博士学位论文全文数据库 基础科学辑》 * |
谭万沛 等: ""暴雨泥石流预报程式"", 《自然灾害学报》 * |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110737931A (zh) * | 2019-08-28 | 2020-01-31 | 中国铁路设计集团有限公司 | 一种基于ArcGIS的铁路桥渡水文关键参数提取方法 |
CN111008259A (zh) * | 2019-12-06 | 2020-04-14 | 中国水利水电科学研究院 | 一种流域降雨相似性搜索方法 |
CN113094916A (zh) * | 2021-04-21 | 2021-07-09 | 东莞理工学院 | 一种通过dem数据精准模拟生成河网的方法 |
CN117171287A (zh) * | 2023-11-02 | 2023-12-05 | 航天宏图信息技术股份有限公司 | 地形流域的提取划分方法、装置、电子设备和存储介质 |
CN117171287B (zh) * | 2023-11-02 | 2024-01-30 | 航天宏图信息技术股份有限公司 | 地形流域的提取划分方法、装置、电子设备和存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109815611B (zh) | 2020-01-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109815611A (zh) | 一种基于数字流域的流域边界生成方法 | |
CN109815305B (zh) | 一种无资料地区场次洪水径流过程反演的方法 | |
CN108874936B (zh) | 一种基于改进新安江模型的适用于山丘区的水文预报方法 | |
CN108399309B (zh) | 一种大尺度复杂地形区分布式水文模型的子流域划分方法 | |
CN107704592A (zh) | 一种基于WebGIS的洪水预报服务构建方法 | |
CN104298841A (zh) | 一种基于历史数据的洪水预报方法和系统 | |
CN106168991A (zh) | 一种基于水动力数值模拟的感潮河网潮位预报方法 | |
CN111898303A (zh) | 基于气象预报及水动力模拟的流域水位及内涝预报方法 | |
CN107330086B (zh) | 一种提高无资料高海拔流域水文过程模拟精度的方法 | |
CN105320827A (zh) | 降雨径流预报系统 | |
CN105893770A (zh) | 一种气候变化和人类活动对流域水资源影响的量化方法 | |
CN109145376A (zh) | 一种基于cad地形图的水库淹没范围确定方法 | |
CN110096751B (zh) | 一种估算无资料地区中小型水库蓄量的方法 | |
CN102902893A (zh) | 一种基于dem的汇水区降雨积水深度的计算方法 | |
CN105303612A (zh) | 一种基于不规则三角网模型的数字河网提取方法 | |
CN104408318A (zh) | 一种河流型水功能区季节性设计流量计算方法 | |
CN112149314A (zh) | 一种基于虚拟库容修正的多沙水库库容冲淤模拟方法 | |
CN109633790B (zh) | 自然流域划分法中确定子流域降雨量时空分布的方法 | |
CN109979172A (zh) | 一种基于新安江模型的动态山洪临界雨量预报方法 | |
Chaube et al. | Synthesis of flow series of tributaries in Upper Betwa basin | |
CN114091277B (zh) | 一种考虑初始状态变量影响的新安江模型参数率定方法 | |
CN108269199A (zh) | 一种面向对象的小水库群时空分布式出流计算方法 | |
CN112904458B (zh) | 一种超渗-蓄满混合产流模式的水文预报方法及系统 | |
CN109636908A (zh) | 一种水电厂流域图绘制与发布方法 | |
CN116305933B (zh) | 一种基于dem数据的坡面产汇流简易计算方法及装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |