一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法与流程

本发明涉及到遥感图像的分割领域，尤其涉及一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法。

背景技术：

图像分割(章毓晋.图像分割.北京：科学出版社,2005.)是图像处理与计算机视觉领域中最为基础和重要的问题之一，也是对图像进行视觉分析和模式识别的基本前提。高分辨遥感图像以其丰富的光谱信息，包含着大量的地物信息特征。多尺度分割算法(＝multiresolutionsegmentation、均值漂移算法)已被广泛应用到高分辨遥感图像分割中。多尺度分割算法中常常需要预先设定一个尺度参数，用来控制分割对象的平均大小，尺度参数的选取对分割效果有着重要影响。因此，尺度参数的选取成为多尺度分割优化的一个重要问题。目前，大部分多尺度分割优化方法集中于从多个尺度的分割结果中选取一个单独的尺度作为最优化的参数。某一尺度下的分割结果的好坏通常由分割对象间的差异性来度量，常用欧式距离、光谱角等评价指标来度量这种差异性，最后通过这个评价指标来选取一个全局最优的尺度参数。然而，这种方法往往难以兼顾高分辨遥感图像中不同地物对分割尺度的需求，极大限制了分割效果。例如，大面积的海域、连接的养殖池塘、规则的农田等不同的地物环境，便不适用同一个分割尺度。此外，全局最优尺度的选取也忽略了其余尺度参数下的分割信息，没有充分利用。

传统的遥感图像分割优化算法仅在多个分割尺度的结果中选取一个单一尺度作为最优分割结果，没有充分利用到多个分割结果中的尺度信息，依然不能保证对遥感图像中不同地物提供适应的分割尺度。

技术实现要素：

针对现有技术的不足，本发明提出一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法，该方法针对多尺度分割算法只能设置一个尺度参数的不足，自适应在不同尺度参数下的分割结果中为相应地物选取最优分割结果，可以更准确地描绘遥感图像中不同环境下的地物边缘。

为达到以上目的，本发明具体技术方案如下：

一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法，其特征在于，该方法包括如下步骤：

(1)设置分割尺度参数范围[1,2,…,l]，使用多尺度分割算法对待分割的高分辨遥感影像进行不同尺度的分割，获得l个不同尺度下的分割结果，作为待优化的初始分割结果；

(2)对每一个尺度下的初始分割结果建立分割对象邻接图rag，用以表示相邻分割对象的空间关系，所述的分割对象为一幅遥感图像分割后划分的子区域，这些子区域组成遥感图像的分割结果；

(3)选择任意一个坐标点，该坐标点所在的分割对象为o，分割尺度为t，t∈[1,2,…,l]，进行如下步骤：

(3.1)计算分割对象的局部光谱差异，包括分割对象间的光谱差异和分割对象内光谱差异，其中，分割对象间的光谱差异角计算如下：

计算分割对象o的光谱均值向量s，并将该分割对象o的相邻分割对象的光谱均值向量组成集合计算分割对象间光谱差异角计算公式如下：

其中，s∈r^b，b是该遥感图像的波段数，每一个xi都是o的相邻分割对象的光谱均值向量，xi∈r^b，u是相邻分割对象的个数；

分割对象内光谱差异角计算如下：

在小于尺度t的分割结果中搜索组成当前分割对象o的分割对象，组成分割对象o的每个分割对象的光谱均值向量的集合为ω，根据下式计算分割对象内光谱差异角

其中，p与q是小于当前尺度t下组成o的任意两个分割对象的光谱均值向量，r是集合ω内的对象数目；若在小于t的尺度下没有找到组成分割对象o的分割对象，则

(3.2)采用下式计算分割对象o的局部光谱差异值

(3.3)计算l个分割结果在当前坐标点的局部光谱差异值eo，形成光谱差异值的曲线，选择光谱差异值最大的尺度对应分割对象作为该坐标下的最优分割对象；

(4)对遥感图像上所有坐标位置重复步骤(3)，获得每个坐标下的最优分割对象；

(5)将最优分割对象组合为最终的优化后分割结果。

进一步地，所述光谱均值向量s计算公式如下：

其中，z为分割对象o内所有像素点的集合，zi∈r^b，m是像素的个数。

进一步地，所述的多尺度分割算法具体为multiresolutionsegmentation。

进一步地，所述的多尺度分割算法具体为均值漂移算法。

与现有技术相比，本发明的有益效果如下：

本发明中，通过对比同一坐标下不同尺度的分割对象的局部光谱值，来选取该坐标下的最优分割对象，可以自适应对不同地物设置相应的分割尺度。相比与传统尺度优化方法，本发明的优化结果能综合利用不同尺度下的分割信息，提高分割的准确率。

附图说明

图1是本发明的基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法的流程框图；

图2是本发明中分割对象光谱间和分割对象内差异值的计算示意图；

图3是本发明的其中一个实施例中的最优分割对象的选取示意图，其中，图3(a)-(e)分别是尺度50、100、150、350、400下该养殖鱼塘的分割结果图，图3(f)是目标养殖鱼塘的参考分割结果图，图3(g)为各个尺度下该位置的局部光谱值e的曲线图；

图4是本发明的其中一个实施例的初始对尺度分割结果与最终分割优化后结果对比图，图4(a)是该遥感图像，图4(b)-(g)分别是尺度70、140、200、250、300、400下的分割结果图，图4(h)是分割优化后的结果图。

具体实施方式

下面根据附图和优选实施例详细描述本发明，本发明的目的和效果将变得更加明白，以下结合附图和实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。

如图1所示，一种基于光谱差异最大化的高分辨遥感图像多尺度分割优化方法，其特征在于，该方法包括如下步骤：

(1)设置分割尺度参数范围[1,2,…,l]，使用多尺度分割算法对待分割的高分辨遥感影像进行不同尺度的分割，获得l个不同尺度下的分割结果，作为待优化的初始分割结果；

(2)对每一个尺度下的初始分割结果建立分割对象邻接图(region-adjacentgraph)rag，用以表示相邻分割对象的空间关系；所述的分割对象为一幅遥感图像分割后划分的子区域，这些子区域组成遥感图像的分割结果；

(3)选择任意一个坐标点，该坐标点所在的分割对象为o，分割尺度为t，t∈[1,2,…,l]，进行如下步骤：

(3.1)计算分割对象的局部光谱差异，包括分割对象间的光谱差异和分割对象内光谱差异，其中，

分割对象间的光谱差异角计算如下：

计算分割对象o的光谱均值向量s，并将该分割对象o的相邻分割对象的光谱均值组成集合计算分割对象间光谱差异角计算公式如下：

其中，s∈r^b，b是该遥感图像的波段数，每一个xi都是o的相邻分割对象的光谱均值向量，xi∈r^b，u是相邻分割对象的个数；如果o与周围相邻对象xi的光谱角越大，表示o与周围是同一种地物的可能性很小；反之，o越可能与周围对象属于同一种地物类别。

分割对象内光谱差异角计算如下：

在小于尺度t的分割结果中搜索组成当前分割对象o的分割对象，组成分割对象o的每个分割对象的光谱均值向量的集合为ω，根据下式计算分割对象内光谱差异角

其中，p与q是小于当前尺度t下组成o的任意两个分割对象的光谱均值向量，r是集合ω内的对象数目；若在小于t的尺度下没有找到组成分割对象o的分割对象，则类似的，公式(2)考虑了分割对象o内部光谱的差异性，如果这种内部的差异性越大，就表示o本身的可分割性比较大。

需要指出的是，在多尺度分割算法中，大尺度下的分割对象都是由小尺度下的分割对象合并而来，因此，在小于t的尺度下搜寻组成o的小分割对象时，不会出现分割对象冲突(相交)的情况，小的分割对象要么可以组合成o，要么没有组成o的对象(o就是该位置的最小分割)。

(3.2)采用下式计算分割对象o的局部光谱差异值

通过计算每一个分割对象的局部光谱差异值e，我们可以用来度量该分割对象所在尺度是否处于合适的分割尺度。综合分割对象间和分割对象内的光谱差异，如果e的值偏小，那么表示该分割对象的可分割性高，反之该分割对象可分割性较低；

(3.3)计算l个分割结果在当前坐标点的局部光谱差异值eo，形成光谱差异值的曲线，选择光谱差异值最大的尺度作为该坐标下的最优分割；

(4)对遥感图像上所有坐标位置重复步骤(3)，获得每个坐标下的最优分割对象；

(5)将最优分割对象组合为最终的优化后分割结果。

优选地，所述光谱均值向量s计算公式如下：

其中，z为分割对象o内所有像素点的集合，zi∈r^b，m是像素的个数。

优选地，所述的多尺度分割算法具体为multiresolutionsegmentation或为均值漂移算法。

综上所述，从优化后的分割结果来看，不同地物类别下的最优尺度可以通过本发明从不同的分割尺度中挑选出来。不同于传统的选取单个最优尺度作为优化结果，本发明直接综合了最优分割对象来形成最终结果。需要指出的是，本发明适用于各种多尺度分割算法。

图2示意了分割对象间和分割对象内光谱差异值的计算示意图。图中，o表示分割尺度t下的某一个分割对象，xi表示o周围的相邻对象，p和q表示在尺度小于t中找到的组合成o的小分割对象。从图2(a)可以看到，o的分割对象间差异考虑到了周围所有的分割对象，而图2(b)表示了分割对象o内光谱差异的计算方式。

图3示意了高分辨遥感图像内某一位置下的最优分割对象的选取过程图。如图所示，图3(a)-(e)分别是尺度50、100、150、350、400下该养殖鱼塘的分割结果图。图3(f)是手动绘制的参考分割结果，图3(g)示意了各个尺度下该位置的局部光谱值e的曲线图，通过选取该曲线上最高位置，即0.25所在位置，获得该位置上的最优尺度350。与图3(a)-(e)对比可以发现，350尺度下对养殖鱼塘的分割较好。

图4示意了整幅高分辨遥感影像不同尺度下的分割结果与分割优化后的结果图。如图所示，图4(a)是该遥感图像伪彩色图，图4(b)-(g)分别是尺度70、140、200、250、300、400下的分割结果图，图4(h)是分割优化后的结果图。从图中可以看出，图4(h)的分割结果融合了不同尺度下的分割对象，例如，左侧海域选取了尺度400的分割对象，保证了海域的完整性，而在中间的养殖鱼塘区域选取了中等尺度的分割对象，勾勒出各个鱼塘的边缘。

本领域普通技术人员可以理解，以上所述仅为发明的优选实例而已，并不用于限制发明，尽管参照前述实例对发明进行了详细的说明，对于本领域的技术人员来说，其依然可以对前述各实例记载的技术方案进行修改，或者对其中部分技术特征进行等同替换。凡在发明的精神和原则之内，所做的修改、等同替换等均应包含在发明的保护范围之内。