申请/专利权人：黄河水利委员会黄河水利科学研究院

申请日：2021-09-16

公开（公告）日：2024-04-09

公开（公告）号：CN113989653B

主分类号：G06V20/13

分类号：G06V20/13;G06V10/26;G06V10/28;G06V10/30;G06V10/34;G06V10/44

优先权：

专利状态码：有效-授权

法律状态：2024.04.09#授权;2022.02.18#实质审查的生效;2022.01.28#公开

摘要：本发明公开了一种河流平面几何指标与拓扑结构的提取方法，其IPC分类号为E02B100。该方法基于复杂网络分析技术，以卫星图片为输入量，由河道水体的提取、河流中心线骨架的刻画，河流分汊点与端点的识别与检测和河流网络拓扑结构与平面几何指标的提取四个主要部分组成，在考虑了江心滩，沙坝等对于河道平面形态所产生的影响的同时，尽可能保证了所提取的河道的连续性，因此对具有复杂平面形态的辫状河流，畸形河湾等同样具有较强的适用性。本发明不仅节约了现场测量所带来的大量人力物力成本，还兼顾了所获取数据的准确性与多样性，为河流研究提供了技术支持。

主权项：1.一种河流平面几何指标与拓扑结构的提取方法，其特征在于，包括以下步骤：S1、确定输入图像，并率定初始参数，包括：S11、输入的图像格式为卫星图片，首先确定所输入图像的分辨率，以及图像尺寸M×N，其中，M，N分别为图像的像素行数与列数；其次，按照改进的归一化差异水体指数公式：计算并储存每个图像的像素点所对应的MNDWI，其中Green与MIR分别为卫星图像中绿色波段与中红外波段的反射率；S12、率定模型初始参数σ1，σ2与ΔW，其中σ1与σ2分别为去除噪声过程中膨胀与腐蚀算法所使用的卷积核的最大尺寸，而ΔW为河宽提取指标参数；初始参数的最优值均通过对所输出的河道数目的敏感性分析而确定，具体过程为，分别从0开始逐渐增加各参数值，当所输出到的河道数目趋于稳定值后，所对应的参数值即为最优值；S2、提取河道水体并消除噪声，包括：S21、提取河道水体：根据S11计算所得到的MNDWI指数，对图像所有像素点运用OTSU阈值分割方法，找出水体像素与陆地像素的阈值T，对于MNDWI指数大于T的像素定义为水体像素，剩下非水体像素则被当作背景去除，并输出河道水体的二值图像，值为1的像素点为河道水体像素，值为0的像素点为背景像素；S22、噪声消除：通过S21所得到的二值图像中的像素点与卷积核进行卷积运算，其中卷积核的边缘像素值设为1，中间像素值为0；首先，遍历所有水体像素点，根据膨胀算法使其分别与尺寸为n1n1＝1，3，5.....σ1的卷积核进行卷积，若结果为0，则标记为零星水域噪声并予以消除；其次，同样地，遍历所有背景像素点，根据腐蚀算法分别使其与尺寸为n2n2＝1，3，5.....σ2的卷积核进行卷积，若结果为4σ2-4，则标记为对河流平面形态无主要影响的较小的江心滩或沙坝并予以填充，最终得到整个河流系统的掩膜；S3、刻画河道中心线骨架；包括：S31、对于步骤S22得到的河流掩膜进行Zhang-Suen快速并行细化算法，得到河流骨架，记为矩阵S＝{Xs,Ys}，其中Xs与Ys分别为第s个骨架像素点的行、列坐标；S32、根据步骤S31所得到的河流骨架，遍历所有骨架像素点，分别采用消除模版与其进行卷积，所得结果若属于集合{65，5，20，80，13，22，52，133，141，54}，则删除该骨架像素点，以保证所得河流骨架为单分辨率宽度；S4、检测河流的分汊点与端点；包括：S41、记骨架像素点的8个邻域中同样为骨架的像素点个数为Ns，遍历所有河流骨架像素点，并检测其Ns的取值，如果Ns＝3，则判定该骨架像素点为河流分汊点；如果Ns＝1，则判定其为河流端点；S42、记录所有得到的河流分汊点与端点坐标，并按照从河流上游到下游的顺序对这些点进行排序；S5、构建河流的拓扑结构，提取河流的平面几何指标，包括：S51、根据步骤S42所得到的河流分汊点与端点进行追踪，首先，按顺序从河流的一个分汊点或者端点i开始，逐个追踪并记录与其相连的骨架像素点及其坐标，直至另一个河流分汊点或端点j结束；然后，按照相同的方式追踪另外一条与i相连的河流骨架；最后，直到所有与i相连的河流骨架均被追踪完成；S52、河流拓扑结构的构建，河流拓扑结构的数学表达通过其河网连通矩阵An,n＝{ai，j}来实现，其中n为河流分汊点与端点的数目，ai，j表示河流中第i个与第j个分汊点或端点的连通关系，ai,j＝1表示两分汊点或端点连通，ai,j＝0则表示不连通；首先，假设连通矩阵An,n为零矩阵，设在骨架追踪过程中，若检测到一条河流骨架的两端分别为点i与j并且i与j不同时为河流端点，则ai,j＝1,否则ai,j＝0，最后，直至所有的河流分汊点与端点均已追踪完成，则所输出的连通矩阵即为该河流的网络拓扑结构表示；S53、河流平面几何指标的提取，在从分汊点或端点i到j的河流骨架追踪的过程中，通过对于骨架元素点坐标运算来提取河流的河道长度与曲率，根据表达式与分别求出该骨架所代表的河道的曲线长度直线长度以及曲率其中与分别为连接点i和j的河流骨架中第k个像素点的横纵坐标；而XiYi与XjYj分别是河流分汊点或端点i与j的横纵坐标；K为骨架中所追踪到的像素点总个数；求出河道的宽度：首先，在追踪过程中，根据表达式来计算连接点i与j的河流骨架中的第k个像素点的局部坡度，垂直于该坡度在河流掩膜中向该像素点左右两侧做垂线得到左右横截面，若其左右横截面差值的绝对值小于ΔW，则该横截面被判定为有效横截面，其长度被记录为像素点k的局部河宽然后，当一条河流骨架被追踪完成后，计算所有被记录的该骨架上的有效横截面宽度的平均值，即为该骨架所对应的河道的河宽。

全文数据：

权利要求：

百度查询： 黄河水利委员会黄河水利科学研究院 一种河流平面几何指标与拓扑结构的提取方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD