一种基于分水岭的山脊线提取方法

摘要

本发明的地形图像山脊线提取方法，它包括找出为分割线，由该分割线在地形图像中分割出区域为第一分割区域、第二分割区域、第三分割区域……；标注出第一分割区域、第二分割区域、第三分割区域……中每个像素的灰度值，并且将上述每个区域中灰度值相同的点就近连接起来，分别画出上述第一分割区域、第二分割区域、第三分割区域……的等高线折线；然后拟合出并画出每根等高线折线的等高线曲线，再计算出每根等高线曲线上各个点的曲率；将区域内最大等高线曲率对应的点按照就近的原则连接成折线，该折线即为分割区域的山脊线。

说明书

技术领域

本发明涉及了一种地形图像山脊线的提取，特别是用于一种利用图像处理的方法进行地形图像的山脊线提取的处理方法。

背景技术

在用图像处理的方法对地形图像进行信息提取时，山脊线作为一个重要的地形特征，获取地形图像的山脊线信息极其最要。

目前已有的提取算法包括等高线骨架法、等高线曲率最大判别法等，等高线骨架法视地形特征线两侧的地形为对称变化，存在的问题是等高线有噪声时提取错误；等高线曲率最大判断法使用全局阈值进行判断，全局阈值过大时会丢失某些方向的山脊线的点；全局阈值过小时会增加噪点。该方案中就是在分水岭算法的基础上利用等高线曲率最大判别法对山脊线进行提取。

发明内容

本申请的发明目的在于提供一种地形图像山脊线的提取方法，该方法在地形图像中进行提取对山脊进行判断分割,然后根据分割后的山脊线计算其等高线以及等高线的曲率,根据等高线曲率判断该点是否属于山脊线，最终得到较好的山谷、山脊线。

为了完成本申请的发明目的，本申请采用以下技术方案：

本发明的一种地形图像山脊线的提取方法，该地形图像包括至少一个连通区域，其中：它包括以下步骤：

（1）对拍摄的地形图像中每个像素的灰度值进行排序，找出灰度极小值，并且统计地形图像的灰度值极小值个数n和每个连通区域的灰度极小值个数m₁、m₂、m₃……,n=>1+m₂+m₃……；

（2）如果n<2，则对地形图像的灰度极小值加1，返回步骤1；

（3）如果n>2，但是m₁、m₂、m₃……的个数﹤2，则对地形图像的灰度极小值加1，返回步骤1；

（4）如果n>2且2≤m₁、m₂、m₃……之一﹤n，说明n个灰度极小值不在一个连通区域内，则对灰度极小值所在的连通区域的边界均匀地进行膨胀，使得多个连通区域彼此连通，则膨胀前边界的像素点的连线即为分割线，由该分割线在地形图像中分割出区域为第一分割区域、第二分割区域、第三分割区域……；

（5）标注出第一分割区域、第二分割区域、第三分割区域……中每个像素的灰度值，并且将上述每个区域中灰度值相同的点就近连接起来，分别画出上述第一分割区域、第二分割区域、第三分割区域……的等高线折线；

（6）根据第一分割区域、第二分割区域、第三分割区域……中的每根等高线折线上的像素坐标拟合出每根等高线折线的等高线曲线方程，根据上述每个等高线曲线方程重新画出第一分割区域、第二分割区域、第三分割区域……中的每根等高线曲线，然后，根据每根等高线曲线y，由公式计算出每根等高线曲线上各个像素坐标点的曲率；

（7）分别将第一分割区域、第二分割区域、第三分割区域……中每个区域内的等高线曲率由大到小进行排序，找出前k*2位大曲率所对应于第一分割区域、第二分割区域、第三分割区域……中特征点，其中k分别为第一分割区域、第二分割区域、第三分割区域……中等高线曲线的个数；

（8）分别将第一分割区域、第二分割区域、第三分割区域……中的特征点按照就近的原则连接成折线，该折线即为第一分割区域、第二分割区域、第三分割区域……山脊线；

本发明的地形图像山脊线的提取方法，其中：对灰度极小值所在的连通区域的边界均匀地进行膨胀是指出膨胀后的边界与膨胀前的边界之间等距；

本发明的地形图像山脊线的提取方法，其中：在步骤（6）中，根据至少5个像素坐标拟合出每根等高线折线的等高线曲线方程；

本发明的地形图像山脊线的提取方法，其中：在步骤（6）中，等高线折线的等高线曲线方程由以下步骤拟合出来：

（a）. 设拟合多项式为：

（b） 找出等高线折线上的各点到这条曲线的距离之和，即偏差平方和如下：

（c）为了求得符合条件的a值，对等式右边求a_i偏导数，得到了：

（d）将等式左边进行一下化简，然后应该可以得到下面的等式：

(f) 把这些等式表示成矩阵的形式，就可以得到下面的矩阵：

（e）将这个范德蒙得矩阵化简后可得到:

(g)上式可以表示为X*A=Y，那么A=(X*Y^T)^-1X^TY，便得到了系数矩阵A，同时，得到了拟合的等高线曲线方程。

本发明的地形图像山脊线的提取方法与现有的方法相比，现有的等高线骨架法视地形特征线两侧的地形为对称变化，存在的问题是等高线有噪声时提取错误，而现有的等高线曲率最大判断法使用全局阈值进行判断，全局阈值过大时会丢失某些方向的山脊线的点；全局阈值过小时会增加噪点。本发明方法是在分水岭算法的基础上利用等高线曲率最大判别法对山脊线进行提取，消除了误提取、噪点增加等问题，提高了提取准确度。

附图说明

图1为在本发明的地形图像山脊线的提取方法中，找出分割线、第一分割区域和第二分割区域的示意图；

图2为在本发明的地形图像山脊线的提取方法中，找出第一分割区域的等高折线的放大示意图；

图3为在本发明的地形图像山脊线的提取方法中，找出第一分割区域的等高曲线的放大示意图。

在图1至图3中，标号1为分割线；标号2为第一分割区域；标号3为第二分割区域；标号4为等高线折线；标号5为等高线曲线；标号6为特征点；标号7为山脊线。

具体实施方式

一种地形图像山脊线的提取方法，该地形图像包括至少一个连通区域，它包括以下步骤：

（1）如图1所示，对拍摄的地形图像中每个像素的灰度值进行排序，找出灰度极小值，并且统计地形图像的灰度值极小值个数n和每个连通区域的灰度极小值个数m₁、m₂、m₃……,n=m₁+m₂+m₃……；

（2）如果n<2，则对地形图像的灰度极小值加1，返回步骤1；

（3）如果n>2，但是m₁、m₂、m₃……的个数﹤2，则对地形图像的灰度极小值加1，返回步骤1；

（4）如果n>2且2≤m₁、m₂、m₃……之一﹤n，最后找到n=7,m₁=3,m₂=4说明n个灰度极小值不在一个连通区域内，则对灰度极小值所在的连通区域的边界均匀地进行膨胀，对灰度极小值所在的连通区域的边界均匀地进行膨胀是指出膨胀后的边界与膨胀前的边界之间等距，使得多个连通区域彼此连通，则膨胀前边界的像素点的连线即为分割线1，由该分割线在地形图像中分割出区域为第一分割区2域和第二分割区域3；

（5）如图2所示，标注出第一分割区域2和第二分割区域3中每个像素的灰度值，并且将上述每个区域中灰度值相同的点就近连接起来，分别画出上述第一分割区域和第二分割区域等高线折线4，在图2中仅画出了第一分割区域2的等高线折线4；

（6）根据第一分割区域、第二分割区域、第三分割区域……中的每根等高线折线4上的像素坐标拟合出每根等高线折线的等高线曲线方程，根据上述每个等高线曲线方程重新画出第一分割区域和第二分割区域中的每根等高线曲线5，然后，根据每根等高线曲线y，由公式计算出每根等高线上各个像素坐标点的曲率；

（7）如图3所示，分别将第一分割区域和二分割区域中每个区域内的等高线曲率由大到小进行排序，找出前k*2位大曲率所对应于第一分割区域和第二分割区域中特征点6，其中k分别为第一分割区域和第二分割区域中等高线曲线的个数，在图3仅画出了3根等高线曲线5；

（8）分别将第一分割区域和第二分割区域中的特征点6按照就近的原则连接成折线，该折线即为第一分割区域好和第二分割区域山脊线7。

在步骤（6）中，根据至少取5个像素坐标拟合出每根等高线折线的等高线曲线方程，等高线折线的等高线曲线方程由以下步骤拟合出来：

（a）设拟合多项式为：

（b） 找出等高线折线上的各点到这条曲线的距离之和，即偏差平方和如下：

（c）为了求得符合条件的a值，对等式右边求a_i偏导数，得到了：

（d）将等式左边进行一下化简，然后应该可以得到下面的等式：

(f) 把这些等式表示成矩阵的形式，就可以得到下面的矩阵：

（e）将这个范德蒙得矩阵化简后可得到:

(g)上式可以表示为X*A=Y，那么A=(X*X^T)^-1X^TY，便得到了系数矩阵A，同时，得到了拟合的等高线曲线方程。特别需要说明的是：上述的多项式拟合公式来源于由哈尔滨工程大学出版社出版的《误差理论与曲线拟合》。

以上描述是对本发明的解释，不是对发明的限定，本发明所限定的范围参见权利要求，在不违背本发明的精神的情况下，本发明可以作任何形式的修改。