基于显著性引导的sar图像变化检测方法
【技术领域】
[0001] 本发明属于雷达图像处理技术领域，涉及SAR图像的变化检测，可用于农业调查， 自然灾害和国土资源监测，军事目标监测与战场打击效果评估。
【背景技术】
[0002] 变化检测方法主要是分析多时图像中的变化信息从而对其进行理解与解译，该方 法近些年已引起了广泛的研宄和关注。特别是在遥感领域，由于合成孔径雷达SAR本身所 具有的全天时全天候的工作特性，使得SAR图像的变化检测方法已广泛应用于农业调查， 自然灾害和国土资源监测，军事目标监测与战场打击效果评估等诸多方面。
[0003] 针对同一地区不同时段获取的SAR图像，传统的SAR图像变化检测方法主要分为 两大类：有监督方法和无监督方法。有监督方法由于其需要有标记的样本来进行模型的训 练，而这类有标记样本恰是不易获取的，从而使得这类方法的使用受到局限。无监督方法不 需要标记样本，直接在构造差分图的基础上进行变化信息的分析和解译，从而使得该类方 法得到广泛应用。
[0004] 比较经典的无监督SAR变化检测方法主要有主成分分析和K均值聚类结合的 方法PCA-K、加入局部信息的模糊聚类算法FLICM、基于马尔科夫随机场的模糊聚类算法 MRFFCM。PCA-K方法在由对数比值算子构造的差分图的基础上，首先进行非重叠分块，利用 主分量分析方法获得特征空间，然后对差分图进行重叠分块，将所有向量化的块投影到该 特征空间以获得特征向量，最后对这些特征向量通过K均值算法进行聚类来获得最终的变 化检测图。FLICM和MRFFCM方法主要是在进行变化检测的过程中考虑到了数据的局部邻域 信息。尽管这些方法都能够有效的对变化信息进行检测，但仍存在检测精度不高和易受噪 声干扰的问题。
【发明内容】
[0005] 本发明的目的在于克服现有技术检测精度不高且易受噪声影响的缺点，提出一种 基于显著性引导的SAR图像变化检测方法，以提高SAR图像变化检测的精度。
[0006] 实现本发明目的的技术方案是：
[0007] -.技术原理
[0008] 在计算机视觉领域，显著图描述的是在视觉上存在强局部对比度的区域，这些区 域能被人能首先并且长时间关注到。这些强对比度来自局部区域的纹理、灰度值、形状、颜 色等等。这使得我们考虑在给定一个初始差分图的情况下SAR图像的变化检测问题，显著 性和变化检测在理论和视觉上存在着共性。显著性主要是用来提取明显区别与局部和全局 区域的那些区域，而在给定初始差分图的情况下，SAR图像的变化检测问题可以看作是寻找 与其他区域具有精确区分的区域，从该角度来看，显著图的提取和SAR图像的变化检测问 题在本质上是一致的，而且在视觉效果上来说，多时SAR图像的变化区域正好也对应着初 始差分图的显著区域。
[0009] 二.实现方案
[0010] 本发明基于显著性引导的SAR图像变化检测方法，包括如下步骤：
[0011] (1)对待检测的两幅SAR图像XjPX2进行对数比值操作，得到一个初始的差分图 Dl：
【主权项】
1. 一种基于显著性引导的SAR图像变化检测方法，包括如下步骤： (1) 对待检测的两幅SAR图像Xi和X2进行对数比值操作，得到一个初始的差分图D^
(2) 对初始差分图町提取显著图S; (3) 对显著图S进行阔值化操作，得到具有精确边缘的二值化显著图S' (4) 根据二值化显著图S'E分别提取两幅SAR图像X1和X2中对应的区域，对该区域使 用窗口大小为3X3的均值滤波器进行降斑，得到降斑后的图像Ii和12,对降斑后的图像进 行对数比值操作，得到差分图町S:
(5) 对差分图町S通过主分量分析提取特征； (6) 用K均值算法对得到的特征矩阵进行聚类，得到最终的变化检测结果。
2. 根据权利要求1所述的基于显著性引导的SAR图像变化检测方法，其中所述步骤 (2)提取初始差分图町的显著图S，按如下步骤进行： (2a)按大小为7X7的窗口W50%的重叠率对初始差分图町进行块划分，并对每一个 块进行向量化操作，得到相对应的向量并构造向量矩阵X，其中X的每一列为一个向量，X的 大小为49XLL为块的个数； (2b)计算X中第i列向量Xi和第j列向量Xj之间的不相似性；
其中，dy(Xi，Xj)为对向量Xi和向量Xj的值进行归一化后的欧式距离，dp(Xi，Xj)为向量Xi和向量XJ之间的欧式距离，C为参数取为3 ; (2c)计算向量Xi在尺度r下的显著值；
其中，K为与向量Xi最为相似的向量的个数，其值取64 ;r为尺度比例，^〇<，為)为在尺 度r下第i列向量為和与向量Xi最为相似的向量集中第h列向量為之间的不相似性，exp 为W自然对数e为底数的指数函数； (2d)计算向量Xi在各尺度下的显著值：
其中，M为尺度的个数，取值为4 ;R为尺度比例集，其包括100 %，80 %，50 %和30 %四 种不同比例；考（i)为第i个列向量和与第i列向量最为接近的向量在尺度r下归一化的 欧氏距离，□是一个归一化算子，用于将不同尺度下得到的显著图插值到与原始SAR图像 具有同样的大小。
3. 根据权利要求1所述的基于显著性引导的SAR图像变化检测方法，其中所述步骤 (3)对显著图S进行阔值化操作，按如下步骤进行： (3a)构造一个与显著图S大小相同的全零矩阵S,，再根据显著值的直方图分布设置阔 值T; (3b)将显著图S中像素点的值与阔值T进行比较；如果显著图S中像素点的值大于 阔值T，则将S冲与该像素点位置对应的像素点的值改为1，否则，仍为0,最终得到由1和 0组成的二值化显著图s'C。
4. 根据权利要求1所述的基于显著性引导的SAR图像变化检测方法，其中步骤（4)所 述的根据二值化显著图S'C分别提取两幅SAR图像X1和X2中对应的区域，是将SAR图像X1 和X2分别与二值化显著图S'E进行对应点的乘积，得到与二值化显著图S'E相对应的区域 Isi和IS2; /、, = 乂,0兩，''=1，2， 其中@为点乘算子。
5. 根据权利要求1所述的基于显著性引导的SAR图像变化检测方法，其中所述步骤 (5)对差分图Du通过主分量分析提取特征，按如下步骤进行： 巧a)将差分图町S划分为大小为3X3的不重叠块，并将每一个块进行向量化操作，即将 每一个块改写成一个向量，构造向量矩阵C，其中C的每一行为一个向量，C的大小为PX9， P为块的个数； 巧b)通过主分量分析对矩阵C生成一个特征空间； 巧C)将差分图町S划分为大小为3X3的重叠块，对每一个块进行向量化操作并投影到 特征空间W生成特征向量，特征向量的维度取为1 ; 巧d)把所有的特征向量组成一个大小为1X化XJ)的特征矩阵F，其中K为差分图町S的宽度值，J为差分图町S的长度值。
【专利摘要】本发明公开了一种基于显著性引导的SAR图像变化检测方法，主要解决现有方法对SAR图像变化区域不能精确检测以及易受斑点噪声影响的问题。其实现步骤是：1.对待检测的两幅SAR图像进行对数比值操作，得到初始差分图；2.对初始差分图提取显著图并进行阈值化操作，得到二值化显著图；3.根据二值化显著图分别提取两幅SAR图像中对应的区域，对该区域使用窗口大小为3×3的均值滤波器降斑，再通过对数比值操作得到差分图；4.对差分图通过主分量分析提取特征；5.用K均值算法对得到的特征矩阵进行聚类，得到最终的变化检测结果。本发明能对变化区域进行准确定位并精确检测，可用于公共安全，雷达图像处理和视频监控领域。
【IPC分类】G06T7-00
【公开号】CN104794729
【申请号】CN201510223998
【发明人】焦李成, 郑耀国, 马文萍, 屈嵘, 侯彪, 王爽, 杨淑媛, 马晶晶
【申请人】西安电子科技大学
【公开日】2015年7月22日
【申请日】2015年5月5日