一种基于sae的遥感图像变化检测方法
【技术领域】
[0001]本发明属于遥感图像变化检测方法，具体涉及一种基于SAE的遥感图像变化检测 方法。
【背景技术】
[0002] 遥感图像变化检测的通常技术是先将两图像进行差分，再将差分图以像素点为单 位进行二分类;或者先对图像进行分割，再以分割后的对象作为基本单位对分割后的图像 进行二分类。在现有技术中，提取哪些特征来进行分类，通常由人为设定，往往会忽略图像 本身含有的许多信息，导致检测结果不够精确。
[0003] 深度神经网络通过逐层训练的方式自主提取图像特征，有效地解决了人为提取特 征过程中信息遗漏的问题。但深度神经网络是一种有监督的学习模型，而遥感图像变化检 测问题是一个无监督的图像分类问题，没有训练数据，也没有类别标签，因此这两种技术通 常被认为很难结合使用。
【发明内容】
[0004] 要解决的技术问题
[0005] 为了避免现有技术的不足之处，本发明提出一种基于SAE的遥感图像变化检测方 法，把深度神经网络技术运用到遥感图像变化检测方面以解决检测精确度不高的问题。
[0006] 技术方案
[0007] -种基于SAE的遥感图像变化检测方法，其特征在于步骤如下：
[0008] 步骤1:对不同时刻同一地点获取的两幅待检测遥感图像进行配准、辐射校正及几 何校正;对校正后的图像数据进行原始特征提取，提取方法是:提取以每个坐标点为中心3* 3的图像块，组成一个18维的向量，并进行归一化;再加入2维的图像点归一化的坐标信息， 拓展成为一个20维向量，作为一个检测样本点的原始特征；
[0009] 步骤2:以提取到的原始特征作为输入对堆叠自编码器SAE进行无监督的训练；
[0010] 步骤3:计算两幅待检测图像的差分图像，对差分图像使用基于同质区域的马尔科 夫随机场算法得到粗变化检测结果图；
[0011]如果待检测图像是光谱图像，则对两幅待检测图像进行绝对值差分，得到差分图 像；
[0012]如果是合成孔径雷达SAR图像，则使用取对数差分，得到差分图像；
[0013] 步骤4:将步骤1获得的原始特征与步骤3中获得的粗变化监测结果图以坐标进行 配对，对于在粗变化检测结果图中与自身8邻域结果均不相同的输出点，去掉与它们配对对 应的样本点；
[0014] 步骤5:在堆叠自编码器SAE后构造逻辑回归分类器形成基于SAE的检测网络，然后 将筛选后的样本点数据，作为基于SAE的检测网络的输入样本，对基于SAE的检测网络进行 训练；
[0015] 训练方法为:将步骤4精炼后的变化检测结果作为输入样本的教师信号，采用误差 的反向传播BP算法微调SAE的参数。迭代次数通常为200到300次；
[0016] 步骤6:训练终止后，利用训练好的基于SAE的检测网络对待测样本进行检测，得到 最终变化检测结果图。
[0017] 所述步骤2的堆叠自编码器SAE的隐层层数为2到9,每层中节点个数取5到35。
[0018] 所述步骤2的训练方法为堆叠自编码器SAE进行无监督训练的常规训练方法。
[0019] 有益效果
[0020] 本发明提出的一种基于SAE的遥感图像变化检测方法，使用堆叠自编码器 (Stacked AutoEncoder，SAE)的训练方式，先训练好一个SAE，然后利用它对原始数据进行 自主式特征提取;再用无监督的变化检测方法对两幅原始图像做变化检测，得到的粗变化 检测结果;在去除掉一些疑似噪音的样本点后，将其作为教师数据对分类器进行训练;在训 练分类器的同时也需要对SAE的参数进行微调。训练结束后，就可得到一种基于SAE的检测 网络，从而利用该检测网络实现遥感图像变化检测。
[0021] 而且，本发明对输入数据的原始特征进行了改进。在现有技术中，通常采用的为差 分图中每个坐标点为中心的3*3区域拉抻成9维向量作为原始特征。本发明中不采用差分 图，更改为:采用两幅图像相同坐标点的3*3区域，拉伸成18维向量，再加入归一化的坐标点 位置信息，拓展成20维向量。
[0022]本发明将深度神经网络这种有监督的学习模型，与遥感图像变化检测这种无监督 的图像分类问题巧妙地结合起来，并改进了特征和训练样本，提高了检测精度和检测的鲁 棒性。
【附图说明】
[0023] 图1:基于SAE的检测网络的训练流程图；
[0024] 图2:利用训练好的基于SAE的检测网络获得变化检测结果图流程图
【具体实施方式】
[0025]现结合实施例、附图对本发明作进一步描述：
[0026]深度神经网络通常包括三个模型，分别为:堆叠自编码器、深度置信网络、卷积神 经网络。本发明使用堆叠自编码器(Stacked AutoEncoder，SAE)的训练方式。
[0027]实施方式的具体步骤如下：
[0028] (1)设时刻ti，t2为两个不同时刻，设Xi，X2为在ti，t 2时刻，同一遥感传感器及同一 波段上对同一地点获取的两幅遥感图像，且Xi，X2已经经过配准和辐射校正及几何校正。对 每个坐标点（i，j)，提取两幅图像中以此坐标点为中心的3*3的像素，并将其拉成一列，作为 一个18维的向量，并对其进行归一化，再加入归一化后的坐标信息获得一个20维 训练样本Tr(i，j)。
[0029] (2)SAE是由多个自动编码器(AutoEncoder，AE)组成。AE由三层节点组成，分别为 输入层、隐藏层和重建层。AE具体训练方式如下：
[0030] 设置输入层到隐藏层的权重为Wy，偏置为by;隐藏层到输出层权重和偏置分别为 Wz，bz，我们对Wy和Wz定义如下的约束：
[0031] ffy=ffzT=ff；
[0032] 定义代价函数：
[0033]
[0034] 其中d代表输出单元的个数，m表示训练数据的个数，Xlk表示第i个输入的第k个单 元，zlk表示第i个输出的第k个单元，分别计算代价函数c对权重W及偏置b y，bz的偏导数，然 后利用BP算法：
[0035]
[0036]来更新网络参数，最小化AE的代价函数。