技术领域本发明涉及一种黄河主溜信息提取方法，特别是涉及一种基于多光谱遥感图像的黄河主溜信息协同提取方法。
背景技术：
遥感图像已经广泛应用于水体识别、河流提取、水质检测、洪灾检测、水域变化检测等领域，这些研究采用的方法基本上是根据水体在图像上所表现出的特性并使用传统的一些分类方法进行检测。黄河主溜信息的提取在传统上都是通过人工绘制，该方法不仅费时、费力，而且容易受天气等自然条件的影响，更重要的是，难以在汛期及时掌握主溜变化情况。遥感技术和遥感图像处理的不断发展，使得运用遥感图像进行黄河主溜线检测成为可能。河道主溜在河势多变的淤积沙床质河道中表现的问题突出，由于国际上对多沙河流河道河势研究较少。文献“Main-streamFeatureExtractionoftheYellowRiverBasedonRegionalSpectralUn-mixing,4thInternationalYellowRiverForum,2009,pp212-219”提出了一中基于局部光谱解混的主溜特征提取方法，利用光谱向量在高维空间中的几何属性进行解混，求取黄河主溜区域的端元，从而提取出典型的主溜信息。然而该方法未能充分利用遥感图像的空间信息，且单种方法很难获得多个方面的鲁棒特性，因而难以提取出连续性好、准确率高的主溜信息。
技术实现要素：
为了克服现有黄河主溜信息提取方法准确率差的不足，本发明提供一种基于多光谱遥感图像的黄河主溜信息协同提取方法。该方法首先由偏度分析方法和基于光谱相似性与空间连续性的方法分别获得一组主溜点位置；然后将偏度分析方法得到的主溜点作为先验信息输入到基于光谱相似性与空间连续性的提取方法中，用于获取矫正后的主溜点；此后，将所获得的主溜点再作为偏度分析方法的输入，进一步矫正主溜点；如此迭代，直至获得的主溜点位置稳定；最后，将得到的主溜点进行连接，形成主溜线。本发明公开的方法在黄河中下游的LandsatTM遥感图像上的实验结果表明，其检测率达到95％时，虚警率仅为1％，且获得的主溜线具有很好的连续性。本发明解决其技术问题所采用的技术方案是：一种基于多光谱遥感图像的黄河主溜信息协同提取方法，其特点是包括以下步骤：(1)河流分割。输入一幅LandsatTM多光谱图像，首先进行河流粗分割。采用光谱分类和匹配方法进行光谱图像分类，并根据黄河河段的形状特点进行分类后处理，采用两种样本分别对图像进行分类，然后合成图像。由于图像中有桥梁和滩涂的影响。(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)中的像元数小于Y(0)中的像元数，偏度分析窗口大小减小1；反之，偏度分析窗口大小增大1；由调整参数后的偏度分析方法，得到其输出的主溜点集合Y(1)。同理，将Y(0)用于指导光谱相似性与空间连续性方法的参数调整，若Z(0)与Y(0)的一致性L(Y(0),Z(0))小于给定阈值，则结束协同过程；否则，若Z(0)中的像元数小于Y(0)中的像元数，th1增加2％；反之，th1减小2％；由调整参数后的偏度分析方法，得到其输出的主溜点集合Z(1)。重复上述协同过程，直至算法因L(Y(t),Z(t))小于给定阈值th2而结束，t为迭代代数。本发明的有益效果是：该方法首先由偏度分析方法和基于光谱相似性与空间连续性的方法分别获得一组主溜点位置；然后将偏度分析方法得到的主溜点作为先验信息输入到基于光谱相似性与空间连续性的提取方法中，用于获取矫正后的主溜点；此后，将所获得的主溜点再作为偏度分析方法的输入，进一步矫正主溜点；如此迭代，直至获得的主溜点位置稳定；最后，将得到的主溜点进行连接，形成主溜线。本发明公开的方法在黄河中下游的LandsatTM遥感图像上的实验结果表明，其检测率达到95％时，虚警率仅为1％，且获得的主溜线具有很好的连续性。下面结合具体实施方式对本发明作详细说明。具体实施方式本发明基于多光谱遥感图像的黄河主溜信息协同提取方法具体步骤如下：(7)河流分割。输入一幅LandsatTM多光谱图像，首先进行河流粗分割。采用光谱分类和匹配技术，如：光谱向量匹配、马氏距离分割和GaussMarkov等方法进行光谱图像分类，并根据黄河河段的形状特点进行分类后处理，如：区域合并等。分类时采用了监督分类方法，由于图像受天气等因素的影响，图像中黄河上游和下游水体的光谱有较大差异，若采用单一样本进行光谱角分类会有很大误差，因此采用两种样本分别对图像进行分类，然后合成图像。由于图像中有桥梁和滩涂的影响，不便于南北两岸水边线的提取，在边缘精度要求适合进一步研究的情况下对图像进行了形态学处理，即图像进行膨胀和腐蚀操作，消除了较小的滩涂和桥梁。(8)河岸线提取。采用Canny算子对边缘进行检测，并对检测出的边缘进行跟踪连接，分类出河流南北岸的水边线，为下一步河流分段做好准备。为了便于提取南北两岸，需对所得到的边缘线进行连接跟踪；根据邻域法，利用一个边缘跟踪算法可到得到一组连续的线段，根据长度去除一些较短的干扰线段，同时可根据滩涂的形状特性，去除滩涂，从而得到一组有用线段。根据黄河水体的分布特征，观察它具有一定的横向和竖向特征，据此可设计算法对每段线段进行南北岸的判断。基于算法效率考虑只需对每段中某些有代表性的点进行南北岸的判断即可得到这整段水边线所属哪一岸的信息。根据统计特性，每段中所判断的点属于某一岸的数目大于某一阈值即可判断这段是属于南岸或北岸，反之，则属于另外一岸，从而得到南北岸图像。根据需要把南北两岸的线段按一定顺序分别进行有序存储到矩阵中，从而得到两岸完整的水边线，以利于下一步的主溜线提取。(9)河流分段。就黄河主溜线问题将黄河分成典型河段和非典型河段。典型河段又分为顺直微弯、弯曲和分汊三类。每一类河段的形状都有很大差别，描述差别的主要方式一是利用弯曲系数和曲率来描述河段的弯曲程度及方向，二是通过南北水边线之间的河流分布来确定是否存在分汊。河段与河段之间的空间联系也是判定的黄河主溜线依据之一，因此，对于河段之间的空间分布关系计算也是非常重要的，本发明对空间分布关系定义为连续弯相接、弯道归顺等河段，利用河段之间的连续变化排列从空间予以分别。具体分段方法如下：a)通过开窗的曲率计算将空间的坝岸变换到曲率域，所得曲率序列的极大值点位置，代表了弯道弯顶的位置所在。b)两个连续极值点之间的极小值点位置，代表两连续河湾间的过渡衔接点位置。(10)利用偏度分析提取主溜初始位置。利用偏度分析提取主溜信息，获取一组初始主溜点Y(0)。偏度分析算法具体步骤如下：Step1获取p×n观测矩阵数据X＝[X1,X2,…,Xn]，其中每一列Xi代表一个观测样本向量，每一行代表一个观测属性；两类样本数据S1,S2；Step2求取X的平均偏差形式B：令则即将坐标轴平移到原数据的中心位置，求S1,S2的样本均值M1，M2；Step3求类间散布矩阵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列的纵坐标平均，作为主溜点位置。(11)利用光谱相似性与空间连续性方法提取主溜初始位置。利用偏度分析提取主溜信息，获取一组初始主溜点Z(0)。光谱相似性与空间连续性方法具体步骤如下：首先通过对人工实地勘测到的主溜进行大量学习，建立主溜光谱库。选取河流入口处与主溜光谱库中最相似的点，作为初始主溜样本。由初始主溜点开始，计算其方向线，并得到其长度和位置，具体步骤如下：Stepl定义方向线为穿过中心像元的一系列线段，它们的长度各不相同，其长度由相邻像元间的光谱相似性测度和阈值来确定。相似性测度为：sim(xi,x′)=-12||xi-x′||Σs2---(b)]]>其中，x′是xi的相邻像元，表示马氏距离。Step2方向线按照特定的规则从中心像元出发朝两边扩展，同时将被选入的新像元加入结果结合。方向线扩展的条件是当前像元的相似性测度小于阈值th1。这里，th1为光谱相似性度阈值，它与同一形状区域内像元灰度的变化程度相关。由于水流流动具有一定的方向性，实际的方向线扩展不可能向着所有方向进行，而是根据河流方向，将与河流方向背离的方向去除，再进行扩展。Step3遍历整条河流，按照Stepl和Step2两步可以跟踪得到初始像元的方向线，同时获得形成该方向线时逐步形成的像素点集合Z(0)。(12)迭代实现主溜信息协同提取。将步骤(5)中得到的主溜点集合Z(0)作为指导信息输入偏度分析方法中，对其内部参数进行调整，若Z(0)与Y(0)的一致性L(Y(0),Z(0))小于给定阈值th2，则结束协同过程；否则，若Z(0)中的像元数小于Y(0)中的像元数，偏度分析窗口大小减小1；反之，偏度分析窗口大小增大1；由调整参数后的偏度分析方法，得到其输出的主溜点集合Y(1)。同理，将Y(0)用于指导光谱相似性与空间连续性方法的参数调整，若Z(0)与Y(0)的一致性L(Y(0),Z(0))小于给定阈值，则结束协同过程；否则，若Z(0)中的像元数小于Y(0)中的像元数，th1增加2％；反之，th1减小2％；由调整参数后的偏度分析方法，得到其输出的主溜点集合Z(1)。重复上述协同过程，直至算法因L(Y(t),Z(t))小于给定阈值th2而结束(t为迭代代数)，此时获得的检测结果为两种算法经协同后趋向一致的结果，具有较高可靠性。