基于多光谱遥感图像的黄河主溜信息协同提取方法技术

本发明专利技术公开了一种基于多光谱遥感图像的黄河主溜信息协同提取方法，用于解决现有黄河主溜信息提取方法准确率差的技术问题。技术方案是首先由偏度分析方法和基于光谱相似性与空间连续性的方法分别获得一组主溜点位置；然后将偏度分析方法得到的主溜点作为先验信息输入到基于光谱相似性与空间连续性的提取方法中，用于获取矫正后的主溜点；此后，将所获得的主溜点再作为偏度分析方法的输入，进一步矫正主溜点；如此迭代，直至获得的主溜点位置稳定；最后，将得到的主溜点进行连接，形成主溜线。本发明专利技术方法在黄河中下游的Landsat TM遥感图像上的实验结果表明，其检测率达到95％时，虚警率仅为1％，且获得的主溜线具有很好的连续性。

【技术实现步骤摘要】

本专利技术涉及一种黄河主溜信息提取方法，特别是涉及一种基于多光谱遥感图像的黄河主溜信息协同提取方法。
技术介绍
遥感图像已经广泛应用于水体识别、河流提取、水质检测、洪灾检测、水域变化检测等领域，这些研究采用的方法基本上是根据水体在图像上所表现出的特性并使用传统的一些分类方法进行检测。黄河主溜信息的提取在传统上都是通过人工绘制，该方法不仅费时、费力，而且容易受天气等自然条件的影响，更重要的是，难以在汛期及时掌握主溜变化情况。遥感技术和遥感图像处理的不断发展，使得运用遥感图像进行黄河主溜线检测成为可能。河道主溜在河势多变的淤积沙床质河道中表现的问题突出，由于国际上对多沙河流河道河势研究较少。文献“Main-stream Feature Extraction of the Yellow River Based on Regional Spectral Un-mixing,4th International Yellow River Forum,2009,pp212-219”提出了一中基于局部光谱解混的主溜特征提取方法，利用光谱向量在高维空间中的几何属性进行解混，求取黄河主溜区域的端元，从而提取出典型的主溜信息。然而该方法未能充分利用遥感图像的空间信息，且单种方法很难获得多个方面的鲁棒特性，因而难以提取出连续性好、准确率高的主溜信息。
技术实现思路
为了克服现有黄河主溜信息提取方法准确率差的不足，本专利技术提供一种基于多光谱遥感图像的黄河主溜信息协同提取方法。该方法首先由偏度分析方法和基于光谱相似性与空间连续性的方法分别获得一组主溜点位置；然后将偏度分析方法得到的主溜点作为先验信息输入到基于光谱相似性与空间连续性的提取方法中，用于获取矫正后的主溜点；此后，将所获得的主溜点再作为偏度分析方法的输入，进一步矫正主溜点；如此迭代，直至获得的主溜点位置稳定；最后，将得到的主溜点进行连接，形成主溜线。本专利技术公开的方法在黄河中下游的Landsat TM遥感图像上的实验结果表明，其检测率达到95％时，虚警率仅为1％，且获得的主溜线具有很好的连续性。本专利技术解决其技术问题所采用的技术方案是：一种基于多光谱遥感图像的黄河主溜信息协同提取方法，其特点是包括以下步骤：(1)河流分割。输入一幅Landsat TM多光谱图像，首先进行河流粗分割。采用光谱分类和匹配方法进行光谱图像分类，并根据黄河河段的形状特点进行分类后处理，采用两种样本分别对图像进行分类，然后合成图像。由于图像中有桥梁和滩涂的影响。(2)河岸线提取。采用Canny算子对边缘进行检测，并对检测出的边缘进行跟踪连接，分类出河流南北岸的水边线。对所得到的边缘线进行连接跟踪；根据邻域法，利用一个边缘跟踪算法得到一组连续的线段，根据长度去除一些较短的干扰线段，同时根据滩涂的形状特性，去除滩涂，得到一组有用线段。根据黄河水体的分布特征，对每段线段进行南北岸的判断。基于算法效率考虑只需对每段中某些有代表性的点进行南北岸的判断即得到这整段水边线所属哪一岸的信息。根据统计特性，每段中所判断的点属于某一岸的数目大于某一阈值即判断这段是属于南岸或北岸，反之，则属于另外一岸，从而得到南北岸图像。根据需要把南北两岸的线段分别进行有序存储到矩阵中，从而得到两岸完整的水边线，进行主溜线提取。(3)河流分段。将黄河分成典型河段和非典型河段。典型河段分为顺直微弯、弯曲和分汊三类。利用弯曲系数和曲率描述河段的弯曲程度及方向，通过南北水边线之间的河流分布确定是否存在分汊。对空间分布关系定义为连续弯相接河段和弯道归顺河段，利用河段之间的连续变化排列从空间予以分别。具体分段方法如下：a)通过开窗的曲率计算将空间的坝岸变换到曲率域，所得曲率序列的极大值点位置，代表了弯道弯顶的位置所在。b)两个连续极值点之间的极小值点位置，代表两连续河湾间的过渡衔接点位置。(4)利用偏度分析提取主溜初始位置。利用偏度分析提取主溜信息，获取一组初始主溜点Y(0)。偏度分析算法具体步骤如下：Step1获取p×n观测矩阵数据X＝[X1,X2,…,Xn]，其中每一列Xi代表一个观测样本向量，每一行代表一个观测属性；两类样本数据S1,S2；Step2求取X的平均偏差形式B：令则即将坐标轴平移到原数据的中心位置，求S1,S2的样本均值M1，M2；Step3求类间散布矩阵Gb，类间散布矩阵Gb是p×p的半正定矩阵，定义为Gb＝(M1-M2)(M1-M2)T；Step4求取类间散布矩阵Gb的各个特征值evali与特征向量eigi，i∈[1,n]；Step5按特征值从大到小选择及与其对应得特征向量，构成变换矩阵T＝(eig1,eig2,…eigm)，m≤n；Step6生成新坐标系中的数据集H：H＝TT·X。对变换后的第一分量进行偏度分析： S k e w n e s s = 1 n - 1 Σ i = 1 n ( h i - h ‾ ) 3 / SD 3 - - - ( a ) ]]>其中，SD为标准差。Skewness＝0说明分布形态与正态偏度相同；Skewness>0，正偏，峰值在左；Skewness<0为负偏，峰值在右。|Skewness|越大，分布形态偏移程度越大。在某一断面的直方图集合中，取偏度系数最大点的位置作为当前断面上主溜线的位置。即用3×3窗口计算偏度系数，并记录影像每列偏度系数最大的点的纵坐标，然后将记录中相邻3列的纵坐标平均，作为主溜点位置。(5)利用光谱相似性与空间连续性方法提取主溜初始位置。利用偏度分析提取主溜信息，获取一组初始主溜点Z(0)。光谱相似性与空间连续性方法具体步骤如下：首先通过对人工实地勘测到的主溜进行大量学习，建立主溜光谱库。选取河流入口处与主溜光谱库中最相似的点，作为初始主溜样本。由初始主溜点开始，计算其方向线，并得到其长度和位置，具体步骤如下：Stepl定义方向线为穿过中心像元的一系列线段，它们的长度各不相同，其长度由相邻像元间的光谱相似性测度和阈值来确定。相似性测度为： s i m ( x i , x ′ ) = - 1 2 | | x i - x ′ 本文档来自技高网...

【技术保护点】
一种基于多光谱遥感图像的黄河主溜信息协同提取方法，其特征在于包括以下步骤：(1)河流分割。输入一幅Landsat TM多光谱图像，首先进行河流粗分割。采用光谱分类和匹配方法进行光谱图像分类，并根据黄河河段的形状特点进行分类后处理，采用两种样本分别对图像进行分类，然后合成图像。由于图像中有桥梁和滩涂的影响。(2)河岸线提取。采用Canny算子对边缘进行检测，并对检测出的边缘进行跟踪连接，分类出河流南北岸的水边线。对所得到的边缘线进行连接跟踪；根据邻域法，利用一个边缘跟踪算法得到一组连续的线段，根据长度去除一些较短的干扰线段，同时根据滩涂的形状特性，去除滩涂，得到一组有用线段。根据黄河水体的分布特征，对每段线段进行南北岸的判断。基于算法效率考虑只需对每段中某些有代表性的点进行南北岸的判断即得到这整段水边线所属哪一岸的信息。根据统计特性，每段中所判断的点属于某一岸的数目大于某一阈值即判断这段是属于南岸或北岸，反之，则属于另外一岸，从而得到南北岸图像。根据需要把南北两岸的线段分别进行有序存储到矩阵中，从而得到两岸完整的水边线，进行主溜线提取。(3)河流分段。将黄河分成典型河段和非典型河段。典型河段分为顺直微弯、弯曲和分汊三类。利用弯曲系数和曲率描述河段的弯曲程度及方向，通过南北水边线之间的河流分布确定是否存在分汊。对空间分布关系定义为连续弯相接河段和弯道归顺河段，利用河段之间的连续变化排列从空间予以分别。具体分段方法如下：a)通过开窗的曲率计算将空间的坝岸变换到曲率域，所得曲率序列的极大值点位置，代表了弯道弯顶的位置所在。b)两个连续极值点之间的极小值点位置，代表两连续河湾间的过渡衔接点位置。(4)利用偏度分析提取主溜初始位置。利用偏度分析提取主溜信息，获取一组初始主溜点Y(0)。偏度分析算法具体步骤如下：Step1获取p×n观测矩阵数据X＝[X1,X2,…,Xn]，其中每一列Xi代表一个观测样本向量，每一行代表一个观测属性；两类样本数据S1,S2；Step2求取X的平均偏差形式B：令则即将坐标轴平移到原数据的中心位置，求S1,S2的样本均值M1，M2；Step3求类间散布矩阵Gb，类间散布矩阵Gb是p×p的半正定矩阵，定义为Gb＝(M1‑M2)(M1‑M2)T；Step4求取类间散布矩阵Gb的各个特征值evali与特征向量eigi，i∈[1,n]；Step5按特征值从大到小选择及与其对应得特征向量，构成变换矩阵T＝(eig1,eig2,…eigm)，m≤n；Step6生成新坐标系中的数据集H：H＝TT·X；对变换后的第一分量进行偏度分析：Skewness=1n-1Σi=1n(hi-h‾)3/SD3---(a)]]>其中，SD为标准差；Skewness＝0说明分布形态与正态偏度相同；Skewness>0，正偏，峰值在左；Skewness<0为负偏，峰值在右；|Skewness|越大，分布形态偏移程度越大；在某一断面的直方图集合中，取偏度系数最大点的位置作为当前断面上主溜线的位置；即用3×3窗口计算偏度系数，并记录影像每列偏度系数最大的点的纵坐标，然后将记录中相邻3列的纵坐标平均，作为主溜点位置；(5)利用光谱相似性与空间连续性方法提取主溜初始位置；利用偏度分析提取主溜信息，获取一组初始主溜点Z(0)；光谱相似性与空间连续性方法具体步骤如下：首先通过对人工实地勘测到的主溜进行大量学习，建立主溜光谱库；选取河流入口处与主溜光谱库中最相似的点，作为初始主溜样本；由初始主溜点开始，计算其方向线，并得到其长度和位置，具体步骤如下：Stepl定义方向线为穿过中心像元的一系列线段，它们的长度各不相同，其长度由相邻像元间的光谱相似性测度和阈值来确定；相似性测度为：sim(xi,x′)=-12||xi-x′||Σs2---(b)]]>其中，x′是xi的相邻像元，表示马氏距离；Step2方向线按照特定的规则从中心像元出发朝两边扩展，同时将被选入的新像元加入结果结合；方向线扩展的条件是当前像元的相似性测度小于阈值th1；这里，th1为光谱相似性度阈值，它与同一形状区域内像元灰度的变化程度相关；由于水流流动具有一定的方向性，实际的方向线扩展不可能向着所有方向进行，而是根据河流方向，将与河流方向背离的方向去除，再进行扩展；Step3遍历整条河流，按照Stepl和Step2两步跟踪得到初始像元的方向线，同时获得形成该方向线时逐步形成的像素点集合Z(0)；(6)迭代实现主溜信息协同提取；将步骤(5)中得到的主溜点集合Z(0)作为指导信息输入偏度分析方法中，对其内部参数进行调整，若Z(0)与Y(0)的一致性L(Y(0),Z(0))小于给定阈值th2，则结束协同过程；否则，若Z(0)中的像元数...

【技术特征摘要】
1.一种基于多光谱遥感图像的黄河主溜信息协同提取方法，其特征在于包括以下步骤：(1)河流分割。输入一幅Landsat TM多光谱图像，首先进行河流粗分割。采用光谱分类和匹配方法进行光谱图像分类，并根据黄河河段的形状特点进行分类后处理，采用两种样本分别对图像进行分类，然后合成图像。由于图像中有桥梁和滩涂的影响。(2)河岸线提取。采用Canny算子对边缘进行检测，并对检测出的边缘进行跟踪连接，分类出河流南北岸的水边线。对所得到的边缘线进行连接跟踪；根据邻域法，利用一个边缘跟踪算法得到一组连续的线段，根据长度去除一些较短的干扰线段，同时根据滩涂的形状特性，去除滩涂，得到一组有用线段。根据黄河水体的分布特征，对每段线段进行南北岸的判断。基于算法效率考虑只需对每段中某些有代表性的点进行南北岸的判断即得到这整段水边线所属哪一岸的信息。根据统计特性，每段中所判断的点属于某一岸的数目大于某一阈值即判断这段是属于南岸或北岸，反之，则属于另外一岸，从而得到南北岸图像。根据需要把南北两岸的线段分别进行有序存储到矩阵中，从而得到两岸完整的水边线，进行主溜线提取。(3)河流分段。将黄河分成典型河段和非典型河段。典型河段分为顺直微弯、弯曲和分汊三类。利用弯曲系数和曲率描述河段的弯曲程度及方向，通过南北水边线之间的河流分布确定是否存在分汊。对空间分布关系定义为连续弯相接河段和弯道归顺河段，利用河段之间的连续变化排列从空间予以分别。具体分段方法如下：a)通过开窗的曲率计算将空间的坝岸变换到曲率域，所得曲率序列的极大值点位置，代表了弯道弯顶的位置所在。b)两个连续极值点之间的极小值点位置，代表两连续河湾间的过渡衔接点位置。(4)利用偏度分析提取主溜初始位置。利用偏度分析提取主溜信息，获取一组初始主溜点Y(0)。偏度分析算法具体步骤如下：Step1获取p×n观测矩阵数据X＝[X1,X2,…,Xn]，其中每一列Xi代表一个观测样本
\t向量，每一行代表一个观测属性；两类样本数据S1,S2；Step2求取X的平均偏差形式B：令则即将坐标轴平移到原数据的中心位置，求S1,S2的样本均值M1，M2；Step3求类间散布矩阵Gb，类间散布矩阵Gb是p×p的半正定矩阵，定义为Gb＝(M1-M2)(M1-M2)T；Step4求取类间散布矩阵Gb的各个特征值evali与特征向量eigi，i∈[1,n]；Step5按特征值从大到小选择及与其对应得特征向量，构成变换矩阵T＝(eig1,eig2,…eigm)，m≤n；Step6生成新坐标系中的数据集H：H＝TT·X；对变换后的第一分量进行偏度分析： S k e w n e s s = 1 n - 1 Σ i = 1 n ( h i - h ‾ ) 3 / SD 3 - ...

【专利技术属性】
技术研发人员：张艳宁，佘红伟，严杭琦，张磊，魏巍，刘学工，季万才，
申请(专利权)人：西北工业大学，
类型：发明
国别省市：陕西;61