基于非局部cv模型的sar图像变化检测方法
【专利摘要】本发明公开了一种基于非局部CV模型的SAR图像变化检测方法，实现步骤包括：(1)输入同一区域不同时间的两幅SAR图像构造差异图；(2)在差异图上做矩形轮廓曲线，初始化水平集函数；(3)计算内部区域的灰度均值和外部区域的灰度均值；(4)构建基于非邻域和基于区域信息的总能量函数E；(5)通过梯度下降法最小化E更新水平集函数φ，获得新的水平集函数φn+1，用其对差异图X进行分割；(6)判断当前的φn+1能否使总能量函数E达到最小值，如果能，停止迭代，输出最终的变化检测结果；否则，返回步骤3。本发明的方法对于差异位置更加敏感，对于初始化和图像质量的鲁棒性有明显增强。
【专利说明】
基于非局部CV模型的SAR图像变化检测方法
技术领域
[0001] 本发明属于图像领域，主要解决遥感图像的变化检测问题，具体是一种基于非邻 域统计学主动轮廓的变化检测模型，实现对SAR图像的变化检测。
【背景技术】
[0002] 随着一系列光学和雷达卫星的运行服务的发展，不断积累了大量的遥感数据，相 继建立了各种基于遥感图像数据的空间数据库。通过对这些图像数据进行进一步的分析、 处理及解译，就可以将经遥感技术获得的图像数据转化为有用的信息，进而服务于国民经 济的各个领域。
[0003] 遥感图像的变化检测技术是遥感图像研究的重要组成部分，它是对同一地域不同 时刻拍摄的多幅遥感影像，采用图像处理和模式识别等手段进行对比分析，根据图像之间 的差异来得到人们所需要的地物或目标随时间发生的变化信息。变化检测技术可以根据不 同时期获取图像的灰度值或局部纹理之间的变化，获得感兴趣的目标在数量、位置、形状及 其它属性的变化情况。这些变化或许是由图像场景的真实变化造成的，也或许是由入射角、 传感器精度、地面湿度、大气条件等变化引起的。遥感图像的变化检测可以应用于军事研究 领域，遥感气象服务，
[0004] 随着遥感技术的快速发展，遥感图像变化检测技术得到了快速的发展。Singh.A将 变化检测方法分为直接比较法(差值法、比值法等)和分类比较法，直接比较法的算法相对 简单，容易实现，但难以克服大气条件和传感器噪声等对检测结果带来的误差，受图像的预 处理精度的影响严重;而分类比较法的变化检测结果的精度容易受分类精度的影响，初始 分类的错误将造成难以弥补的误检或漏检现象。根据操作的目标对象可以把变化检测方法 分为基于特征级的变化检测、基于像素级的变化检测和基于目标级的变化检测。按照是否 需要训练样本这一准则，变化检测方法可以分为有监督的变化检测和无监督的变化检测。
[0005] 主动轮廓模型(active contour model)在图像处理中得到了广泛的研究和应用， 它是基于曲面演化和几何流形理论的模型。由Osher和Sethian提出的水平集方法[33]常常 被用来求解主动轮廓模型。一般来讲，主动轮廓模型可以分为两类，基于边界的和基于区域 的主动轮廓模型。每一类各有其缺点和优点。
[0006] 基于边界的主动轮廓模型把图像梯度作为曲线收敛到理想边界的附加约束，通常 会使用一个边界检测函数(或称为边界停止函数)来实现这一边界约束。最早的，也是最为 经典的就是Kass等提出的Snakes模型[34]。之后，Snakes模型得到了发展，更为一般的、只 依赖于曲线几何特征的GAC模型[35]被提出，其成功之处在于对曲线参数的独立性。为增强 抓取边界的能力，常数膨胀项(控制曲线收缩或膨胀)往往被加入到GAC模型的演化方程中。 然而，常数膨胀项虽然减少了迭代时间，但是合适的选值却成为它的瓶颈。过小或过大的膨 胀项都会导致结果的不精确。早期主动轮廓模型[34,35,36,37]的能量函数都是驱使曲线 向图像局部特征演化，最常用的就是边界特征，这也是此类模型被称为基于边界的模型的 原因。由于使用了局部特征，基于边界的模型往往陷入局部最优解，并且分割结果很大程度 上受到噪声、纹理和混杂特征的影响，因此，它们对复杂图像的分割结果差强人意。
[0007] 基于区域的模型[38]往往使用图像的区域信息作为约束条件，由于此类模型的优 化相对于局部的边界特征更接近全局解，因此和基于边界的模型相比具有一定的优势：1、 不使用图像梯度等局部信息，对边界模糊甚至没有明显边界的目标分割结果也很理想；2、 优化的全局解使得初始曲线的位置没有非常严格的要求(在基于边界的模型中，往往要求 初始曲线要贴近目标边界）。早期的基于区域的模型往往是基于图像模型[39,40,41，42， 43]或含参数的图像统计量的[44](如，前景区域和背景区域的均值和方差等）。
[0008] 目前的SAR变化检测算法大部分为中低分辨率、象元级、单极化的SAR变化检测。随 着SAR技术的日益成熟，图像质量逐渐提高，分辨率逐渐增强，获取数据的能力和精度越来 越高，图像获取越来越便捷，SAR变化检测的前景会更加广泛。
【发明内容】
[0009] 本发明的目的是提供一种新的基于非局部CV模型的SAR图像变化检测方法，该方 法是将原有的CV模型和非邻域模型相结合，充分利用了感兴趣的目标像素点附近的非局部 信息。该方法避免了水平集模型在变化检测中出现的局部优化问题，同时更好的利用了边 界和区域信息。
[0010] 本发明的技术关键是利用Chan-Vese模型和非邻域模型，定义一条闭合轮廓曲线 来表示目标的边界轮廓，并定义一个包含轮廓曲线信息的能量泛函，将图像变化检测过程 转变成为求解所定义的能量泛函最小值的过程。
[0011] 本发明的技术方案是如下几个步骤：
[0012] 1、输入同一区域不同时间的两幅SAR图像，构造差异图。
[0013] 2、在差异图上做个矩形轮廓曲线，初始化水平集函数。
[0014] 3、计算内部区域Q1"的灰度均值C1和外部区域Ω 2η的灰度均值C2。
[0015] 4、构建基于非邻域和基于区域信息（即CV模型）的总能量函数E。
[0016] 5、通过梯度下降法最小化总能量函数E更新水平集函数Φ，获得新的水平集函数 Φ n+1，用该水平集函数Φ n+1对差异图X进行分割。
[0017] 6、判断当前的水平集函数巾"+1能否使总能量函数E达到最小值，如果能，停止迭 代，否则，返回步骤3。
[0018] 本发明与现有技术相比具有如下优点：
[0019] 1、本发明的方法对于差异位置更加敏感，避免了原有的Chan-Vese模型在变化检 测中出现的局部优化问题；
[0020] 2、更好的利用了边界和区域信息，反映出来的结果更好；
[0021] 3、对于初始化和图像质量的鲁棒性有明显增强。
【附图说明】
[0022]图1为本发明的实现流程图；
[0023] 图2为水平集方法应用于SAR图像变化检测的整体流程图；
[0024] 图3为反映黄河的两幅前后变化图；
[0025] 图4为原有的Chan-Vese模型与本发明在图3变化检测中对比结果图；
[0026] 图5为反映伯尔尼城市的两幅前后变化图；
[0027] 图6为原有的Chan-Vese模型与本发明在图5变化检测中对比结果图。
【具体实施方式】
[0028]参照图1，本发明的实现步骤如下：
[0029]步骤1，对同一地域不同时间获取的两幅遥感图像心和办，构造差异图X。
[0030]步骤2,初始化水平集函数Φ为符号距离函数形式，根据水平集函数值的正负，根 据水平集函数的正负，将差异图X的整个图像区域分为内部区域〇:和外部区域ω2,在差异 图X上做个矩形轮廓曲线，初始化水平集函数Φ = ±D。当Φ >0,表示矩形轮廓曲线的内部 区域Ω :，当φ <〇时，表示矩形轮廓曲线的内部区域Ω 2,当Φ =0时，表示差异图中的像素点 在矩形轮廓曲线上。其中，D表示差异图中的像素点到矩形轮廓曲线的欧氏距离。
[0031]步骤3,设循环次数η的初始值为n = 0,设第η次循环时的水平集函数Φη=Φ根据 当前水平集函数Φ "将差异图X分成内部区域Ω f和外部区域Ω /，计算内部区域Ω f的灰度 均值Cl和外部区域Ω 2n的灰度均值C2。
[0032] 步骤4,结合步骤2和步骤3构建基于区域信息的能量函数Ecv:
[0033]
[0034] 其中X(x，y)表示差异图X中像素点的坐标，Ω表示整个差异图区域，C1为内部区域 Ω ιη的灰度均值，C2为外部区域Ω 2n的灰度均值，表示对水平集函数求梯度，为梯度 范式，μ为调节长度约束]
的权重参数，并且Η( Φ )为Heaviside函数，δ〇( φ ) 为Η(Φ)的导数。
[0035] 击骤Fi,构律某干韭郃域信息的能量函教Rm，
[0036]
[0037] 其中，<
X表示像素点X的邻域点集合，σ是高斯函数的标准差，d (px，py)表示像素点x，y的邻域点集合px和py之间的距离。
[0038] 步骤6,基于步骤4和步骤5,添加可取代重新初始化过程的距离正则项：
，该距离正则项中的函I
。构建总能量函数E为：
[0039] E = Ecv+Enl+D
[0040] 步骤7,通过梯度下降法最小化总能量函数E更新水平集函数Φ，获得新的水平集 函数Φ n+1，用该水平集函数Φ n+1对差异图X进行分割。
[0041] 步骤8,判断当前的水平集函数巾"+1能否使总能量函数E达到最小值，如果能，停止 迭代，此时的水平集函数Φ η+1对应的零水平集就是变化区域的边界，将变化区域的像素灰 度值赋值I，非变化区域的像素灰度值赋值0,得到变化检测结果图，否则，返回步骤(3)，用 当前的水平集函数Φ n+1替代步骤(3)中的Φ "继续循环迭代。
[0042]本发明的效果可以通过以下仿真进一步说明：
[0043] 1.仿真条件
[0044] 本实例在AMD Athlon(tm)P340Dual_Core Processor 2 · 20GHz Windows 8 · 1系统 下，Mat lab 2014a运行平台上，完成本发明仿真实验。
[0045] 2.仿真实验内容
[0046] 本实验内容包括两部分：
[0047] 1、本发明的算法模型在SAR图像中的变化检测；
[0048] 2、已有的Chan-Vese模型在SAR图像中的变化检测。
[0049] 选择了两个例子，第一组图是反映黄河的前后变化图，是Radarsat_2W卫星拍摄于 2008年6月和2009年6月。
[0050]第二组图是反映城市伯尔尼附近区域的前后变化图，是欧洲的遥感卫星拍摄于 1999年4月和5月。
[0051 ] 两组图的大小均为301*301。
[0052] 3.仿真实验结果及分析
[0053]图4和图6给出了Chan-Vese模型与本发明的结果图，图（a)是表示Chan-Vese模型 在变化检测中的结果图，图（b)是表示本发明在变化检测中的结果图。比较前后两幅图，可 见本发明相比Chan-Vese模型在边界上更加平滑，细节保持的清晰完整，变化检测的效果相 对比较好。使用Chan-Vese模型与加入了非邻域信息的CV模型的效果指标如表1所示。
[0054]表1CV模型与NL-CV模型进行变化检测的效果指标
[0056]从表1中可以看出，本发明的方法的变化检测漏警率和误警率相对原有的CV模型 较低，PCC指数总体效果较好些，检测的结果优良。
[0057]根据以上结果可知，本发明的方法对于差异位置更加敏感，避免了原有的Chan-Vese模型在变化检测中出现的局部优化问题，更好的利用了边界和区域信息，反映出来的 结果更好，对于初始化和图像质量的鲁棒性有明显增强。
[0058]总体来说，本发明方法明显提高了 SAR图像变化检测的效果。
【主权项】
1. 一种基于非局部CV模型的SAR图像变化检测方法，其特征在于所述方法包括下列步 骤： (1) 对同一地域不同时间获取的两幅遥感图像Xi和拉，构造差异图X; (2) 初始化水平集函数Φ为符号距离函数形式； (3) 设循环次数η的初始值为n = 0,设第η次循环时的水平集函数φη= Φ，根据当前水平 集函数Φ η将差异图X分成内部区域Ω 1"和外部区域Ω 2";计算内部区域Ω 1"的灰度均值C1和 外部区域Ω 2"的灰度均值C2; (4) 结合步骤2和步骤3构建基于区域信息的能量函数Ecv:其中X(x，y)表示差异图X中像素点的坐标，Ω表示整个差异图区域，Cl为内部区域Ωι" 的灰度均值，C2为外部区域Ω 2"的灰度均值，表示对水平集函数求梯度，为梯度范式， μ为调节长度约束巧的权重参数，并且Η( Φ )为化aviside函数，δ〇( Φ )为Η (Φ)的导数； (5) 构建基于非邻域信息的能量函数Enl:其中，Ρχ表示像素点X的邻域点集合，0是高斯函数的标准差，d (Px， Py)表示像素点X，y的邻域点集合Ρχ和Py之间的距离； (6) 基于步骤4和步骤5,添加可取代重新初始化过程的距离正则项：，该距离正则项中的函数Ρ(|νψ;构建总能量函数E为： E = Ecv+Enl+D (7) 通过梯度下降法最小化总能量函数E更新水平集函数Φ，获得新的水平集函数 Φ n+l，用该水平集函数Φ n+l对差异图X进行分割； (8) 判断当前的水平集函数φη+1能否使总能量函数E达到最小值，如果能，停止迭代，否 贝1J，返回步骤(3)，用当前的水平集函数φη+1替代步骤(3)中的φη继续循环迭代。2. 根据权利要求书1所述的基于非局部CV模型的SAR图像变化检测方法，其特征是:所 述步骤(2)包括下列步骤： (2a)根据水平集函数值的正负，将差异图X的整个图像区域分为内部区域Ωι和外部区 域Ω2; (2b)在差异图X上做个矩形轮廓曲线，初始化水平集函数Φ=±0;当Φ>0,表示矩形 轮廓曲线的内部区域Ωι，当Φ<0时，表示矩形轮廓曲线的内部区域Ω 2,当Φ=0时，表示差 异图中的像素点在矩形轮廓曲线上;其中，D表示差异图中的像素点到矩形轮廓曲线的欧氏 距离。3.根据权利要求书1所述的基于非局部CV模型的SAR图像变化检测方法，其特征是:所 述步骤(8)包括下列步骤： (3a)判断当前的水平集函数φη+1能否使总能量函数E达到最小值； (3b)否则，返回步骤(3)，用当前的水平集函数φη+1替代步骤（3)中的φη继续循环迭 代； (3c)如果能，停止迭代，此时的水平集函数φη+1对应的零水平集就是变化区域的边界， 将变化区域的像素灰度值赋值1，非变化区域的像素灰度值赋值0,得到变化检测结果图； (3d)结束变化检测过程。
【文档编号】G06T7/00GK105844637SQ201610168000
【公开日】2016年8月10日
【申请日】2016年3月23日
【发明人】公茂果, 王倩, 刘嘉, 马文萍, 张普照, 赵秋楠, 马晶晶, 胡天妤
【申请人】西安电子科技大学