基于sne流形学习的河道主溜线检测方法
【技术领域】
[0001] 本发明设及一种河道主溜线检测方法，特别设及一种基于S肥流形学习的河道主 溜线检测方法。
【背景技术】
[0002] 主溜线是河道横断面上垂线流速最大的点的连线，是对河道主溜的概化表示，被 用来描述河势的变化情况。传统获取方法是通过人工坐船查河确定主溜位置，然后手工在 地形图上勾绘、转绘而成，该方法往往需要较长时间、安全性差，且难W获取洪水期的主溜 变化情况。当前，国外还没有基于遥感影像解译河道主溜线的应用研究成果报道，国内的研 究也处于刚刚起步阶段.
[0003] 文南犬('Main-stream of the Yellow River detection based on spectral fe曰ture 曰nd the Dynsmic Tr曰nsmission Model,Internstionsl Conference on Esrth Observation Data Processing and Analysis(ICEODPA),WuHan,China.Proc.of SPIE Vol.7285,72854Ul-8"公开了一种基于光谱相似性描述的河道主溜线检测方法。该方法是 将单个像素点每个波段的灰度值作为特征向量，利用光谱相似性度量和设置阔值的方法， 通过在上游主溜区选择一个已知主溜像素点来寻求其下游的主溜点位置，并对算法的检测 结果进行了主观评价。文献所述方法将图像处理技术和水流原理相结合，在一定程度上实 现了主溜线检测，但是由于河道主溜自身的特征受河道来水来沙、边界条件W及河型等多 因素影响，因此主溜自身的特征，包括主溜带深度、宽度、泥沙含量W及表面波浪的起伏度 等都存在很大的不确定性，仅利用单像素点的光谱值很难准确描述不同情形下的河道主溜 特征，且检测过程中往往还会受到噪声的影响，降低了主溜线的检测精度，在一些宽、浅、 散、乱，河道主溜特征不十分明显的河段检测算法失效。另外，算法中设置了相似性度量的 阔值，该阔值对整个检测过程有较大的影响，当阔值设置较大时，断点增多；当阔值设置较 小时，检测出的主溜点可能会偏离主溜区，但要确定合理的阔值却十分困难，算法的稳健性 较差。
【发明内容】
[0004] 为了克服现有河道主溜线检测方法检测精度差的不足，本发明提供一种基于S肥 流形学习的河道主溜线检测方法。该方法W单像素点为中屯、构建3X3大小的窗口，利用该 窗口中屯、点像素的光谱值、窗口灰度的统计值W及窗口梯度构建一个24维的特征向量来描 述主溜特征。利用流形学习将24维特征空间降至2维空间，然后在2维空间中寻求与已知主 溜像素点距离最近的点，将该点标记为下游的主溜位置，W此类推直至整个河段结束。该方 法利用窗口构建主溜的高维特征向量，囊括了主溜的光谱特征、纹理特征和几何特征;利用 流形学习算法进行维数约简，能够在构建的高维特征空间中寻求与已知主溜点距离最近的 点的低维实现。克服了主溜线检测中特征的不确定性运一技术难题，提高了主溜的检测精 度和稳健性。
[0005] 本发明解决其技术问题所采用的技术方案:一种基于S肥流形学习的河道主溜线 检测方法，其特点是包括W下步骤：
[0006] 步骤一、分割河流，生成河流分割影像。
[0007] 在一景TM遥感影像上，选择对水体最敏感的bands进行河流粗分割。利用区域生长 法，通过选择种子点像素和定义相似性度量值对河流进行增长，生成与TM遥感影像相同大 小的二值图像，在该二值图像中河流区域的值为1，非河流区域的值为0。然后利用生成的二 值图像对TM遥感影像进行掩膜处理，得到分割后的遥感影像，该遥感影像中河流部分的像 素值为TM遥感影像的光谱值，非河流部分全部为0。
[000引步骤二、提取河道中屯、线。
[0009] 在步骤一生成的二值图像基础上，利用形态学骨架化方法提取河流骨架，去除河 流骨架多余毛刺，得到一个完整的只有一个方向的河道骨架，即河道中屯、线。
[0010] 步骤=、选取初始主溜样本点。
[0011] W河流分割的多光谱遥感影像为基础，在河道上游的主溜区内人工选择一个像素 点作为初始主溜样本点，并将该初始主溜样本点标识为主溜点。
[001^ 步骤四、根据主溜样本点特征向量计算方法，计算主溜点的特征向量Xix24dZ = [S1,S2, ...S6,]il,li2,…，Ji6,〇2i,〇22,…，〇32,gi,g2 …，g6...〇其中，Si, i = l ,2, ...6 表示窗口中屯、 点在i波段的灰度值;Ui，i = 1，2，- '6表示窗口在i波段的灰度均值；巧2,1 = 1，2，- '6表示窗 口在i波段的灰度方差;gi，i = l，2,…，6表示在i波段上的窗口梯度。生成的主溜样本特征 向量T共24维。各参数的计算方法如下：
[OOU] 窗口的灰度均值ii:
[0014] 窗口的灰度方差 [001引窗口梯度g:
[0016] 式中，m ? n表不窗口的大小，设定为3 X 3
u(x,y)表不窗 口中屯、点在(x，y)的m ? n窗口的灰度均值。
[0017] 步骤五、计算主溜样本点的水流方向，利用动态演进模型，选择主溜样本点下游的 主溜候选点，并利用主溜样本特征计算方法，计算各候选点的特征向量YiX24,其中i表示主 溜候选点的个数，24表示特征向量的维数。
[0018] a.计算主溜样本点的水流方向。
[0019] 将河道中屯、线上与主溜点具有相同横坐标的点的切线方向作为该主溜点的水流 方向。其方法是在河道中屯、线上，利用距离该主溜点下游点和上游点3个像素的点位置进行 计算，用公式表示为：
[0020] d = 1:an( (y2-yi)/(X2-xi))
[002。 式中，d表示水流方向，y2、yi分别表示中屯、线上下游点、上游点的纵坐标，X2、xi表 示中屯、线上下游点、上游点的横坐标。
[0022] b.确定下游的主溜样本候选点。
[0023] 建立河道水流演进模型，该演进模型将水流方向简化为(0子斗f.-f.孚，-孚，-死)8 个方向。当上述计算的d值落在两个方向区间时，选择距离近的方向作为水流方向，当 ^ €脚言]时，令(1 = 0;当d e吃，言]时，令(1 =言。确定了主溜样本点的水流方向后，根据水 流演进模型确定主溜样本候选点。
[0024] C.计算主溜样本候选点特征向量Yi。
[0025] 主溜样本候选点特征向量Yi的计算方法与步骤四主溜样本点特征向量的计算方 法相同，计算时是W主溜候选点为中屯、像素点进行的计算。
[0026] 步骤六、利用S肥流形学习算法，对主溜样本点和主溜样本候选点同时进行降维处 理，得到主溜样本点和各候选点的2维特征向量，分别表示为X和yi，i表示第i个主溜候选 点。
[0027] 将主溜样本点和主溜样本候选点的特征向量合并成一个向量I，I = [X，Yi](w)x24， 然后将向量I(i+l)X24输入至S肥流形学习算法中。S肥流形学习算法确定的参数包括混乱度 Perp、迭代次数T、学习速率II、初始动量系数a(t)、最终动量系数丫（t)和输出维数。
[0028] 步骤屯、在主溜候选点中确定下游的主溜样本点。
[0029] 在降维后的2维特征空间中，计算主溜样本点X与各候选点Yi之间的距离，将距离 最小的候选点标识为主溜点，并将其作为新的主溜样本点。
[0030]
[0031] 式中X表示主溜样本点，Yi表示主溜样本候选点，X表示降维后主溜样本点的二维 特征向量，yi表示降维后主溜样本候选点的二维特征向量，n是候选点的个数，n等于5。
[0032] 步骤八、重复步骤四至步骤屯，直至整个河道检测结束，并将所有标记为主溜的点 显示出来，完成整个河段的主溜线检测。
[0033] 本发明的有益效果是:该方法W单像素点为中屯、构建3X3大小的窗口，利用该窗 口中屯、点像素的光谱值、窗口灰度的统计值W及窗口梯度构建一个24维的特征向量来描述 主溜特征。利用流形学习将24维特征空间降至2维空间，然后在2维空间中寻求与已知主溜 像素点距离最近的点，将该点标记为下游的主溜位置，W此类推直至整个河段结束。该方法 利用窗口构建主溜的高维特征向量，囊括了主溜的光谱特征、纹理特征和几何特征;利用流 形学习算法进行维数约简，能够在构建的高维特征空间中寻求与已知主溜点距离最近的点 的低维实现。克服了主溜线检测中特征不确定性运一技术难题，提高了主溜的检测精度和 稳健性。
[0034] 下面结合附图和【具体实施方式】对本发明作详细说明。
【附图说明】
[0035] 图1是本发明基于S肥流形学习的河道主溜线检测方法的流程图。
[0036] 图2是本发明方法建立的河道水流演进模型。
[0037] 图中，^表示水流方向，表示主溜样本点，醒隶示主溜样本候选点。
【具体实施方式】
[0038] 参照图1-2。本发明基于S肥流形学习的河道主溜线检测方法具体步骤如下：
[0039] 1.分割河流，生成河流分割影像。
[0040] 在一景TM遥感影像上，选择对水体最敏感的bands进行河流粗分割。利用区域生长 法，通过选择种子点像素和定义相似性度量值对河流进行增长，生成与原始影像相同大小 的二值图像，在该图像中河流区域的值为1，非河流区域的值为0。然后利用生成的二值图像 对原始影像进行掩膜处理，得到分割后的遥感影像，该影像中河流部分的像素值为原始影 像的光谱值，非河流部分全部为0。
[0041 ] 2.提取河道中屯、线。