一种图像变化检测方法
【专利摘要】本发明公开了一种图像变化检测方法，应用于同一地区不同时期的遥感图像变化区域检测问题，所述方法包括：通过对所述图像的对数比值差异图进行双边滤波得到初始差异图；通过对上述初始差异图进行频率调谐显著性区域检测得到显著图；通过对所述显著图进行局部熵计算得到加熵显著图；根据小波变换法融合初始差异图和加熵显著图得到融合图；通过最大类间方差法对融合图进行分割得到最终的变化检测结果图。本发明将显著性检测模型引入到变化检测问题中，对差异图和显著图进行融合，在保留图像细节的基础上最大程度地突出变化区域，利于进一步提取和分析。实验结果表明，本发明方法能够有效地降低误检率，提高了变化检测结果的精度。
【专利说明】
-种图像变化检测方法
技术领域
[0001] 本发明属于图像处理技术领域，尤其设及一种图像变化检测方法。
【背景技术】
[0002] 遥感图像变化检测技术是通过提取、分析和判别多时相遥感图像中的同一目标或 区域，来获取其中的变化信息。遥感图像变化检测技术从人类获得多时相的航拍照片开始， 就成为遥感研究的重要方向，从早期的人工目视解译、胶片重叠等方法，到计算机自动进行 分析处理，变化检测技术取得了很大进步。随着各国高分辨率遥感卫星的次第发射W及对 全球环境变化的关注，遥感图像变化检测技术的研究成为遥感技术研究的热点领域之一。
[0003] 在利用计算机进行遥感图像变化检测之前，主要通过人工的方法对遥感图像进行 观察、比较和分析、标注，尽管效率低，但能够保证检测结果的精度。随着计算机技术的发 展，遥感图像自动变化检测的相关研究逐渐增多，提出了多种实现方法，但是相对于人工判 另IJ，计算机始终存在局限性，检测结果的精度收到限制。因此，研究高效可行的图像变化检 测算法具有十分重要的理论和研究价值。
[0004] 在计算机视觉和模式识别领域，视觉注意力模型的研究得到了很大的发展。使用 计算机进行人类视觉注意力系统模拟的模型，能够将人眼所能观察到的感兴趣区域从一幅 图像中提取出来，而对于变化检测问题，通常变化区域的面积小于非变化区域的面积，因 此，变化区域能够被视觉注意力模型当作显著区域被进一步突出。近年来，已有学者从不同 的角度对视觉注意力模型应用于图像变化检测进行了探讨，典型方法有W下两种：（1)基于 ITTI模型的遥感图像变化检测方法，该方法通过模拟人眼视觉注意机制，来进行变化检测 研究；（2)基于视觉注意和上下文的非监督变化检测方法，能够有效降低噪声对变化检测的 影响。但是，W上两种方法都是直接对显著图进行分割得到检测结果，并且由于采用的ITTI 视觉注意力模型要做下采样处理，故降低了分辨率，显著图中变化区域轮廓不明显，所W检 测精度受到限制。
【发明内容】
[0005] 有鉴于上述问题，本发明期望提供一种图像变化检测方法，W降低误检率，提高变 化检测结果的精度。
[0006] 本发明实施例的技术方案是运样实现的：
[0007] -种图像变化检测方法，应用于同一地区不同时期的遥感图像，其特征在于，包括 W下步骤：
[000引通过对所述遥感图像的对数比值差异图进行双边滤波得到初始差异图；
[0009] 通过对所述初始差异图进行频率调谐显著性区域检测得到显著图；
[0010] 通过对所述显著图进行局部赌计算得到加赌显著图；
[0011] 根据小波变换法融合所述初始差异图和加赌显著图得到融合图；
[0012] 通过最大类间方差法对所述融合图进行分割得到最终的变化检测结果图。
[0013] 上述方案中，所述通过对数比值差异图进行双边滤波得到初始差异图的步骤，过 程为：
[0014] 对已配准的两幅遥感图像进行计算，得到对数比值差异图；
[0015] 对所述对数比值差异图进行双边滤波，得到初始差异图。
[0016] 上述方案中，所述对数比值差异图的计算公式为：
[0017]
[001引其中，&和拉分别为已配准的两幅遥感图像，Dl为所求得的对数比值差异图。
[0019] 上述方案中，所述双边滤波的公式为：
[0020]
[0021] 其中，Dl(X)为像素点X的对数比值差异图，Di(X)为像素点X的初始差异图；
[0022] 丑:中.归一化累擲HxH+當公古九，
[0023] 其中，基于空间距离的高斯权重cU，x)的计算公式天
[0024] 巧中，某干像素间巧似忡度的高斯权軍S(化化），化(X))的计算公式为：
[0025] 其中，da, X)和S(化化），化(X))分别表示图像像素点之间的空间距离和灰度距离， Od和Or分别表示控制空间距离的参数和控制灰度距离的参数。
[0026] 上述方案中，所述通过对所述初始差异图进行频率调谐显著性区域检测得到显著 图的步骤，包括：
[0027] 步骤1:使用高斯滤波器对所述初始差异图进行平滑处理；
[002引步骤2:将平滑处理后的图像从RGB颜色空间转换到Lab颜色空间，W获得亮度特征 和颜色特征；
[0029] 步骤3:计算各亮度特征和颜色特征在整幅图像上的均值；
[0030] 步骤4:在各亮度特征和颜色特征下，计算每个像素点与整幅图像均值之间的距 离；
[0031] 步骤5:融合所有亮度特征和颜色特征，得到每个像素点的最终显著性值，从而得 到显著图。
[0032] 上述方案中，所述通过对所述显著图进行局部赌计算得到加赌显著图，求取过程 为：
[0033] 确定局部窗口的范围，定义为MXN;
[0034] 计算所述局部窗口M X N内的灰度分布，记为Pij;
[0035] 利用所述灰度分布PU计算图像中像素点（i，j)的局部赌，公式为：
[0036] 得到全图局部赌，记为加赌显著图。
[0037] 上述方案中，所述根据小波变换法融合所述初始差异图和加赌显著图得到融合 图，方法如下：
[0038] 选择化ar小波基对所述初始差异图和加赌显著图分别进行小波变换，分解层数为 两层；
[0039] 使用加权平均法对小波变换后的图像进行融合；
[0040] 对融合后的图像进行反变换，得到最终融合图。
[0041 ]上述方案中，所述加权平均法的融合规则为：
[0042] 低频系数公式为：Cj+l,F(m，n)=0.75XCj+l,I(m，n)+0.25XCj+l,s(m，n)，其中，Cj+l,F (m，n)，Cj+i,i(m，n)和Cj+i,s(m，n)分别表示融合图化、初始差异图化和加赌显著图化在像素点 (m，n)处的第j+1层的小波变换近似系数；
[0043] 高频系数为加赌显著图经小波变换后的高频系数。
[0044] 有益效果:本发明的方法基于显著性检测和小波融合，将基于频率调谐的显著性 检测模型引入到变化检测问题中，该模型不用进行下采样的计算，输出的显著图和原始图 像大小一致，同时，该方法将差异图与显著图进行融合，不仅保留了变化信息的细节，同时 突出了变化区域，有利于下一步的提取。实验结果表明，本发明方法能够有效地降低误检 率，提高了变化检测结果的精度。
【附图说明】
[0045] 图1为本发明实施例提供的图像变化检测方法的流程图；
[0046] 图2为具体实施方案中在化tawa真实遥感图像数据集中五种方法变化检测结果对 比图，（a)~(e)分别表示SOBO方法、LOGO方法、NOENO方法、ITTI方法和本发明方法。
[0047] 为了能明确实现本发明的实施例的结构，在图中标注了特定的尺寸、结构和器件， 但运仅为示意需要，并非意图将本发明限定在该特定尺寸、结构、器件和环境中，根据具体 需要，本领域的普通技术人员可W将运些器件和环境进行调整或者修改，所进行的调整或 者修改仍然包括在后附的权利要求的范围中。
【具体实施方式】
[0048] 在W下的描述中，将描述本发明的多个不同的方面，然而，对于本领域内的普通技 术人员而言，可W仅仅利用本发明的一些或者全部结构或者流程来实施本发明。为了解释 的明确性而言，阐述了特定的数目、配置和顺序，但是很明显，在没有运些特定细节的情况 下也可W实施本发明。在其他情况下，为了不混淆本发明，对于一些众所周知的特征将不再 进行详细阐述。
[0049] -、发明方案
[0050] 为了解决现有图像变化检测方法存在的技术问题，本发明实施例期望提供一种图 像的变化检测方法，应用于同一地区不同时期的两幅遥感图像，如图1所示，所述方法包括 W下步骤：
[0051] (1)通过对两幅遥感图像的对数比值差异图进行双边滤波得到初始差异图；
[0052] 本步骤中，首先对已配准的两幅遥感图像进行计算，得到对数比值差异图；然后对 上述对数比值差异图进行双边滤波，在有效去除噪声的基础上，保留了图像边缘信息，从而 得到初始差异图。
[0053] (2)通过对上述初始差异图进行频率调谐显著性区域检测得到显著图；
[0054] 在初始差异图中，变化区域像素具有较大值，非变化区域像素值接近于零值，并且 通常变化区域在整幅图像中所占面积比较小，因此，变化区域能够被视觉注意力模型当作 显著区域而检测出来。
[0055] 基于此，本发明实施例将频率调谐显著性区域检测方法引入到变化检测领域，来 模拟人类视觉对变化区域的感知，W弥补计算机相对于人类大脑存在的局限性，提高检测 结果的精度。该方法具体实施步骤为:首先对图像进行高斯平滑处理，然后将Lab颜色空间 作为图像特征，计算每个颜色通道中某个像素与整幅图像的平均色差之间的欧几里得距 离，最后将运3个通道值相加，得到的结果作为该像素的显著性值。
[0056] (3)通过对所述显著图进行局部赌计算得到加赌显著图；
[0057] 局部赌反映了图像灰度的离散程度，而与图像的灰度均值无关，依照局部赌值，可 W从图像中把相对均匀的目标分割出来。
[0058] 在生成的显著图中，变化区域进一步地突出，但是一些小目标和边缘信息随之减 弱，为了弥补该缺点造成的结果不准确，本发明实施例对显著图进行局部赌计算，得到加赌 显著图，具体实现方式为:确定局部窗口的范围，计算该局部窗口内的灰度分布，利用上述 灰度分布计算图像上所有点的局部赌。
[0059] (4)根据小波变换法融合初始差异图和加赌显著图得到融合图；
[0060] 小波变换是一种多尺度的分析方法，具有时频局部化和多分辨率特性，因而它可 W同时进行时域和频域的分析。
[0061] 为了得到既能突出变化区域，又可W保留细节信息及边缘信息的差异图，本发明 实施例使用基于小波变换的融合方法来融合初始差异图和加赌显著图，该融合方法减少了 层间的相关性，能够获得所需的融合效果，具体实现方式为:首先选择化ar小波基对图像进 行小波变换，再使用加权平均法对上述小波变换后的图像进行融合，最后对融合图像进行 反变换，得到融合图。
[0062] (5)通过最大类间方差法对融合图进行分割得到最终的变化检测结果图。
[0063] 具体地，步骤(1)包括：
[0064] (Ia)对已配准的两幅遥感图像进行计算，得到对数比值差异图；
[0065] Xi和拉为已配准的两幅遥感图像，分别代表不同时间、同一地理位置所获得的遥感 图像。按照W下公式得到对数比值差异图化：
[0066]
[0067] (化)对上述对数比值差异图进行双边滤波，得到初始差异图。
[0068] 双边滤波算法除了使用像素间几何上的靠近程度外，还考虑了像素之间像素值的 差异，因此，双边滤波不仅能够有效地去除图像噪声，同时还保留了边缘信息。同时，由于图 像噪声对变化检测的结果有很大的影响，而传统的低通滤波器在滤波的过程中又会不可避 免地丢失边缘信息，因此本发明选用双边滤波对对数比值差异图进行降噪，得到初始差异 图。
[0070]
[00例已知，Dl(X)为像素点X的对数比值差异图，Di(X)为像素点X的初始差异图，双边滤 波过程表示为：
[0073] 其中，基于巧间距离的高斯权重cU，x)的计算公式为：
[0071]
[0072]
[0074]
[0075] 其中，基于像素间相似性度的高斯权重S(化化），化(X))的计算公式为：
[0077」 具甲，像奈间置间化尚卿像奈化尚护」订昇公式为：
[0076；
[007引
[0079] 其中，Od为控制空间距离的参数，Or为控制灰度距离的参数。
[0080] 具体地，步骤(2)包括：
[0081 ] (2a)伸巧高斯滤妮器对輪入图像I进行平滑处理，得到输出图像Ig,公式为：
[0082]
[0083] 式中，（i，j )表示像素坐标，符号敏表示卷积运算，G表示高斯滤波器滤波算子，其 尺度大小一般选择为3 X 3;
[0084] (2b)将图像Ig从RGB颜色空间转换到Lab颜色空间，从而获得亮度特征L和颜色特 征a, b，特别地，对于待灰度遥感图像来说，仅需使用一维的灰度特征作为亮度特征来进行 计算；
[0085] (2c)计算每个特征在整幅图像上的均值，公式为：
[0089]其中，彩色图像I的大小为HXW，其中H为高度，W为宽度。
[0086]
[0087]
[008引
[0090] (2d)在各特征下，计算每个像素点与整幅图像均值之间的距离，公式为：
[0091] SlQ, j) = (X(i, j)-UL)2
[0092] Sa(i，j) = (a(i，j)-Ua)2
[0093] Sb(i，j) = (b(i，j)-叫）2
[0094] (2e)融合上述所有特征，得到每个像素点的最终显著性值，公式为：
[0095] S(i，j)=SL(i，j)+Sa(i，j)+Sb(i，j)
[0096] 最后得到显著图。
[0097] 具体地，步骤(3)包括：
[0098] 设DsQ，j)为显著图中像素点（i，j)处的灰度，对于图像中一个MXN大小的局部窗 口，定或其局部'協为，
[0099]
[0100] I：为：
[0101]
[0102] 推广得到全图局部赌，记为加赌显著图。
[0103] 具体地，步骤(4)包括：
[0104] (4a)选择化ar小波基对图像进行小波变换，分解层数为两层，采用二维Mallat算 偏？ _要完加1下.
[0105]
[0106] 其中，Hr、Gr分别表示镜像共辆滤波器H和G在行上的作用分量，出、Gc分别表示镜像 共辆滤波器H和G在列上的作用分量，Cj为第j层的小波系数，Cw为第j+1层的小波变换近似 系数，为计算得到的小波变换细节系数，分别表示行方向、列方向和对 角线方向的细节。
[0107] 进行二维小波变换时，首先进行行方向的一维小波变换，得到左右两张子图，再对 两张子图分别进行列方向的一维小波变换，共得到四张变换后图像，分别代表低频分量、水 平方向高频分量、垂直方向高频分量W及对角线方向高频分量。
[0108] (4b)使用加权平均法对上述小波变换后图像进行融合，融合规则为：
[0109] 低频系数公式：
[0110] Cj+i'F(m，n) = 0.75 XCj+1'i(m，n)+0.化 XCj+i's(m，n)
[0111 ] 其中，Cj+i,F(m，n)，Cj+i,i(m，n)和Cj+i,s(m，n)分别表不融合图Df、初始差异图Di和加 赌显著图化在像素点(m，n)处第j+1层的小波变换近似系数；
[0112] 高频系数即为加赌显著图经小波变换后的高频系数。
[0113] M一对融!会図條袖:軒亦拍L _值至Il县激融!会図舌於]管油
[0114]
[011引其中办麻为H，G的共辆转置矩阵。
[0116] 具体地，步骤(5)的实现方法为：
[0117] 由于最大类间方差法在图像处理领域的广泛应用及稳健表现，本发明采用该方法 来对融合后的图像进行分割，进而提取出变化区域。
[011引二、实验验证
[0119] (1)实验数据说明
[0120] 本实验采用五组遥感图像数据集，包括模拟数据集及基准变化图和真实遥感图像 数据集。其中，真实遥感图像数据集包括Mexico数据集、Sardinia数据集、Bern数据集和 Ottawa数据集和它们的基准变化图。
[0121] (2)验证方法说明
[0122] 对W下五种遥感图像变化检测方法进行对比实验：
[0123] a)为了体现本发明方法良好的差异图构造能力，将由图像差值法获得的差异图直 接进行大津阔值分割得到的结果与本发明方法的结果进行比较。图像差值图进行大津阔值 分割的方法记为SUBO方法。
[0124] b)为了说明本发明方法改进了传统对数比值差异图构造方法，将由对数比值法获 得的差异图直接进行大津阔值分割得到的结果与本发明方法的结果进行比较。对数比值法 进行大津阔值分割的方法记为LOGO方法。
[0125] C)为了体现本发明方法对显著图计算局部赌的作用，将不加赌法获取的差异图直 接进行大津阔值分割得到的检测结果与本发明方法的结果进行对比，该方法记为NOENO方 法。
[0126] d)为了说明本发明方法进行变化检测的效果，将基于视觉注意机制ITTI模型的方 法与本发明方法进行对比，该方法记为ITTI方法。
[0127] e)本发明所述的基于显著性和小波融合的变化检测方法。
[0128] (3)实验结果及分析
[0129] 分别使用上述五种方法在五组遥感图像数据集上进行测试。图2所示为使用五种 方法在化tawa真实SAR图像数据集上进行变化检测得到的结果图像，从图中可W直观地看 出：SUBO方法的虚警过多检测效果较差;L0G0、N0EN0方法的效果较为接近，存在一定的虚警 和噪声点；ITTI方法虚警数较少，但存在较多漏检的情况，丢失了很多变化区域的细节信 息;而本发明所述方法在背景噪声的抑制和变化区域轮廓及内部细节的检测等方面都要优 于其他四种方法。
[0130] 夫1节种方巧对节纽违胎掛据隹的结罢忡能评价夫
[0131]
[0132]
[0133] 从数据角度分析，如表1所示，对于模拟遥感图像数据集，五种方法的检测结果正 确率都比较高，其中ITTI方法和本发明所述方法表现最优;对于真实遥感图像数据集，本发 明的方法体现出较强的优越性：由实验结果数据可W看出，在保证漏检数在可接受范围内 的同时，本发明所述方法的虚警数较少，在上述四个真实遥感图像数据集中总错误数最少、 正确率表现最优。
[0134] 表2五种方法在五组实验数据集上的平均正确率表 「01351
[0136] 再由表2可知，在五个实验数据集上，本发明方法的平均正确率分别比洲BO方法、 LOGO方法、NOENO方法和口TI方法高5.92%、0.36%、0.43和0.56%，说明了本发明方法能够 有效地提高变化检测结果的精确度。
[0137] 由上述的实验结果分析可知，本发明提出的基于显著性检测和小波融合的遥感图 像变化检测方法可W减少漏检数，能够有效地检测到变化信息，并且在真实遥感图像数据 集上表现出更强的优越性，适宜于实际应用。
[0138] S、方案及实验总结
[0139] 本发明将基于频率调谐的显著性检测模型引入到遥感图像的变化检测问题中，提 出了一种基于显著性检测和小波融合的遥感图像变化检测方法，运种方法将对数比值差异 图和加赌显著图进行融合，从而得到一幅既能有效地突出变化区域，又能保留细节及边缘 信息的差异图，为后续的变化区域提取提供了有利的条件。实验结果表明，本发明方法能够 有效地降低误检率，提高了变化检测结果的精度。
[0140] 在本申请所提供的具体实施方案中，应该理解到，所掲露的方法，可W通过其它的 方式实现。
[0141] 本领域普通技术人员可W理解:实现上述方法实施例的全部或部分步骤可W通过 程序指令相关的硬件来完成，前述的程序可W存储于一计算机可读取存储介质中，该程序 在执行时，执行包括上述方法实施例的步骤;而前述的存储介质包括:移动存储设备、只读 存储器（Read-Only Memoir ,ROM)、随机存取存储器（Random Access Memory ,RAM)、磁碟或 者光盘等各种可W存储程序代码的介质。
[0142] W上所述，仅为本发明的【具体实施方式】，但本发明的保护范围并不局限于此，任何 熟悉本技术领域的技术人员在本发明掲露的技术范围内，可轻易想到变化或替换，都应涵 盖在本发明的保护范围之内。因此，本发明的保护范围应W所述权利要求的保护范围为准。
【主权项】
1. 一种图像变化检测方法，应用于同一地区不同时期的遥感图像，其特征在于，包括以 下步骤： 通过对所述遥感图像的对数比值差异图进行双边滤波得到初始差异图； 通过对所述初始差异图进行频率调谐显著性区域检测得到显著图； 通过对所述显著图进行局部熵计算得到加熵显著图； 根据小波变换法融合所述初始差异图和加熵显著图得到融合图； 通过最大类间方差法对所述融合图进行分割得到最终的变化检测结果图。2. 根据权利要求1所述的方法，其特征在于，所述通过对数比值差异图进行双边滤波得 到初始差异图的步骤，过程为： 对已配准的两幅遥感图像进行计算，得到对数比值差异图； 对所述对数比值差异图进行双边滤波，得到初始差异图。3. 根据权利要求2所述的方法，其特征在于，所述对数比值差异图的计算公式为：其中，心和办分别为已配准的两幅遥感图像，Dl为所求得的对数比值差异图。4. 根据权利要求2所述的方法，其特征在于，所述双边滤波的公式为：其中，DUx)为像素点X的对数比值差异图，D:(x)为像素点X的初始差异图； 其中，归一化系数k(x)计算公式为其中，基于空间距离的高斯权重cU，x)的计算公式为其中，基于像素间相似性度的高斯权重s(DlU)，Dl(x))的计算公式为：其中，(Κξ，χ)和S(DlU)，Dl(x))分别表示图像像素点之间的空间距离和灰度距离，〇d和 分别表示控制空间距离的参数和控制灰度距离的参数。5. 根据权利要求1所述的方法，其特征在于，所述通过对所述初始差异图进行频率调谐 显著性区域检测得到显著图的步骤，包括： 步骤1:使用高斯滤波器对所述初始差异图进行平滑处理； 步骤2:将平滑处理后的图像从RGB颜色空间转换到Lab颜色空间，以获得亮度特征和颜 色特征； 步骤3:计算各亮度特征和颜色特征在整幅图像上的均值； 步骤4:在各亮度特征和颜色特征下，计算每个像素点与整幅图像均值之间的距离； 步骤5:融合所有亮度特征和颜色特征，得到每个像素点的最终显著性值，从而得到显 著图。6. 根据权利要求1所述的方法，其特征在于，所述通过对所述显著图进行局部熵计算得 到加熵显著图，求取过程为： 确定局部窗口的范围，定义为MXN; 计算所述局部窗口 Μ X N内的灰度分布，记为pij; 利用所述灰度分布PU计算图像中像素点（i，j)的局部熵，公式为：得到全图局部熵，记为加熵显著图。7. 根据权利要求1所述的方法，其特征在于，所述根据小波变换法融合所述初始差异图 和加熵显著图得到融合图，方法如下： 选择Haar小波基对所述初始差异图和加熵显著图分别进行小波变换，分解层数为两 层； 使用加权平均法对小波变换后的图像进行融合； 对融合后的图像进行反变换，得到最终融合图。8. 根据权利要求7所述的方法，其特征在于，所述加权平均法的融合规则为： 低频系数公式为：(^+1,卩(111，11)=0.75\(^+1,1(111，11)+0.25\(^ +1,5(111，11)，其中，(^+1,卩(111， n)，Cj+i,i(m，n)和Cj+i,s(m，n)分别表不融合图Df、初始差异图Di和加熵显著图Ds在像素点(m， η)处的第j+1层的小波变换近似系数； 高频系数为加熵显著图经小波变换后的高频系数。
【文档编号】G06T7/00GK105957054SQ201610245590
【公开日】2016年9月21日
【申请日】2016年4月20日
【发明人】张浩鹏, 姜志国, 黄洁, 史振威, 谢凤英, 赵丹培, 尹继豪, 罗晓燕
【申请人】北京航空航天大学