一种基于协同进化的高光谱图像波段选择方法及系统

技术领域

本发明涉及遥感图像处理技术领域，特别是涉及一种基于协同进化的高光谱图像波段选择方法及系统。

背景技术

近年来，随着遥感技术的发展，高光谱图像(HSI)得到了越来越广泛的应用。高光谱分辨率遥感(Hyperspectral Remote Sensing)是用窄而连续的光谱通道对地物持续遥感成像的技术。在可见光到短波红外波段，其光谱分辨率可达纳米(nm)数量级，通常具有波段多的特点。高光谱数据可以获得几乎连续的地物波谱，在获取信息的同时，给高光谱图像分类带来了巨大的计算量和冗余，在某些具体实际应用中，数据量的增加并没有增加原有的信息量。对高光谱图像进行波段选择是一种降低数据冗余，提高高光谱遥感影像分类的精度和速度的方法。

然而，目前为止还没有一种能够高效率的实现高光谱图像波段选择的方法。

发明内容

本发明的目的是提供一种基于协同进化的高光谱图像波段选择方法及系统，能够对高光谱遥感图像进行波段选择，提高图像分类与地物检测等实际应用的效率与精度，有效解决算法使用时的缺陷。

为实现上述目的，本发明提供了如下方案：

一种基于协同进化的高光谱图像波段选择方法，具体包括：

对协同进化水稻育种优化算法的种群进行初始化；

执行协同进化水稻育种优化算法对种群中每一种子的位置进行更新；

对每一种子的种子位置进行二进制编码，得到该种子的二进制序列；

利用所述二进制序列对高光谱图像进行波段选择，得到若干波段子集；

利用目标函数评价若干所述波段子集的适应度，得到当前次最优波段子集，并将所述当前次最优波段子集的适应度与预设阈值进行比较，若所述当前次最优波段子集的适应度满足预设阈值，则将所述当前次最优波段子集作为高光谱图像波段选择结果；若所述当前次最优波段子集的适应度不满足预设阈值，则返回所述“执行协同进化水稻育种优化算法对种群中每一种子的位置进行更新”步骤。

一种基于协同进化的高光谱图像波段选择系统，系统包括：

初始化模块，用于对协同进化水稻育种优化算法的种群进行初始化；

协同进化模块，用于执行协同进化水稻育种优化算法对种群中每一种子的位置进行更新；

波段选择模块，用于对每一种子的种子位置进行二进制编码，得到该种子的二进制序列；利用所述二进制序列对高光谱图像进行波段选择，得到若干波段子集；

终止判断模块，用于利用目标函数评价若干所述波段子集的适应度，得到当前次最优波段子集，并将所述当前次最优波段子集的适应度与预设阈值进行比较，若所述当前次最优波段子集的适应度满足预设阈值，则将所述当前次最优波段子集作为高光谱图像波段选择结果；若所述当前次最优波段子集的适应度不满足预设阈值，则返回所述“执行协同进化水稻育种优化算法对种群中每一种子的位置进行更新”步骤。

根据本发明提供的具体实施例，本发明公开了以下技术效果：

本发明提供了一种基于协同进化的高光谱图像波段选择方法及系统，将协同进化水稻育种优化算法与高光谱图像的波段选择相结合，能够在不损失重要信息的情况下从原始数据中尽可能的保留地物光谱辐射信息，剔除冗余和有噪声的波段影像、降低数据维度，提高图像分类与地物检测等实际应用的效率与精度，为遥感图像技术的发展提供了有力支撑。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例提供的一种基于协同进化的高光谱图像波段选择方法流程图；

图2为本发明实施例提供的协同进化过程中全局最优值交流原理图；

图3为本发明实施例提供的一种基于协同进化的高光谱图像波段选择系统框图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

由于波段选择问题是一个NP难题，到目前为止还不能完全解决。本发明的目的是提供一种基于协同进化的高光谱图像波段选择方法及系统，基于新近提出的一种性能优良的协同进化水稻育种优化算法，通过对种群进行协同进化，运用协同进化水稻育种优化算法并对最优个体进行柯西变异，使其性能能够进一步提升。

为使本发明的上述目的、特征和优点能够更加明显易懂，下面结合附图和具体实施方式对本发明作进一步详细的说明。

实施例1

如图1所示，本实施例提供了一种基于协同进化的高光谱图像波段选择方法，具体包括：

步骤101：对协同进化水稻育种优化算法的种群进行初始化；

步骤102：执行协同进化水稻育种优化算法对种群中每一种子的位置进行更新；

步骤103：对每一种子的种子位置进行二进制编码，得到该种子的二进制序列；利用所述二进制序列对高光谱图像进行波段选择，得到若干波段子集；

利用目标函数评价若干所述波段子集的适应度，得到当前次最优波段子集，并将所述当前次最优波段子集的适应度与预设阈值进行比较，若所述当前次最优波段子集的适应度满足预设阈值，则将所述当前次最优波段子集作为高光谱图像波段选择结果；若所述当前次最优波段子集的适应度不满足预设阈值，则返回所述“执行协同进化水稻育种优化算法对种群中每一种子的位置进行更新”步骤。

由此，本实施例提供的高光谱图像波段选择方法，结合协同进化水稻育种优化算法，能在不损失重要信息的情况下从原始数据中尽可能的保留地物光谱辐射信息，剔除冗余和有噪声的波段影像、降低数据维度，提高图像分类与地物检测等实际应用的效率与精度，有效解决算法使用时的缺陷。

协同进化思想为：将群体中的个体划分为多个子群，每个子群代表解空间中的一个子空间，所有子群共享初始的全局最优种子gbest，每隔R代，将M个子群中当前的gbest

基于协同进化思想，首先对协同进化水稻育种优化算法的种群规模、最大自交次数、最大迭代次数以及种群初始位置等进行初始化。具体包括：

其中，X为含有d维变量的水稻个体，X

为了便于种群划分，保证各个子种群的个体数目一致，从而提高协同进化的精确度，可将种群规模设置为子种群个数的整数倍，本实施例中具体设置种群大小为3的倍数。

初始化后，接着对种群进行划分，将种群划分为保持系子种群、不育系子种群和恢复系子种群。如图2所示，本实施例通过将种群整体划分为多个子种群，所有子种群并行展开搜索，将搜索到的较好个体在不同子种群间进行迁徙，对于第i个子种群，其全局最优位置gbest

为了提高算法的搜索效率，本实施例中，在对种群进行划分时，先根据个体的适应度值将种群中所有个体进行排序：

X

X

其中X

然后利用杂交操作对不育系子种群进行更新；利用自交操作对恢复系子种群进行更新。

对于不育系子种群的更新，本实施例采用杂交操作。每一次育种进化中，进行杂交的次数与不育系子种群中个体的数量相等，杂交的雄性亲本或雌性亲本都可以分别从不育系子种群和保持系子种群中选择，由父本和母本的基因重组产生具有新基因的新个体。选择的方式包括随机选择或对应选择。具体的：

初始化次数变量为零；

根据公式

其中，

将

将次数变量与所述不育系子种群中的个体数量进行比较，若所述次数变量小于所述不育系子种群中的个体数量，则返回“根据公式

对于恢复系子种群的更新，本实施例选择自交操作。依次遍历每一个恢复系子种群的个体，通过其与全局最优水稻个体及随机挑选的一恢复系个体的差值的结合，利用梯度下降法进行最优值搜索。

其中，对所述恢复系子种群中的某一个体进行更新包括：

根据X

其中，X

将X

为了避免恢复系子种群的更新陷入局部最优，对于恢复系子种群的每个个体设置了自交次数变量。当自交次数变量的值累加到超过预设自交阈值时，则对该个体执行重置操作并将该个体的自交次数置零。

重置操作如下：

X

这里X

在一次更新完成后，将种子位置进行二进制编码，利用编码得到的二进制序列对高光谱图像进行波段选择，得到若干波段子集。

在编码时，将种子的位置值转化为[0,1]之间的一个概率值。一般情况下，使用sigmoid函数将种子的基因值转换为概率值，依据概率值大小判定当前种子的位置值。但本实施例采用双曲正切函数代替sigmoid函数对位置值进行转化，获得了更高的适应度值。此时，每个种子的位置采用公式进行转化，每一种子得到一二进制序列：

其中，tanh()为双曲正切函数，t代表当前迭代次数。

然后利用二进制序列选择高光谱图像的波段子集，使所述二进制序列中的每一位分别对应所述高光谱图像的波段，选择与所述二进制序列中值为1的位对应的波段组成波段子集。例如，一副高光谱图像共包含10个波段，则可以采用“0100101010”形式进行编码，用对10个波段的光谱特征进行标记。由该编码形式可以看出，在所有波段中只有编号为2、5、7和9四个波段“被选择”组成波段子集，采用波段子集替代原有数据集，降低数据维度，避免“维数灾难”问题。

具体的，根据波段子集，通过支持向量机对高光谱图像数据进行分类，基于支持向量机的特征选择，每个种子的位置是一个二进制向量，其中每个数字对应一个原始波段。0表示对应的波段不被选择，1表示对应的波段被选择。得到与波段子集对应的高光谱图像数据组成波段子集图像数据。

为了进一步降低运算量，在对高光谱图像进行处理前，还包括对高光谱图像预处理的步骤，通过将高光谱图像数据归一化处理，将数据限制在所需要的一定范围内，去除冗杂数据，便于后续数据处理。

在得到各种子位置对应的波段子集后，根据以下适应度函数对波段子集进行评估：

其中fitness是适应度函数；Error是分类错误率，可以通过替代训练集或原始训练集来测量，

评估后，得到当前次迭代的最优波段子集，将与当前次迭代的最优波段子集对应的当前次最优种子位置通过全局多个子种群之间的最优值交流，使得各个子种群中的最优种子位置与当前次最优种子位置一致，然后对该当前次自由种子位置进行柯西变异：

p′

X

s(t+1)＝s(t)-(1/Iter

Iter

柯西分布具有较高的两翼概率特性，容易产生一个远离原点具有更宽的分布范围的随机数，本实施例通过柯西变异，使算法可以在更宽的范围进行搜索，同时降低种子陷入局部最优的危险。

对种子的最优位置进行柯西变异后，通过目标函数评价种子变异后的位置，将较优的位置作为全局最优种子位置赋给种子，让种子群体拥有一个较好的领导种子，可以有效提高算法的收敛精度。

在得到全局最优种子位置后，将与该全局最优种子位置对应的波段子集作为当前次最优波段子集，将当前次最优波段子集的适应度与预设阈值进行比较，若当前次最优波段子集的适应度满足预设阈值，则输出当前次最优波段子集及其对应的适应度值，将当前次最优波段子集作为高光谱图像波段选择结果；若当前次最优波段子集的适应度不满足预设阈值，则继续下一次更新，直至得到符合条件的波段子集或达到最大迭代次数。

实施例2

本实施例用于提供一种基于协同进化的高光谱图像波段选择系统，如图3所示，所述系统包括：

初始化模块M1，用于对协同进化水稻育种优化算法的种群进行初始化；

协同进化模块M2，用于执行协同进化水稻育种优化算法对种群中每一种子的位置进行更新；

波段选择模块M3，用于对每一种子的种子位置进行二进制编码，得到该种子的二进制序列；利用所述二进制序列对高光谱图像进行波段选择，得到若干波段子集；

终止判断模块M4，用于利用目标函数评价若干所述波段子集的适应度，得到当前次最优波段子集，并将所述当前次最优波段子集的适应度与预设阈值进行比较，若所述当前次最优波段子集的适应度满足预设阈值，则将所述当前次最优波段子集作为高光谱图像波段选择结果；若所述当前次最优波段子集的适应度不满足预设阈值，则返回所述“执行协同进化水稻育种优化算法对种群中每一种子的位置进行更新”步骤。

本说明书中每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想；同时，对于本领域的一般技术人员，依据本发明的思想，在具体实施方式及应用范围上均会有改变之处。综上所述，本说明书内容不应理解为对本发明的限制。