本发明涉及多时相遥感影像自动检测领域，尤其是涉及一种面向大尺度复杂场景下遥感影像的自动变化检测方法。
背景技术：
变化检测是利用不同时期的遥感影像，定量地分析和确定地物变化的特征和过程。近年来，随着陆表的加速变化，准确自动地识别多时相遥感影像中的土地覆盖变化已成为一个日益迫切的任务。在过去的几十年里，各国学者不断地提出许多新颖的变化检测技术，一些较为先进的变化检测技术可以实现遥感影像精细尺度上的稳健检测，并在不同的遥感应用中得到有效利用如农业、林业、城市以及灾害监测等。
但在实际的变化检测应用中，遥感变化检测仍然面临着许多问题和挑战。一项突出的问题是如何在大尺度复杂场景影像上有效使用变化检测技术，实现自动化、高准确性的变化识别。一方面，当前大多数的变化检测技术主要被开发并集中应用于小范围区域的图像对中。影像的覆盖范围相对较小，变化类型相对简单，并且易于通过遵循一定的统计分布(例如，高斯分布)进行建模。一旦所考虑的变化检测问题扩展到大范围，则可能存在更复杂的条件，如更多不确定类型的变化。一些非监督的变化检测算法(例如，期望最大化)可能会由于无法转化为局部最优解而失败。在这种情况下，监督的方法可以提供更好的变化检测解决方案，但是大量训练样本的选择始终是一项耗时的任务。尤其是在大尺度图像场景中，现实中几乎不可能拥有较为全面的地面真实数据或先验知识。另一方面，传统变化检测方法的设计高度依赖于精确的多时相遥感影像配准过程。此过程通常是手动执行，配准残差可能导致同一土地覆盖物错位进而引起非相关变化的误检。
技术实现要素：
本发明的目的就是为了克服上述现有技术存在的缺陷而提供一种识别精度高的面向大尺度复杂场景下遥感影像的自动变化检测方法。
本发明的目的可以通过以下技术方案来实现：
一种面向大尺度复杂场景下遥感影像的自动变化检测方法，包括以下步骤：
s1：获取前后时相的遥感影像数据对；
s2：从遥感影像数据对中提取特征点，并进行影像配准；
s3：基于配准后的遥感影像数据对，通过差值法获取差值影像；
s4：提取差值影像的显著度，生成变化类伪训练样本和不变化类伪训练样本；
s5：将变化类伪训练样本和不变化类伪训练样本输入分类器中，对步骤s3获取的差值影像进行二值分类，获取关于变化和不变化两种类别的二值检测结果。
进一步地，所述步骤s2具体包括以下步骤：
s201：利用sift算法对遥感影像数据对进行粗匹配，提取特征点；
s202：根据特征点，生成sift特征向量，并进行sift特征向量间相似性判定度量的计算；
s203：根据sift特征向量间相似性判定度量的计算结果，进行特征点匹配，获取特征点对；
s204：去除特征点对中的误匹配点；
s205：根据去除误匹配点后的特征点对，估算仿射模型参数，对步骤s1中获取的遥感影像数据对进行重采样和仿射变换，获取配准影像。
进一步地，步骤s201中，利用sift算法进行粗匹配包括，采用高斯差分dog算子检测特征点，通过拟合三维二次函数确定特征点的位置和尺度，采用图像梯度的方法为每个特征点分配基准方向。
具体地，粗匹配包括尺度空间极值检测、关键点定位、方向分配以及关键点描述四个部分。首先使用dog算子在尺度空间检测稳定的关键点。由于检测到的极值点是离散空间的极值点，通过拟合三维二次函数精确确定关键点的位置和尺度，同时去除低对比度的关键点和不稳定的边缘响应点。为了使描述符具有旋转不变性，需要利用图像的局部特征为给每一个关键点分配一个基准方向，因此使用图像梯度的方法求取局部结构的稳定方向。最后将4×4个样本子区域上的所有方向直方图级联生成128维的sift特征向量。
进一步地，步骤s202中，生成sift特征向量具体为，将4×4个样本子区域上的所有方向直方图级联生成128维的sift特征向量。
进一步地，步骤s202中，当两张影像的sift特征向量生成后，采用关键点特征向量的欧式距离作为两幅影像中关键点的相似性判定度量，所述相似性判定度量的计算表达式为：
式中，df为相似性判定度量，a为不同维度，f1,i为第a个维度中第1个sift特征向量，f2,i为第a个维度中第2个sift特征向量。
进一步地，步骤s204中，为了消除初始候选中的错误匹配，采用了两种离群值去除策略。一种有效的度量是最近邻距离比，它表示最近邻距离与次近邻距离之比。另一种是利用随机样本一致性ransac方法，通过鲁棒性的几何约束去除误匹配点。
进一步地，步骤s3中，提取差值影像的显著度具体为，首先利用dog滤波和均值滤波函数分别对差值影像进行滤波，并将dog滤波后的图像从rgb颜色空间转换到lab颜色空间，然后计算差值影像中每个像素的显著度，每个像素的显著性定义为均值滤波得到的平均矢量fμ与转换到lab颜色空间的像素矢量fωhc之间的欧式距离，所述每个像素的显著度的计算表达式为：
sm(i,j)＝||fμ-fωhc(i,j)||
式中，sm(i,j)为第i行第j列像素的显著度，fμ为均值滤波后的差值影像的平均矢量，fωhc(i,j)为差值影像由dog滤波转换到lab颜色空间的图像中第i行第j列像素的像素矢量，||·||为范数计算符。
进一步地，步骤s3中，所述不变化类伪训练样本包括第一部分像素和第二部分像素，所述第一部分像素为特征点及其5×5领域中显著度小于第一阈值的像素，所述第二部分像素为从显著度小于所述第一阈值的像素中随机选取5％的像素，所述第一阈值为计算差值影像显著度获取的视觉显著图中像素值的均值。
进一步地，步骤s3中，所述变化类伪训练样本包括从显著度大于第二阈值的像素中，随机选取5％的像素，所述第二阈值为计算差值影像显著度获取的视觉显著图中像素值的二倍均值。
进一步地，步骤s4中，所述分类器为线性支持向量机。
与现有技术相比，本发明具有以下优点：
(1)本发明面向大尺度复杂场景下遥感影像的自动变化检测方法结合了影像的自动配准和检测过程，自动生成关于变化和不变化两类的伪训练样本，然后载入分类器对差值影像分类，与常规的非监督变化检测方法相比，本发明在实现较高自动化程度的基础上，可获得更高的检测识别精度且兼顾效率，且可应用于大尺度复杂场景下遥感影像的变化检测。
(2)本发明采用lsvm分类器对差值影像分类，保证了分类的精度和运算效率。
(3)本发明为消除特征点对中的误匹配点，提供了两种离群值去除策略，一种基于最近邻距离比进行度量；另一种是利用随机样本一致性ransac方法，通过鲁棒性的几何约束去除误匹配点。
(4)经过试验对比，证明本发明方法具有更高的变化检测性能，而且适用于检测覆盖范围较大、地类较为复杂的遥感影像场景。
附图说明
图1为本发明方法的流程示意图；
图2为前后时相的遥感影像数据对的变化强度图；
图3(a)到图3(e)为不同方法得到的二值变化检测图，其中，图3(a)为ki二值检测图，图3(b)为otsu二值检测图，图3(c)为fcm二值检测图，图3(d)为本发明基于变化强度图分类得到的二值检测图，图3(e)为本发明基于差值影像分类得到的二值检测图；
图4(a)到图4(f)为火灾局部区域a的对比图，其中，图4(a)为局部变化强度图，图4(b)为ki的局部二值检测图，图4(c)为otsu的局部二值检测图，图4(d)为fcm的局部二值检测图，图4(e)为本发明基于变化强度图分类得到的局部二值检测图，图4(f)为本发明基于差值影像分类得到的局部二值检测图；
图5(a)到图5(f)为水体局部区域b的对比图，其中，图5(a)为局部变化强度图，图5(b)为ki的局部二值检测图，图5(c)为otsu的局部二值检测图，图5(d)为fcm的局部二值检测图，图5(e)为本发明基于变化强度图分类得到的局部二值检测图，图5(f)为本发明基于差值影像分类得到的局部二值检测图。
具体实施方式
下面结合附图和具体实施例对本发明进行详细说明。本实施例以本发明技术方案为前提进行实施，给出了详细的实施方式和具体的操作过程，但本发明的保护范围不限于下述的实施例。
实施例1
如图1所示，本实施例提供一种面向大尺度复杂场景下遥感影像的自动变化检测方法，本方法主要由以下三个步骤组成：
1)获取大区域场景下的前后时相遥感影像数据对，并进行自动配准；
2)根据匹配提取的特征点以及配准后差值影像提取的显著图，自动生成变化和不变化两种伪训练样本；
3)将生成的伪训练样本输入到lsvm对差值影像执行二值分类，生成关于变化和不变化两种类别的检测结果，从而实现大尺度复杂场景下遥感影像的自动变化检测。
下面对本实施例自动变化检测方法关键部分和试验对比进行详细描述。
1、基于sift的影像自动配准
利用sift算法实现粗匹配并生成初始候选点。sift粗匹配包括尺度空间极值检测、关键点定位、方向分配以及关键点描述四个部分。首先使用dog算子在尺度空间检测稳定的关键点。由于检测到的极值点是离散空间的极值点，通过拟合三维二次函数精确确定关键点的位置和尺度，同时去除低对比度的关键点和不稳定的边缘响应点。为了使描述符具有旋转不变性，需要利用图像的局部特征为给每一个关键点分配一个基准方向，因此使用图像梯度的方法求取局部结构的稳定方向。最后将4×4个样本子区域上的所有方向直方图级联生成128维的sift特征向量。
当两张影像的sift特征向量生成后，采用关键点特征向量的欧式距离作为两幅影像中关键点的相似性判定度量：
其中，df是两幅影像上相应特征向量f1和f2之间的特征距离，i代表不同维度。
初始候选点中的误匹配点去除。为了消除初始候选中的错误匹配，采用了两种离群值去除策略。一种有效的度量是最近邻距离比，它表示最近邻距离与次近邻距离之比。另一种是利用随机样本一致性ransac方法，通过鲁棒性的几何约束去除误匹配点。
选择仿射模型进行几何变换和影像变形。根据去除误匹配后的特征点对估算仿射模型参数，然后对输入影像进行重采样和仿射变换，最终生成用于变化检测的配准影像。
2、生成伪训练样本
针对配准后的影像，首先通过差值法得到前后时相影像的差值影像，得到差值影像xd：
xd＝x2-x1(2)
其中，x1为配准的前一时相影像，x2为配准的后一时相影像，xd代表前后时相影像的差值影像，b代表差值影像的波段数，ρ代表差值影像的变化强度图。
提取差值影像的视觉显著图。首先利用dog滤波和均值滤波函数分别对xd进行滤波，并将dog滤波后的图像从rgb颜色空间转换到lab颜色空间，此时获得视觉显著图sm。每个像素的显著性定义为均值滤波得到的平均矢量fμ与转换到lab颜色空间的像素矢量fωhc之间的欧式距离：
sm(i,j)＝||fμ-fωhc(i,j)||(4)
其中，fμ代表xd均值滤波得到的均值图像，fωhc代表经dog滤波转换到lab颜色空间的图像，sm(i,j)代表对应像素(i,j)处的显著性值。
最后利用匹配提取的特征点以及差值影像提取的视觉显著图，自动生成变化和不变化两种伪训练样本。不变类伪训练样本是由匹配提取的特征点和差值影像提取的显著图两部分生成。第一部分来自特征匹配提取的特征点及其5×5邻域显著性小于给定阈值tnc(tnc：定义为sm图像素值的均值)的像素。第二部分是从显著性小于给定阈值tnc的像素中随机生成5％的像素。两部分像素共同作为不变类伪训练样本集合。变化类伪训练样本是基于差值影像提取的显著图生成。从显著性大于给定阈值tc(tc：定义为sm图像素值的2倍均值)的像素中随机选择5％的像素作为变化类伪训练样本。
3、执行lsvm分类，实现大尺度复杂场景下遥感影像的二值变化检测
选用主流的svm作为本实施例的分类器。对比svm中线性核(l)和高斯核(rbf)两种核函数的性能，发现二者对大尺度遥感影像分类的精度相差不大，但lsvm(线性支持向量机)的运算时间远低于rbfsvm(高斯核支持向量机)。因此将变化和不变化两种伪训练样本输入到效率较高的lsvm中对xd执行二值分类，生成关于变化和不变化两种类别的二值检测结果，从而实现大尺度复杂场景下遥感影像的自动变化检测。
4、试验对比
为综合对比不同变化检测方法的性能并验证本实施例所提自动变化检测方法的有效性和优势，选取kittler-illingworth阈值法(ki)、大津阈值法(otsu)以及模糊c均值聚类法(fcm)对变化强度图进行自动阈值分割，得到二值变化检测结果。为了保证对比实验的一致性，将本实施例分别基于变化强度图以及差值影像执行二值分类，并与其它三种方法一同进行对比分析。
4.1、具体实施方式
试验数据采用中等分辨率landsat-8卫星遥感1级(l1)产品数据，地图投影为utm-wgs84南极洲极地投影，源于美国地质调查局(usgs)官网。本实验使用了覆盖美国门多西诺2018年火灾区的整景影像数据，火灾前后时相影像的获取时间分别是2018年7月17日和2018年8月11日，影像大小为7750×7870像素，地表覆盖类型较为复杂。选择landsat-8影像数据中显示变化效果最好的短波红外2(swir2)、短波红外1(swir1)和近红外(nir)三个波段作为变化检测的输入影像。
4.2、试验结果
对比分析不同方法得到的二值变化检测结果
手动选取10000像素的真实训练样本(变化类和不变化类各5000像素)，通过计算二值变化检测图与真实训练样本的混淆矩阵，评价ki、otsu、fcm以及本实施例(基于变化强度图和差值影像分类两种方式)4种变化检测方法的精度指标。选取总体精度(overallaccuracy，简称oa)，kappa系数(kappacoefficient，简称k)，错分误差(commissionerror，简称ce)，漏分误差(omissionerror，简称oe)，总误差(totalerror，简称te)五类评价指标。评价结果见表1，然后对比分析4种方法的二值变化检测图，同时选取两处局部区域(a和b)深入对比分析。
表1精度评价表
从精度评价结果来看，本实施例提出的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，无论是基于变化强度图还是基于差值影像分类，总精度均远高于其它三种方法。由于本实施例基于差值影像分类方式输入的信息更加全面，精度要高于基于变化强度图分类的结果。因此，本实施例最终采用的是基于差值影像分类的方式。
整景影像的变化可从图2给出的变化强度图看出，图中白点表示变化类像素，黑点表示不变化类像素。强度图的左侧区域基本都是海水，由于前后时相影像仅相差一个月，大量的变化其实是来自于陆地上的火灾区域。图3给出了4种方法的全局二值检测图。与变化强度图相比，可以明显地看出ki(图(3a))、otsu(图(3b))以及fcm(图(3c))三种方法的二值检测图上存在大量的错分像素。本实施例的两种方法都能很好地抑制非变化类的背景干扰。图4和图5分别是局部区域a和b的放大对比图。其中，a主要反映火灾区域的变化，b主要反映水体区域的变化。与图(4a)给出的强度图相比，效果最好的是本实施例基于差值影像分类的结果(图(4f))，图(4e)给出的基于变化强度图分类结果虽然也能较好地抑制背景像素的错分，但火灾区域像素存在一些漏检现象。其余三种方法均将大量的背景像素错分为变化类。对比图(5a)给出的强度图，ki(图(5b))、otsu(图(5c))以及fcm(图(5d))三种方法对水体区域的检测效果较差，效果最好的依旧是本实施例方法(图(5f))。因此，本实施例提出的自动检测方法在用于大尺度复杂场景遥感影像变化检测时效果最好。
综上，经过一系列的对比分析，证明所提出的一种面向大尺度复杂场景下遥感影像的自动变化检测方法与所参照的其它方法对比来看，具有更高的变化检测性能，而且适用于检测覆盖范围较大、地类较为复杂的遥感影像场景。
以上详细描述了本发明的较佳具体实施例。应当理解，本领域的普通技术人员无需创造性劳动就可以根据本发明的构思作出诸多修改和变化。因此，凡本技术领域中技术人员依本发明的构思在现有技术的基础上通过逻辑分析、推理或者有限的实验可以得到的技术方案，皆应在由权利要求书所确定的保护范围内。
1.一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，包括以下步骤：
s1：获取前后时相的遥感影像数据对；
s2：从遥感影像数据对中提取特征点，并进行影像配准；
s3：基于配准后的遥感影像数据对，通过差值法获取差值影像；
s4：提取差值影像的显著度，生成变化类伪训练样本和不变化类伪训练样本；
s5：将变化类伪训练样本和不变化类伪训练样本输入分类器中，对步骤s3获取的差值影像进行二值分类，获取关于变化和不变化两种类别的二值检测结果。
2.根据权利要求1所述的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，所述步骤s2具体包括以下步骤：
s201：利用sift算法对遥感影像数据对进行粗匹配，提取特征点；
s202：根据特征点，生成sift特征向量，并进行sift特征向量间相似性判定度量的计算；
s203：根据sift特征向量间相似性判定度量的计算结果，进行特征点匹配，获取特征点对；
s204：去除特征点对中的误匹配点；
s205：根据去除误匹配点后的特征点对，估算仿射模型参数，对步骤s1中获取的遥感影像数据对进行重采样和仿射变换，获取配准影像。
3.根据权利要求2所述的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，步骤s201中，利用sift算法进行粗匹配具体为，采用高斯差分dog算子检测特征点，通过拟合三维二次函数确定特征点的位置和尺度，采用图像梯度的方法为每个特征点分配基准方向。
4.根据权利要求3所述的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，步骤s202中，生成sift特征向量具体为，将4×4个样本子区域上的所有方向直方图级联生成128维的sift特征向量。
5.根据权利要求4所述的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，步骤s202中，所述相似性判定度量的计算表达式为：
式中，df为相似性判定度量，a为不同维度，f1,i为第a个维度中第1个sift特征向量，f2,i为第a个维度中第2个sift特征向量。
6.根据权利要求2所述的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，步骤s204中，采用基于最近邻距离比的离群值去除策略，去除特征点对中的误匹配点。
7.根据权利要求1所述的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，步骤s3中，提取差值影像的显著度具体为，首先利用dog滤波和均值滤波函数分别对差值影像进行滤波，并将dog滤波后的图像从rgb颜色空间转换到lab颜色空间，然后计算差值影像中每个像素的显著度，所述每个像素的显著度的计算表达式为：
sm(i,j)＝||fμ-fωhc(i,j)||
式中，sm(i,j)为第i行第j列像素的显著度，fμ为均值滤波后的差值影像的平均矢量，fωhc(i,j)为差值影像由dog滤波转换到lab颜色空间的图像中第i行第j列像素的像素矢量，||·||为范数计算符。
8.根据权利要求1所述的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，步骤s3中，所述不变化类伪训练样本包括第一部分像素和第二部分像素，所述第一部分像素为特征点及其5×5领域中显著度小于第一阈值的像素，所述第二部分像素为从显著度小于所述第一阈值的像素中随机选取5％的像素，所述第一阈值为计算差值影像显著度获取的视觉显著图中像素值的均值。
9.根据权利要求1所述的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，步骤s3中，所述变化类伪训练样本包括从显著度大于第二阈值的像素中，随机选取5％的像素，所述第二阈值为计算差值影像显著度获取的视觉显著图中像素值的二倍均值。
10.根据权利要求1所述的一种面向大尺度复杂场景下遥感影像的自动变化检测方法，其特征在于，步骤s4中，所述分类器为线性支持向量机。
技术总结