基于生成式特征表示网络的多元遥感影像变化检测方法

技术领域

本发明属网络计算及表示学习技术领域，具体涉及一种基于生成式特征表示网络的多元遥感影像变化检测方法。

背景技术

变化检测是根据对同一物体或现象不同时间的观测来确定其不同的处理过程。实时掌握地物变化信息，有利于人类对现实存在问题的认知，以及更好的预测未来发展。变化检测可以分为以下两类：1)二元变化检测；2)多元变化检测。二元方法把各种变化看作一个变化类，因此，他们的目标是在考虑的特征空间中找出变化的和没有变化的像素，而忽略了可能发生的不同变化的语义。多元变化检测在检测出变化区域的基础上，还能对变化的类别进行甄别，更能满足现实的需求。

文献“F.Bovolo,S.Marchesi,and L.Bruzzone,“A framework for automaticand unsupervised detection of multiple changes in multitemporal images,”IEEETransactions on Geoscience and Remote Sensing,vol.50,no.6,pp.2196–2212,2012.”记载了变化向量分析(Compressed change vector analysis，C

发明内容

为了克服现有技术的不足，本发明提供一种基于生成式特征表示网络的多元遥感影像变化检测方法。设计了一种深度特征学习框架，该框架由超聚类和自适应循环学习机制驱动，学习不同类型变化的显著表示和分布。在整体框架中，采用了一种学习鉴别性差异特征的循环训练策略，该策略以自适应循环学习机制为指导，通过自适应方式优化变化类型的数量，并最终收敛到其稳定状态；引入变化向量分析算法(C

一种基于生成式特征表示网络的多元遥感影像变化检测方法，其特征在于步骤如下：

步骤1：给定两幅遥感影像I

步骤2：使用变化向量分析算法对训练样本进行分类处理，得到每个样本的初始标签；

步骤3：设置分类数目和最大迭代次数，初始迭代次数k＝1，按以下过程进行网络训练：

步骤a：分别将影像I

步骤b：以特征提取网络输出的遥感影像的高级特征作为噪声数据，以样本标签作为真实数据，输入到生成对抗网络GAN对其生成器和判别器进行训练，得到训练好的生成对抗网络GAN；特征提取网络和生成对抗网络共同构成生成特征表示网络；

步骤c：使用变化向量分析算法对特征提取网络输出的遥感影像的高级特征进行分类处理，以分类结果作为新的样本标签，更新训练数据；

步骤d：更新分类数目；

步骤e：令迭代次数k＝k+1，判断k是否达到最大迭代次数，如果不满足条件，则返回步骤a，进行下一次迭代训练，否则，此时的生成特征表示网络即为最优网络，其中，生成对抗网络的生成器输出即为最终的变化检测结果。

进一步地，步骤1中所述的以像素点为中心展开得到向量的具体过程为：设定以像素点为中心的大小为w×w的邻域，将此邻域中的像素点按行展开，得到由像素点值构成的w

进一步地，步骤3的步骤a中所述的特征提取网络的代价函数设置如下：

其中，L

训练过程中网络参数θ

其中，θ

进一步地，步骤3的步骤b中所述的对生成对抗网络GAN对其生成器和判别器进行训练的具体过程如下：

首先，设定训练参数，包括训练次数k、学习率α、截断阈值c和批数据规模m，然后，分别输入噪声样本和数据样本至生成对抗网络GAN进行训练；所述的噪声样本为从噪声数据中随机采样的m个样本；所述的数据样本为从真实数据中随机采样的m个样本；

训练过程中按以下公式更新生成对抗网络GAN判别器参数：

其中，

按以下公式更新生成器参数：

其中，

网络训练时按上述过程对判别器和生成器参数进行迭代更新，直至迭代次数达到设定的训练次数k，完成生成对抗网络GAN的生成器和判别器的优化，得到训练好的GAN网络。

进一步地，步骤3的步骤d所述的更新分类数目的具体过程如下：

首先，按照下式计算得到提升误差e：

其中，e

然后，按照下式计算得到更新后的分类数目M：

其中，λ∈[0,1]是自定义常数，M′表示更新前的分类数目，

本发明的有益效果是：(1)利用深度学习善于挖掘原始数据的潜在分布、能提取原始像素信息的关键性特征的优点，针对直接在像素空间进行差异分析而无法应对像素强度异质的弊端，先用深度学习提取像素的高维特征，然后在高维特征中实现差异分析，最后通过分类器将差异向量归类，实现最终的变化检测。相比于现有方法，这样一种的新的检测范式，一方面避免了原始像素值的直接比较，另一方面打破了传统变化检测框架三步走的壁垒，使得最终的变化检测结果不再依赖于差异图的生成质量。(2)在变化检测过程中，为进一步解决传统预检测算法无法提供相对精确的多类标签、深度网络的训练过程艰难的问题，引入变化向量分析算法(C

附图说明

图1是本发明的基于生成式特征表示网络的多元遥感影像变化检测方法框图；

图2是双时相多光谱遥感影像数据集Xi'an-2中的图像；

图中，(a)-拍摄于2013年8月19日的图像，(b)-拍摄于2015年8月29日的图像，(c)-多类别检测参考图；

图3是Xi'an-2数据集差异信息的极坐标表示；

图4是采用本发明方法对Xi'an-2数据集图像进行多类变化检测结果图像；

图中，(a)-20类检测结果图，(b)-16类检测结果图，(c)-10类检测结果图；(d)-5类检测结果图，(e)-2类检测结果图，(f)-检测参考图；

图5是采用本发明方法对Xi'an-2数据集进行处理的优化曲线；

图6是采用本发明方法对Xi'an-5数据集图像进行迁移变化检测结果图像；

图中，(a)-5类检测结果图，(b)-变化强度分布图，(c)-二分类变化检测结果图。

具体实施方式

下面结合附图和实施例对本发明进一步说明，本发明包括但不仅限于下述实施例。

本发明提供了一种基于生成式特征表示网络的多元遥感影像变化检测方法，如图1所示，其基本实现过程如下：

1、构建样本集

给定两幅遥感影像I

假设影像的大小为r×c，M

2、使用C

文献“Bovolo F,Marchesi S,Bruzzone L.A framework for automatic andunsupervised detection of multiple changes in multitemporal images[J].IEEETransactions on Geoscience and Remote Sensing,2012,50(6):2196-2212.”中记载了变化向量分析(Compressed change vector analysis，C

本发明采用C

3、网络训练

设置初始分类数目M

(1)优化特征提取网络

由于C

分别将影像I

为了实现特征提取的目标，设置网络的代价函数如下：

其中，L

以上代价函数可微，可进行优化。训练过程中网络参数θ

其中，θ

特征提取网络预训练好后，输入步骤2中的训练样本数据，在网络输出端得到影像数据高级特征。

(2)优化生成对抗网络GAN

生成对抗网络GAN由判别器和生成器两个部分构成，具体优化过程如下：

首先，设定训练参数，包括训练次数k、学习率α、截断阈值c和批数据规模m，如本实施例设定k＝5，α＝0.00005，c＝0.01，m＝64。输入训练数据集X，然后以特征提取网络的输出作为噪声数据P(z)，从中采样m个噪声样本{z

其中，

按以下公式更新生成器参数：

其中，

网络训练时按上述过程对判别器和生成器参数进行迭代更新，直至迭代次数达到设定的训练次数k，完成生成对抗网络GAN的生成器和判别器的优化，得到训练好的GAN网络。

(3)更新样本标签

使用C

(4)自适应更新分类数目

在训练过程中，初始分类数目M

首先，按照下式计算得到提升误差e：

其中，e

然后，按照下式计算得到更新后的分类数目M：

其中，λ∈[0,1]是自定义常数，M′表示更新前的分类数目，

(5)迭代

令迭代次数k＝k+1，判断k是否达到最大迭代次数E，如果不是，则返回步骤(1)，进行下一次迭代训练，否则，此时的特征提取网络即为最优网络，网络输出即为最终的变化检测结果。

为了验证算法的性能，在Intel(R)Core(TM)i5-3210M CPU 2.5GHz Windows 10系统下和Python3.6运行平台上进行实验验证。实验中采用双时相多光谱遥感影像数据集Xi'an-2，如图2所示，是从用GF-1卫星拍摄的大图(分辨率为2M)中截取的小图，分别拍摄于2013年8月19日(图(a))和2015年8月29日(图(b))，两幅截取的影像是配准好的。真实变化检测结果图是用两幅原始影像，图(c)为基于先验信息和摄像机判读技术得到的真实变化检测图(人工标注)，用于评估算法的检测结果。

图3给出了采用本发明方法对数据集Xi'an-2进行变化向量分析时得到的极坐标图，可以清晰地看出有5个突出的部分(如图中虚线圈所示)，表明该数据集具有5类变化。图4给出了设定不同聚类数目时，采用本发明方法进行变化检测的结果图像，当设定的聚类数目逐渐减少时，相似的类别会逐渐合成为同一类，当聚类数目接近真实变化数目时，多类变化检测的结果也更加接近真实分布，变化区域的轮廓也愈发清晰，证明了本发明自适应聚类数目调整策略的有效性。图5给出了采用本发明方法对数据集Xi'an-2进行变化检测时的优化过程，每幅图中，横轴表示迭代周期，纵轴表示代价函数值。由于整个优化过程采用迭代优化方式，每次更新分类数目后，优化在前一次迭代的基础上开始。但是由于训练标签改变，因此目标函数的值会出现跳跃。总的来说，在每个迭代周期内，目标函数值是逐步减少的。图6给出了数据迁移实验的结果，在此实验中将在数据集Xi'an-2上优化好的本发明模型应用到Xi'an-5数据集上，其中图(a)是5类变化检测结果图，图(b)是变化强度分布图，颜色越深表示越有可能是变化区域，图(c)则是二分类变化检测结果图。从结果可以看出，迁移结果能辨别并显示变化的区域和类型，由于Xi'an-5数据集仅有3类变化，而数据集Xi'an-2是5类变化，所以迁移后的检测结果是检测了5类的变化，因此还存在一些检测的错误，不过整体效果在一定程度上验证了本发明学习数据分布的有效性。

实验选择正确分类率(Percentage correct classification,PCC)评估算法性能好坏，其计算公式如式(21)所示；同时，为了更加精确地评估生成的变化检测结果图与真实参考图的接近程度，还引入Kappa系数，其计算公式如式(22)所示。

其中，M

分别采用C

表1