基于多粒度级联森林模型的SAR图像变化检测方法与流程

本发明属于图像处理
技术领域：
，尤其涉及一种sar图像变化检测方法，可用于环境检测、农业调查、自然灾害检测、森林资源监测等。
背景技术：
：合成孔径雷达sar具有高分辨率，可全天时、全天候工作的优势，被广泛应用于民用和军事领域。目前，自然环境的恶劣变化，城市的迅速发展，使得图像变化检测技术迅速崛起。变化检测是用于检测同一地点在一段时间内所发生的变化情况，传统的变化检测方法采用bruzzone等人在2002年提出的经典三步流程范式：1)预处理；2)生成差异图；3)分析差异图。由于sar图像深受相干斑的影响，这种传统的方法的检测结果很大程度上依赖于差异图的性能，往往当差异图的性能不好时，产生的精度会很低。为了克服这个缺点，可利用具有强大的抽象表达能力的深度学习来进行变化检测。随着近年来深度学习的火热，深度学习的应用越来越频繁。理论分析证明，深度学习能够学习表达抽象的复杂函数，深度结构由多层非线性结构组成。很多学者都利用深度学习来解决变化检测问题，例如gong等人提出基于深度神经网络的sar图像变化检测，利用堆叠的限制性波尔兹曼机，通过逐层的预训练及对整个网络的精调，取得了较高的精度；如gao等人提出的基于pcanet的sar图像变化检测，利用一种简化的cnn模型pcanet从多时相图像中提取鲁棒的更易于分类的特征来实现变化检测。上述大多数利用深度学习进行变化检测的方法是基于深度神经网络的，由于深度神经网络需要大量的训练样本，但是实际中可用于训练的具有标签的样本非常有限，这会严重影响此方法的性能；其次深度神经网络的超参数很多，调参困难，导致效率较低。技术实现要素：本发明的目的在于针对传统方法的不足，提出一种基于多粒度级联森林模型的sar图像变化检测方法，以避免检测结果对差异图的强依赖性，减少训练所需要的样本数及可调节超参数的个数，提高检测精度和效率。为实现上述目的，本发明基于多粒度级联森林模型的sar图像变化检测方法，包括如下：1.一种基于多粒度级联森林模型的sar图像变化检测方法，其特征在于，包括如下步骤：(1)给定两幅配准的具有相同尺寸的多时相sar图像i1和i2，并对i1和i2使用对数比算子，生成归一化的差异图id；(2)以归一化的差异图id的每个像素点为中心按大小为21×21取块，每个图像块作为一个样本，也就是每个样本对应一个像素点，然后对每个样本提取类哈尔特征；(3)随机选取少量样本作为训练样本，利用训练样本的类哈尔特征对多粒度级联森林模型进行训练，得到训练好的多粒度级联森林模型；(4)利用训练好的多粒度级联森林模型，对所有样本进行测试，得到预测结果(5)生成两个概率图i1和i0将所有样本的类哈尔特征输入到已训练好的多粒度级联森林模型，得到每个样本属于变化类的概率p1和属于非变化类的概率p0，其中，对每个样本有p1+p0＝1，用所有样本中属于变化类的概率p1构成第一概率图i1，用所有样本中属于非变化类的概率p0构成第二概率图i0；(6)利用差异图id、第一概率图i1和第二概率图i0构成新特征，重新训练多粒度级联森林模型，得到新预测结果并将新预测结果的kappa系数与上一次预测结果的kappa系数进行比较，选择kappa系数较高的预测结果作为最终的变化检测结果。本发明与现有技术相比具有以下优点：第一，由于本发明提取了输入sar图像的差异图的类哈尔特征，并利用该特征对多粒度级联森林模型进行训练，两者的结合克服了传统方法对差异图的强依赖性，提高了sar图像变化检测的精度。第二，由于本发明采用了多粒度级联森林进行变化检测，在训练样本数量较少的情况下也能保持很好的检测效果；且多粒度级联森林的超参数较少，便于调节。第三，由于本发明采用的多粒度级联森林是一种基于决策树的集成方法，便于对模型进行并行处理，可同时保证较高的检测准确率和时间效率。附图说明图1为本发明的实现流程图；图2为用本发明对黄河农田c地区sar图像的变化检测的仿真实验结果图；图3为用本发明对黄河农田d地区sar图像的变化检测的仿真实验结果图；图4为用本发明对bern地区sar图像的变化检测的仿真实验结果图。具体实施方法下面结合图1对本发明的实施例和效果进行详细的描述。参照图1，本发明的实现步骤如下：步骤1：生成归一化的差异图。(1a)给定两幅配准的具有相同的尺寸的多时相sar图像i1和i2，对i1和i2使用对数比算子，生成差异图：其中，log代表取自然对数操作，|·|表示取绝对值操作；(1b)对差异图id1进行归一化，得到归一化后的差异图id：其中，min(·)表示取最小值操作，max(·)表示取最大值操作。步骤2：提取类哈尔特征。(2a)对差异图id以每个像素点为中心取大小为w×h图像块；(2b)计算图像块的积分图，积分图s定义为此图像块从左上角(0,0)到右下角(x,y)范围内所有像素值的和：其中，1≤x≤w,1≤y≤h，(x',y')表示为左上角(0,0)到右下角(x,y)范围内任意像素点的位置，i代表图像块；(2c)将每个图像块内的任意垂直矩形表示为：ri＝(a,b,c,d)，其中i∈{1,...,n}，n表示图像块内的矩形个数，(a,b)表示矩形左上角像素点的位置，1≤a≤w,1≤b≤h，c表示矩形ri的宽度，d表示矩形ri的长度，c,d>0；(2d)将矩形ri中所有像素值的和表示为rs(ri)，通过步骤(2b)得到的积分图s，计算矩形ri的像素值之和rs(ri)：rs(ri)＝s(a-1,b-1)+s(a+c-1,b+d-1)-s(a-1,b+d-1)-s(a+c-1,b-1)；(2e)根据矩形ri的像素值之和rs(ri)，计算得到每个图像块的类哈尔特征fi:其中，ωi表示第i个矩形的权重，ri表示第i个矩形，n表示矩形的个数。步骤3：训练多粒度级联森林模型。(3a)多粒度扫描：利用三个尺度分别为15×1，30×1和60×1的滑动窗口来扫描每一个类哈尔特征：假设有一个m维的类哈尔特征，利用一个15×1的滑动窗口进行扫描，产生(m-14)个的15×1维的特征；再利用一个30×1的滑动窗口进行扫描，产生(m-29)个的30×1维的特征；再利用一个60×1的滑动窗口进行扫描，产生(m-59)个的60×1维的特征；通过多粒度扫描，一个m维的类哈尔特征总共能产生(3m–102)个新特征，此操作使得特征更加丰富；(3b)训练级联森林级联森林由多个层级构成，它的每一层均是一个集成的决策树森林，设置每一层包括2个随机森林和2个完全随机森林，每个森林都包含有20棵决策树；利用多粒度扫描产生的特征对级联森林进行训练，得到训练好的多粒度级联森林模型。步骤4：将所有样本输入到训练好的多粒度级联森林模型，对所有样本进行测试，得到预测结果步骤5：生成两个概率图i1和i0。将所有样本的类哈尔特征输入到已训练好的多粒度级联森林模型，得到每个样本属于变化类的概率p1和属于非变化类的概率p0，且对每个样本其p1+p0＝1；用所有样本中属于变化类的概率p1构成第一概率图i1；用所有样本中属于非变化类的概率p0构成第二概率图i0。步骤6：生成最终的变化检测结果。(6a)利用差异图id、概率图i1和概率图i0构成三通道图像；(6b)对三通道图像的每个通道提取类哈尔特征，将三个通道对应位置提取的类哈尔特征首尾相接构成一个新的特征；(6c)利用新特征重新训练多粒度级联森林模型，得到新的多粒度级联森林模型和新的预测结果(6d)设kappa系数为衡量变化检测结果图与对应参考图一致性的参数，其中参考图为人工标记的真实变化图像；计算第i+1次变化检测结果的kappa系数kappai+1和第i次变化检测结果的kappa系数kappai；(6d1)计算第i+1次变化检测结果的kappa系数kappai+1；首先，计算第i+1次被正确检测为变化类的像素点数：tpi+1＝nc-fni+1，其中，nc表示参考图中变化类的像素点数，fni+1表示第i+1次将变化类错检为未变化类的像素点数；接着，计算第i+1次被正确检测为非变化类的像素点数：tni+1＝nu-fpi+1，其中，nu表示参考图中非变化类的像素点数，fpi+1表示第i+1次将非变化类错测为变化类的像素点数；接着，利用tpi+1和tni+1计算第i+1次的实际一致率：其中，n表示图像像素点的总数；接着，利用所得的tpi+1，tni+1，fpi+1，fni+1，nu，nc和n计算第i+1次的理论一致率：最后，利用prai+1和prei+1计算得到第i+1的kappa系数：(6d2)计算第i次变化检测结果的kappa系数kappai：首先，计算第i次被正确检测为变化类的像素点数：tpi＝nc-fni，其中，nc表示参考图中变化类的像素点数，fni表示第i次将变化类错检为未变化类的像素点数；接着，计算第i次被正确检测为非变化类的像素点数tni＝nu-fpi，其中，nu表示参考图中非变化类的像素点数，fpi表示第i次将非变化类错测为变化类的像素点数；接着，利用tpi和tni计算第i次的实际一致率：其中，n表示图像中像素点总数；接着，利用所得的tpi，tni，fpi，fni，nu，nc和n计算第i次的理论一致率最后，利用prai和prei计算得到第i的kappa系数：(6d3)将kappai+1和kappai进行比较，选择kappa系数较高的预测结果作为最终的变化检测结果。下面结合仿真实验对本发明的效果做进一步的说明。1.仿真实验数据本发明的仿真实验使用了三组真实的sar图像数据及相应的参考图，其中：第一组真实的sar图像数据及相应的参考图是黄河农田c地区的sar图像，如图2所示，图像大小为257×289，其中图2(a)是2008年6月黄河农田c地区的sar图像，图2(b)是2009年6月黄河农田c地区的sar图像，图2(c)是黄河农田c地区相应的参考图。第二组真实的sar图像数据及相应的参考图是黄河农田d地区的sar图像，如图3所示，图像大小为306×291，其中图3(a)是2008年6月黄河农田d地区的sar图像，图3(b)是2009年6月黄河农田d地区的sar图像，图3(c)是黄河农田d地区相应的参考图。第三组真实的sar图像数据及相应的参考图是bern地区的sar图像，如图4所示，图像大小为301×301，其中图4(a)是1999年4月bern地区的sar图像，图4(b)是1999年5月bern地区的sar图像，图4(c)是bern地区相应的参考图。2.仿真实验的评价指标漏检fn：将变化类错检为未变化类的像素点数；错检fp：将未变化类错检为变化类的像素点数；oe：错检和漏检的总像素点数；kappa系数：衡量变化检测图与参考图一致性的参数。3.仿真内容与结果使用本发明在上述三组真实sar图像数据集上进行测试，得到的变化检测图如图2(d)，3(d)，4(d)所示，得到的变化检测评价指标如表1所示。表1数据集fpfnoekappa黄河农田c70355312560.9432黄河农田d1951383330.9666bern1061212270.8965从表1中可以看出，三组数据均得到较高的kappa系数。将本发明在黄河农田c地区产生的变化检测图2(d)与对应的参考图2(c)进行对比，将本发明在黄河农田d产生的变化检测图3(d)和对应的参考图3(c)进行对比，将本发明在bern地区产生的变化检测图4(d)和对应的参考图4(c)进行对比，可见本发明采用的方法能有效抑制差异图中的噪声对最终结果的影响，有效地降低了错检fp和漏检fn，提高了变化检测的kappa系数。当前第1页12